double_pendulum.al 427 B

12345678910111213141516171819202122232425
  1. MOTIONVARIABLES' Q{2}', U{2}'
  2. CONSTANTS L,M,G
  3. NEWTONIAN N
  4. FRAMES A,B
  5. SIMPROT(N, A, 3, Q1)
  6. SIMPROT(N, B, 3, Q2)
  7. W_A_N>=U1*N3>
  8. W_B_N>=U2*N3>
  9. POINT O
  10. PARTICLES P,R
  11. P_O_P> = L*A1>
  12. P_P_R> = L*B1>
  13. V_O_N> = 0>
  14. V2PTS(N, A, O, P)
  15. V2PTS(N, B, P, R)
  16. MASS P=M, R=M
  17. Q1' = U1
  18. Q2' = U2
  19. GRAVITY(G*N1>)
  20. ZERO = FR() + FRSTAR()
  21. KANE()
  22. INPUT M=1,G=9.81,L=1
  23. INPUT Q1=.1,Q2=.2,U1=0,U2=0
  24. INPUT TFINAL=10, INTEGSTP=.01
  25. CODE DYNAMICS() some_filename.c