solve_pbtk_lifestage function

Solve the pbtk_lifestage model, which has time-dependent parameters