func9.asv

来自「数值分析最常用的四十种算法」· ASV 代码 · 共 11 行

x = linspace(0,10,6)
y=[4,2.5,1.5,0.5,-0.5,4.0]
xx = 0:.25:10;
pp=csapi(x,y)        %求三次样条插值多项式,in not-a-knot end condition
yy=ppval(pp,xx)
disp(xx)
plot(x,y,'o',xx,yy)

x1=5
yd=polyder(pp)
y1=polyval(