$\%$ please solve this with ode$45,$ and explain how $.$ I need the matlab plots, please provide me plots screen shots and matlab codes.

clc;

clear;

clear all;

$\%$ System numerical values

m$\_$crank $=1$;

I$\_$crank $=0\mathrm{.}5$;

m$\_$rod $=2$;

m$\_$slider $=3$

l$\_$crank $=1$;

l$\_$rod $=2$;

g $=9\mathrm{.}81$;

$\%$ Symbolic equations

syms theta x

$\%$ Angulat velocity $=$ omega $/$ Slider velocity $=$ v$\_$slider

omega $=$ diff$($theta$)$;

v$\_$slider $=$ diff$($x$)$;

$\%$ C will be used in T$\_$rod

c $=$ sqrt$($l$\_$rod$^2-($l$\_$crank$^2*$sin$($theta$)^2))$;

$\%*************************************************************************$

$\%$ Kinetic energy formulas

$\%$ Crank kinetic energy:

T$\_$crank $=(1/2)*(1/12)*$m$\_$crank$*$l$\_$crank$^2*$omega$^2$;

$\%$ Connected rod kinetic energy:

T$\_$rod $=$ m$\_$rod$*$omega$^2*(1/6*($l$\_$rod$*$l$\_$crank$*$cos$($theta$)/$c$)^2+...$

$1/2*($l$\_$crank$*$sin$($theta$))^2+(1/2)*($l$\_$crank$^3/$c$)*$cos$($theta$)*$sin$($theta$)^2)$;

$\%$ Slider kinetic energy:

T$\_$slider $=0\mathrm{.}5*$m$\_$slider$*$v$\_$slider$^2$;

$\%*************************************************************************$

$\%$ Potantial energy formula

$\%$ Connected rod potential Energy

V$\_$rod $=0\mathrm{.}5*$m$\_$rod$*$g$*$l$\_$crank$*$sin$($theta$)$;

$\%$ Lagrange equation

T $=$ T$\_$crank $+$ T$\_$rod $+$ T$\_$slider;

V $=$ V$\_$rod;

L $=$ T $-$ V

$\%$ Lagrange diff

L$\_$theta $=$ diff$($L$,$ theta$)$

L$\_$theta$\_$dot $=$ diff$($L$\_$theta, omega$)$

Q$\_$theta $=$ sym$(\text{'}$tau$\text{'})$; $\%$ Applied torque to crank

eq$1=$ L$\_$theta$\_$dot $-$ L$\_$theta $==$ Q$\_$theta

$\%$ Initial conditions

theta$\_0=$ pi$/4$; $\%$Initial angle

theta$\_$dot$\_0=0$; $\%$ Initial angular velocity