

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



function dy=g100(x,y)

dy=(.2438*y^2+.2563*y)/(.2563-.2438*y);
%tol=10^(-4);
%[x,yn]=ode23('g100',0,6.252*10^(-4)-.00001,1,tol,1);
%plot(x,yn,'o');grid
%hold on
%plot(x,dy,'.r')