

%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
%xf=sc-10^(-9)

function dy=claus(x,y,m,n,p)
%z1=['dy'];
%z2=['=(m*y^2+n*y)/(n-p*y);'];
%w=[z1 z2];
%eval(w)
dy=(m*y^2+n*y)/(n-p*y);

