Claiming the smallest remaining mass is best was misleading.
Plot rho as well.
Print final L and rho as well.
Also check if M <= 0.