由对角特征值生成的3d曲面图（Matlab错误）

我有以下循环：我已经求解了特征值，但希望为循环的每个连续迭代生成e，r和c（I，j）的曲面图。我知道每个输入都是一个标量，但我能生成一个3×3的对角矩阵D。我怎么做曲面图？ 这是我的代码：（我不知道如何上传matlab代码，所以这是图片）

---

将代码粘贴为文本而不是图像：无论循环如何，如何使用矩阵D中指定的特征值以及数量e和r生成曲面图？检查

surf

所需的输入尺寸。您正试图根据x和y标量绘制一个3x3 Z矩阵。如果您将执行

surf（ee，er，c）

它将绘制曲面（不确定是否为所需曲面）。将代码粘贴为文本而不是图像：无论循环如何，如何使用矩阵D中指定的特征值和数量e和r生成曲面图？检查

surf

所需的输入尺寸。您正试图根据x和y标量绘制一个3x3 Z矩阵。如果您将执行

surf（ee，er，c）

操作，它将绘制曲面（不确定是否为所需曲面）。

L=1.25;            
k=1.5;             
m=0.2;              
M=45;                          
p=0.625;            
Ig=(1/12)*M*L^(2);

for i=1:3

for j=1:3

ee(i)=0.1 + i.*0.1;
er(j)=0.1+j.*0.1;

e=ee(i);r=er(j);
Q=(((m*(r.^(2))*(L^(2))*k)./e)-M.*p.*k-m.*k.*L.*r.^(2));
N=(M*(p^(2))+Ig+(m.*(r.^(2)).*(L^(2)))/4*e.^(2));
J=-(m.*(r.^(2))*(k^(2))*L)./e; 
y=-k./e;
x=(L./e)-1;
A=[0,1,0; 0,Q./N,J./N;0,x,y];
[V,D] = eig(A);
c(i,j)=D(2,2);
u(i,j)=D(3,3)

end
end
surf(e,r,c)