function dotq = damped_vibration_global (t, q)
% 微分方程式（大域変数でパラメータの値を受け渡す）
    global b;
    x = q(1);
    v = q(2);
    dotx = v;
    dotv = -b*v - 4*x;
    dotq = [dotx; dotv];
end
