function discrep=ckf(x,xl,a,alph,bet,gam,delt) %function discrep=ckf(x,xl,a,alph,bet,gam,delt) discrep=zeros(2,1); C=x(1);K=x(2);CL=xl(1);KL=xl(2); discrep(1)=C+K-exp(a)*KL^alph-delt*KL; discrep(2)=CL^(-gam)-(exp(a)*KL^(alph-1)*alph+delt)*bet*C^(-gam);