make_matrices_dens function

Build a standardised P-Spline design matrix and the associated P-Spline penalty matrix