%   määritellään haluttu funktio:
f = @(x,y) sin(3*y-x.^2+1)+cos(2*y.^2-2*x);

%   haluttu xy-alue, jako 20 osaan:
x = linspace(-2,2,20);
y = linspace(-1,1,20);

%   meshgridin toiminnan ymmärtää parhaiten esim. tutkimalla 
%   [X,Y]=meshgrid(0:4) palautusarvoja
[X,Y]=meshgrid(x,y);
Z=f(X,Y);

%   3D-pinta:
figure;
surf(x,y,Z);

%   ja sama tasa-arvokäyrinä:
figure;
contour(x,y,Z);
