Consider the following two mass-spring-damper system: The equations of motions for the system shown in Figure 1 are: m_{1} \bar{x}_{1}+c_{1} \dot{x}_{1}+k_{1} x_{1}-c_{1} \dot{x}_{2}-k_{1} x_{2}=f m_{2} \bar{x}_{2}+\left(c_{1}+c_{2}\right) \dot{x}_{2}+\left(k_{1}+k_{2}\right) x_{2}-c_{1} \dot{x}_{1}-k_{1} x_{1}=0

a) Implement the system of equations above in Simulink using the following parameters: m1 = 10; m2 = 100; c1 = 100;% I c2 = 1000; k1 = 1e4; k2 = 1e5; Define the model parameters in a separate .m file and use the ode45 Solver inside of Simulink. Make sure to decrease the maximum step size if the plots are not smooth. b) Simulate the response of the system assuming that f(t) is a step function of magnitude 5N. Plot the response of the systems (the two positions x1(t) and x2(t)) in two separate figures. c) Simulate the response of the system assuming that f (t) is a sinusoidal function:flt)=3 sin (10t). Plot the response of the systems (the two positions x1(t) and x,(t)) in two separate figures.