急求，为何运行后surf不显示图像

gh784478335

syms x  z q   m t   c   
+ x- ]' p, i! Z* M" Lm=0.5;
R=40;
/ l# R( m  _$ U, ^, Vq=200;
e=1;f=1;
Z= []
for X=[-55:2:55];
+ L/ u; s* _8 S; h9 ?  for Y=[-55:2:55]                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
syms r y
if ((X.*X)./3600+(Y.*Y)./400)<1;
$ @% _3 p2 t& {  n% C# Y
  j3 z0 R- }& D3 r' q# bZ(e,f)=nan;
; l6 u, k$ Y% n7 C0 b7 x
f=f+1;
. ^- c# X: g( v% scontinue
. k$ Z$ G" D4 k* J2 I1 a      else
1 U7 {/ ^0 S3 O      a=X.*X
2 F; L- O3 M8 ?3 C     b=Y.*Y
     
     end
      [r]=solve((a./(R.^2*(1./r+m.*r).^2)+b./(R.^2.*(1./r-m.*r).^2))==1,r);
      [y]=solve((a./(4.*R.^2.*m.*(cos(y)).^2)-b./(4.*R.^2.*m.*(sin(y)).^2))==1,y);
5 w1 l# r6 n3 _4 J4 Z  W  ^0 q   
% p. j  m6 q2 ]1 S9 G       y=vpa(y);
# D' K5 Q# x7 s; q  bJ=y(y>0);
K=J(J<2*pi);
7 @( Z: N" N& tL=K(imag(K)==0);

r=vpa(r);
3 h$ F# ~" L& JG=r(r>0);
H=G(G<1);
I=H(imag(H)==0);
0 _# V* j, n; a* n/ |% G

! i* N( b$ r% P$ `9 s* @
x=0.5*pi;
/ y4 r1 ]' y3 \2 f. |8 A+ V  _
, u9 n7 ~5 q0 a/ ~0 @0 DA=q*real(((2.*exp(2.*i.*x)-m).*(I.*exp(i.*L)).^2-1)./(m.*(I.*exp(i.*L)).^2-1));
' g; E3 M& F- ~& U. TB=((q*(m*I^4+(I*exp(i*L))^2)*(I*exp(i*L))^2)./(I^4*(m-(I*exp(i*L))^2./I^4)*(m*(I*exp(i*L))^2-1))*((2*exp(2*i*x)-m+m*(1+m*(I*exp(i*L))^2-2*exp(2*i*x)*(I*exp(i*L))^2)/(m*(I*exp(i*L))^2-1)))+q/(I^2*(m-(I*exp(i*L))^2/I^4))*(exp(-2*i*x)-((3*exp(2*i*x)*(I*exp(i*L))^2+m*exp(2*i*x)-m^2-1)*(I*exp(i*L))^2)/(m*(I*exp(i*L))^2-1)+((exp(2*i*x)*(I*exp(i*L))^2+m*exp(2*i*x)-m^2-1)*(I*exp(i*L))^4)*2*m/(m*(I*exp(i*L))^2-1)^2));

3 C1 ?/ b2 v2 A. \O=B./((I.*exp(i.*L)).^2.*(R.*(m-1./(I.*exp(i.*L)).^2))).*(I.^2.*R.*(m-(I.*exp(i.*L)).^2./I.^4));

6 O: }$ \& C; k0 K7 |D=(A+real(O))./2;
& s8 o/ }+ w7 }& X! @* LE=(A-real(O))./2;
5 y8 `$ g/ b- ]/ `F=sqrt(D.^2+E.^2);
2 U0 I: Y: k: h/ b, d$ wZ(e,f)=E

( }' Z6 t3 `! V7 gf=f+1;
$ _. W7 q7 @' ?, H# {2 W/ j  end
- q/ J2 V) z5 p2 l  e=e+1;
  f=1;
end
[d]=[-55:2:55];
7 ]5 N7 P0 l$ u* L! ^2 l2 Z; a* W' l9 H[c]=[-55:2:55];
%idy =d.^2./3600+c.^2./400< 1;
%d(idy) = nan;
%c(idy) = nan;
# n$ R) t6 F1 g( x
figure(3)
' W( ~9 h5 K8 tsuRF(d,c,Z)
%surf(d,c,Z)
  P0 J! O9 M* J: Y- z6 D%view([0,1,0]);
shading interp;  
7 P" g6 h5 M) c/ _0 @colorbar;  
8 V/ o! I, t2 a; q6 W- ]%axis equal;
2 j0 a; u/ p5 w) K! b$ S$ U这个函数运行后可以产生Z的矩阵值，但是surf不显示图像，但我将xy的范围定义在-45到45的时候会显示图像，请问是什么问题

0 n" C5 J. r% e
* c0 h+ C& S7 w# _

Liberallh21

没看懂