Double Pendulum

  • #1
Hello all,
I have written a python code implementing the Runge-Kutta fourth order method for higher orders to approximate the motion of a double pendulum. The problem I am having is that my plot of theta2 (angle of second rod) Vs time looks a little off and I am curious if I have it correct.
Attached is a graph of both rod angles Vs time. The blue is the first rod and the green is the second rod. I figure the opinion of some physicists may be insightful. Thanks :)
-thespaceman
 

Attachments

  • DoublePendulum_N-1000000_1stIC.png
    DoublePendulum_N-1000000_1stIC.png
    20.7 KB · Views: 529

Answers and Replies

  • #2
vanhees71
Science Advisor
Insights Author
Gold Member
2021 Award
19,540
10,298
It's not so easy to see just from the plot, whether there are numerical problems or not. A good check is, whether the total energy of the system remains conserved in the numerical simulation.
 

Related Threads on Double Pendulum

  • Last Post
Replies
2
Views
5K
  • Last Post
Replies
1
Views
5K
Replies
7
Views
1K
  • Last Post
Replies
4
Views
811
Replies
2
Views
2K
  • Last Post
Replies
3
Views
1K
  • Last Post
Replies
4
Views
849
Replies
6
Views
651
Replies
1
Views
3K
  • Last Post
Replies
9
Views
2K
Top