import matplotlib.pyplot as plt
import modes as ms

wavelengths = np.linspace(1.3, 1.6, 10)
nsi = [ms.materials.si(w) for w in wavelengths]

plt.plot(wavelengths, nsi)
plt.xlabel('wavelength (nm)')
plt.ylabel('Refractive index')
plt.title('Silicon refractive index')