本篇将主要讲解程序的框架部分。
该程序主要分为三个部分，首先是初始化部分，主要为设置离子质荷比、初始位置、速度。
其次为求解轨迹部分，通过离子位置获取对应位置的电场，并经由空间电荷效应修改电场后，通过数值求解方法获取新的位置，最后在经由碰撞模型获取最终的离子位置。
最后为终止条件部分，主要检测离子是否脱离束缚区域，如果脱离束缚区域，记录条件并退出求解轨迹；否则，将继续运行求解轨迹，直至满足终止条件或超过最大运行时间。