ϕ[L_,n_,x_]:=(Sqrt[(2/L)])Sin[((n π)/L)x];

(* Grafica de la delta de Dirac centrada en 1/2, L=1*)
δ[x_,xp_]=Sum[ϕ[1,i,x]ϕ[1,i,xp],{i,1,30}];
Plot[δ[x,1/2],{x,0,1},PlotRange->All]

(* Grafica de la delta de Dirac centrada en 1/3, L=1 *)δ[x_,xp_]=Sum[ϕ[1,i,x]ϕ[1,i,xp],{i,1,30}];
Plot[δ[x,1/3],{x,0,1},PlotRange->All]

(* Grafica de la delta de Dirac centrada en 1/3, L=1 pero con 300 terminos en la suma *)
δ[x_,xp_]=Sum[ϕ[1,i,x]ϕ[1,i,xp],{i,1,300}];
Plot[δ[x,1/3],{x,0,1},PlotRange->All]