[x,y]=meshgrid([-3:0.1:3],[-3:0.1:3]);
S=x.*y;
mesh(x,y,S);
[c,h]=contour(x,y,S,[-5 -2 0 3 8]); clabel(c,h);
diary off