Another way to add a smoothness constraint is to add the - norm of the derivative to the objective.
or, in matrix form,
The objective function becomes
norm of diff(h) added to the objective function ( ):
Six times the norm of diff(h) added to the objective function ( ):