Defined Q's and summed up the ones I assume(?) have the same rates. eps has been defined in the loop. rho is now being updated in the loop with the equation of state. Replaced exp with np.exp. TODO: It seems r is a function and is also assumed to be a radius array??? Fix this.