Newmark-beta法解微分方程的MATLAB程序。 function [q,v,a]=newmark(M,C,K,F,q0,v0,dt,nt) % newmark-beta method % [q,v,a]=newmark(M,C,K,F,q0,v0,dt,nt) % obtain the response of the dynamic system % M - mass matrix % K -