ẋ = (z-b)x - dy
ẏ = dx + (z-b)y
ż = c + az - z³/3
- (x²+y²)(1+ez) + fz·x³

Lyapunov: computing...