get_2ndDerivPenalty function

Compute Integrated Squared Second Derivative Penalty Matrix for Smoothing Splines