import h_transport_materials as htm

import matplotlib.pyplot as plt

u = htm.ureg

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

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

htm.plotting.plot(D_1, label="diffusivity 1")
htm.plotting.plot(D_2, label="diffusivity 2")

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