%% Pintapiirros
% luodaan piirtomuuttujat
t1 = linspace(-2,2,30);
t2 = linspace(-1,1,20);

[x y] = meshgrid(t1,t2);
%% lasketaan ja piirretään pinta
z = sin(3*y-x.^2+1)+cos(2*y.^2+2*x);
figure(1);
surf(x,y,z)
figure(2)
%% Korkeuskäyrät 
contour(x,y,z);
[px,py] = gradient(z,0.2);
hold on 
quiver(x,y,px,py)
%%