Please write code in PYTHON to solve the following question,

ok to use numpy, scipy etc

with explanations to steps. The code needs to be able to run AND solved in PYTHON language. no other languages please!!

assuming:

$q=10kg\frac{f}{m}m$

$$ Titanium alloy $E=100$GPa;

Second area moment $J=1440m{m}^4.$ Which is equal to rectangle with

sides $10mm$ and $12mm,$ and force acts along longer side;

Circular freq. $=\frac{0\mathrm{.}11}{s}$;

$$ Length $1=100mm$;

$$ Rigid constraint on left bound $y(0)=0,y^{\text{'}}(0)=0$;

$$ Bending force on the left bound $y^{\text{'}\text{'}\text{'}}(0)=\frac{q*l}{E*J}$;

$$ Bending moment on the left bound $y^{\text{'}\text{'}}(0)=-\frac{q*l^2}{2*E*J}$;

Solve this equation assuming:

q$=10$ kgf$/$mm;

Titanium alloy E$=100$ GPa;

Second area moment J$=1440$ mm$^4.$ Which is equal to rectangle with sides $10$ mm and $12$ mm$,$ and force acts along longer side;

Circular freq. $\backslash$omega $=0\mathrm{.}11/$s;

Length l$=100$ mm;

Rigid constraint on left bound y$(0)=0,$ y$(0)=0$;

Bending force on the left bound y$\text{'}\text{'}\text{'}(0)=($q$$l$)/($E$$J$)$;

Bending moment on the left bound:

y$\text{'}\text{'}(0)=-($q$($l$^2))/(2$E$$J$)$;