[sigmoid10]

Mass scales like density times radius cubed. The first M would technically be a function of position. So you would integrate over slices of varying M not just because of varying density.