Cyklus byl chyba. Odstranil jsem přičítání konstanty při výpočtu horní meze upper_bound a podmínil jsem volání toms748_solve nezápornou a nenulovou total_mass.