Coefficients

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

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

Maximum fit

a= 8.1107566

b= 24.770667

c= 0.12386137

d= 8.8783877

e= 0.018050805

f= 4.1693711

Average fit

a= 7.0401813

b= 24.577236

c= 0.09226137

d= 6.8409626

e= 0.0014101056

f= 2.8620123

Minimum fit

a= 6.8154628

b= 24.840721

c= 0.08961893

d= 6.9984775

e= 0.0000945208

f= 1.3446629