3D-сюжет в Matlab

Question

Я пытаюсь нарисовать фигуру трехмерного полярного графика, преобразовав ее из полярных в декартовые координаты и используя surf() для ее построения.

Проблема: Я получаю ошибку ??? CData must be an M-by-N matrix or M-by-N-by-3 array. Где я неправ? Я совершенно новичок в MATLAB и не понимаю, что происходит.

Код MATLAB

очистить все; закрыть все;

N=50;
%define matrices for ploting
x=zeros(N,N,N);
y=zeros(N,N,N);
z=zeros(N,N,N);
f=zeros(N,N,N);

%define basic input variables
r=linspace(0,2,N);
theta=linspace(0,pi,N);
phi=linspace(0,2.*pi,N);

for ii=(1:N) %use ii, jj to avoid confusion with the imaginary units
    for jj=(1:N)
        for kk=(1:N)
            x(ii,jj,kk)=r(ii).*sin(theta(jj)).*cos(phi(kk)); %%not using 
            %%for loop probably can work too, test later
            y(ii,jj,kk)=r(ii).*sin(theta(jj)).*sin(phi(kk));

            z(ii,jj,kk)=r(ii).*cos(theta(jj));

            f(ii,jj,kk)=r(ii).*exp(-r(ii)).*cos(theta(jj));
        end
    end
end

figure;
surf(x,y,z,f);
colormap([1,1,1]);

Answer 1

surf строит поверхность: каждая точка (x, y) имеет высоту, z.

Например,

N=50;
[X Y] = meshgrid(linspace(0,2,N));
Z = zeros(size(X));

for ii=(1:N)
    for jj=(1:N)
        x = X(ii,jj);
        y = Y(ii,jj);
        [theta, r] = cart2pol(x,y);
         Z(ii,jj) =r * exp(-r) * cos(theta);
    end
end

figure;
surf(X,Y,Z);