function [f] = frictionfactor(Re)
r = .00001;
%f = 0.3164./(Re.^.25);       
f = 4*(log10(r/3.7+(5.74./(Re.^.9)))).^2;
f = 1./f;
