function?gm（x0）;
clc
format?long;
if?length（x0（:1））==1????%對輸入矩陣進(jìn)行判斷，如不是一維列矩陣，進(jìn)行轉(zhuǎn)置變換
????x0=x0‘;
end
n=length（x0）;
%由數(shù)列x0生成數(shù)列x1;
s=0;
for?i=1:n
????s=s+x0（i）;
????x1（i）=s;
end
%由數(shù)列x1生成矩陣C和A；?
for?i=1:（n-1）;?
A（i）=-（x1（i）+x1（i+1））/2;?
end?
A=[A‘?ones（n-11）];?
for?i=1:（n-1）;?
C（i）=x0（i+1）;?
end?
C=C‘;
beta=inv（A‘*A）*A‘*C;
%計(jì)算出a和u;
a=beta（1）;