---

T_g, T_k, eta_g, cp, Hu, T0 = 1133, 600, 0.99, 1005, 42900000, 293
q_T = ((4e-5 * T_g**2 + 1.0795 * T_g - 59.843) - (4e-5 * T_k**2 + 1.0795 * T_k - 59.843)) / (Hu * eta_g - (0.0004 * T_g**2 + 2.1436 * T_g - 281.82) + cp * T0)

print(f"Относительный расход топлива q_T = {q_T:.5f}")
print(f"Удельная тяга P_уд = {(1 + q_T * (1 - 0.05)) * 600 - 250:.2f} Н·с/кг")