function [ ff ] = colebrook( f,Re,roughness )

ff = 1.14-2*log10(roughness+9.35/(Re*sqrt(f)))-1/sqrt(f);

end

