Please write code in PYTHON to solve the following question,

$($

$!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!$ NO USE OF INBUILT FUNCTIONS $!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!$

$)$

$,$

$$with explanations to steps. The code needs to be able to run in PYTHON language. no other languages please!! Please don't let my question go to waste...

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$)$;