文档介绍：%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%functionp=newpoly(x,y,n,x)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%p=newpoly(x,y,n)%(x,y)为节点坐标%p为所求的的差值多项式的系数向量%n为牛顿插值次数,n=时为线性插值,n=时为二次插值...N=length(x);d=[x;y];ifnargin==fori=:N%根据插值误差估计式(根据距x的远近)重新排列插值节点forj=i+:Nifabs(d(,i)-x)>abs(d(,j)-x)t=d(:,i);d(:,i)=d(:,j);d(:,j)=t;endendendenddq=zeros(n+,n+);dq(:,)=d(,:n+)';dq(:,)=d(,:n+)';fori=:n+forj=:i+dq(i,j)=(dq(i,j-)-dq(i-,j-))/(dq(i,)-dq(i-j+,));%差商表endends=diag(dq(:end,:end))';p=dq(,);t=;fori=:nt=conv(t,[-dq(i,)]);p=polysum(p,s(i).*t);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function[x,k,flag]=SOR(A,b,delta,w,step)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%函数格式:[x,k,flag]=SOR(A,b,delta,w,step)%A为方程组的系数矩阵%b为方程组右端项%delta为精度要求%step为最大迭代次数,缺省值为%w为松弛因子,w<时为低松弛法,w>时为超松弛法,w=时为高斯-塞德尔迭代%x为方程组的解%k为迭代次数%flag为收敛标志ifnargin<step=;endifnargin<w=;endifnargin<ep=e-;endn=length(A);k=;x=zeros(n,);flag=;whiley=x;fori=:nz=b(i);forj=:nifj~=iz=z-A(i,j)*x(j);endendifabs(A(i,i))<e-||k==stepflag=;return;endz=z/A(i,i);x(i)=(-w)*x(i)+w*z;endifnorm(y-x,inf)<deltabreak;endk=k+;end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%functiony=GaussLegendre(func,a,b,n)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%y=GaussLegendre(func,a,b,n)%func:被积函数表达式,例如'x.^.*cos(x)'%a:积分下限%b:积分上限%n:节点数xTable=[....NaN;.-.;.-..;.-..-.;.-..-..;.-..-..-.;];ATable=[....;..;...;....;.....;......;];f=inline(func);x=xTable(n,:);A=ATable(n,:);T=zeros(,n);T=(a+b)/+(b-a)/*x;y=(b-a)/*sum(A.*feval(f,T));%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function[y,Ck,Ak]=NewtonCotes(func,a,b,n)%%%%%%%%%%%%%%%%%