h=stepplot(sys_cl,t);
s_info=stepinfo(sys_cl);
grid
title('Step Response with Text')
h.showCharacteristic('PeakResponse')
h.showCharacteristic('SettlingTime')
h.showCharacteristic('SteadyState')
text(s_info.PeakTime+0.05,h.Responses.Data.Amplitude(end)*0.6,sprintf('PeakTime (s) = %.3f', s_info.PeakTime))
text(s_info.PeakTime+0.05,s_info.Peak,sprintf('Overshoot (%%) = %.3f', s_info.Overshoot))
text(s_info.SettlingTime+0.05,h.Responses.Data.Amplitude(end)*0.4,sprintf('SettlingTime (s) = %.3f', s_info.SettlingTime))
text(h.Responses.Data.Time(end)*0.5,h.Responses.Data.Amplitude(end)*1.05,sprintf('Final Value (s) = %.3f', h.Responses.Data.Amplitude(end)))