from pylab import *
from IndiceVerre import *
n1 = IndiceVerre(1.51,60)
n2 = IndiceVerre(1.60,38)
Lambda = arange(400,700,1)
n1a = n1.indiceArray(Lambda)
n2a = n2.indiceArray(Lambda)
figure()
plot(Lambda,n1a,color='r',label='Crown')
plot(Lambda,n2a,color='b',label='Dense flint')
xlabel('lambda (nm)')
ylabel('n')
legend()
grid(True)