%定义第八曲面的相关数据。 9~/k25P  
C=22.334133;                                                               %透镜曲率 7o_1PwKS6  
r=(-2.881923:0.001:2.881923);                                              %Semi-meter     '[Ap/:/UY  
k=-1000;                                                                   %圆锥系数 ?Z9C}t]  
a1=0.26015;                                                               %2nd Order Term [H<![Z1*r  
a2=-0.045105;                                                              %4nd Order Term |`B*\\ 1  
a3=0.003017477;                                                            %6nd Order Term Up*1j:_O  
a4=-0.0001295027;                                                           %8nd Order Term K%h9'}pq>1  
z=C.*r.^2./(1+sqrt(1-(1+k)*C^2.*r.^2))+a1.*r.^2+a2.*r.^4+a3.*r.^6+a4.*r.^8; FrT.<3  
plot(z,r,'r'),grid on,axis  2.)xWCG  
 _4Z|O]  
从你的公式来看,你的c应该取曲率,而不是半径,没有验证过,你可以试试。