We compare trans2E3 (2 E3) with a simple transmission.

from exojax.spec.rtransfer import trans2E3

def simple_trans(dtau, mu):
    return jnp.exp(-dtau/mu)

import jax.numpy as jnp
dtau_arr=jnp.logspace(-3,1,100)

import matplotlib.pyplot as plt
plt.plot(dtau_arr,trans2E3(dtau_arr),color="black",label="$\\mathcal{T}= 2E_3 = 2 \\int_0^1 d \\mu \\, \\mu \\, e^{-\\Delta \\tau/\\mu}$")
plt.plot(dtau_arr,simple_trans(dtau_arr,1),ls="dotted",color="gray", label="$\\mathcal{T}= e^{-\\Delta \\tau/\\mu} \\, \\, (\\mu=1)$")
plt.plot(dtau_arr,simple_trans(dtau_arr,2.0/3.0),ls="dashed",color="gray", label="$\\mathcal{T}= e^{-\\Delta \\tau/\\mu} \\, \\, (\\mu=2/3)$")
plt.plot(dtau_arr,simple_trans(dtau_arr,0.3),ls="dashdot",color="gray", label="$\\mathcal{T}= e^{-\\Delta \\tau/\\mu} \\, \\, (\\mu=0.3)$")
#plt.yscale("log")
plt.legend()
plt.xscale("log")
plt.xlabel("$\\Delta \\tau$")
plt.ylabel("transmission")
plt.savefig("transrt.png")
plt.savefig("transrt.pdf")
plt.show()