Astronomical Applications Department, U.S. Naval Observatory Precession Maple Page 24
K
K
K
0
2
a
( )
sin
b
( )
sin
P
( )
cos
=
-
( )
sin
V
-
+
-
( )
sin
(
)
-
Q
2 A
C
( )
cos
(
)
+
( )
sin
( )
cos
V
2
( )
cos
U
( )
cos
A
C
-
Q
( )
sin
( )
cos
V
d
d
Evaluating the integral in
, we find
map int
(
)
map
,
,
collect
(
)
evalm
(
)
op
,
[
]
, ,
2 1 1
(
)
subs
,
=
Prho torque_cone_integral
,
,
=
..
a
( )
sin
b
( )
sin
factormat
(
)
subs
,
,
=
Prho
=
V
Q U
(
)
map
,
expand
(
)
evalm %
collect
,
,
(
[
]
,
,
,
,
,
,
( )
cos
A
C
Q U
factor
,
)
(
)
subsop
,
=
[
]
,
2 1
% %%%
:=
torque_cone
%
K
K
K
0
2
P
( )
cos
(
)
-
a
b
=
1
2
(
)
+
a
b U Q
( )
sin
-
(
)
+
a
b
( )
cos
U
( )
sin
2
2
3
+
+
a
2
a b
b
2
( )
sin
2
A
C
( )
cos
1
2
(
)
+
a
b U Q
( )
sin
-
-
+
1
2
(
)
+
a
b
( )
cos
U
( )
sin
2
1
3
+
+
a
2
a b
b
2
( )
sin
2
Q
+
-
1
2
(
)
+
a
b
( )
cos
U
( )
sin
2
1
3
+
+
a
2
a b
b
2
( )
sin
2
Q
d
Page 24
