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)
plt.plot(T, pippard.j_0_t(T, *args), "-")
plt.xlabel("$T / T_{c}$")
plt.ylabel("$J(0,T)$")
plt.show()