function?[Xcerror1error2]=verr（X1k）
%?建立函數[Xcerror1error2]=example9_3_2_3（X1k）
%?其中X1為輸入序列，k為預測長度，
%?X為預測輸出序列，c為后驗差檢驗數，error1為殘差，error2為相對誤差
format?long;
n=length（X1）;
X0=[];
X0（1）=X1（1）;
for?i=2:n
????X0（i）=X1（i）-X1（i-1）;?????????????%計算累減生成序列
end
for?i=1:n-1
B（i1）=-0.5*（X1（i）+X1（i+1））;???%計算B，Yn
????B（i2）=0.25*（X1（i）+X1（i+1））^2;
????Y（i）=X0（i+1）;
end
alpha=（B‘*B）^（-1）*B‘*Y‘;????????????%做最小二乘估計
a=alpha（11）;
b=alpha（21）;
c=a*X1（1）;???????????????????????