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

q = np.logspace(-4, 4, 200)
args = (0.0, 10.0, 1.43e-3, 30.0, 300.0, 40.0)
k = np.array([pippard.K_Pippard(qq, *args) for qq in q])
plt.plot(q, k, "-")
plt.xlabel("$q$ (nm$^{-1}$)")
plt.ylabel(r"$K_{\mathrm{Pippard}}(q)$ (nm$^{-2}$)")
plt.xscale("log")
plt.yscale("log")
plt.show()