how to find transfer function of a second order system using matlab commands? can anyone help me with it ?

Consider the system shown in following figure, where damping ratio is 0.6 and natural undamped frequency is 5 rad/sec. Obtain the rise time tr, peak time tp, maximum overshoot Mp, and settling time 2% and 5% criterion ts when the system is subjected to a unit-step input. What are the commands to introduce num and den , since i get an error if i use num = [wn^2] den = [s^2+2*zeta*wn*s] sys = tf(num, den) and how to use commands to find tr, ts, mp and to plot in graph

[Image Analyst]

What figure? You didn't insert or attach anything.