Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Solving system of ODE

  1. Jan 9, 2009 #1
    Dear all,

    I have been trying to solve the following system of first-order ordinary differential equations for a week:
    x' = y * (a1*x + a2*y + c1),
    y' = y * (a3*x + a4*y + c2),
    where x and y are functions of t, and ai and ci are constants. This system seems not very complex, but I have not found the its exact solutions.

    Can anyone tell me the exact solution for x(t)? Many thanks!

    Best regards,
    Last edited: Jan 9, 2009
  2. jcsd
  3. Jan 9, 2009 #2
    I assume you mean a4 in the 2nd equation, can't you just look at the quotient of the two?
  4. Jan 9, 2009 #3
    Hi 'NoMoreExams', thank you very much. The typo has been corrected.

    Could you tell me how you solve this problem by making use of the quotient? Many thanks for your patience and time.
  5. Jan 9, 2009 #4
    Well what does y'/x' look like? Note that you are dropping a solution.
  6. Jan 9, 2009 #5
    I tried this way... y'/x' = (a3*x + a4*y + c2)/(a1*x + a2*y + c1), and we could find some relationship between x*y', y^2, x'*y, x^2, x', and y'. But what follows? Sorry, I am confused.
  7. Jan 9, 2009 #6
    No... if y' = dy/dt and x' = dx/dt then y'/x' = ?
  8. Jan 9, 2009 #7
    That's cool. But does that mean y MUST BE A FUNCTION of x? Originally this does not necessarily hold. (For example, can we use this method if x = sin(t) and y = t^2?)
  9. Jan 9, 2009 #8
    ... why doesn't it work in that case? x = sin(sqrt(y)) or x = sin(-sqrt(y))
  10. Jan 9, 2009 #9
    You are very smart :)
    I am not sure if I am right now, because it becomes more and more complicated as I continue this procedure... Could you tell me the result, if it is simple to you? Thanks once again.
  11. Jan 9, 2009 #10
    Well it looks like it can be made into an exact ODE so I'd go that route probably
  12. Jan 10, 2009 #11
    You should also try to go about this in a different way, for your own education. Find the fixed points, figure out their stability and try to draw the phase portrait yourself.
  13. Jan 10, 2009 #12
    Hi Kenneth,

    For the system of ODEs as you have them in full generality, one cannot hope to write the nontrivial solutions in terms of elementary functions (a computer package-generated solution involving the Lambert W function et al instructs me as much). Investigating the behaviour of solutions by way of phase portraits, as NoMoreExams suggested, would seem more appropriate. Indeed, performing this kind of analysis for various values of the constants hints at the impossibility of constructing an exact solution.
  14. Jan 11, 2009 #13
    Hi Unco and 'NoMoreExams',

    Thank you very much for your help. Actually I am just interested in the case where c2 = -c3. I have obtained the numerical solution using MATLAB but it is not enough since I need to find the property of all possible x or \int x dt. Now I am trying the way of phase portraits :)
  15. Jan 12, 2009 #14
    Hi Unco,

    Could you please tell me what package you used to generate the solution involving the Lambert W function (and the solution you found)? It is OK that the solution could not be written in terms of elementary function. I just need to analyze the behavior of the solution.

    Many thanks!

Share this great discussion with others via Reddit, Google+, Twitter, or Facebook