困扰我很久的一个问题，关于线性振动方面的

对于一个线性结构振动问题，我想用龙格库塔求解，所以我就找了个简单算例试了一下，但是结果发散。我看了好久都不明白什么地方不对，是不是龙格库塔方法只能用在振型叠加法，不能用在直接积分法上面啊。下面是我做的简单算例，大家可以看看。