func1p2.m

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

%func1p2
%求自然边界条件的三次样条
x = [27.7,28,29,30];
y = [4.1,4.3,4.1,3.0];
cs = spline(x,[0 y 0]);
xx = linspace(27.7,30,101);
plot(x,y,'o',xx,ppval(cs,xx),'-');
%样条函数请见图像