矢量化蛮力重力N体模拟 我们利用numpy的广播功能来矢量化简单的重力蛮力N体模拟。 然后,我们在天体物理系统的研究中采用这种方法。 代码应该简单易读。 下面的所有应用程序都是代码的应用程序,用于为给定的文章生成数据。 init_cond中的所有内容都会生成初始条件的样本,并且集成器会以给定的时间间隔和给定的时间步长对系统进行集成。