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([