%% H2T6R.m
%   M"a"aritell"a"an haluttu funktio:

f = @(x,y) sin(3*y-x.^2+1)+cos(2*y.^2-2*x);

% Funktion f on toimittava samankokoisille matriisiargumenteille x,y
% pisteitt"ain ("array smart").

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

%   meshgridin toiminnan ymm"art"a"a parhaiten esim. tutkimalla komennon
[X,Y]=meshgrid(0:3,-3:-1) % palautusarvoja X ja Y

% Perusteellisemmin skriptiss"a:  meshgridi.m 

% Palataan teht"av"a"an, x ja y muodostettiin edell"a:

[X,Y]=meshgrid(x,y);
Z=f(X,Y);

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

%   ja sama tasa-arvok"ayrin"a:
figure;
contour(x,y,Z);
