function df=dphi(sigma);
sigma=sigma.*(sigma<0);
df=0.5*(-1+cosh(sigma));
