数值分析,任玉杰教程课后题..docVIP

杨颖 应用数学131802班 201318050228 6.2 1.已知函数在[1,7]上具有二阶连续导数，，且满足条件f(x)=1,f(7)=12,.求线性插值多项式和函数值f(3.5)，并估计其误差. 解： 输入程序 X=[1,7];Y=[1,12]; l01= poly(X(2))/( X(1)- X(2)), l11= poly(X(1))/( X(2)- X(1)), l0=poly2sym (l01),l1=poly2sym (l11), P = l01* Y(1)+ l11* Y(2), L=poly2sym (P),x=3.5; Y = polyval(P,x) 运行后输出基函数l0和l1及其插值多项式的系数向量P（略）、插值多项式L和插值Y为 l01 = -0.1667 1.1667 l11 = 0.1667 -0.1667 l0 =7/6 - x/6 l1 =x/6 - 1/6 P =1.8333 -0.8333 L =(11*x)/6 - 5/6 Y =5.5833 6.3 5.给出节点数据 ，，,作五阶牛顿插值多项式和差商，并写出其估计误差的公式. 解 （1）保存名为newpoly.m的M文件. function [A,C,L,wcgs,Cw]= newploy(X,Y) n=length(