# Path: proptools/isentropic.py
# def stag_temperature_ratio(M, gamma):
# def stag_pressure_ratio(M, gamma):
# def stag_density_ratio(M, gamma):
# def velocity(v_1, p_1, T_1, p_2, gamma, m_molar):    # pylint: disable=too-many-arguments


import numpy as np
from matplotlib import pyplot as plt
from proptools import isentropic
"""Plot isentropic relations."""

M = np.logspace(-1, 1)
gamma = 1.2

plt.loglog(M, 1 / isentropic.stag_temperature_ratio(M, gamma), label='$T / T_0$')

# Path: proptools/isentropic.py
# def stag_temperature_ratio(M, gamma):
# def stag_pressure_ratio(M, gamma):
# def stag_density_ratio(M, gamma):
# def velocity(v_1, p_1, T_1, p_2, gamma, m_molar):    # pylint: disable=too-many-arguments


# Path: docs/source/examples/plots/isentropic_relations.py
import numpy as np
from matplotlib import pyplot as plt
from proptools import isentropic

"""Plot isentropic relations."""

M = np.logspace(-1, 1)
gamma = 1.2

plt.loglog(M, 1 / isentropic.stag_temperature_ratio(M, gamma), label='$T / T_0$')

plt.loglog(M, 1 / isentropic.stag_pressure_ratio(M, gamma), label='$p / p_0$')