Numerically solving the equation of motion for inflation in cosmology

So I want to solve the equation of inflation involving a scalar field numerically using python libraries such as odient or scipy. The equation is:
$$defphi{varphi}
phi”(t)+3Hphi'(t)+V'(phi(t))=0$$
with given condition: $3Mp^2H^2 = V(phi)$.
Can anyone help me with the code and advice on how to progress on this?