Penalty.matrix function

Penalty Matrix for Higher Order Differences