s=%s //G=k/(s*(s+a)) and T=G/(1+G) Mr=1.04 omegaR=11.55 disp("for zeta>0.7 , Mr does not exist, so neglect higher value") zeta=0.6021 disp(zeta) omegaN=omegaR/(sqrt(1-2*(zeta)^2)) disp(omegaN,"omegaN = ") k=omegaN^2 disp(k,"k=") disp(2*zeta*omegaN," a = ") disp(omegaN*sqrt(1-2*zeta^2 + sqrt(2-4*zeta^2+4*zeta^4))," B.W. = ") disp(4/(zeta*omegaN)," Ts = ")