clear all
close all
% spiegazione di meshgrid
a=[1 2 3]
b=[11 12 13]
[c,d]=meshgrid(a,b)

% maniera alternativa di definire un vettore da
%  a a b con elementi spaziati di d: vettore=[a:d:b]
s1=[-1:0.1:1];
t1=[-1:0.1:1];
[s,t]=meshgrid(s1,t1);
x=s;
y=t;
z=s.^2+t.^2;
surf(x,y,z)

% apre una nuova figura (quindi con num 2 in questo caso)
figure

x=3*s;
y=t;
z=s.^2+t.^2;
surf(x,y,z)
% forza fgli assi a stare tutti nella stessa scala
axis equal

% apre una nuova figura con numero 5
figure(5)
axis equal

% disegna cilindro
[s,phi]=meshgrid(linspace(0,2,20),linspace(0,2*pi,40));
x=cos(phi);
y=sin(phi);
z=s;
surf(x,y,z)
axis equal
hold on

% disegna emisfero sopra cilindro
[theta,phi]=meshgrid(linspace(0,pi/2,20),linspace(0,2*pi,40));
x=sin(theta).*cos(phi);
y=sin(theta).*sin(phi);
z=2+cos(theta);
surf(x,y,z)
