Astronomical Applications Department, U.S. Naval Observatory Precession Maple Page 37
:=
L
(
)
location
,
eqs
(
)
op %
subsop
,
=
%
(
)
op %%
(
)
sin
eqs
select
,
,
has %
t
:=
L
(
)
location
,
%%
(
)
op %
subsop
,
=
%
(
)
op %%
(
)
sin
%%%
:=
odesubs
%
(
)
mat %
t
(
)
-
(
)
-
1
(
)
cos
(
)
+
1
=
(
)
-
+
(
)
cos
Z
Y
(
)
sin
( )
cos
X
(
)
sin
( )
sin
G (
)
+
X
( )
cos
Y
( )
sin
I
xy
+
(
)
sin
/
t
+
(
)
- +
1
2
(
)
cos
(
)
sin
(
+
=
(
)
-
+
(
)
cos
Z
Y
(
)
sin
( )
cos
X
(
)
sin
( )
sin
G
(
)
-
-
+
(
)
sin
Z
(
)
cos
Y
( )
cos
( )
sin
(
)
cos
X
I
xy
)
t
(
)
+
(
)
+
(
)
cos
2
1
(
)
- +
1
(
)
cos
=
(
)
-
+
(
)
cos
Z
Y
(
)
sin
( )
cos
X
(
)
sin
( )
sin
G (
)
+
X
( )
cos
Y
( )
sin
(
)
cos
I
xy
-
/
(
)
sin
Differentiate the equations of motion.
(
)
map
,
,
diff eqs t
:=
eqs2
%
(
)
mat %
+
(
)
cos
t
t
(
)
sin
2
t
2
=
Page 37
