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

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] ;
: b+ b) d- W! S9 f: n1 Ot=30+1;   %城市的数目是30个
) U( g. V8 G- u- a" s& Ys=10;   %样本的数目是1400个

$ p, j- R( {' o/ \8 ^$ lx=cities(1,;
2 ~1 m( M5 K# @" r8 Xy=cities(2,;
. t6 k1 L& D. l3 S0 h! g
5 |: x$ c) J  G: i# P5 opop=zeros(s,t);  %得初始的pop矩阵，矩阵的最后一列表示所在行的样本的路径距离
8 `# v) c7 B( ^for i=1:s
   pop(i,1:t-1)=randperm(t-1);   %随机产生1—（t-1）的t-1个打乱的随机数
end
! z+ m' Q2 a; c1 f  q/ tM=100;
; s+ w0 `9 m% d; ~2 p  Zfor l=1:1:2for i=1:2:s-1
%随机选择一个交叉点
0 k6 c8 d8 {8 l2 ]' A, O% j& U   m=randperm(t-3)+1;
   crosspoint=min(m(1),m(2));
, v9 G  C, ?) }# a%任意两行交叉   
( h; E2 ^; Q$ t& f) M" ]   x1=pop(i,;
   x2=pop(i+1,;
6 p: P  B) u+ @1 H5 U) {+ c5 a   k=0;
  for K=1:10
8 i5 P0 l- E0 G) \; {1 T  W     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
            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)];
; A# {  @- @5 g            %检查x1左边的重复性并得到x1的左边
          for j=1:crosspoint
           while find(x1(crosspoint+1:t)==x1(j))
              zhi=find(x1(crosspoint+1:t)==x1(j)); %确定重复位置
+ s8 T" G" ]* j( [4 J              temp=x2(crosspoint+zhi);
              x1(j)=temp;
2 H! R4 ]  V/ w+ |4 n* R" [           end
0 D0 _* @. T/ t# e; K4 U          end
         for j=1:crosspoint
" U, }0 q2 e5 S% B           while find(x2(crosspoint+1:t)==x2(j))
! f+ `" v( D0 y7 v              zhi=find(x2(crosspoint+1:t)==x2(j)); %确定重复位置
, V& A& A" X8 R$ g( l, H              temp=x1(crosspoint+zhi);
% K; |- a9 n5 E! F& T" a              x2(j)=temp;
3 t8 |0 R- c9 H7 U  Y: |' f! O! r           end
+ z+ o  K; Y3 z4 h: I8 O          end
8 _/ |* s' S  L- X( I0 Z* d            k=1;
       else
+ j2 x% X, I- ?            m=randperm(t-3)+1;
& |5 l9 A; I$ J, t5 L5 v% d; \            crosspoint=min(m(1),m(2));
            k=0;
# F9 K% @7 N) N* K       end
     end
- P" `: M  Z( n  y+ w. b   end
end
, p) k$ U6 B0 S9 s* s3 G. kendbestL=min(pop(:,t))
J=pop(:,t);
fi=1./J;
0 F! r8 F7 D* J% j7 w4 ^
[Oderfi,Indexfi]=sort(fi);   %对于fi进行排序
) ?4 D4 F6 E1 r+ H) a: dBestS=pop(Indexfi(s),;     %得到最短路
, `8 J  j$ T& g+ N% i3 H$ r# _
- G7 G0 p& d' o. b. O& V  r' wI=BestS;
2 q" e  d5 t, ?  q* Y& E: P3 t8 L
+ [4 U, i* ~5 w& P9 T* Yfor i=1:1:t-1
    x1(i)=x(I(i));
    y1(i)=y(I(i));
end
) a4 {. ^7 E( G' kx1(t)=x(I(1));
y1(t)=y(I(1));

cities_new=[x1;y1];
disp('Best Route is:');disp(cities_new);
8 X9 [: W/ h0 K$ A- ]$ x  }, opos=[cities_new cities_new(:,1)];
' \& ~* p3 l: _0 h/ S2 o
) T: q4 l) x9 v' Z7 [, Elentemp=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);
( x2 n, V' C+ }- d  Q% {    lentemp=lentemp+temp;
/ T, b( f) a) i; O4 Fend
disp('Shortest Length is:');disp(lentemp);

figure(1);
subplot(1,2,1);       %窗口分割的左边部分
  x1 N: N( q) [( ~) v, tx(t)=x(1);y(t)=y(1);
; i6 W% T0 Y) ?. S* Gplot(x,y,'-or');
xlabel('X axis'), ylabel('Y axis'), title('原始路径');
/ r- ~; |6 |: saxis([0,1,0,1]);
& l- @$ y1 j5 b+ b& _/ S/ Z& }axis([0,100,0,100]);
axis on
4 M0 J. a0 Z* }7 Ihold on;
5 l: s) O( C* S' K0 D8 [  vsubplot(1,2,2);        %窗口分割的右边部分
1 V2 z4 h) S/ l7 o" a7 T$ l" _plot(x1,y1,'-or');
xlabel('X axis'), ylabel('Y axis'), title('最新的路径');
axis([0,1,0,1]);
axis([0,100,0,100]);
axis on

) q0 d1 Z4 h' W3 P$ ^这只是遗传算法的交叉部分，现在有点问题，就是为什么去掉红色字体就可以运行出来，而加上后程序就没反应了呢？
$ D1 u2 d( Z7 V* f那位高人能找出原因，麻烦看看啊 ，不胜感激啊！