Coefficients

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

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

Maximum fit

a= 11.525870

b= 28.642040

c= 0.36249083

d= 11.591293

e= 0.003646159

f= 2.7521489

Average fit

a= 12.808760

b= 29.327155

c= 0.30927066

d= 11.178678

e= 0.0026380049

f= 2.6686171

Minimum fit

a= 23.276952

b= 33.970953

c= 0.41030425

d= 11.915613

e= 0.001870912

f= 2.5824408