Variational Characterization of Cubic Smoothing Splines