Coefficients

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

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

Maximum fit

a= 6.7996021

b= 23.063099

c= 0.05884421

d= 5.9419158

e= 0.0029401137

f= 1.8816427

Average fit

a= 7.2873359

b= 24.837850

c= 0.09743608

d= 6.9263302

e= 0.0010123872

f= 2.5846868

Minimum fit

a= 5.2797486

b= 23.005729

c= 0.07580663

d= 6.8415790

e= -1.4407086

f= 322.04722