Mathematica code:
S[n_, t_] := Sin[n* Pi/25 + t]
ListAnimate[
Show[
Table[
Plot[
100 - n + Sum[(8*S[n,t] + .02)*Exp[-(x - k*5*S[n, 0])^2/Abs[.3*S[n,t]]], {k,{-1,1}}],
{x, -10, 10},
PlotStyle -> Directive[Black, Thick], PlotRange -> {{-7, 7}, {0, 100.5}},
Filling -> Axis, FillingStyle -> White, Axes -> False, AspectRatio -> Full,
ImageSize -> {500, 750}],
{n, -10, 100, 1}]],
{t, .001, 2 Pi + .001, (2 Pi + .001)/30}, AnimationRunning->False]