Write a Python code.

Write a python code for the problem below.

For turbulent flow, the Colebrook equation can be used to calculate the friction factor,

$\frac{1}{\sqrt[\phantom{2}]{f}}=-2log(\frac{}{3\mathrm{.}7D}+\frac{2\mathrm{.}51}{Re*\sqrt[\phantom{2}]{f}})$

Where $f$ is the friction factor, $$ is the relative roughness, $D$ is the pipe diameter, and Re is the Reynolds number which is defined as:

$Re=\frac{*V*D}{}$

Use newton's method to find the friction factor when density $()$ is $3\mathrm{.}23k\frac{g}{m^3},$ velocity $(V)$ is $20\frac{m}{s},=0\mathrm{.}0015mm,D=0\mathrm{.}001m,$ and the viscosity $()$ is $1\mathrm{.}79**{10}^{-5}{N}^{**}\frac{s}{m^2}.$ Make sure your solution converges to $5$ decimal places.