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)
isp = np.array([pippard.integrand_specular_profile(qq, *args) for qq in q])
plt.plot(q, isp, "-")
plt.xlabel("$q$ (nm$^{-1}$)")
plt.ylabel("$I(q)$ (nm)")
plt.xscale("log")
plt.yscale("log")
plt.show()