Please do this math, if possible do this math by hand.

$m_1{m}_6=2\mathrm{.}5$;$k_1{k}_3=3200,k_4=2400,k_5=2000$ and $k_6=1600$; Assume ${}_r=0\mathrm{.}02$ for all modes

a$).$ Determine the mass, damping, and stiffness matrices and loading factor vector.

b$).$ Use 'eig' command determine mode shape and ${}_n$

c$).$ Determine EOM in modal domain

d$).$

${}_r(t)={\displaystyle \underset{0}{\overset{t}{}}}p()h(t-)d$ for $t=T_p$

Where, $h(t)=\frac{1}{M_r{}_d}{e}^{-{}_{n}t}sin({}_{d}t)$

$u(t)=**$

Obtain the solution using Matlab$/$Simulink$.$