#************************************************************************* # 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);