%% H3teht4
%{
Olkoon
f(x,y)=sin(3y - x^2+1)+cos(2y^2-2x). Piirrä pintakuva ja korkeuskäyräpiirros, 
jälkimmäinen sekä contour että ezcontour-funktioilla. Tässä on mahdollisuus 
kokeilla yllä esiteltyjä korkeuskäyrien valitsemistapoja, myös clabel. 
Ota alueeksi vaikka  [-2 2 -1 1] 
%}
%% 1) ezcontour
close all
ezcontour('sin(3*y - x^2+1)+cos(2*y^2-2*x)',[-2 2 -1 1])
% Onpas helppoa! 
%% 2) contour  % Lisää joustavuutta ja mahdollisuuksia.

figure
x=-2:.01:2; y=-1:.01:1;
[X,Y]=meshgrid(x,y);
Z=sin(3*Y - X.^2+1)+cos(2*Y.^2-2*X);
contour(x,y,Z,20)    % 20 korkeuskäyrää.

% Antaa seuraavan pyyhkiä edellinen.

Z_vali=[min(min(Z)) max(Z(:))]  % Tahallaan kahdella tavalla
korkeudet=linspace(Z_vali(1),Z_vali(2),30);
contour(x,y,Z,korkeudet)   % Korkeuskäyrät näille korkeuksille.

[C,h]=contour(x,y,Z,korkeudet);
clabel(C,h,korkeudet(1:6:30))   % Osalle käyristä numeroarvot

clabel(C,h,'manual')    % Lisätään arvoja manuaalisesti
