Astronomical Applications Department, U.S. Naval Observatory Sensitivity Integral Page 28
1
9240 a
4
(
)
+
+
2
3
3 (
)
+
+
+
+
+
+
6
6
5
15
4
21
3
18
2
9
3
1
8
/
10395 a
3
(
)
+
2 (
)
+
+
2
2
2 (
)
+
+
+
+
4
4
3
6
2
4
2
1
7
+
11880 a
2
(
)
+
+
+
+
+
+
6
7
5
21
4
35
3
35
2
21
7
1
6
+
13860 a
1
(
)
+
2 (
)
+
+
2
3
3 (
)
+ +
2
1
1
5
+
16632 a
0
(
)
+
+
+
+
5
10
10
2
5
3
4
1
4
+
:=
tmp
(
)
subs
,
,
(
)
seq
,
=
a
k
.
a k
=
k
..
0
4
=
1
1
%
(
)
cost
(
)
optimize
,
tmp tryhard
+
+
+
31 additions
47 multiplications
divisions
10 assignments
(Remember this has to be evaluated only once per integration over the sky.)
:=
G
(
)
optimize/makeproc
,
[
]
(
)
optimize
,
tmp tryhard
=
parameters
[
]
,
,
1
(
)
seq
,
.
a i
=
i
..
0
4
Now create the geometry function Maple procedure:
(
)
select
,
,
has
(
)
rhs error_expr
s
2520 (
)
+
s
2
r
2
(
)
+ +
g
2
g
1
a
4
-
(
)
+
+
+
+
+
+
+
+ +
+
10
462
5
55
8
330
3
330
6
55
165
7
11
9
11
165
2
462
4
1
10
2772 (
)
+
s
2
r
2
(
)
+ +
g
2
g
1
a
3
(
)
+
2 (
)
+
+
+
+
4
3
3
4
2
2
1
-
(
)
+
+
+
+
5
10
10
2
5
3
4
1
9
3080
-
(
)
+
s
2
r
2
(
)
+ +
g
2
g
1
a
2
(
)
+
+
2
3
3 (
)
+
+
+
+
+
+
6
6
5
15
4
21
3
18
2
9
3
1
8
3465 (
)
+
s
2
r
2
(
)
+ +
g
2
g
1
a
1
(
)
+
2 (
)
+
+
2
2
2 (
)
+
+
+
+
4
4
3
6
2
4
2
1
7
-
3960 (
)
+
s
2
r
2
(
)
+ +
g
2
g
1
a
0
(
)
+
+
+
+
+
+
6
7
5
21
4
35
3
35
2
21
7
1
6
-
:=
tmp
(
)
subs
,
,
(
)
seq
,
=
a
k
.
a k
=
k
..
0
4
=
1
1
%
H :=
(
)
optimize/makeproc
,
[
]
(
)
optimize
,
% tryhard
=
parameters
[
]
,
,
, , ,
1
(
)
seq
,
.
a i
=
i
..
0
4
s r g
Finally, convert the Maple functions into fortran subroutines.
fortran G optimized
=
mode
double
=
precision
double
,
,
,
,
(
=
filename
"d:/FAMEStuff/PhotonSensitivity/error_freq.f" )
Page 28
