附录六  球面的计算过程


    这是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中。