Solving 4-Coupled ODEs with Mathematica

  • Thread starter Thread starter Qyzren
  • Start date Start date
  • Tags Tags
    Mathematica Odes
Click For Summary
SUMMARY

The discussion focuses on solving a system of four coupled ordinary differential equations (ODEs) using Mathematica. The equations are defined as follows: (1) dy1/dt = y2y3 - μy1, (2) dy2/dt = y1y4 - μy2, (3) dy3/dt = 1 - y1y2, and (4) dy4/dt = 1 - y1y2. The steady state solutions are established as y1 = ±k, y2 = ±k⁻¹, y3 = μk², and y4 = μk⁻¹, where μ(k² - k⁻²) = A, a constant. The user initially struggled with using DSolve in Mathematica but resolved the issue by recognizing that steady state solutions can be found by setting the derivatives to zero, simplifying the problem to algebra.

PREREQUISITES
  • Understanding of ordinary differential equations (ODEs)
  • Familiarity with Mathematica version 12.3 or later
  • Knowledge of steady state analysis in dynamical systems
  • Basic algebraic manipulation skills
NEXT STEPS
  • Explore the use of Mathematica's DSolve function for solving ODEs
  • Learn about steady state solutions in nonlinear systems
  • Investigate numerical methods for solving coupled ODEs
  • Study the implications of parameter μ in dynamical systems
USEFUL FOR

This discussion is beneficial for mathematicians, physicists, and engineers working with dynamical systems, particularly those dealing with coupled ordinary differential equations and seeking to utilize Mathematica for analysis and solutions.

Qyzren
Messages
41
Reaction score
0
Hi guys, i have 4-coupled ode's that are giving trouble

(1) \frac{dy_1}{dt}=y_2y_3-\mu y_1, \hspace{1cm} \\(2) \frac{dy_2}{dt}=y_1y_4-\mu y_2, \hspace{1cm} \\(3) \frac{dy_3}{dt}=1-y_1y_2, \hspace{1cm} \\(4) \frac{dy_4}{dt}=1-y_1y_2
I need to show that the steady state solutions are
y_1=\pm k, y_2=\pm k^{-1}, y_3=\mu k^2, y_4 = \mu k^{-1} where \mu (k^2-k^{-2})=A a const.
now in an earlier part of the question, I was able to show that y_3-y_4=A.
But trying to solve these coupled ODE's is giving trouble. I tried solving this in mathematica using DSolve as well, and mathematica doesn't seem to know how to do it either.

mathematica code:
\text{DSolve}\left[\left\{\text{y1}'[t]==-\mu \text{y1}[t]+\text{y2}[t] \text{y3}[t],\text{y2}'[t]==-\mu \text{y2}[t]+\text{y1}[t] \text{y4}[t],\text{y3}'[t]==1-\text{y1}[t] \text{y2}[t],\text{y4}'[t]==1-\text{y1}[t] \text{y2}[t]\right\},\{\text{y1}[t],\text{y2}[t],\text{y3}[t],\text{y4}[t]\},t\right]

PS: how do i get each of the differential equations on a new line? \\ and \newline didn't work.
each of the ODE's starting with the derivative should be on a new line
 
Last edited:
Physics news on Phys.org
You don't have to solve the ODE's to find steady state solutions, do you? Just put all of the derivatives equal to zero. Now it's just an algebra problem.
 
thanks, problem solved
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 3 ·
Replies
3
Views
3K
Replies
6
Views
2K
  • · Replies 2 ·
Replies
2
Views
5K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 6 ·
Replies
6
Views
3K
Replies
2
Views
2K
  • · Replies 6 ·
Replies
6
Views
3K