附录六 球面的计算过程
这是1985年的计算。使用袖珍程序电算器,型号为EC-512,共128步存贮容量,即
只能设置128步程序;有9个数据存贮器。这在当时已是最好的电算器,因为那时只有
286的电脑。
1. 把公式数值化
在计算公式中代入表1中的单位,例如
类似的换算使公式(33),(34)和(36)成为
把上面 a
的算式代入 Zv
和 Ze
的算式中,消去 a
,然后用公式(37)进行计算。
由于电算器只允许有128步程序,所以必须利用数据存贮器,特作如下设置:
STO1=27.85918952
STO2=0.240582497
STO7=H
STO9=M
在存贮公式之前,先把常数 STO1
、STO2
、STO7
、STO9
键入数据存储器。
设 m
为自变量,按如下的计算顺序输入电算器:
2.
H
、M
、m
的计算方法:
1) 设H=H1
,M=M1
,当给 m
一个数值时,计算器将显示一个Z
值,直到输入的
m
使 Z=2.3
为止。这时求出一组数:
( H1 M1
m1
)
2) 设H=H1
,M=M2
,适当选取 m
之后,得到一组数:
( H1 M2
m2
)
3) 设H=H1
,M=M3
,适当选取 m
之后,得到一组数:
( H1 M3
m3
)
把所设的M
计算完之后,再设H=H2
,并且重复上述计算步骤,得到:
( H2 M1
m1
)
( H2
M2
m2
)
( H2
M3
m3
)
.....................
3
. 球面宇宙参数的计算
1)在上面的计算中已经求出 a
,于是用公式(8)算出 a
0
;
2)有了 a
0
,就可用公式(7)算出 L
0
;
3)把 a
0
和 L
0
代入公式(10)算出 L
;
4)用公式(15)计算 Hmax
;
5)用公式(17)计算 r
;
6)计算 Zmax
要用球外公式,参见附录8,与 Zmax
对应的就是 rs
。
最外层的计算结果见表2(附录12);a
,L
0
和 L
的计算结果列入附录15中。