Halvorsen Chaotic Attractor

A Python code in order to illustrate and animate the Halvorsen chaotic attractor, which is described by a particular system of 3 non-linear ODEs.

$$ \begin{align} \begin{cases} \frac{dx}{dt}=-ax-4y-4z-y^2\\ \frac{dy}{dt}=-ay-4z-z^2\\ \frac{dz}{dt}=-az-4x-4y-x^2 \end{cases} \end{align} $$

Here the link for the Animated Gif