import matplotlib.pyplot as plt
from matplotlib import cm
import numpy as np

#Taille du tableau NxN
N=100

#Déclaration et Initialisation du tableau
Ecl=[[0 for j in range(N)] for i in range(N)]

# Définition des paramètres physiques de l'expérience
lambda0=589E-9
alpha=3.1e-4


# Définition de la zone de calcul
R=12.5e-3
pas=R/(N-1)


# Calcul de l'éclairement à l'écran
for i in range(N):
    for j in range(N):
        x=(i-N/2)*pas
        y=(j-N/2)*pas
        
        #Différence de marche
        delta=2*alpha*x
        
        #Eclairement au point (x,y)
        Ecl[j][i]=2*(1+np.cos(2*np.pi*delta/lambda0))


#Tracé du tableau

plt.imshow(Ecl,cmap=cm.gist_gray)
plt.colorbar()
plt.axis('off')
plt.show()