function f = friction_coef(flow,Re)
