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

gh784478335

syms x  z q   m t   c   
% V% \1 H! S$ e) Vm=0.5;
% h1 U/ n1 W4 _+ u* \& ]% a: w8 cR=40;
- C3 U& Y( w' i7 a3 Qq=200;
e=1;f=1;
. e# y( f" W; {3 G! ~6 FZ= []
0 ^  T+ Q. F% F0 l6 c) ?4 ffor X=[-55:2:55];
4 N. ^* o* W8 ]0 s7 g3 |  for Y=[-55:2:55]                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
6 z: l5 E- S4 \# Zsyms r y
if ((X.*X)./3600+(Y.*Y)./400)<1;
: n& s8 n- O! k6 I# V0 |
8 Y/ m, v4 f5 iZ(e,f)=nan;
2 B( ?! m5 H! `. r$ w
' z) P" m% l$ q+ if=f+1;
continue
9 w1 [# g; O2 K9 w# |9 g      else
      a=X.*X
/ i( [* a/ Z" z* X0 `+ d+ @     b=Y.*Y
     
     end
+ \, u( c( |9 ]2 r& c      [r]=solve((a./(R.^2*(1./r+m.*r).^2)+b./(R.^2.*(1./r-m.*r).^2))==1,r);
" U" z" \$ Q2 v5 G# p- {0 Z: ^5 @      [y]=solve((a./(4.*R.^2.*m.*(cos(y)).^2)-b./(4.*R.^2.*m.*(sin(y)).^2))==1,y);
   
9 l! }  E# z. ?  O8 b+ J0 P$ O       y=vpa(y);
! G" o& J, @# _  dJ=y(y>0);
) q  Q2 [' L5 O: J& f  mK=J(J<2*pi);
8 I/ @, s1 D3 S% }5 c. h+ r: ]; uL=K(imag(K)==0);

0 [" y# b. ~/ N# c7 [r=vpa(r);
G=r(r>0);
H=G(G<1);
I=H(imag(H)==0);



. U+ g+ g) {( t5 k6 X$ M$ vx=0.5*pi;
* f9 C0 g5 h0 r1 t" y4 y& x; Z
A=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));

# y, s/ k9 P7 n( j* p' R- x! |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));

$ P: V  W9 ~' iD=(A+real(O))./2;
E=(A-real(O))./2;
- x. c& h$ m" H# o( bF=sqrt(D.^2+E.^2);
- V1 v4 z2 [% N) cZ(e,f)=E
8 W: ]- v# m/ L# h
1 S, `' [$ j2 P( f5 H; ~5 S0 tf=f+1;
; s. T6 W2 ^" ?6 J1 L  end
: r2 F/ h$ E! }$ ]8 Z8 N9 j* d  e=e+1;
  f=1;
end
[d]=[-55:2:55];
8 V/ [! y( a! V6 E; O- f1 c3 w[c]=[-55:2:55];
. m7 V# T9 t9 [5 X3 i' e%idy =d.^2./3600+c.^2./400< 1;
%d(idy) = nan;
%c(idy) = nan;

! _' k* C( s7 d% k  d( Vfigure(3)
1 Z/ ?0 A8 r2 t" P$ A! W9 t" nsuRF(d,c,Z)
, p& Y4 [. G' E& Q$ v8 ?% v9 H%surf(d,c,Z)
+ V$ O/ i: ^! ?%view([0,1,0]);
/ p, g& Y. d  U- R! z1 G  Qshading interp;  
/ A7 W( p8 B9 d* g* R' ?1 c9 Scolorbar;  
%axis equal;
这个函数运行后可以产生Z的矩阵值，但是surf不显示图像，但我将xy的范围定义在-45到45的时候会显示图像，请问是什么问题

Liberallh21

没看懂