I try to find the friction factor for Moody diagram (f) by this formula using sympy
from sympy import * f = Symbol (\'f\') Re = 3000 relative_roughness=0.1 so