clear;clc;clf;x1=02 04 06 08 10;y1=098 092 081 064 038;n=lengthy1;c=y1;for j=2n 求差商 for i=n-1j ci=ci-ci-1x1i-x1i-j+1; endendsyms x df d;df1=1;d1=y11;for i=2n 求牛顿差值多项式 dfi=dfi-1x-x1i-1;

本代码实现了基于五个插值点的4次牛顿插值多项式和三次样条插值函数，并将其在图像中进行比较。

首先，定义了插值点的横纵坐标$x1,y1$，并计算出差商$c$，用于求解牛顿插值多项式。

然后，定义了符号变量$x$，并利用差商和$x1$求解出牛顿插值多项式$P4$。

接着，调用了MATLAB自带的三次样条插值函数csape，并将其结果存储在变量pp中。

最后，利用ezplot函数将牛顿插值多项式和三次样条插值函数的图像进行比较。同时，还绘制了样条插值点的图像，并用绿色的圆点标注出来。