difsm function

Smoothing with a finite difference penalty