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

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] ;
0 E/ A8 a" T0 D+ P7 Wt=30+1;   %城市的数目是30个
6 \, c  x+ N  y, [3 hs=10;   %样本的数目是1400个

x=cities(1,;
y=cities(2,;

pop=zeros(s,t);  %得初始的pop矩阵，矩阵的最后一列表示所在行的样本的路径距离
for i=1:s
   pop(i,1:t-1)=randperm(t-1);   %随机产生1—（t-1）的t-1个打乱的随机数
end
0 Y. e/ N! G" x2 L8 CM=100;
; ?1 l7 P  B8 h* S" b& sfor l=1:1:2for i=1:2:s-1
%随机选择一个交叉点
2 L) e, E6 @; F. [9 Q& Y   m=randperm(t-3)+1;
   crosspoint=min(m(1),m(2));
" p2 O9 K/ n6 K% B%任意两行交叉   
4 E6 a3 x# F, _" b1 W) u   x1=pop(i,;
$ w8 F& d* w3 q. w4 b   x2=pop(i+1,;
   k=0;
  for K=1:10
4 l  D  i5 ~+ E3 Y  F" z) n! N) }     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
' p4 x  V7 l* d0 `% m5 [            pop(i,=[pop(i,1:crosspoint),pop(i+1,crosspoint+1:t)];
, L5 Y; E6 l# L- o& x            pop(i+1,=[pop(i+1,1:crosspoint),pop(i,crosspoint+1:t)];
            %检查x1左边的重复性并得到x1的左边
          for j=1:crosspoint
           while find(x1(crosspoint+1:t)==x1(j))
) g7 @; J% V1 D2 h6 A% {              zhi=find(x1(crosspoint+1:t)==x1(j)); %确定重复位置
              temp=x2(crosspoint+zhi);
              x1(j)=temp;
: g2 _0 j" d2 Y' G           end
" I0 E- k1 y- [( M( x          end
         for j=1:crosspoint
           while find(x2(crosspoint+1:t)==x2(j))
              zhi=find(x2(crosspoint+1:t)==x2(j)); %确定重复位置
5 l3 [- }7 `9 ^- Z              temp=x1(crosspoint+zhi);
              x2(j)=temp;
1 B' ]$ \1 D; w           end
) _% Z/ ~2 R- s8 L0 i/ ~          end
            k=1;
       else
& G) K7 G3 y8 O/ [            m=randperm(t-3)+1;
. r* @% N) Z8 O            crosspoint=min(m(1),m(2));
            k=0;
$ l- u: t" r6 Z% {       end
     end
2 P0 ]) p+ h- q  r9 H$ }* g8 \" a- o; G   end
' P9 q$ l4 C1 o2 q  L  S( C; i- g4 uend
# e; E5 b$ `; f1 f& O, V. j+ q7 Q. OendbestL=min(pop(:,t))
J=pop(:,t);
9 F. i0 ~) a) w9 q4 E# C* Ufi=1./J;
  {$ b$ I1 i- A7 p( `+ q$ B2 O5 I- w
[Oderfi,Indexfi]=sort(fi);   %对于fi进行排序
) P) i2 u5 t4 o. yBestS=pop(Indexfi(s),;     %得到最短路

I=BestS;
# i# G# p  Y* ^& |/ ^0 N
8 G# i  ?0 m# W& d+ U) I5 Sfor i=1:1:t-1
    x1(i)=x(I(i));
    y1(i)=y(I(i));
7 T0 C/ L) E! K  ^end
  _  y" X; {1 @# q& j: b! H$ L6 Y0 jx1(t)=x(I(1));
- d) u2 _. u0 K. ?3 |3 `8 e5 [" xy1(t)=y(I(1));

% \8 H$ k9 H" R% @" acities_new=[x1;y1];
- R9 [( c2 B6 L' s7 Ldisp('Best Route is:');disp(cities_new);
pos=[cities_new cities_new(:,1)];

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);
    lentemp=lentemp+temp;
end
disp('Shortest Length is:');disp(lentemp);
$ f& F" S' t+ q/ A. S) U# V
figure(1);
subplot(1,2,1);       %窗口分割的左边部分
9 P' w/ W/ B- ~4 t( Z' Gx(t)=x(1);y(t)=y(1);
0 e8 J4 g5 Q, F" s: L2 L, e9 f" cplot(x,y,'-or');
xlabel('X axis'), ylabel('Y axis'), title('原始路径');
axis([0,1,0,1]);
axis([0,100,0,100]);
axis on
; n  p3 H5 `$ Y( I1 [5 O$ M. Ohold on;
. q0 a, V) l4 ~+ `% \0 M) g" N7 gsubplot(1,2,2);        %窗口分割的右边部分
plot(x1,y1,'-or');
xlabel('X axis'), ylabel('Y axis'), title('最新的路径');
axis([0,1,0,1]);
axis([0,100,0,100]);
axis on

这只是遗传算法的交叉部分，现在有点问题，就是为什么去掉红色字体就可以运行出来，而加上后程序就没反应了呢？
那位高人能找出原因，麻烦看看啊 ，不胜感激啊！