Colebrook-White Equation

MATLAB code to compute the friction factor in pipes for given values of the Reynolds number (Re) and the relative roughness coefficient (epsilon).

Syntax:
f = colebrook(Re,epsilon)

Example 1: Single Re, single epsilon
Re = 1e5;
epsilon = 1e-4;
f = colebrook(Re,epsilon)

Example 2: Multiple Re, single epsilon
Re = 5000:1000:100000;
epsilon = 1e-4;
f = colebrook(Re,epsilon);
plot(Re,f)

Example 3: Single Re, multiple epsilon
Re = 1e5;
epsilon = linspace(1e-4,1e-1,100);
f = colebrook(Re,epsilon);
plot(epsilon,f)

Example 4: Multiple Re, multiple epsilon
Re = logspace(4,8,100);
epsilon = linspace(1e-4,1e-1,100);
[RE,EPSILON] = meshgrid(Re,epsilon);
F = colebrook(RE,EPSILON);
surf(RE,EPSILON,F)

References:
[1] Colebrook, C. F., & White, C. M. (1937). Experiments with fluid
friction in roughened pipes. Proceedings of the Royal Society of
London. Series A - Mathematical and Physical Sciences, 161(906),
367-381.
[2] Colebrook, C. (1939). Turbulent Flow in Pipes, with Particular
Reference to the Transition Region between the Smooth and Rough
Pipe Laws. Journal of the Institution of Civil Engineers, 11(4),
133-156.