function out=dphi_dlambda(lambda,parameter1,parameter2)



out=(lambda+1)^2-parameter1-parameter2+xstar(lambda,parameter1,parameter2);
