function ts = damped_vibration_setting_time ( b )
% Œ¸�Š�U“®‚Ì�®’èŽžŠÔ
    xinit = 2.00;
    ratio = 0.01;
    interval = 0:0.001:10;
    qinit = [xinit;0.00];
    vibration = @(t,q) damped_vibration_param (t,q,b);
    [time,q] = ode45(vibration,interval,qinit);
    plot(time,q(:,1));
    [ts,ks] = setting_time( time, q(:,1), 0.00, ratio*xinit);
end
