% 整定時間が最小となる粘性の値を求める

binit = 5;
[bopt,topt] = fminsearch(@damped_vibration_setting_time,binit);
fprintf("bの初期値 %f  bの最適値 %f  最小整定時間 %f\n", binit, bopt, topt);
