Solve the level occupation with Newton iteration