MATLAB code to compute the Darcy-Weisbach friction factor for turbulent flow in pipes. It includes a subroutine to estimate the relative roughness coefficient from values of the friction factor and the Reynolds number measured in different operating points.