x1=[-5:.5:5]';
x2=x1;
[X1,X2]=meshgrid(x1,x2);
s=X1.^2+X1.*(X2.^2);
mesh(X1,X2,s);
xlabel('x1'); ylabel('x2'); zlabel('S(x)');
% cinvert
figure;contour(X1,X2,s,50);
xlabel('x1'); ylabel('x2');
