This project involves analyzing the stationary solutions of a system of differential equations depending on the parameter ( p4 ). The primary goal is to plot the dependence of the stationary solutions ( x1(p4) ), ( x2(p4) ), and ( x3(p4) ) on ( p4 ), and mark stable and unstable stationary points, as well as points of real and complex bifurcation if they exist.
The differential equations describing the system are:
Where:
$p1 = 8.4 \times 10^{-6}$ $p2 = 6.6667 \times 10^{-4} $ -
$p4$ is the varying parameter $p5 = 2 $ -
$p6 = 10$ or$p6 = 100$
-
Install dependencies: Ensure you have the necessary Python packages installed. You can do this via pip:
pip install numpy matplotlib
-
Run the script: Execute the main script to perform the bifurcation analysis.
python bifurcation_analysis.py
-
Output: The script will generate and display several plots:
- Plots for ( x1 ), ( x2 ), and ( x3 ) as functions of ( p4 ).
- Stability of stationary points will be indicated with different colors (e.g., green for stable, red for unstable).
Calculates the Jacobian matrix for the system at given values of ( x1 ), ( x2 ), ( x3 ), and ( p4 ).
Checks if the provided stationary points satisfy the system of equations.
Main function to compute and plot the stationary solutions and their stability. It iterates over a range of ( p4 ) values, computes the stationary points and their stability, and then generates the plots.
The results of the analysis are saved in text files and plotted using matplotlib. The plots show how the stationary solutions ( x1 ), ( x2 ), and ( x3 ) vary with ( p4 ), and they also indicate the stability of these solutions.