func1p3-fail.m

数值分析最常用的四十种算法

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),'-');