Interior Penalty Function Method