山东师范大学数学科学学院实验报告
2) 选取插值节点，计算插值节点画出cosx 在[0,1.2]上图像,如图3-1，并标出插值算法所计算的点。
>> X=[0 0.2 0.4 0.6 0.8 1.0 1.2;1 0.980067 0.921061 0.825336 0.696707 0.540302 0.362358]; >> x1=0:0.1:1.2; >> format long
>> [A,Pn]=Newton(X,x1) >> x2=0:0.01:1.2; >> plot(x2,cos(x2)) >> hold on
>> plot(x1,Pn,'*')
>> legend('cosx','Newton 插值');
用Newton 插值法计算出x1=0:0.1:1.2处的值为表3-2：
表3-2 x
0 0.1 0.2
0.3 0.4 0.5 0.6 Pn 1 0.995005 0.980067 0.955337 0.921061 0.877583 0.825336
x 0.7 0.8 0.9 1 1.1 1.2
Pn 0.764843 0.696707 0.62161
0.540302 0.453596 0.362358
计算插值节点画出cosx 在[0,1.2]上图像，并标出插值算法所计算的点。
图3-1
总结
（1）使用Newton 插值法编程，在计算差商表时应注意把),,,(10k x x x f 到
)
,,,(1n k n k n x x x f +--),,2,1（n k =都计算出来； （2）Newton 插值法的误差大约为-7101⨯。
