四届龙格库塔法求解线性、非线性单自由度振动系统微分方程
1