上次采用过Mathematica自带的微分方程求解器求解，虽然求解的精度非常的不错，但是求解的速度非常慢。为了提高地震波反应谱的计算速度，并且能够同时计算不同阻尼比下的多个反应谱，需要一种新的计算方法...
Continue reading...