%not a function file
1;

% incorrected numbered axis
R = 1;
H = 1.5;
fx = @(th,rho) rho.*cos(th);
fy = @(th,rho) rho.*sin(th);
fz = @(th,rho) rho.*H/R;

figure;
ezsurf(fx,fy,fz,[0,1],[0,2*pi])

% Corrected numbered axis
rhovec = linspace(0,R,40);
thvec = linspace(0,2*pi,40);
[Rho,Th] = meshgrid(rhovec,thvec);
xx = Rho.*cos(Th);
yy = Rho.*sin(Th);
zz = Rho.*H/R;

figure;
surf(xx,yy,zz)


