penalty.matrix function

Calculating the penalty matrix P(lambda)