function [g] = G(epsilonstar,epsilonmin,gamma2)
g=(1-(epsilonstar./epsilonmin).^(-gamma2)).*(epsilonstar>=epsilonmin)+0.*(epsilonstar<epsilonmin);