#*****************************************************
# programme de visualisation d'interférences de Young
# Dominique Lefebvre  mars 2013
# www.tangenteX.com
#*****************************************************

# chargement des packages
with(plots):

# initialisation des variables
e := 0.5E-4:      # largeur des fentes
a := 0.2E-3:      # distance entre les fentes
d := 2:           # distance fentes - écran
lambda := 0.7E-6: # longueur d'onde

# définition de la fonction sinus cardinal
sinc := x -> sin(Pi*x)/(Pi*x):

# fonction de calcul de l'intensité lumineuse 
f := x -> 4*cos(Pi*a*x/(lambda*d))^2*sinc(e*x/(lambda*d))^2:

# tracé de la figure d'interférences
densityplot(f(x),x=-0.1 .. 0.1,y = -1..1, grid=[400,50], axes = none, 
title ="Interférences de Young",style =patchnogrid);