Calculates the variance matrix of the random effects for a natural cubic smoothing spline
Useful links