http://octave-online.net/
[X,Y] = meshgrid(-2:.2:2, -2:.2:2);
Z = X .* exp(-X.^2 - Y.^2);
mesh(X,Y,Z)
[X,Y] = meshgrid(-5:.5:5, -5:.5:5);
Z1 = X.^2 + Y.^2 -1;
mesh(X,Y,Z1);
f=@(X,Y) X + Y + 1;
[X Y] = meshgrid(-5:.5:5, -5:.5:5);
Z2 = f(X, Y)
mesh(X, Y, Z2, 'edgecolor', 'none', 'facecolor', 'r', 'facealpha', 0.3);
[X,Y] = meshgrid(-2:.2:2, -2:.2:2);
Z = X .* exp(-X.^2 - Y.^2);
mesh(X,Y,Z)
[X,Y] = meshgrid(-5:.5:5, -5:.5:5);
Z1 = X.^2 + Y.^2 -1;
mesh(X,Y,Z1);
f=@(X,Y) X + Y + 1;
[X Y] = meshgrid(-5:.5:5, -5:.5:5);
Z2 = f(X, Y)
mesh(X, Y, Z2, 'edgecolor', 'none', 'facecolor', 'r', 'facealpha', 0.3);