function DF=newton_df(x)
pi2 = 0.5* pi;
DF(1,1) = 2*x(1);
DF(1,2) = 2*x(2);
DF(2,1) = pi2*cos(pi2*x(1));
DF(2,2) = 3*x(2)^2;
return