LU 分解：[L U］=lu(A)；2
追赶法function x=zhuiganfa（A，b）［n,n]=size（A); for i=1：n if（i==1） l(i）=A(i,i）; y（i)=b（i）/l(i）； else l（i）=A（i，i）—A（i，i-1）＊u(i—1）； y（i）=(b(i)—y（i-1)＊A(i，i—1))/l（i); end if(i>n=101；c1=ones（1，n—1)；a1=diag（c1,—1)； 这个-1 说明行位置—1c2=12＊ones（1,n);a2=diag(c2）; c3=ones(1，n—1）;a3=diag（c3，1); a=a1+a2+a3;3
拉格朗日插值function yh=lage(x,y,xh）n=length（x)；m=length（xh)；yh=zeros（1,m）；c1=ones（n—1，1）；c2=ones(1，m）;for i=1：nxp=x(［1：i—1 i+1：n]);yh=yh+y（i)＊prod((c1*xh—xp'＊c2）
/（x(i）—xp'＊c2））；endend>〉 x=［11,12］;>〉 y=[2，4]；>〉 xh=［11
75]；〉〉 lage（x，y，xh）ans =3
50004 最小二乘法1
最小二乘的 xi 和 yi 为：xi19531964198219902000yi5
66要拟合的函数为：y=a+bx—cxy 注意不是多项式2
编程函数为:function z = erchen（x，y）x1=ones（5,1）； A=［x1,x，-x
＊y]; 注意点乘z=A\y； 注意左除a=z（1）；b=z（2）；c=z(3）；end 输入：≻≻ x=[1953 1964 1982 1990 2000
