Velocity-Verlet Double-Well Algorithm Python

Question

I am implemententing the Verlet algorithm for a double well potential V(x) = x^4-20x^2, as to create a simple phase portrait.The generated phase portrait has an au