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

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个
5 r1 v: o/ ^5 |
% L" v; ?3 Y* D/ ?# t9 Xx=cities(1,;
5 P0 y7 a8 E% W+ X- B* Qy=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
( s/ \" J0 x' {, w" S$ mM=100;
) u& G+ ]/ q( `+ \: W# G. ]for l=1:1:2for i=1:2:s-1
%随机选择一个交叉点
   m=randperm(t-3)+1;
$ t( O5 I, [& v* H8 @2 }   crosspoint=min(m(1),m(2));
$ Z9 u  g  \- e/ N: P%任意两行交叉   
   x1=pop(i,;
   x2=pop(i+1,;
   k=0;
  for K=1:10
7 i% v% W! t% Z     while k==0
( a6 {! M0 C- J* A: q) X       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
" I! F$ A( X# g4 ]9 l- b: }            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的左边
1 T* e6 g8 K, f+ `          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
          end
         for j=1:crosspoint
           while find(x2(crosspoint+1:t)==x2(j))
              zhi=find(x2(crosspoint+1:t)==x2(j)); %确定重复位置
              temp=x1(crosspoint+zhi);
, {' `7 c- d7 w" ]( m              x2(j)=temp;
; ^: N6 K5 b) p           end
% `" k* f$ Y' o8 t          end
            k=1;
% |# {8 i4 F5 a; Q# Q* [$ u       else
            m=randperm(t-3)+1;
, {8 C6 m, v* ]3 }6 h9 b            crosspoint=min(m(1),m(2));
            k=0;
       end
     end
   end
end
endbestL=min(pop(:,t))
+ |. F# C3 l* o8 h6 ~J=pop(:,t);
fi=1./J;

[Oderfi,Indexfi]=sort(fi);   %对于fi进行排序
# U6 n: ^+ \3 b# D; d, t8 XBestS=pop(Indexfi(s),;     %得到最短路
1 v! v$ J/ p4 s: u( m8 E1 a
& P' y1 o$ U$ V/ B# qI=BestS;
( {% t) J' y: P
for i=1:1:t-1
: D) h( l  b7 G+ E! w! Q    x1(i)=x(I(i));
    y1(i)=y(I(i));
0 h  Y. |$ Q; G( Z: U8 P- Uend
9 p. F1 o& Y5 u: N- q! g) jx1(t)=x(I(1));
y1(t)=y(I(1));
% q* F2 k9 z3 w" `, K& G
7 J( \: h/ j8 \4 m! y! Vcities_new=[x1;y1];
disp('Best Route is:');disp(cities_new);
pos=[cities_new cities_new(:,1)];

lentemp=0;
for i=1:1:t-1
" h% C4 h, ]8 T' \    temp=sqrt((pos(1,i)-pos(1,i+1))^2+(pos(2,i)-pos(2,i+1))^2);
    lentemp=lentemp+temp;
end
, j8 d) {! s7 E2 v9 P2 pdisp('Shortest Length is:');disp(lentemp);

figure(1);
subplot(1,2,1);       %窗口分割的左边部分
x(t)=x(1);y(t)=y(1);
& {1 A% [$ A$ y! k7 s* F$ U3 {plot(x,y,'-or');
xlabel('X axis'), ylabel('Y axis'), title('原始路径');
axis([0,1,0,1]);
3 k" L0 t; e5 a6 r- k: b* ]4 D6 Jaxis([0,100,0,100]);
8 T2 B' L. m% L. O) }# Xaxis on
4 z0 p7 m) i  X( w) ehold on;
subplot(1,2,2);        %窗口分割的右边部分
; {3 ^+ s% e" h& p$ e! F3 n; kplot(x1,y1,'-or');
# e6 b9 o$ V5 t5 }& }) ~xlabel('X axis'), ylabel('Y axis'), title('最新的路径');
3 z/ f. J' a/ ^0 V$ \2 e5 Zaxis([0,1,0,1]);
axis([0,100,0,100]);
axis on

: j; N7 O' K4 }7 g4 G5 L这只是遗传算法的交叉部分，现在有点问题，就是为什么去掉红色字体就可以运行出来，而加上后程序就没反应了呢？
) R2 ]% G0 L4 N, d2 t那位高人能找出原因，麻烦看看啊 ，不胜感激啊！