用于求解简谐激励下单自由度系统稳态响应的fortran程序,采用机械阻抗法
1