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

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] ;
! X$ y  I! N; R; X( F$ g' |t=30+1;   %城市的数目是30个
s=10;   %样本的数目是1400个

! D0 r6 a8 X4 ~- k# `x=cities(1,;
% V/ [2 i% [5 m; C& My=cities(2,;
' y) S% Z0 ]  I: X/ ~
pop=zeros(s,t);  %得初始的pop矩阵，矩阵的最后一列表示所在行的样本的路径距离
" S' d( j9 T1 yfor i=1:s
, k% I6 x& y# l, h   pop(i,1:t-1)=randperm(t-1);   %随机产生1—（t-1）的t-1个打乱的随机数
end
$ D' M5 L: i0 ~M=100;
for l=1:1:2for i=1:2:s-1
%随机选择一个交叉点
" i; q$ a4 l; t   m=randperm(t-3)+1;
   crosspoint=min(m(1),m(2));
" G& ~+ g4 c, }0 F) U. K%任意两行交叉   
& R1 p' A$ K: l, U   x1=pop(i,;
   x2=pop(i+1,;
: e, d: ~0 h- Q9 h   k=0;
6 h$ [% j0 ~% @$ g+ V  for K=1:10
* g+ z$ Y. O, H% g% U7 t     while k==0
% @7 a; z8 E  v+ x6 y1 M: O% \       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)];
            %检查x1左边的重复性并得到x1的左边
6 R# K# w+ b8 R          for j=1:crosspoint
           while find(x1(crosspoint+1:t)==x1(j))
              zhi=find(x1(crosspoint+1:t)==x1(j)); %确定重复位置
              temp=x2(crosspoint+zhi);
              x1(j)=temp;
           end
) U7 `, ^- C+ |! `          end
- S) J& c2 _+ M( z! h7 d9 T  h0 ?         for j=1:crosspoint
! A9 b% l( j6 b7 N0 ?5 ^# P9 L* |           while find(x2(crosspoint+1:t)==x2(j))
) d" N% J2 d, Q/ z) d# u" T: E              zhi=find(x2(crosspoint+1:t)==x2(j)); %确定重复位置
/ W( |" R% k2 ^2 ^) X/ B9 m3 h              temp=x1(crosspoint+zhi);
              x2(j)=temp;
           end
          end
, q. i3 O/ o1 ]- F6 s$ ~            k=1;
       else
  A4 _. G9 m, }' f            m=randperm(t-3)+1;
. \  Z7 f9 Q8 _2 P8 E  k. d            crosspoint=min(m(1),m(2));
            k=0;
$ y  L4 L) S& f$ q" \; u6 D* ?/ T       end
$ Y& V  L( k3 h# w+ y5 u6 {' G     end
% `- l9 Y4 N  r: n& U1 J   end
, a+ c; j. N' {1 ]9 Xend
endbestL=min(pop(:,t))
J=pop(:,t);
fi=1./J;
, Q+ g( z' I# x. f
[Oderfi,Indexfi]=sort(fi);   %对于fi进行排序
9 }5 p* [: K! q: uBestS=pop(Indexfi(s),;     %得到最短路

3 f' `9 Q( g9 H5 `I=BestS;

. a4 f4 G. |5 j8 {" H7 d8 ^, T, Pfor i=1:1:t-1
    x1(i)=x(I(i));
4 \, q, D% Q, ^& b) O. u. v    y1(i)=y(I(i));
end
x1(t)=x(I(1));
- D. O% [: N: s0 G- K+ b! N5 ky1(t)=y(I(1));

cities_new=[x1;y1];
disp('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;
7 Z8 _" \5 I: b' P3 w! A! h! Mend
disp('Shortest Length is:');disp(lentemp);
- b" l7 l4 N  w/ F& ?6 N! P0 H
figure(1);
8 ~1 i( ^1 s9 dsubplot(1,2,1);       %窗口分割的左边部分
0 _! d* l' R4 w) f8 S. {; l( rx(t)=x(1);y(t)=y(1);
plot(x,y,'-or');
xlabel('X axis'), ylabel('Y axis'), title('原始路径');
( L0 @* z1 B4 @% E+ D( U0 R* Yaxis([0,1,0,1]);
axis([0,100,0,100]);
' k4 A; R2 g! p: ?1 iaxis on
hold on;
1 S8 q9 A  x' @7 s/ U6 n( \, Fsubplot(1,2,2);        %窗口分割的右边部分
6 ]; Q4 a2 a7 s5 K5 Gplot(x1,y1,'-or');
% v# G' G! v, z+ Zxlabel('X axis'), ylabel('Y axis'), title('最新的路径');
axis([0,1,0,1]);
3 e* S3 ]7 ]+ t$ n- Eaxis([0,100,0,100]);
- ]. v( [/ J7 q9 }! l* zaxis on

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