If I have a 3D differential system of 3 variables and considerable number of parameters, and I on PURPOSE would like this system to exhibit chaos, how would I go about ensuring that the system has positive Lyapunov exponents analytically (rather than using the hit-and-try method)? My first thought is to ensure that the divergence of the system is less than zero AND use the Routh Hurwitz criterion on the solution of eigenvalues to make sure that the system is stable. But,all this still doesn't give analytical proof that the Lyapunov exponents are positive (unless I run it in a program and check).

Similar questions and discussions