我们提出了一种计算框架，使用高阶隐式龙格-库塔物理信息神经网络（IRK-PINNs）方案来获取非线性耦合微分方程组的多维相空间解。基于最初解决依赖于坐标的场的微分方程的基础工作 [J. Comput. Phys. 378, 686 (2019)]，我们调整了该方案以适应将坐标视为函数的上下文。这种修改使我们能够有效地求解粒子在外部场中的运动方程。我们的方案对于显式时间无关和周期性场特别有用。我们将此方法应用于成功求解放置在中心力场中的质量粒子和周期性电场中的带电粒子的运动方程。