The function returns 4f(x2 - x1)/D, and not f(x2 - x1)/D. If the latter is used as a fanno parameter anywhere else, it is probably wrong.