Astronomical Applications Department, U.S. Naval Observatory Precession Maple Page 26
=
( )
sin
( )
cos
( )
sin
( )
sin
-
( )
cos
-
( )
sin
( )
cos
0
( )
cos
( )
cos
( )
cos
( )
sin
( )
sin
&*
(
)
R
,
,
&*
X
Y
Z
(
)
factormat
,
,
(
)
evalm %
collect [
]
,
( )
sin
( )
cos
=
( )
sin
(
)
-
-
( )
sin
( )
cos
( )
cos
(
)
cos
( )
sin
X
(
[
( )
sin
(
)
-
+
( )
sin
( )
sin
( )
cos
(
)
cos
( )
cos
Y
( )
sin
( )
cos
(
)
sin
Z
+
+
)
( )
sin
( )
sin
(
)
-
( )
cos
( )
cos
( )
sin
(
)
cos
( )
sin
X
(
+
( )
sin
(
)
+
( )
cos
( )
sin
( )
sin
(
)
cos
( )
cos
Y
( )
sin
( )
sin
(
)
sin
Z
+
+
)
( )
cos
X
( )
cos
(
)
sin
( )
sin
Y
( )
cos
(
)
sin
( )
cos
Z
( )
cos
(
)
cos
-
+
-
]
(
)
-
+
( )
cos
( )
cos
( )
sin
(
)
cos
( )
sin
X
(
[
(
)
-
-
( )
cos
( )
sin
( )
sin
(
)
cos
( )
cos
Y
( )
sin
(
)
sin
Z
+
-
)
( )
sin
(
+
(
)
-
-
( )
sin
( )
cos
( )
cos
(
)
cos
( )
sin
X
(
)
-
+
( )
sin
( )
sin
( )
cos
(
)
cos
( )
cos
Y
Z
( )
cos
(
)
sin
+
+
)
( )
cos
]
( )
cos
(
)
-
-
( )
sin
( )
cos
( )
cos
(
)
cos
( )
sin
X
(
[
( )
cos
(
)
-
+
( )
sin
( )
sin
( )
cos
(
)
cos
( )
cos
Y
( )
cos
( )
cos
(
)
sin
Z
+
+
)
( )
sin
( )
cos
(
)
-
( )
cos
( )
cos
( )
sin
(
)
cos
( )
sin
X
(
+
( )
cos
(
)
+
( )
cos
( )
sin
( )
sin
(
)
cos
( )
cos
Y
( )
cos
( )
sin
(
)
sin
Z
+
+
)
( )
cos
X
( )
sin
(
)
sin
( )
sin
Y
( )
sin
(
)
sin
( )
cos
Z
( )
sin
(
)
cos
+
-
+
]
:=
Psubs
[
]
(
)
seq
,
=
(
)
lhs %
,
k 1
(
)
rhs %
,
k 1
=
k
..
1
3
We must transform the torque components from the conical frame to the body cartesian
frame.
torque_cone
Page 26