

%k=.5;
%g=10;
%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);
%yc=n/p;
%sc=.0101
function dy=g1a(x,y)

dy=(m*y^2+n*y)/(n-p*y);
%[x,yn]=ode23('g1a',0,.0101-.0001,1);
%plot(x,yn,'o');grid