说明：四阶龙格库塔算法，用来求解二阶偏微分方程，在转子系统的动力学方程数值求解中普遍应用

四阶龙格库塔