用matlab做灰色预测，第一次运行程序是好的，可是再运行就不对了，求大神帮忙，急急急！！！

%程序中的变量定义:alpha是包含α、μ值的矩阵;ago是预测后累加值矩阵;var是预测值矩阵;error是残差矩阵; c是后验差比值
 function gm(x);          %定义函数gm(x)
 clc;                     %清屏,以使计算结果独立显示
 format long;             %设置计算精度
 if length(x(:,1))==1     %对输入矩阵进行判断,如不是一维列矩阵,进行转置变换
     x=x';
 end;
 n=length(x);             %取输入数据的样本量
 z=0;
 for i=1:n                %计算累加值,并将值赋与矩阵be
     z=z+x(i,:);
     be(i,:)=z;
 end
 for i=2:n                %对原始数列平行移位
     y(i-1,:)=x(i,:);
 end
 for i=1:n-1              %计算数据矩阵B的第一列数据
     c(i,:)=-0.5*(be(i,:)+be(i+1,:));
 end
 for j=1:n-1              %计算数据矩阵B的第二列数据
     e(j,:)=1;
 end
 for i=1:n-1              %构造数据矩阵B
     B(i,1)=c(i,:);
     B(i,2)=e(i,:);
 end
 alpha=inv(B.'*B)*B.'*y;  %计算参数α、μ矩阵
 for i=1:n+3              %计算数据估计值的累加数列,如改n+1为n+m可预测后m-1个值
     ago(i,:)=(x(1,:)-alpha(2,:)/alpha(1,:))*exp(-alpha(1,:)*(i-1))+alpha(2,:)/alpha(1,:);
 end
 var(1,:)=ago(1,:)
  for i=1:n+2  %如改n为n+m-1,可预测后m-1个值
      var(i+1,:)=ago(i+1,:)-ago(i,:);    %估计值的累加数列的还原,并计算出下一预测值
              end
          for i=1:n
       error(i,:)=var(i,:)-x(i,:);     %计算残差
   end
     c=std(error)/std(x);      %调用统计工具箱的标准差函数计算后验差的比值c
     ago             %显示输出预测值的累加数列
     alpha           %显示输出参数α、μ数列
     var             %显示输出预测值
     error           %显示输出误差
     c               %显示后验差的比值c
x=[19519,19578,19637,19695,16602,25723,30379,34473,38485,40514,42400,48337];gm(x)
一直说第五行有问题，求大神帮忙！！！跪求！
	           

答案：var =
       19519
 
ago =
  1.0e+005 *
   0.19519000000000
   0.36476693858308
   0.55334597565305
   0.76305640868351
   0.99626613551813
   1.25560839095269
   1.54401147930921
   1.86473183872077
   2.22139081046496
   2.61801552851843
   3.05908439102787
   3.54957762712807
   4.09503353007087
   4.70161099160872
   5.37615904372558
 
alpha =
  1.0e+004 *
  -0.00001062104750
   1.39999674117304
 
var =
  1.0e+004 *
   1.95190000000000
   1.69576938583078
   1.88579037069969
   2.09710433030461
   2.33209726834619
   2.59342255434559
   2.88403088356519
   3.20720359411566
   3.56658971744186
   3.96624718053473
   4.41068862509443
   4.90493236100196
   5.45455902942797
   6.06577461537853
   6.74548052116857
 
error =
  1.0e+003 *
                  0
  -2.62030614169218
  -0.77909629300307
   1.27604330304606
   6.71897268346191
   0.21122554345592
  -1.53869116434810
  -2.40096405884344
  -2.81910282558145
  -0.85152819465270
   1.70688625094428
   0.71232361001964
 
c =
   0.23750894298831
 
 
 
运行了两次没有出现问题？

上一个：谁有matlab的入门的简单程序？还有数学建模常用的模型
下一个：我用matlab编的程序，另存了一个文件，结果一个能运算出来，一个就不行。都是放在同一个work文件里面的。