Astronomical Applications Department, U.S. Naval Observatory 3 D plot fits Page 1
restart
:=
S
(
)
readdata
,
"E:/TNOs/aeV-R.txt" 3
:=
F
(
)
, ,
x y z
+
+
+
+
+
+ + + +
x
2
b y
2
c z
2
2 d x y
2 e x z
2 f y z
g x
h y
i z
j
:=
eqns
(
)
map
,
@
F
op
(
)
convert
,
S set
(
)
linalg
leastsqrs
,
eqns {
}
, , , , , , , ,
b c d e f g h i j
=
c
-1594.827695
=
d
-36.73945508
=
e
7.008045672
=
f
1734.181721
=
b
-6339.295891
,
,
,
,
,
{
=
g
-58.86482353
=
h
3616.689156
=
i
715.1995081
=
j
311.6372556
,
,
,
}
(
)
subs
,
%
(
)
F
, ,
x y z
:=
fit
(
)
fn
, , ,
% x y z
fit
(
)
, ,
x y z
x
2
6339.295891 y
2
1594.827695 z
2
73.47891016 x y
14.01609134 x z
-
-
-
+
:=
3468.363442 y z
58.86482353 x
3616.689156 y
715.1995081 z
311.6372556
+
-
+
+
+
ContourPlanePlot
::
F procedure
::
v1min {
}
,
integer float
::
v1max {
}
,
integer float
,
,
,
proc(
:=
::
v2min {
}
,
integer float
::
v2max {
}
,
integer float
::
v0 {
}
,
integer float
::
v3min {
}
,
integer float
::
plane string
,
,
,
,
,
::
Ncontours posint
::
cgrid list
::
contour_colors list
,
,
)
local
;
, ,
p f tmp
=
plane
"xy"
if
then
p
[
]
plots contourplot
(
)
F
, ,
x y v0
=
x
..
v1min
v1max
=
y
..
v2min
v2max
=
contours
Ncontours
,
,
,
,
(
:=
=
filled
true
=
coloring
contour_colors
=
grid
cgrid
,
,
);
:=
f
(
)
[
]
plottools transform
(
)
,
x y
[
]
, ,
x y v3min
=
plane
"yz"
elif
then
p
[
]
plots contourplot
(
)
F
, ,
v0 y z
=
y
..
v1min
v1max
=
z
..
v2min
v2max
=
contours
Ncontours
,
,
,
,
(
:=
=
filled
true
=
coloring
contour_colors
=
grid
cgrid
,
,
);
:=
f
(
)
[
]
plottools transform
(
)
,
y z
[
]
, ,
v3min y z
else
p
[
]
plots contourplot
(
)
F
,
,
x v0 z
=
x
..
v1min
v1max
=
z
..
v2min
v2max
=
contours
Ncontours
,
,
,
,
(
:=
=
filled
true
=
coloring
contour_colors
=
grid
cgrid
,
,
);
:=
f
(
)
[
]
plottools transform
(
)
,
x z
[
]
,
,
x v3min z
fi;
( )
f p
end
Warning, `v3min` is a lexically scoped parameter
Warning, `v3min` is a lexically scoped parameter
Warning, `v3min` is a lexically scoped parameter
:=
VRmin
.2
:=
VRmax
1.2
:=
amin
17
:=
amax
50
:=
emin
0
:=
emax
.7
Page 1
