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