matlab仿真代码总是出问题。提示说for循环那行有问题，找不出来

three

U=1000;
  V. Z6 ]/ U5 T) U* Ba=1e-6;
b=1e-6;
5 q: s- V( J! R7 B' L7 P1 |D=0.25;
  M, N* {8 U+ ]' E: t; X; T% rlamda=1.225e-9/sqrt(U*(1+0.9783e-6*U));
4 s' \6 T9 a* b. p3 W& aA=pi*a/lamda*D;
" v% R, j" B' @( y8 OB=pi*(a+b)/lamda*D;
for(i=0;i<=200;i++)
4 U, _6 Q! h' }- c8 a) E{x=5e-5*(2*rand-1);
y=4e-5*(2*rand-1);
! }, d! Q0 U( x& rPw=(sin(A*x))^2*(cos(B*x))^2/(A*x)^2<=1;
; Y0 ^& g8 \4 D4 Y9 tM=rand;
7 q) ?. o& X3 e6 F0 E# mif(Pw>=M)
3 Y, s4 v+ o" V6 _5 rplot(x,y,'.r','markersize',10);
pause(200);}
$ {8 l4 E) I8 E: b' Z; v) L1 a# Z' ]

land

U=1000;
: A8 r  F) a# i" [a=1e-6;
b=1e-6;
D=0.25;
lamda=1.225e-9/sqrt(U*(1+0.9783e-6*U));
, U8 p6 ~- B& t8 x) i) g# f( t" fA=pi*a/lamda*D;
B=pi*(a+b)/lamda*D;
hold on
for i=0:200
% i' t$ M( o7 g! m/ Ax=5e-5*(2*rand-1);
y=4e-5*(2*rand-1);
Pw=(sin(A*x))^2*(cos(B*x))^2/(A*x)^2<=1;
M=rand;
" P% Q! }* t5 j/ e( Y! Q, Q3 gif Pw>=M
0 |' R/ T! t8 W, F3 Bplot(x,y,'.r','markersize',10)
5 W  K9 w; x( c- j5 p( dend
pause(.02)
" w3 S. _! M2 T; u6 \7 _7 I- lend
+ S5 y# e# x7 l6 z  U

rergr

再看看别人是怎么说的