#*************************************************************************
# programme de tracé d'une figure de diffraction
# Dominique Lefebvre  novembre 2010
# www.tangenteX.com
#*************************************************************************
#
# Tracé d'une figure de diffraction par une fente rectangulaire sous Maple
with(plots):
#
# Définition des paramètres
lx := 1E-4:            # largeur de la fente en m
ly := 2E-4:            # hauteur de la fente en m
lambda := 0.65E-6:     # longeur d'onde de la lumière incidente
L := 1:                # distance fente-écran
K := Pi/(L*lambda):
# Définition de l'étendue de calcul
XMIN := -0.02:
XMAX := 0.02:
YMIN := -0.02:
YMAX := 0.02:
#
# Définition de la fonction sinc(t) = sin(t)/t
sinc := t->sin(t)/(t):
#
# Calcul de l'intensité lumineuse
IL := (x,y)->(sinc(K*lx*x))^2*(sinc(K*ly*y))^2:
#
# tracé de la figure
densityplot(IL(x,y),x=XMIN..XMAX,y=YMIN..YMAX,grid=[200,200],
style=patchnogrid, color=RED, contrast=0.8);