Calculates the design matrix for fitting the random component of a natural cubic smoothing spline
Useful links