% Generic 2nd order system
dur=5;
%b=3; c=43;             % existing system
b=4.9; c=12.25;         % desired system

num=[c];
den=[1 b c];
gen=tf(num,den);
step(gen,dur); axis([0 dur 0 1.5]); grid on;

ome_n=sqrt(c)
zeta=b/(2*sqrt(c))
sigma=zeta*ome_n
ome_d=ome_n*sqrt(1-zeta^2)
phi=atan2(ome_d,sigma)

rise=1.8/ome_n
PO=exp(-zeta*pi/sqrt(1-zeta^2))
settling=4.6/(zeta*ome_n)