function v=dg1(x)
v(1,1)=exp(x(1))*cos(x(2));
v(1,2)=exp(x(1))*sin(x(2));
v(2,1)=-exp(x(1))*sin(x(2));
v(2,2)=exp(x(1))*cos(x(2));
