求助 哪有问题啊！！！！！

zyli

cities= [18 87 74 71 25 58 4 13 18 24 71 64 68 83 58 54 51 37 41 2 7 22 25 62 87 91 83 41 45 44;54 76 78 71 38 35 50 40 40 40 42 44 60 58 69 69 62 67 84 94 99 64 60 62 32 38 46 26 21 35] ;
t=30+1;   %城市的数目是30个
* z/ ]5 T! T" U) Xs=10;   %样本的数目是1400个
/ o5 h8 W/ v7 u0 B
x=cities(1,;
y=cities(2,;

/ }# I5 h0 p9 G9 b2 I" D- K9 Vpop=zeros(s,t);  %得初始的pop矩阵，矩阵的最后一列表示所在行的样本的路径距离
+ \% m* b) U- b# }for i=1:s
9 C" ^- K) W% L/ a( `/ L8 k: a; W   pop(i,1:t-1)=randperm(t-1);   %随机产生1—（t-1）的t-1个打乱的随机数
end
6 L4 D9 [! S# m& v$ e) \M=100;
for l=1:1:2for i=1:2:s-1
: K. S; [9 O7 R; L: b%随机选择一个交叉点
0 @1 F7 }* h$ F' I' g  j7 p   m=randperm(t-3)+1;
, J5 l* S; ?* O7 e( t) c   crosspoint=min(m(1),m(2));
%任意两行交叉   
; ]9 A  a; f: \* p   x1=pop(i,;
   x2=pop(i+1,;
, m' p( D$ l$ q; e, |   k=0;
/ Z! U3 G% a9 D& h5 G  for K=1:10
     while k==0
       if sqrt((x(x1(crosspoint))-x(x2(crosspoint)))^2+(y(x1(crosspoint))-y(x2(crosspoint)))^2)<M|sqrt((x(x2(crosspoint))-x(x1(crosspoint)))^2+(y(x2(crosspoint))-y(x1(crosspoint)))^2)<M
8 q, e; E& N7 I" `& N' M! l9 f8 A            pop(i,=[pop(i,1:crosspoint),pop(i+1,crosspoint+1:t)];
            pop(i+1,=[pop(i+1,1:crosspoint),pop(i,crosspoint+1:t)];
            %检查x1左边的重复性并得到x1的左边
- g2 n" A$ P: u" l          for j=1:crosspoint
           while find(x1(crosspoint+1:t)==x1(j))
+ Z) `9 u- l' m/ p" A2 k1 X2 x              zhi=find(x1(crosspoint+1:t)==x1(j)); %确定重复位置
8 J. p1 t6 {& |: ]) P8 ~9 X% n              temp=x2(crosspoint+zhi);
, l) i8 I# h3 J, `              x1(j)=temp;
0 L( ]5 w+ D. q/ g3 b9 _           end
          end
" K* o# c* e$ a( c3 W6 P- r% A         for j=1:crosspoint
' W$ N$ N! u# R8 }  z3 H           while find(x2(crosspoint+1:t)==x2(j))
              zhi=find(x2(crosspoint+1:t)==x2(j)); %确定重复位置
. L: T( v8 P! j, q( k7 L; E0 s              temp=x1(crosspoint+zhi);
              x2(j)=temp;
* y/ d9 |: X* r( ?. B           end
          end
* H2 I, O/ |! N            k=1;
       else
            m=randperm(t-3)+1;
            crosspoint=min(m(1),m(2));
            k=0;
) V; G* @% O/ O& m0 [& O0 m       end
     end
   end
end
endbestL=min(pop(:,t))
3 B* l( W: M1 f7 y0 J2 kJ=pop(:,t);
, x- v8 @$ _0 \1 i& F  Pfi=1./J;
5 Y6 [" t/ {* |; \% l% y5 {! t* ]
! J0 }# W; m& a[Oderfi,Indexfi]=sort(fi);   %对于fi进行排序
BestS=pop(Indexfi(s),;     %得到最短路
6 H# p$ O* j" W
% \; A* r2 d; T  z, ZI=BestS;

for i=1:1:t-1
    x1(i)=x(I(i));
% b  U1 y3 ]% M0 Y- @# ^    y1(i)=y(I(i));
end
x1(t)=x(I(1));
y1(t)=y(I(1));
6 z/ y- @" `! \9 L' x4 i9 U9 t
1 Q. u0 a/ U9 J5 }cities_new=[x1;y1];
disp('Best Route is:');disp(cities_new);
pos=[cities_new cities_new(:,1)];
' W. a. |* u/ P: L  n
lentemp=0;
for i=1:1:t-1
    temp=sqrt((pos(1,i)-pos(1,i+1))^2+(pos(2,i)-pos(2,i+1))^2);
1 D$ z  K) U$ o9 ^- O4 R* S    lentemp=lentemp+temp;
: ~; j6 i9 ^6 L0 O6 @0 E- Rend
disp('Shortest Length is:');disp(lentemp);

figure(1);
- y4 L2 G+ A5 E- z' ~subplot(1,2,1);       %窗口分割的左边部分
x(t)=x(1);y(t)=y(1);
plot(x,y,'-or');
xlabel('X axis'), ylabel('Y axis'), title('原始路径');
- ~' e1 G0 }& v' Iaxis([0,1,0,1]);
axis([0,100,0,100]);
axis on
8 u/ ]5 Q; @+ b( Z) nhold on;
subplot(1,2,2);        %窗口分割的右边部分
, g2 @6 m# Y. c" D$ Hplot(x1,y1,'-or');
% S+ ^1 A. X9 H4 n, I# Dxlabel('X axis'), ylabel('Y axis'), title('最新的路径');
axis([0,1,0,1]);
axis([0,100,0,100]);
% a3 W0 H. c" z, Z& `axis on
# z) D  ^# o. Z! K+ d
这只是遗传算法的交叉部分，现在有点问题，就是为什么去掉红色字体就可以运行出来，而加上后程序就没反应了呢？
那位高人能找出原因，麻烦看看啊 ，不胜感激啊！