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

prop1 = htm.Diffusivity(1, 0.1)
prop2 = htm.Diffusivity(2, 0.2)

group = htm.PropertiesGroup([prop1, prop2])
mean_property = group.mean()

plot(group, alpha=0.5)
plot(mean_property, color="black")

plt.annotate("mean property", (0.0025, 2e-2))
plt.yscale("log")
plt.xlabel("1/T (K$^{-1}$)")
plt.show()