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

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] ;
- K3 i( z9 |, u7 L. |% }t=30+1;   %城市的数目是30个
) _! F- o; s* A2 i# x2 X1 A3 d9 C: Zs=10;   %样本的数目是1400个

: E& ^4 k1 q: `9 f2 r8 ux=cities(1,;
) @. L& `" I7 \0 M! Z& A# zy=cities(2,;

$ i5 O' m7 e2 [( [0 c+ ]3 Epop=zeros(s,t);  %得初始的pop矩阵，矩阵的最后一列表示所在行的样本的路径距离
8 S, W$ ]! j8 s  t/ Z$ Afor i=1:s
   pop(i,1:t-1)=randperm(t-1);   %随机产生1—（t-1）的t-1个打乱的随机数
' i7 G! N' B. ~- Cend
  V% M" A1 ^9 m: S/ j5 i3 gM=100;
for l=1:1:2for i=1:2:s-1
2 T: K  Y0 k. L$ B/ u( j%随机选择一个交叉点
   m=randperm(t-3)+1;
* @3 S' [- ]: m7 d   crosspoint=min(m(1),m(2));
%任意两行交叉   
   x1=pop(i,;
   x2=pop(i+1,;
) S( R3 O- \8 f. l' s  Q8 \6 T   k=0;
" m, m" k/ n# ?; }7 d  for K=1:10
8 G0 g: [8 C  L     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)];
, S7 Z5 l% T" x  j            %检查x1左边的重复性并得到x1的左边
          for j=1:crosspoint
) L; N+ F* O/ F# U           while find(x1(crosspoint+1:t)==x1(j))
6 Q' H& k; s: ~8 c& r" m              zhi=find(x1(crosspoint+1:t)==x1(j)); %确定重复位置
              temp=x2(crosspoint+zhi);
' q$ L2 C2 a9 Y7 y% z2 [1 R              x1(j)=temp;
$ v; q0 S8 m/ e           end
" ?, `- W& S! V          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);
/ Y) F. v* `* t2 d9 \              x2(j)=temp;
           end
          end
6 f1 R/ t7 }* p9 h" d" W2 a, Q            k=1;
       else
            m=randperm(t-3)+1;
            crosspoint=min(m(1),m(2));
            k=0;
9 v3 A8 D' ?1 S5 f$ \( G; j       end
1 J2 g- W8 q( _  |7 s     end
   end
end
1 E& F4 p' X7 N5 f: ~endbestL=min(pop(:,t))
! \* p- ?" o% e' E) V  |J=pop(:,t);
fi=1./J;

[Oderfi,Indexfi]=sort(fi);   %对于fi进行排序
7 Q* b3 I; I, v9 D# @0 b1 [1 z% ]. _BestS=pop(Indexfi(s),;     %得到最短路
3 D9 W( A9 O8 z! d7 V4 f7 `
I=BestS;

for i=1:1:t-1
9 F9 ]' k8 A' H' t) i0 {+ a, ^( l. z    x1(i)=x(I(i));
! S* Y) E2 U4 u. e9 ?    y1(i)=y(I(i));
end
' r/ w0 ]6 Y( N7 O" Mx1(t)=x(I(1));
1 n* Z2 m% M0 ]7 W2 E" ?y1(t)=y(I(1));

" c+ |. [7 p/ `* y+ ^0 t7 P$ p4 [3 Ccities_new=[x1;y1];
disp('Best Route is:');disp(cities_new);
/ ^. `: g" z5 u6 ]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);
- C! i2 t* h) O$ R2 @% X0 b    lentemp=lentemp+temp;
: ^! A) X4 x! K4 E# X5 C* Z8 Q( yend
disp('Shortest Length is:');disp(lentemp);
: L1 F4 s2 r- V4 y
figure(1);
2 v1 Z# w: _' j5 |* Ssubplot(1,2,1);       %窗口分割的左边部分
2 h3 {) t* d0 X4 P: ^/ r1 Px(t)=x(1);y(t)=y(1);
plot(x,y,'-or');
xlabel('X axis'), ylabel('Y axis'), title('原始路径');
axis([0,1,0,1]);
axis([0,100,0,100]);
' _% X4 v% s: O' ~5 d! baxis on
/ k# `# D: n# ]hold on;
subplot(1,2,2);        %窗口分割的右边部分
plot(x1,y1,'-or');
xlabel('X axis'), ylabel('Y axis'), title('最新的路径');
3 K- \5 `9 V. I1 l) v  [& Maxis([0,1,0,1]);
" L- o( V# [9 M3 e: ?axis([0,100,0,100]);
7 x5 G% [7 S% Maxis on
6 O$ c1 {# b$ v; a( t
这只是遗传算法的交叉部分，现在有点问题，就是为什么去掉红色字体就可以运行出来，而加上后程序就没反应了呢？
6 W0 P( o# u) I* o3 W( \3 @那位高人能找出原因，麻烦看看啊 ，不胜感激啊！