|
楼主 |
发表于 2011-12-26 13:56:09
|
显示全部楼层
现在我假设我们这个宇宙周围有六个宇宙(这个假定可能根本不成立!),计算一下引力,(远处的更多宇宙(如果假设都在同一个平面均匀分布)引力情况应该一样(增加或减少)),
前面(大前天)写的‘程序一’改为:
程序1A:
double precision a,s,x,ss
! 单位圆(半径或许可能是一千亿光年)上均匀分布六个单位质量的物体
write(*,*) ' 序号 离开中心的距离 引力的和'
m=100
do ia=0,50
a=ia/(m+0.0d0) ! 距离中心的距离
n=3
s=0
do i=1-n,n
x=cos((n-i)*3.1415926535897932384626d0/3.0d0)
! y**2=1-x**2, 六个单位质量的物体之一的坐标是x,y
! r=sqrt(a*a+1+2*a*x) ! =sqrt((a+x)**2+y**2)
! cos((seta)=(-a-x)/r
ss=(-a-x)/sqrt((a*a+1+2*a*x)**3) ! =(1/r**2)*cos((seta)
s=s+ss ! 引力主要看增加减少,所以没有考虑引力的单位
!write(*,'(i5,3f10.5)') i,x,s,ss
enddo
write(*,'(i6,f9.2,f16.7)') ia,a,s
enddo
end
这个Fortran程序运行的结果:
序号 离开中心的距离 引力的和
0 0.00 0.0000000
1 0.01 0.0300034
2 0.02 0.0600271
3 0.03 0.0900916
4 0.04 0.1202180
5 0.05 0.1504281
6 0.06 0.1807444
7 0.07 0.2111909
8 0.08 0.2417930
9 0.09 0.2725777
10 0.10 0.3035740
11 0.11 0.3348130
12 0.12 0.3663286
13 0.13 0.3981572
14 0.14 0.4303383
15 0.15 0.4629148
16 0.16 0.4959333
17 0.17 0.5294443
18 0.18 0.5635026
19 0.19 0.5981678
20 0.20 0.6335043
21 0.21 0.6695821
22 0.22 0.7064770
23 0.23 0.7442709
24 0.24 0.7830527
25 0.25 0.8229182
26 0.26 0.8639711
27 0.27 0.9063233
28 0.28 0.9500957
29 0.29 0.9954186
30 0.30 1.0424325
31 0.31 1.0912892
32 0.32 1.1421520
33 0.33 1.1951972
34 0.34 1.2506147
35 0.35 1.3086095
36 0.36 1.3694025
37 0.37 1.4332322
38 0.38 1.5003559
39 0.39 1.5710518
40 0.40 1.6456202
41 0.41 1.7243865
42 0.42 1.8077027
43 0.43 1.8959506
44 0.44 1.9895444
45 0.45 2.0889345
46 0.46 2.1946111
47 0.47 2.3071082
48 0.48 2.4270093
49 0.49 2.5549524
50 0.50 2.6916368
这个计算结果说明:离开中心的距离越远,周围六个宇宙引力的和越大。也就是说我们的宇宙越膨胀则速度越快。
以上只是个朴素计算,假定了有很多宇宙按普通方式存在,这个假定可能根本不成立! |
|