added rho via solution; added ode solving