import numpy as np
import matplotlib.pyplot as plt
import aero.Rayleigh as ray
npoints = 200
gam = 1.4
Mmin = 0.1
Mmax = 4.
Mach = np.log10(np.logspace(Mmin, Mmax, npoints + 1))
Tparam = ray.maxTiratio_Mach(Mach, gam)
Ts = ray.Ts_Tscri(Mach, gam)
Ti = ray.Ti_Ticri(Mach, gam)
Ps = ray.Ps_Pscri(Mach, gam)
Pi = ray.Pi_Picri(Mach, gam)
V = ray.V_Vcri(Mach, gam)
dS = ray.NormdS(Mach, gam)
fig = plt.figure(1, figsize=(10, 8))
fig.suptitle('Ratio to critical state, $\gamma = %.1f$' % gam,
fontsize=12,
y=0.93)
#plt.plot(Mach, Tparam, 'k--')
plt.plot(Mach, Ti, '-', color='#ff0000')
plt.plot(Mach, Ts, '-', color='#882222')
plt.plot(Mach, Pi, '-', color='#0000ff')
plt.plot(Mach, Ps, '-', color='#222288')
plt.plot(Mach, V, '-', color='#009999')
plt.legend([
