Coefficients

Ne_ave = No*(a*(r)^-b + c*(r)^-d + e*(r)^-f)

No = 3.14e8 (Rsun/(G Msun mp) )

Maximum fit

a= 7.4522137

b= 23.861773

c= 0.08522234

d= 6.6028767

e= 0.0016693167

f= 2.0853350

Average fit

a= 7.2184916

b= 24.769254

c= 0.09614883

d= 6.9220274

e= 0.001174982

f= 2.5166538

Minimum fit

a= 5.1334363

b= 22.808642

c= 0.05404971

d= 6.1909661

e= 0.008304779

f= 7.0457379