function [y,dy]=fu3(x)
y(1)=x(1)^2+x(2)^2-4;
y(2)=x(1)*x(2)-1;
dy(1,1)=2*x(1);
dy(1,2)=2*x(2);
dy(2,1)=x(2);
dy(2,2)=x(1);