This is 2D case only?
Is there parametrized interpolation like evalAt(t=0...1)?
It fails in general, for messy points x=[0; 6 ; 3 ; 12; 8] and y=[0;-3; 6; 12; 0]?
If I am not wrong, it works similar as function provided with Matlab...
If You could extend such possibilities in 3D - it would be great, because Matlab is lacking with CAD features.
Richard has sent me an email with this solution for his problems which stem from using a different compiler (VC++ on Windows rather than gcc on Linux):
"Under Windows it is important that all files have the ending .cpp instead of .c so that the right compiler is taken. After correcting that (*.c à *.cpp) and stepping through all calls of .c files in the code which also have to be corrected, of course, everything works perfect!"