aquastar
/ gist:861f0eed23ad77e03050982a26544d38
Created
July 23, 2024 23:17
The code simulates the chaotic behavior of the Lorenz system for two slightly different initial conditions, calculates the trajectory differences and Lyapunov exponent over time, plots the results in three subfigures, and saves the output as a PDF file with enlarged font size for better readability. In the results, the first subfigure shows the …
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| import numpy as np | |
| import matplotlib.pyplot as plt | |
| from scipy.integrate import solve_ivp | |
| # Set the font size | |
| plt.rcParams.update({'font.size': 16}) | |
| # Lorenz system parameters | |
| sigma = 10.0 | |
| rho = 28.0 |