solve dV/dt= a(z(t))*V^2 + b(z(t))*V + c(z(t)) for V