Documentation/Maxima/E4C_delta_omega.max

   1 Momega : MATRIX([cos(omega), 0, -sin(omega)], [0, 1, 0], [sin(omega), 0, cos(omega)]);
   2 Mchi : MATRIX([1, 0, 0], [0, cos(chi), -sin(chi)], [0, sin(chi), cos(chi)]);
   3 Mphi : MATRIX([cos(phi), 0, -sin(phi)], [0, 1, 0], [sin(phi), 0, cos(phi)]);
   4 M2theta : MATRIX([cos(2*theta), 0, -sin(2*theta)],[0, 1, 0],[sin(2*theta), 0, cos(2*theta)]);
   5
   6 R : Momega.Mchi.Mphi;
   7
   8 Q : M2theta.[ki, 0, 0]-[ki, 0, 0];
   9
  10 INVERT(R).Q, TRIGEXPAND = TRUE, EXPAND, omega = theta + C;
  11
  12 MM : RATSIMP(TRIGSIMP(%));