---

# ==========================================
# 6. РАСЧЕТ ВРЕМЕНИ ПОЛНОГО ОСТЫВАНИЯ
# ==========================================
print("\n--- Анализ времени остывания ---")

# 1. Вывод характерного времени (tau) по методичке
a2_theory = 2.3e-5
tau = (R**2) / a2_theory
print(f"1. Характерное время (tau = R^2/a^2): {tau:.1f} сек")

# 2. Поиск времени, когда температура упадет до 1 градуса (условное "полное остывание")
# Для этого будем перебирать время с шагом в 1 секунду, пока температура не станет < 1.0 °C
t_cool = 0
while get_temperature(t_cool, a2_theory) > 1.0:
    t_cool += 1
print(f"2. Время остывания до 1°C (по аналитике): {t_cool} сек")

# 3. Время остывания для реального ЛОГОСа (с подогнанным a2)
t_cool_logos = 0
while get_temperature(t_cool_logos, a2_fitted) > 1.0:
    t_cool_logos += 1
print(f"3. Время остывания до 1°C (для условий ЛОГОСа): {t_cool_logos} сек")