This indicates we have a good model for the fit is almost perfect.

We optimize the objective function - no energy costs.

Input := 


P[t_] = -.5 asol[t] + 3.5 bsol[t] + 

			.25 csol[t]/.solk[[2]];

Input := 


Plot[P[t],{t,0,200}]

Output =


-Graphics-

Input := 


solp = FindRoot[P'[t]==0,{t,40}]

Output =


{t -> 29.5001}

Input := 


P[t]/.solp[[1]]

Output =


1.99135