%線性最小二乘逼近的求法
x=input（‘輸入x觀測值：‘）;
y=input（‘輸入y觀測值：‘）;
n=input（‘輸入階次n：‘）;
m=length（x）;
b=0;
A=0;
if?n>m
????disp（‘階次應小于觀測值個數(shù)‘）;
else
for?i=1:m
????for?j=0:n
????????A（ij+1）=x（i）^j;
????end
end
G=A‘*A;
for?j=0:n
????d=0;
????for?i=1:m
????????d=d+y（i）*x（i）^j;
????end
????b（j+1）=d;
end
?