theta = linspace(0,2*pi,25);
r = linspace(0,5,25);
[TH,R] = meshgrid(theta,r);
Z = sinc(R);  % Use matrix from MESHGRID
[x,y,z] = pol2cart(TH,R,Z);
surf(x,y,z)