% plotfunction

X=linspace(-3,3,50)';             %X wordt lineair verdeeld over [-3,3] in 39 stukken
Y = X;                            % ook Y=[-3:.1:3]'; werkt beter bij bekende stapgrootte
[x,y]=meshgrid(X,Y);              %maak rooster van x- en y-waarden
[x,y]=meshgrid(-3:.2:3,-3:.2:3);  %kan natuurlijk ook z=peaks(x,y);
z=y.*y+x.*y;

mesh(x,y,z)                       %3-D plaatje
axis tight                        %soms zijn assen aan de ruime kant
figure(1)                         %brengt figuur window naar voren (i.p.v. klikken met muis)
pause(2)  % pause 2 seconden')


meshc(x,y,z); figure(2)
%pause(2)

contour(x,y,z); figure(2)
%pause(2)

%[px,py]=gradient(z,.2,.2);
%contour(x,y,z) 

%hold on  % houd grafiek vast

%quiver(x,y,px,py)
%figure(3)
%hold off