Survspline function

Royston/Parmar spline survival distribution