[x,y] = meshgrid(0:0.1:3);
f=x.*exp(y);
[c,h] = contour(x,y,f,1:1:10);
clabel(c,h);