penalty_uni function

Penalty approximation of unimodality constraints for univariates smooths