ruletest2.py 820 B

12345678910111213141516171819202122
  1. import sympy.physics.mechanics as _me
  2. import sympy as _sm
  3. import math as m
  4. import numpy as _np
  5. x1, x2 = _me.dynamicsymbols('x1 x2')
  6. f1 = x1*x2+3*x1**2
  7. f2 = x1*_me.dynamicsymbols._t+x2*_me.dynamicsymbols._t**2
  8. x, y = _me.dynamicsymbols('x y')
  9. x_d, y_d = _me.dynamicsymbols('x_ y_', 1)
  10. y_dd = _me.dynamicsymbols('y_', 2)
  11. q1, q2, q3, u1, u2 = _me.dynamicsymbols('q1 q2 q3 u1 u2')
  12. p1, p2 = _me.dynamicsymbols('p1 p2')
  13. p1_d, p2_d = _me.dynamicsymbols('p1_ p2_', 1)
  14. w1, w2, w3, r1, r2 = _me.dynamicsymbols('w1 w2 w3 r1 r2')
  15. w1_d, w2_d, w3_d, r1_d, r2_d = _me.dynamicsymbols('w1_ w2_ w3_ r1_ r2_', 1)
  16. r1_dd, r2_dd = _me.dynamicsymbols('r1_ r2_', 2)
  17. c11, c12, c21, c22 = _me.dynamicsymbols('c11 c12 c21 c22')
  18. d11, d12, d13 = _me.dynamicsymbols('d11 d12 d13')
  19. j1, j2 = _me.dynamicsymbols('j1 j2')
  20. n = _sm.symbols('n')
  21. n = _sm.I