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

gh784478335

syms x  z q   m t   c   
m=0.5;
R=40;
( N5 B# Q1 o. ]) k0 ~6 G* ?& oq=200;
' q8 B4 [; {. b8 j+ L. we=1;f=1;
Z= []
for X=[-55:2:55];
1 L& l* y9 T# U( {7 L  for Y=[-55:2:55]                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
+ Z+ n0 @2 k7 ?' Psyms r y
+ r0 K  P' J1 u% P. Y  }% ? if ((X.*X)./3600+(Y.*Y)./400)<1;
; j5 b5 _+ q! `! w- w& z
; o" l* i: ^- \# KZ(e,f)=nan;
: Y2 H( V/ L9 `6 f% ~
f=f+1;
continue
      else
      a=X.*X
     b=Y.*Y
3 r& F0 W' p8 v( x0 T     
" O' S, e/ x- X( x5 p! U" o! v% g     end
. |! Q0 N$ [( }6 o$ T      [r]=solve((a./(R.^2*(1./r+m.*r).^2)+b./(R.^2.*(1./r-m.*r).^2))==1,r);
: j  s- v$ k  T0 K. p! f      [y]=solve((a./(4.*R.^2.*m.*(cos(y)).^2)-b./(4.*R.^2.*m.*(sin(y)).^2))==1,y);
0 R2 a( H+ t# t( v4 N# O   
0 z! P4 Z6 y1 q# J       y=vpa(y);
9 l) n) D" f4 t6 I! |+ B  dJ=y(y>0);
9 D7 s$ p& H, n5 Y' CK=J(J<2*pi);
+ o! K( {5 G& }# ^" `; a" uL=K(imag(K)==0);

r=vpa(r);
G=r(r>0);
) _3 T/ r' j& WH=G(G<1);
9 a( \, U7 ~- ]! U$ lI=H(imag(H)==0);

! ^6 R- z4 H3 h- q* Q& D1 I

9 \" ]7 _" l% J# `( M, f6 Zx=0.5*pi;

# e6 ^2 f/ \6 b1 s' V, Q9 J/ nA=q*real(((2.*exp(2.*i.*x)-m).*(I.*exp(i.*L)).^2-1)./(m.*(I.*exp(i.*L)).^2-1));
B=((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));
, Z! y! Z, S4 x6 {4 U
/ n3 Q) v! f; \4 [* C$ vO=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));
4 W( l9 c6 D  f0 x- ]' H
3 @5 f' ?; ^5 s; q: RD=(A+real(O))./2;
E=(A-real(O))./2;
$ v# n' \* l& ?: m" z2 `% rF=sqrt(D.^2+E.^2);
& _: g1 a' ]/ s. Y0 G% pZ(e,f)=E
  Y! Y3 S" {. d$ P* g! ]+ [
f=f+1;
  end
& J0 ^) B1 |) A/ {) x5 R3 R& n  e=e+1;
/ [/ P" }! N8 D  w4 U  f=1;
/ k1 C4 v0 f" g$ \end
* h9 v( D; C/ g1 e/ S[d]=[-55:2:55];
) o7 n5 X# H. ~9 L7 d; y- ~5 N[c]=[-55:2:55];
7 E0 a1 X  j$ ~5 B) t; p! \%idy =d.^2./3600+c.^2./400< 1;
%d(idy) = nan;
%c(idy) = nan;

) G1 ]# }( \, [- J  X9 r0 g+ m" xfigure(3)
+ L) O2 C, ~! {1 Z$ k# tsuRF(d,c,Z)
2 d3 P! b6 \: d# L%surf(d,c,Z)
%view([0,1,0]);
shading interp;  
, o4 h7 P) \( v7 o" ?colorbar;  
) g. J2 G8 ]1 n' D8 l%axis equal;
& @) \$ z; N% y* y这个函数运行后可以产生Z的矩阵值，但是surf不显示图像，但我将xy的范围定义在-45到45的时候会显示图像，请问是什么问题

Liberallh21

没看懂