

%k=.04;
%g=90;
%a=k+(1-k)/g
%m=(1-k)*(1-a)
%p=k*(1-a)
%sc=log((1-k)/(1-a))-(k/(1-k))*log(a/k)
%yc=n/p
%se lanseaza ultimile doua randuri !function corespunde la k=.04 & g=90


function dy=g11(x,y)

dy=(.9114*y^2+.0486*y)/(.0486-.0380*y);
%[x,yn]=ode23('g11',0,.0013-.00001,1);
%plot(x,yn,'-.g');grid