说明：变步长龙格库塔法，用于解决四阶微分方程组求解析解的问题，尤其适用于两自由度间隙碰撞问题的解决。

间隙碰撞 间隙-MATLAB 间隙-碰撞 龙格库塔-间隙 碰撞微分方程