#*************************************************************************
# programme de tracé d'un paquet d'ondes de spectre gaussien
# Dominique Lefebvre  mars 2013
# www.tangenteX.com
#*************************************************************************

# chargement des packages
with(plots):
with(inttrans):

# initialisation des variables
w0 := 10:
v := 2:

# définition du spectre gaussien centré en omega0
g := w -> exp(-(w-w0)^2):
# affichage du spectre gaussien
plot(g(w),w=0..2*w0, numpoints=1000);

# calcul de la fonction d'onde en appliquant la définition
psi := (x,t) -> fourier(exp(I*(w/v)*x)*g(w),w,t):

# tracé du paquet d'ondes
animate(evalc(Re(psi(x,t))),x=0..30,t=0..10, numpoints=200);