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