$$\left\{ \begin{array}{lll}\dot{x}_{1} & = & x_{2}\\\dot{x}_{2} & = & \frac{8}{25}\cdot x_{1}^{5}-\frac{4}{3}\cdot x_{1}^{3}+\frac{4}{5}\cdot x_{1}\end{array}\right.$$
Search space: $[-3, 3]\times[-3, 3.1]\setminus{\left\{\mathbf{x}\mid x_{1}^{2}+x_{2}^{2}\leq 0.5\right\}}$
from pyinvariant import *
# Define the search space
space = IntervalVector([[-3, 3],[-3,3.1]])
# Create the grpah structure
smartSubPaving = SmartSubPaving(space)
# Create the Domain
dom = Domain(smartSubPaving, FULL_DOOR)
dom.set_border_path_in(False)
dom.set_border_path_out(False)
f_sep = Function("x[2]", "(x[0])^2+(x[1])^2-(0.5)^2")
s = SepFwdBwd(f_sep, GEQ) # possible options : LT, LEQ, EQ, GEQ, GT
dom.set_sep(s);
# Create the Dynamics
f = Function("x[2]", "(x[1],8.0/25.0*x[0]^5-4.0/3.0*x[0]^3+4.0/5.0*x[0])")
dyn = DynamicsFunction(f, FWD_BWD)
# Create the Maze associated with the Domain and the dynamics
maze = Maze(dom, dyn)
# Contract the system
for i in range(17): # Number of bisections
print(i)
smartSubPaving.bisect()
maze.contract()
# Visualization
visu = VibesMaze("Zang Invariant", maze)
visu.setProperties(0,0,512,512)
visu.show()