dx = (z-b) * x - d*y

dy = d * x + (z-b) * y

dz = c + a*z - z3 /3 - x2 + f * z * x3