差商

计算⽅法_Newton插值法（差商）

clearall;%函数插值——Newton插值

n=input('请输⼊插值节点个数（正整数）:');

fori=1:n

F(i,1)=input('输⼊节点值:');

F(i,2)=input('输⼊对应节点函数值值:');

end

forj=3:n+1

fori=j-1:n

F(i,j)=(F(i,j-1)-F(i-1,j-1))/(F(i,1)-F(i+2-j,1));%差商

end

end

F%矩阵F即为差商表

运⾏结果：

请输⼊插值节点个数（正整数）:4

输⼊节点值:>>-1.1

输⼊对应节点函数值值:-1.12

输⼊节点值:-0.5

输⼊对应节点函数值值:0

输⼊节点值:0.9

输⼊对应节点函数值值:1.8

输⼊节点值:1.7

输⼊对应节点函数值值:2.2

F=

-1.1000-1.1200000-0.500001.866700

0.90001.80001.2857-0.29050

1.7000

2.20000.5000-0.3571-0.0238