（2） 四个粒子的运动
其中有两个p粒子（红色），两个n粒子（黑色）。由于初始位置的对称，在开始一段时间四个粒子的运动几乎各自表现为简谐振动，但后来由于求解微分方程不断累积误差，使得各个粒子的轨迹不断偏离理论值，最后呈现出杂乱无章的状态。
（3） 六个粒子的运动
其中有两个p粒子（红色），四个n粒子（黑色）。六个粒子的运动状态最初很有规律，但随着时间增加，运动轨迹的复杂度急速增加，最后变成一团毛线。
5 结语
该文档中的程序只处理了平面中的运动情况，而且设计并非完美，只希望起到抛砖引玉的作用，希望读者能根据该内容，做出更加完美的多体运动的matlab动画演示，得到更加准确的数值解，或者能在更大的时间范围内求解数值解。
更多内容见文库内本人的其他matlab 动画演示。