Untitled

def findDefiningTemperature (lam_or_turb: float,T: float, M: float) -> float: # Функция поика температуры восстановления и определяющей температуры
    T_z = 10**6
    T_z_previous = 0
    rl = 0.83
    rt = 0.88
    stepen = 0.5
    r = rl
    if (lam_or_turb == 1):
        stepen = 1/3
        r = rt
    while abs(T_z - T_z_previous) > (T_z * 0.01):
        T_z_previous = T_z
        Tr = T * (1 + r * (k - 1) * 0.5 * M ** 2)
        T_z = 0.5 * (Tst + T) + 0.22 * (Tr - T)
        Cpi = SpecificHeat_Cp(T_z)
        mui = DynamicViscosity_mu(T_z)
        Lambdai = CoeffThermalConductivity_Lambda(T_z)
        Prl = Cpi * mui / Lambdai
        r = Prl ** stepen
    return Tr,T_z