function dotq = beam_supported_dynamic_equation_param (t,q, n,M,B,K,force)
% ƒr�[ƒ€‚Ì‰^“®•û’öŽ®
    un = q(1:n);
    vn = q(n+1:2*n);
    
    mat = M;
    vec = -K*un-B*vn + force(t,n,un,vn);
    sol = mat \ vec;
    dotvn = sol(1:n);
    
    dotq = [ vn; dotvn ];
end
