I keep on getting the error that I'm missing a positional argument even though its been clearly defined

Question

re = np.logspace(3,5,20)        
epsD = 0.001
fguess = 0.0001
def F(f,re,epsD):
    return 1/np.sqrt(f)+2*np.log10(epsD/3.7 + 2.51/re/np.sqrt(f))
f=np.empty(len(re))