import numpy as np
import matplotlib.pyplot as plt
from hyperfine.superconductivity import pippard

T = np.linspace(0.0, 1.0, 100)
args = (1.0, 1.43e-3, 200.0, 50.0)
xi = np.array([pippard.xi_Pippard(tt, *args) for tt in T])
plt.plot(T, xi, "-")
plt.xlabel("$T / T_{c}$")
plt.ylabel(r"$\xi_{0}(T)$ (nm)")
plt.show()