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

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个
s=10;   %样本的数目是1400个
+ ~8 s. W5 V% ^) K" M
! R1 X0 O# W8 f$ b& fx=cities(1,;
y=cities(2,;
7 R" F" b+ o! P9 q% k3 L
" ^7 J/ {* F0 |# r* _pop=zeros(s,t);  %得初始的pop矩阵，矩阵的最后一列表示所在行的样本的路径距离
+ [  A8 O$ x+ }6 G* Yfor i=1:s
: z! z+ c- m* R3 Z   pop(i,1:t-1)=randperm(t-1);   %随机产生1—（t-1）的t-1个打乱的随机数
end
M=100;
for l=1:1:2for i=1:2:s-1
# o1 V8 i- j) y%随机选择一个交叉点
   m=randperm(t-3)+1;
   crosspoint=min(m(1),m(2));
# d! r  x) g7 r2 O7 e8 {7 x6 L$ G%任意两行交叉   
   x1=pop(i,;
% |+ _8 T( I5 c  n   x2=pop(i+1,;
6 s6 y5 D3 X) |   k=0;
  for K=1:10
     while k==0
2 Z2 P8 U$ v/ P/ M/ z% u: k1 A       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
7 G1 U  J& c6 b# o$ J            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)];
, y3 D* a/ g5 z9 V            %检查x1左边的重复性并得到x1的左边
% ?; r3 h# ^0 z! K1 p          for j=1:crosspoint
           while find(x1(crosspoint+1:t)==x1(j))
              zhi=find(x1(crosspoint+1:t)==x1(j)); %确定重复位置
              temp=x2(crosspoint+zhi);
7 r/ R# [* M! k3 W1 s8 p% E              x1(j)=temp;
           end
          end
         for j=1:crosspoint
           while find(x2(crosspoint+1:t)==x2(j))
  ^( B  m3 p' Q              zhi=find(x2(crosspoint+1:t)==x2(j)); %确定重复位置
              temp=x1(crosspoint+zhi);
              x2(j)=temp;
: t1 h6 V+ p& C: ?: \% A           end
9 G! M1 f: s* E* _: P- `% q          end
            k=1;
       else
" f+ M, [6 y5 [  Q: t            m=randperm(t-3)+1;
            crosspoint=min(m(1),m(2));
) f" h. n( a6 a. S1 E8 n4 l0 G1 f            k=0;
$ d0 Q2 q  ?4 X) N. b8 W$ U8 I9 T* N       end
2 [) E/ \, q  o6 f) C, r; F     end
   end
- W  p) P4 A  d% ]& y# B" a1 u4 B: ~* send
) R7 T( A1 Y. nendbestL=min(pop(:,t))
! r' _. z5 K( _! H' b9 d' x' ~/ ]J=pop(:,t);
fi=1./J;
# l3 l( _  s% J9 S( y( m, G3 e6 {6 p
  I: p1 r; S( h8 C: S6 ^+ \[Oderfi,Indexfi]=sort(fi);   %对于fi进行排序
9 p9 j! F7 ?9 @; ?! p# `3 mBestS=pop(Indexfi(s),;     %得到最短路

I=BestS;
1 Z* @' R& Q$ \' J/ y6 H0 C0 C
1 W( K- O* N. P* V, r/ C: wfor i=1:1:t-1
' p6 p9 u) H  C0 D* E+ c5 U- T# m    x1(i)=x(I(i));
0 y" }$ u: y* b6 \# ]6 R    y1(i)=y(I(i));
end
. ^; _' o: }7 Z+ k2 b, J6 U$ X& I0 Qx1(t)=x(I(1));
8 v4 }/ X7 l. l( qy1(t)=y(I(1));

cities_new=[x1;y1];
disp('Best Route is:');disp(cities_new);
pos=[cities_new cities_new(:,1)];
2 `3 B( v3 r2 h+ E) H# A
lentemp=0;
8 U0 i& P! v- A, ^5 W# r+ t. @) bfor i=1:1:t-1
    temp=sqrt((pos(1,i)-pos(1,i+1))^2+(pos(2,i)-pos(2,i+1))^2);
3 p  }" u" p3 _: i) O    lentemp=lentemp+temp;
1 \5 l9 o+ u. N& N. Tend
3 L0 \4 [1 i# M) {; \+ _disp('Shortest Length is:');disp(lentemp);
8 a" g) d. P' h: @7 D' g, N
% _3 y$ z& f; E9 ]" Ufigure(1);
subplot(1,2,1);       %窗口分割的左边部分
x(t)=x(1);y(t)=y(1);
% X/ X( M: J7 |7 Jplot(x,y,'-or');
: ^5 ?. K8 T" X' [9 s( n& u/ \. \xlabel('X axis'), ylabel('Y axis'), title('原始路径');
& t+ W9 A/ ~/ M+ m! Daxis([0,1,0,1]);
axis([0,100,0,100]);
$ s* c4 Z9 q' z& [axis on
hold on;
4 f8 V, I2 v" I* P7 G9 x) `& psubplot(1,2,2);        %窗口分割的右边部分
5 T! L5 v2 t9 |( s: Oplot(x1,y1,'-or');
) v5 B" S0 I& txlabel('X axis'), ylabel('Y axis'), title('最新的路径');
$ f4 L9 h# l& a0 saxis([0,1,0,1]);
axis([0,100,0,100]);
2 Q' C8 O( v6 }& M$ v3 s# Y$ s, Naxis on
! x9 q" n% l0 e6 |8 Y8 C. X
这只是遗传算法的交叉部分，现在有点问题，就是为什么去掉红色字体就可以运行出来，而加上后程序就没反应了呢？
' W1 [: }( o9 V# L+ y! {* }那位高人能找出原因，麻烦看看啊 ，不胜感激啊！