For a synchronized ODE system,I wanted to know if there is a program code in MATLAB available for plotting the conditional Lyapunov exponent. For example, synchronization of identical Rossler system with a coupling term $F$ present in both the response and the driver.The $F$ is not a coefficient of any variable. This is where I have a problem.
- Do we have to include the coupling in the jacobian matrix and if yes how?
- How to plot the Lyapunov exponent as a function of different values of $F$?
Thank you for any help.