Implement solve_period() in Model class

Initial implementation of solve_period() function in Model class. Still
to add convergence diagnostics and count of number of calls per period.