x1=-6:0.1:6;
x2=x1;
[X,Y]=meshgrid(x1,x2);
Z=X.*Y;
mesh(X,Y,Z);
G=X.^2+Y.^2-25;
hold on;
contour(X,Y,G,[0 0]);
