ψ[x_, t_, α_] = 
      Exp[-((x^2)/(2 (1 - I α t)))]/(Sqrt[(Sqrt[π]) (1 + I α t)])

Plot[Abs[ψ[x, 0, 2]]^2, {x, -10, 10}, 
      PlotRange ->  {0,  .6}]

n = 100
Tmax = 10
graficas = 
      Table[Plot[Abs[ψ[x, i Tmax/n, 2]]^2, {x, -10, 10}, 
          PlotRange ->  {0,  .55}], {i, 0, n}]