function g = gavu1(u1)
g=u1*exp(4*(1-u1)/(1+0.2*(1-u1)));