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

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] ;
$ v& |0 F  o; O+ X4 F; i# U3 I- Rt=30+1;   %城市的数目是30个
7 q& Q# i5 g3 \8 j; w" Us=10;   %样本的数目是1400个

x=cities(1,;
' f: S/ z* }' P# M9 my=cities(2,;

( t  T9 ], v! ~; y% D0 h$ epop=zeros(s,t);  %得初始的pop矩阵，矩阵的最后一列表示所在行的样本的路径距离
for i=1:s
   pop(i,1:t-1)=randperm(t-1);   %随机产生1—（t-1）的t-1个打乱的随机数
end
' L- ?& U8 v% m! M5 q8 K  zM=100;
1 S( L- T$ m) G3 I4 G- d/ D: K, e2 N! nfor l=1:1:2for i=1:2:s-1
8 Y3 Y! B. X5 R0 x" _: _: ^- l" D, i9 a%随机选择一个交叉点
   m=randperm(t-3)+1;
   crosspoint=min(m(1),m(2));
) W7 H. W! z0 ^9 A* m%任意两行交叉   
  d* c3 i& r: \$ o% {$ o" y   x1=pop(i,;
2 r, |+ w! O6 B; i- e   x2=pop(i+1,;
   k=0;
  for K=1:10
     while k==0
% k$ ?$ |* [" s       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的左边
          for j=1:crosspoint
  {; {3 L- b- }' ~* R, }3 a           while find(x1(crosspoint+1:t)==x1(j))
2 ?5 E% ?2 c, L! W+ r              zhi=find(x1(crosspoint+1:t)==x1(j)); %确定重复位置
1 Z6 F1 p7 b4 C              temp=x2(crosspoint+zhi);
  M* d" `7 F  h+ _1 s3 n+ K              x1(j)=temp;
           end
          end
) U  m8 Q! R$ G         for j=1:crosspoint
           while find(x2(crosspoint+1:t)==x2(j))
              zhi=find(x2(crosspoint+1:t)==x2(j)); %确定重复位置
  z, R  L/ h5 a" l9 v/ s              temp=x1(crosspoint+zhi);
              x2(j)=temp;
4 |5 W  i8 }8 d' }3 r/ n' k: F' X           end
          end
            k=1;
) ?7 z, C- Q, h7 @. D& E       else
            m=randperm(t-3)+1;
3 ~. a: }( [/ N% X6 T$ p. c            crosspoint=min(m(1),m(2));
% f" M+ d3 w) }' q            k=0;
       end
* A9 m( j/ r! N# C     end
2 ]6 v; L1 x8 f   end
end
4 z  Z# P  ]' @% O5 E# p* U6 v! {5 _endbestL=min(pop(:,t))
; h4 N  L4 _4 z- C, lJ=pop(:,t);
fi=1./J;
2 }: _1 r+ i9 l+ W' W
( u& e4 u# J/ ~3 Q3 p- }[Oderfi,Indexfi]=sort(fi);   %对于fi进行排序
BestS=pop(Indexfi(s),;     %得到最短路

I=BestS;

for i=1:1:t-1
2 W; y8 W# a; p, E    x1(i)=x(I(i));
    y1(i)=y(I(i));
end
. e  V: _+ d/ W9 Xx1(t)=x(I(1));
y1(t)=y(I(1));

cities_new=[x1;y1];
disp('Best Route is:');disp(cities_new);
pos=[cities_new cities_new(:,1)];

& Q, W  d' a1 [  v8 A; Q! ]; {/ Blentemp=0;
for i=1:1:t-1
$ [' H- D8 Y# n6 N% g$ m    temp=sqrt((pos(1,i)-pos(1,i+1))^2+(pos(2,i)-pos(2,i+1))^2);
( e) O: Y7 c& x7 B( X    lentemp=lentemp+temp;
end
disp('Shortest Length is:');disp(lentemp);

figure(1);
2 I9 Z( w7 h6 \/ T! q1 H+ Asubplot(1,2,1);       %窗口分割的左边部分
x(t)=x(1);y(t)=y(1);
plot(x,y,'-or');
/ X% Z" X+ ^( D( h/ Pxlabel('X axis'), ylabel('Y axis'), title('原始路径');
axis([0,1,0,1]);
axis([0,100,0,100]);
0 e8 m1 s& b4 ^# b9 T7 u( Iaxis on
3 \  e+ l/ y- \4 t2 T3 U* d" F* {hold on;
6 M5 T% _# I- _8 u, ^/ Vsubplot(1,2,2);        %窗口分割的右边部分
$ |* \; U4 Z6 N$ q, c% z. Yplot(x1,y1,'-or');
xlabel('X axis'), ylabel('Y axis'), title('最新的路径');
axis([0,1,0,1]);
axis([0,100,0,100]);
axis on
) ?$ R! H- S/ l/ u! c: z% e
这只是遗传算法的交叉部分，现在有点问题，就是为什么去掉红色字体就可以运行出来，而加上后程序就没反应了呢？
那位高人能找出原因，麻烦看看啊 ，不胜感激啊！