Animate[L = 10.0;
hbar = 4.14;
hbarc = 197;
mp = 940;
En[n_] := (n^2)(((hbarc^2)*(π^2))/(2*mp*(((2*L))^2)));
ω[n_]  := ((En[n] )/hbar);
(Subscript[ω, 8]) = ((En[8])/hbar);
bn[n_] := ((Sqrt[2])/π)×(16/((n^2) - 64));
(Subscript[N, start]) = 1;
(Subscript[N, stop]) = 81;
(Subscript[N, step]) = 4;
phi[x_] := Sum[ bn[n]×(1/(Sqrt[L]))×Sin[((n*π*x)/(2*L))]× Exp[-I*ω[n]*time],
		        {n, (Subscript[N, start]), (Subscript[N, stop]), (Subscript[N, step])}] +
	Sum[ -bn[n]×(1/(Sqrt[L]))×Sin[((n*π*x)/(2*L))]× Exp[-I*ω[n]*time],
		        {n, (Subscript[N, start])+2, (Subscript[N, stop]) + 2, (Subscript[N, step])}] +
	(1/(Sqrt[2]))×(1/(Sqrt[L]))×Sin[((8*π*x)/(2*L))]×Exp[-I*(Subscript[ω, 8])*time];
Plot[Evaluate[phi[x]*Conjugate[phi[x]]],{x,0.1,20},PlotRange->{{0,20},{0,0.22}},Frame->True,FrameLabel->{"x(fm)","(SuperStar[ψ])ψ",StringForm["t=`` \[Cross] (10^(-21)) s",NumberForm[time,{4,2}]," "]},
BaseStyle->FontSize->16,ImageSize->7.5*72
],{time,0,10.0,0.005}]