Matlabの3Dプロット

Question

3D極座標の図を極座標から直交座標に変換し、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);