For the circle (disk) do this:
[T,R] = meshgrid(linspace(0,2*pi,64),linspace(0,2,16));
X = R.*cos(T);
Y = R.*sin(T);
Z = X.^2-Y.^2;
surf(X,Y,Z)
For the annulus, change first line to:
[T,R] = meshgrid(linspace(0,2*pi,64),linspace(1,3,16));
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!