smooth_dens_construct function

Build the design and penalty matrices for smooth density estimation