phi=@(x,d)[x-d*(2*cos(x)-exp(x))./(-2*sin(x)-exp(x))]