from math import log, pi M = 67 c = 3.7 rho = 1.038 K = 5.4e-3 Ty = 70 Tw = 100 T0 = 20.0 a = (M**(2/3) * c * rho**(1/3)) / (K * pi**2 * (4 * pi / 3)**(2/3)) b = log ( 0.76 * (T0-Tw) / (Ty-Tw) ) time = a * b minutes = int (time / 60) seconds = int (time % 60) print(f"Starting from {T0} degrees, it takes {minutes} minutes and {seconds} seconds to reach {Ty} degrees.")