Integrate[ (Sin[t]^2 Exp[(x-x0)^2 / 2 / s])^2, {x,0,L}]