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

z = np.linspace(0.0, 200.0, 100)
args = (100.0, 10.0, 50.0, 190.0, 0.0)
plt.plot(z, london.screening_profile_film(z, *args), "-")
plt.xlabel("$z$ (nm)")
plt.ylabel("$B(z)$ (nm)")
plt.show()