We compute the net return as a function of time and energy rate s $/min to run the reaction.

Input := 


net[t_,s_] = - s t -.5 asol[t] + 3.5 bsol[t] +

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