func1p3-spline edi.m

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

x = 0.4:0.1:0.8;
y = [-.916291 -.693147 -.510826 -.356675 -.223144];
cs = spline(x,[0 y 0]);
xx = linspace(0.4,0.8,101);
plot(x,y,'o',xx,ppval(cs,xx),'-');
yyy=ppval(cs,0.54)