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(htm.PropertiesGroup([D_1, D_2]))

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