Qpenalty function

Penalty function for log-likelihood of selm models