尊敬的各位大侠：

    你好，我现在想要用fortran编三次样条插值程序,想要实现和matlab的spline函数一样的效果，这个spline函数边界条件是not-a-knot边界条件，也叫非扭结边界条件。

    我上网搜了很多fortran写的三次样条插值程序，没有一个是按照这种边界条件写的，很多都是用自然边界条件，因为这样能简化计算，结果也就可想而知的不好了。

    哪位大侠能否帮忙用fortran写一个基于非扭结边界条件的三次样条插值程序？