Can someone help me solve this with clear steps?

One$-$Mass model of laryngeal airflow

Consider the $1$ mass model of laryngeal airflow shown above.

The system of equations for the vocal fold vibration, airflow velocity and pressure are given by

$\frac{dx}{dt}=v,\frac{dv}{dt}=P*L(\frac{(d)}{(m)})-{}_0^{2}x-(\frac{b}{(m)})v$

$u_g=\frac{2}{}\sqrt[\phantom{2}]{(\frac{3(L)}{x^2}{)}^2+P_s}-\frac{2}{}(\frac{3(L)}{x^2})$

$P=\frac{P_s-\frac{3}{4}{u}_g^2}{2}$

where x and v are the vocal fold displacement and velocity, respectively, $u_g$ is the airflow velocity, P is the glottal pressure.

Parameters are given by $P_s=1$kPa,$L=2cm,=2e-5$Pasec,$=1\mathrm{.}2k\frac{g}{m^3},f_0=100Hz,{}_0=2f_0,(\frac{d}{m})=10,$ and $(\frac{b}{m})=0\mathrm{.}1.$

a$.$ Find the numerical solution of the system using the RK$4$ method. Recommended time step size $dt=1e-4($ sec $).$ Use the initial condition, $x=0,v=0,u_g=0,P=0\mathrm{.}5P_s,$ and the contact model; $($if $x0$ then $x=0,u_g=0).$ Plot time profiles of $x$ and $ug.$

b$.$ Find the frequency of glottal airflow fluctuation. Is the frequency same with the natural frequency of the vocal fold $(f_0)?$ If not, discuss the possible reasons.