i try to using 100 points and 5000 points using the same interpolation as you did in the tov program,but the m-r outcome is different.
however the nb-M relation is the same,nb-R relation is very different.
and there may be some problem in your interpolation?in:
den = neos[i-1]*np.exp(np.log(P/peos[i-1])
*np.log(neos[i-1]/neos[i])
/np.log(peos[i-1]/peos[i]));
shouldnt it be:
den = neos[i]*np.exp(np.log(P/peos[i])
*np.log(neos[i-1]/neos[i])
/np.log(peos[i-1]/peos[i]));
?
i try to using 100 points and 5000 points using the same interpolation as you did in the tov program,but the m-r outcome is different.
however the nb-M relation is the same,nb-R relation is very different.
and there may be some problem in your interpolation?in:
den = neos[i-1]*np.exp(np.log(P/peos[i-1])
*np.log(neos[i-1]/neos[i])
/np.log(peos[i-1]/peos[i]));
shouldnt it be:
den = neos[i]*np.exp(np.log(P/peos[i])
*np.log(neos[i-1]/neos[i])
/np.log(peos[i-1]/peos[i]));
?