psp function

Tool to build the basis matrix and the penalty matrix of P-splines.