Mathematica code:
S[n_, t_] := Sin[n*3 Pi/50 + t]; Manipulate[ Show[ Table[ Plot[ 100 - n + (20*Abs[S[n,t]] + .02)*Exp[-(x - 4*S[n, t])^2/Abs[S[n,t]]], {x, -10, 10}, PlotStyle -> Directive[White, Thick], PlotRange -> {{-7, 7}, {0, 100}}, Filling -> Axis, FillingStyle -> White, Axes -> False, AspectRatio -> Full, ImageSize -> {500, 750}, Background -> RGBColor[54/255, 70/255, 93/255]], {n, 1, 100, 1}]], {t, .001, 2 Pi + .001, (2 Pi + .001)/30}]













