%% Example: meshgrid, surf, mesh contour
clc,close all
x=linspace(-2,2,30);y=x;
[X,Y]=meshgrid(x,y);
f=@(x,y) (x.*y.^2)./(x.^2+y.^4)
surfc(x,y,f(X,Y)); 
colorbar; shg
figure
contour(x,y,f(X,Y))
contour(x,y,f(X,Y),linspace(-1,1,20))
% tihennetään hilaa.
x=linspace(-2,2,60);y=x;
[X,Y]=meshgrid(x,y);contour(x,y,f(X,Y),linspace(-1,1,20));shg
[c,h] = contour(x,y,f(X,Y)); clabel(c,h), colorbar
%%
x=linspace(-2,2,6), y=linspace(-1,1,3)
[X,Y]=meshgrid(x,y) 

for i=1:6                     % Luodaan solumatriisi, jonka
   for j=1:3                  % alkioina on 2:n pituiset vektorit
      XY{i,j}=[x(i),y(j)];    % [x(i),y(j)];    
   end
end
XY=XY'                        % Transponoidaan, sillä meshgrid:n
cellplot(XY)                  % vaaka-akseli on x-akseli ja pysty y.

%g=inline('x(1).*x(2)','x');

XY{1,:}                       % Katsotaan solumatriisin 1. rivi

XY2(:,:,1)=X; XY2(:,:,2)=Y;   % Kootaan X- ja Y-"tasot" 3-ulotteiseksi
XY2                         
                               
% matriisiksi