GM(1,1)
X0=[29720041.01 28313169.6 29042078.1 29142684.71 29270689.83
32690539.37 33746240.65 30686234.79 30467186 30788868.15 32833784.44
34579255.77 35473440.19 37056010.27 38038330.97 38735790.11 38352027.74
37464074.42 40570777.1 42601876.43 43210774.71
];
n=length(X0);
X1=[];
X1(1)=X0(1);
for i=2:n
X1(i)=X1(i-1)+X0(i);
end
for i=1:n-1
B(i,1)=-0.5*(X1(i)+X1(i+1));
B(i,2)=1;
Y(i)=X0(i+1);
end
alpha=(B'*B)^(-1)*B'*Y';
a=alpha(1,1);
b=alpha(2,1);
d=b/a;
c=X1(1)-d;
X2(1)=X0(1);
X(1)=X0(1);
for i=1:n-1;
X2(i+1)=c*exp(-a*i)+d;
X(i+1)=X2(i+1)-X2(i);
