Untitled

# ПС
def CriticalReynoldsNumber(Tr, M):
    Tst_otn = Tst / Tr  # Относительная температура стенки
    x = (Tst_otn - 1) / M ** 2
    y = 1 - 16 * x - 412.5 * x**2 - 35000 * x**3 - 375000 * x**4
    return criticalReynoldsNumber * y

# Критические числа Рейнольдса для 1 и 2 грани
Re_kr1 = CriticalReynoldsNumber(Trl1, M1)
Re_kr2 = CriticalReynoldsNumber(Trl2, M2)

# Координаты точек перехода
x_kr1 = Re_kr1 * mu1 / (V1 * ro1)
x_kr2 = Re_kr2 * mu2 / (V2 * ro2)