第二章 再次模拟（第3页）

大脑高速运转，大量与数学，计算机相关的知识融会贯通，被沈程理解，固化，想要编写一个程序，相应代码组合，算法应用，随之并现出来，感受非常新奇。

早上。

沉睡中沈程被舍友们叫起，睡意惺惺，满脸的困倦。

“老沈，快点，上课要迟到了！”刘洪涛催促道。

“早上第一节是‘严老虎’的课，你慢慢磨蹭，我们先走啦！”吴志坤提醒了沈程一句，拿着课本先走了。

刘洪涛也不管沈程了，跟着吴志坤，李沐辰，先一步离开宿舍。

见舍友们丢下他，沈程不由急了，脸也没空洗了，穿上鞋，拿了课本，小跑追上舍他们。

早上几节课，沈程都处于走神状态，思考应对小行星撞击地球的事。

自建末日生存堡垒肯定保不了命，将小行星撞击地球的事提前透露给相关部门，与国家合作，想办法阻止小行星的撞击才是合适的方向。

想到这里，沈程霍然醒悟，随即从桌底拿出圆珠笔，草稿纸，按照脑海中的数据记忆，迅速展开推导计算。

地球围绕太阳转动的轨道大体呈偏心率很小的椭圆，半长轴（a）1.496×108千米，半短轴（b）1.4958×108千米；半焦距（c）25×105千米；周长（l）9.4×108千米。

椭圆轨道偏心率（e）和扁率（f）分别为：160和17000……，

轨道倾角在：28.58°与18.28°之间变化

升交点赤经：125.08°

近地点辐角：318.15°

……

今年2027年3月22日17点34分，太阳光线直射赤道时是北半球的春分日，……，末日陨石撞击南美洲，西半球，截撞击地球，介于：西经35-74度，北纬5度-南纬20度之间。

由于在太阳的有心力场中，引力作用在行星上的力矩为0，因而行星的角动量是守恒，……。

设地球在t=T时的r值，误差为e，t=0时的初值r0，a0，v0为已知，将0到T的时间间隔划分为n个dt，即令计算步长dt=Tn，……，步长dt从t=0时的初值点推算到t=n·dt=T时的r值……。

末日小行星侧面撞击地球，小行星轨道与地球轨道重叠，……。

根据奥伯斯方法和高斯方法，小行星轨道接近抛物线﹐计算其轨道将其作为抛物线处理，结合天体在抛物线轨道上时刻向径和弦关系的欧拉方程，……。

……。

“沙沙……。”沈程手上的水性笔在草稿纸上飞快的列写公式。