import h_transport_materials as htm
from h_transport_materials.plotting import plot
import matplotlib.pyplot as plt

D_1 = htm.Diffusivity(
   D_0=1 * htm.ureg.m**2 * htm.ureg.s**-1,
   E_D=0.2* htm.ureg.eV * htm.ureg.particle**-1,
)

D_2 = htm.Diffusivity(
   D_0=2 * htm.ureg.m**2 * htm.ureg.s**-1,
   E_D=0.3* htm.ureg.eV * htm.ureg.particle**-1,
)

plot(D_1, linestyle="-.", color="tab:red", linewidth=3, label="D1")
plot(D_2, linestyle="dashed", color="tab:blue", linewidth=2, label="D2")

plt.yscale("log")
plt.legend()
plt.show()