---

import math

# ==========================================
# 1. ВВОДНЫЕ ДАННЫЕ (Подставлены из Варианта 1)
# ==========================================
H = 7260        # Высота полета, м
Mn = 0.9        # Число Маха (в задании 9.1, но это опечатка, для высоты 7260 м это сверхзвук)
pi_k = 9.1      # Степень повышения давления в компрессоре
T_g = 1133      # Температура газа перед турбиной, К

# Коэффициенты
eta_k = 0.86    # КПД компрессора
sigma_ks = 0.95 # Коэфф. восстановления в КС
eta_g = 0.99    # Полнота сгорания
lambda_v = 0.5  # Приведенная скорость на входе
lambda_k = 0.3  # Приведенная скорость на выходе

# Константы
R = 287
k = 1.4
cp = 1005
Hu = 42900 * 1000
T0 = 293

# ==========================================
# РАСЧЕТ
# ==========================================
# Атмосфера
if H <= 11000:
    T_h = 288.15 - 0.0065 * H
    p_h = 101325 * (T_h / 288.15) ** 5.25588
else:
    T_h = 216.65
    p_h = 22632 * math.exp(-0.0001577 * (H - 11000))

a_h = math.sqrt(k * R * T_h)
V_n = Mn * a_h

# Вход
sigma_vh = 0.97 if Mn <= 1 else 0.97 - 0.1 * (Mn - 1)**1.5
p_v = p_h * (1 + 0.2 * Mn**2)**3.5 * sigma_vh
T_v = T_h * (1 + 0.2 * Mn**2)
a_kp_v = 18.3 * math.sqrt(T_v)
c_v = a_kp_v * lambda_v
T_v_st = T_v - (c_v**2) / (2 * cp)
p_v_st = p_v * (T_v_st / T_v)**(k / (k - 1))

# Компрессор
p_k = p_v * pi_k
T_k = T_v * (1 + (pi_k**((k - 1) / k) - 1) / eta_k)
a_kp_k = 18.3 * math.sqrt(T_k)
c_k = a_kp_k * lambda_k
T_k_st = T_k - (c_k**2) / (2 * cp)
p_k_st = p_k * (T_k_st / T_k)**(k / (k - 1))

# Камера сгорания
p_g = p_k * sigma_ks
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
q_t = (cp_Tg - cp_Tk) / (Hu * eta_g - cpn_Tg + cp * T0)

# ==========================================
# ВЫВОД
# ==========================================
print("="*40)
print("РЕЗУЛЬТАТЫ РАСЧЕТА (Вариант 1)")
print("="*40)
print(f"V_n = {V_n:.2f} м/с")
print(f"p_в* = {p_v:.2f} Па, T_в* = {T_v:.2f} К")
print(f"c_в = {c_v:.2f} м/с")
print(f"T_в = {T_v_st:.2f} К, p_в = {p_v_st:.2f} Па")
print("-" * 40)
print(f"p_к* = {p_k:.2f} Па, T_к* = {T_k:.2f} К")
print(f"c_к = {c_k:.2f} м/с")
print(f"T_к = {T_k_st:.2f} К, p_к = {p_k_st:.2f} Па")
print("-" * 40)
print(f"p_г* = {p_g:.2f} Па")
print(f"q_т = {q_t:.5f}")
print("="*40)