数值分析差商表和牛顿差值的matlab程序

差商表functionA=csb(X,Y)n=length(X);A=zeros(n,n+1);A(:,1)=X';A(:,2)=Y';forj=3:n+1fori=j-1:nA(i,j)=(A(i,j-1)-A(i-1,j-1))./(A(i,1)-A(i-j+2,1));endend牛顿差值function[NR]=ndcz(X,Y,x,M);n=length(X);A=csb(X,Y);N=A(1,2);a=zeros(n,1);fori=1:n-1;a(i)=1;forj=1:i;a(i)=a(i)*(x-X(j));endN=N+a(i)*A(i+1,i+2);endk=1;fori=1:n;k=k*(x-X(i))/i;k=abs(k);endR=M*k;