从对角线特征值生成的3d曲面图（Matlab误差）

Question

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

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)