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
tol=10^(-5);
d=[sc-10^(-9)]
e=[tol];
%a1=['[x,yn]'];
%a2=['=ode23('claus',0.0006152,1,.00001,1);'];
   %b=[a1 a2];
   %eval(b)
[x,yn]=ode23('claus',0,0.0006152,1,.00001,1)
%plot(x,yn,'o');grid
%hold on
%[x,yn]=ode23('claus',0,-0.1,1,.000001,1)
%plot(x,yn,'o');grid


