please write code in Julia

Consider the Colebrook equation, which phenomenologically describes the Darcy friction factor for rough pipes for Reynolds number greater than 4000 . The equation has the implicit form fD1=2log10(3.7/D+RefD2.51) where Re is the Reynolds number, /D is the relative roughness, and fD is the Darcy friction factor. Using fixed point iteration to solve for fD, construct the Moody diagram for the turbulent region (Re from 4000 to 108 )