---

import math

def calculate_q_T(T_g, T_k, eta_g, cp, Hu, T0, delta_otb, c_c, V_n):
    # 1. Расчет теплоемкостей по формулам с картинки
    cp_Tg = 4e-5 * T_g**2 + 1.0795 * T_g - 59.843
    cpn_Tg = 0.0004 * T_g**2 + 2.1436 * T_g - 281.82
    cp_Tk = 4e-5 * T_k**2 + 1.0795 * T_k - 59.843
    
    # 2. Расчет q_T (формула Ильичева)
    q_T = (cp_Tg - cp_Tk) / (Hu * eta_g - cpn_Tg + cp * T0)
    
    # 3. Расчет удельной тяги (твоя формула)
    P_ud = (1 + q_T * (1 - delta_otb)) * c_c - V_n
    
    # Вывод прямо из функции
    print(f"q_T = {q_T:.5f}")
    print(f"P_уд = {P_ud:.2f} Н·с/кг")
    
    return q_T, P_ud

# ===== ЗАПУСК (СЮДА ПИШИ СВОИ ЗНАЧЕНИЯ) =====
calculate_q_T(
    T_g=1133,         # Температура газа перед турбиной
    T_k=600,          # Температура за компрессором
    eta_g=0.99,       # Полнота сгорания
    cp=1005,          # Теплоемкость воздуха
    Hu=42900000,      # Теплота сгорания керосина
    T0=293,           # Стандартная температура
    delta_otb=0.05,   # Отбор воздуха
    c_c=600,          # Скорость истечения
    V_n=250           # Скорость полета
)