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

Found a cool Non-autonomous Coupled System of ODEs

  1. Mar 6, 2014 #1
    By a fortuitous mistake in copy/pasting, I happened across this system. It exhibits very chaotic behavior for some initial values and spiral-like shapes for others. (About a 70%/30% split, respectively.)


    Here's an album of it plotted from t=0 to 1000. The titles of the plot are x0_y0.
    A lot of those have serious errors in their long term trends, but the overall type of pattern they take is, I believe, realistic. I believe so, because, as you can see in the video below, their vector fields are locally linear. From that, I intuitively draw the conclusion that there could be no serious numerical error, like jumping a phase line, but only small compounding error that only affects long term trend. That is to say, the only kind of numerical error present is the kind inherent to numerical methods, not anything specific to this ODE. (Correct me if my intuitive reasoning is incorrect.)

    And, here's a video showing (x0_y0)=(1,1) with surrounding vector field from t=0 to 250:

    If anyone would like, I can post a video that is zoomed closer for the entire video, so as to see the local linearity better in the faster moving parts. I have only uploaded the one below so far because, IMO, it is more fun to watch that the one above.
    Last edited: Mar 6, 2014
  2. jcsd
  3. Mar 6, 2014 #2
    Also, I've reduced it to a system of very complex 2nd order noncoupled ODEs, but that doesn't seem to help in solving it. Mathematica can't solve it either. I'm almost convinced it isn't solvable. So it seems the only way to explore it is with numerical methods.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook