Hi,I found a bug in csape from the package splines v1.0.4 that can be reproduced by doing:
>> x = linspace(0,1,3); >> pp = csape(x, x, 'variational'); >> xs = linspace(0,1,100); >> ys = ppval(pp,xs); >> plot(xs,ys)Unfortunately I don't have time to debug this myself at the moment, but I attach a function implementing cubic splines with conditions on the 2nd derivatives that I had implemented myself some time ago, maybe it can serve as a reference to someone willing to solve the problem in csape.
c.
naturalcubicspline.m
Description: Binary data
------------------------------------------------------------------------- This SF.net email is sponsored by: Microsoft Defy all challenges. Microsoft(R) Visual Studio 2008. http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/
_______________________________________________ Octave-dev mailing list Octave-dev@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/octave-dev
