Astronomical Applications Department, U.S. Naval Observatory Linear Least Squares Page 21
A
+
-
-
-
+
4 T
7
2
T
8
2
4 T
7
T
8
T
5
T
1
4 T
4
T
1
4 T
6
T
1
S
0
=
(
)
-
+
+
-
4 T
7
T
3
4 T
0
T
1
2 T
8
T
3
2 T
2
T
1
S
1
(
)
+
-
-
T
2
T
1
2 T
7
T
3
T
8
T
3
2 T
0
T
1
S
2
+
+
(
)
-
+
+
-
+
+
-
4 T
3
T
6
T
3
T
5
4 T
3
T
4
4 T
0
T
7
2 T
7
T
2
2 T
0
T
8
T
2
T
8
S
3
+
/
B
(
)
-
+
+
-
2 T
7
T
3
2 T
0
T
1
T
8
T
3
T
2
T
1
S
0
-
+
2 T
1
N 2 T
3
2
S
1
-
T
1
N T
3
2
S
2
+
+
=
(
)
-
(
)
- +
T
8
2 T
7
N T
3
(
)
-
2 T
0
T
2
S
3
+
/
C
-
-
+
+
-
+
+
1
2
T
2
T
8
2 T
0
T
7
T
0
T
8
T
7
T
2
2 T
3
T
6
1
2
T
3
T
5
2 T
3
T
4
S
0
=
(
)
-
(
)
- +
T
8
2 T
7
N T
3
(
)
-
2 T
0
T
2
S
1
+
-
1
2
T
8
T
7
N
1
2
T
3
(
)
-
2 T
0
T
2
S
2
+
+
+
-
-
2 T
6
1
2
T
5
2 T
4
N
1
2
(
)
-
2 T
0
T
2
2
S
3
+
/
The following commands produce an optimized procedure, which we will not show here.
:=
tmp
(
)
array
(
)
map
,
x
( )
rhs x
(
)
eval
,
soln _subslist
(
)
optimize
(
)
makeproc
,
tmp
=
parameters [
]
, ,
N Y t
Sinusoid with Exponentially Decaying Component
:=
model
=
Y
+
+
+
A B
(
)
cos 2
1
t
i
C
(
)
sin 2
1
t
i
e
(
)
-
t
i
(
)
+
D
(
)
cos 2
2
t
i
E
(
)
sin 2
2
t
i
(
)
leastsqrs
, , ,
model t i [
]
, , , ,
A B C D E
leastsqrs[4]: normal equations
=
+
+
+
-
+
2 N A (
)
-
+
2 N 4 T
11
B 4 C T
12
(
)
-
4 T
13
2 T
15
D 2 S
3
4 E T
18
0
(
)
-
+
2 N 4 T
11
A (
)
-
+
+
8 T
11
2 N 8 T
3
B (
)
-
8 T
6
4 T
12
C
+
+
(
)
-
-
+
+
4 T
13
4 T
9
2 T
15
8 T
7
D (
)
-
8 T
10
4 T
18
E 2 S
3
4 S
5
+
+
+
-
0
=
=
+
+
+
-
+
4 A T
12
(
)
-
8 T
6
4 T
12
B 8 C T
16
(
)
-
8 T
17
4 T
19
D 4 S
4
8 E T
8
0
(
)
-
4 T
13
2 T
15
A (
)
-
-
+
+
4 T
13
4 T
9
2 T
15
8 T
7
B (
)
-
8 T
17
4 T
19
C
+
+
(
)
-
+
8 T
14
8 T
0
2 T
1
D (
)
-
8 T
2
4 T
4
E 4 S
0
2 S
1
+
+
-
+
0
=
=
+
+
+
-
+
4 A T
18
(
)
-
8 T
10
4 T
18
B 8 C T
8
(
)
-
8 T
2
4 T
4
D 4 S
2
8 E T
5
0
Page 21
