function A=dJack(sigma);
A=[1+dphi(sigma(1)),dphi(sigma(2));dphi(sigma(1)),1+dphi(sigma(2))];
