

%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)
%c=n/p



function dy=g10(x,y)

dy=(.225*y^2+.275*y)/(.275-.225*y);

%[x,yn]=ode23('g1a',0,.0101-.0001,1,.00001,1)
%plot(x,yn,'o');grid
%hold on
%[x,yn]=ode23('g1a',0,-0.1,1,.000001,1)
%plot(x,yn,'o');grid