function [gg,emspp]=mgm(t,x,n1)
cc=size(x);%定義數據的個數.
x0=x;%產生原始數據
nn=cc(1);n=cc(2);
%[a b]'=(A'*A)^(-1)*A'*X
%準備相關數據
B=[];y=[];x1=[];A=[];X=[];
for j=1:
x1(1,1)=x0(1,1);%定義累加第一個數據
%----------計算矩陣A,x1----------------------------------
for i=1:n-1
    x1(1,i+1)=x1(1,i)+x0(1,i+1);   
    B(i,1)=-(x1(1,i+1)+x1(1,i))/2;
    X(i,1)=x0(1,i+1);
end
A=[B,ones(n-1,1)];
auz=(A'*A)^(-1)*A'*X;
%------ 預測值計算---------
g0=[]; ems=[];
g0(1,1)=x(1,1);ems(1,1)=0;a=auz(1);b=auz(2);
for i=2:n1
  g0(1,i)=(x(1,1)-b/a)*(1-exp(a))*exp(-a*i);
end
%for i=1:n
 %   ems(1,i)=abs(g0(1,i)-x(1,i));
%end
%gg=g0';emspp=ems';