function dy=g1a(x,y)
k=.5;
g=40;
a=k+(1-k)/g;
m=(1-k)*(1-a);
n=a*(1-k);
p=k*(1-a);
sc=log((1-k)/(1-a))-(k/(1-k))*log(a/k);
%c=n/p





dy=(m*y^2+n*y)/(n-p*y);