асчёт p̅к, Mк, Cx

#\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\ Расчёт p̅к, Mк, Cx \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\

T_0 = T * (1 + ((k - 1) / 2) * M ** 2)
T_k = T_0 * (1 - ((k - 1) * V_k ** 2 / (2 * k * R * T_0)))
M_k = V_k / np.sqrt(k * R * T_k)
p0_shtrih = p * ((k+1)/(2*k*M**2*(np.sin(θ_c))**2-(k-1)))**(1/(k-1)) * ((k+1)*M**2*(np.sin(θ_c))**2*((k-1)*M**2+2)/(2 * ((k-1)*M**2*(np.sin(θ_c))**2 + 2)))**(k/(k-1)) 
P_k = p0_shtrih * (1 + ((k - 1) * M_k ** 2)/ 2) ** ((- k)/(k - 1))
p̅k = (P_k - p) * 2 / (ro * V ** 2)
#print("T_0",T_0)
#print("T",T)
#print("k",k)
#print("M",M)
#print("T_k",T_k)
#print("M_k",M_k)
#print("P_k",P_k)
print("Коэффициент давления p̅k: " + str(round(p̅k ,3)) )

#\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\ Расчёт по приближенным зависимостям \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\

# Погрешность 5%
θ_c_emp = np.arcsin(np.sqrt((1 / M ** 2) + (k + 1) * (np.sin(β_k) ** 2) / 2))
p̅k_5_proc = 0.002 * (0.8 + (1 / M ** 2)) * np.degrees(β_k) ** 1.7
# Погрешность 3%

print("Угол θ приближенный: " + str(round(θ_c_emp, 3)) + ' рад; ' + str(round(np.degrees(θ_c_emp), 3)) + ' град')
print("Коэффициент давления p̅k (5%): " + str(round(p̅k_5_proc ,10)) )