ruletest10.al 781 B

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758
  1. % ruletest10.al
  2. VARIABLES X,Y
  3. COMPLEX ON
  4. CONSTANTS A,B
  5. E = A*(B*X+Y)^2
  6. M = [E;E]
  7. EXPAND(E)
  8. EXPAND(M)
  9. FACTOR(E,X)
  10. FACTOR(M,X)
  11. EQN[1] = A*X + B*Y
  12. EQN[2] = 2*A*X - 3*B*Y
  13. SOLVE(EQN, X, Y)
  14. RHS_Y = RHS(Y)
  15. E = (X+Y)^2 + 2*X^2
  16. ARRANGE(E, 2, X)
  17. CONSTANTS A,B,C
  18. M = [A,B;C,0]
  19. M2 = EVALUATE(M,A=1,B=2,C=3)
  20. EIG(M2, EIGVALUE, EIGVEC)
  21. NEWTONIAN N
  22. FRAMES A
  23. SIMPROT(N, A, N1>, X)
  24. DEGREES OFF
  25. SIMPROT(N, A, N1>, PI/2)
  26. CONSTANTS C{3}
  27. V> = C1*A1> + C2*A2> + C3*A3>
  28. POINTS O, P
  29. P_P_O> = C1*A1>
  30. EXPRESS(V>,N)
  31. EXPRESS(P_P_O>,N)
  32. W_A_N> = C3*A3>
  33. ANGVEL(A,N)
  34. V2PTS(N,A,O,P)
  35. PARTICLES P{2}
  36. V2PTS(N,A,P1,P2)
  37. A2PTS(N,A,P1,P)
  38. BODIES B{2}
  39. CONSTANT G
  40. GRAVITY(G*N1>)
  41. VARIABLE Z
  42. V> = X*A1> + Y*A3>
  43. P_P_O> = X*A1> + Y*A2>
  44. X = 2*Z
  45. Y = Z
  46. EXPLICIT(V>)
  47. EXPLICIT(P_P_O>)
  48. FORCE(O/P1, X*Y*A1>)
  49. FORCE(P2, X*Y*A1>)