def mass_integral(self, x, axis_ratio):
r = x * axis_ratio<EOL>return <NUM_LIT:2> * np.pi * r * self.convergence_func(x)<EOL>
Routine to integrate an elliptical light profiles - set axis ratio to 1 to compute the luminosity within a \
        circle