Gu function

Tool to build the basis matrix and the penalty matrix of natural cubic splines.