p58_csape.m

来自「数值分析, 同济大学教材<<现代数值数学和计算>>数值分」· M 代码 · 共 14 行

x = [1 2 4 5];
y0 = [2 3 4 2];
y = [0 1 3 4 2 0];
pp = csape(x, y, 'complete')
[breaks, coefs, npolys, ncoefs, dim] = unmkpp(pp)

pp1 = csape(x, y, 'second')
[breaks, coefs, npolys, ncoefs, dim] = unmkpp(pp1)

pp2 = csape(x, y0, 'periodic')
[breaks, coefs, npolys, ncoefs, dim] = unmkpp(pp2)