the graphs of my provided code look approximately like the picture below. How can I plot such gpaph, help me please. Any suggestion appreciated. Thanks a lot.

```
Subscript(C, i)=2.5*10^6
Subscript(k, e)=315
σ=1*10^-9
Subscript(S, e)=1.58*10^-5
g=2.3*10^16
Subscript(C, e)=2.1*10^4
τ=1*10^-15
a=1/τ
Subscript(w, 1)=1
Subscript(s, 1)=y/(Subscript(w, 1)*σ)
Subscript(b, 1)=g/Subscript(C, e)*(1+(Subscript(k, e)*Subscript(s, 1)^2)/g)
Subscript(Δ, 1)=Sqrt(Subscript(b, 1)^2-4*Subscript(k, e)*Subscript(s, 1)^2*g/(Subscript(C, i)*Subscript(C, e)))
Subscript(p, 11)=(-Subscript(b, 1)+Subscript(Δ, 1))/2
Subscript(p, 12)=(-Subscript(b, 1)-Subscript(Δ, 1))/2
Subscript(T, i)=(Subscript(S, e)*g)/(2*π*τ*Subscript(C, i)*Subscript(C, e))*NIntegrate(BesselJ(0,y)*Exp(-((σ^2*Subscript(s, 1)^2)/4))*(Exp(-a*t)/((a+Subscript(p, 11))*(a+Subscript(p, 12)))+1/(Subscript(p, 11)-Subscript(p, 12))*(Exp(Subscript(p, 11)*t)/(Subscript(p, 11)+a)-Exp(Subscript(p, 12)*t)/(Subscript(p, 12)+a)))*y/(σ*Subscript(w, 1))^2,{y,0,100})
Subscript(T, e)=Subscript(T, i)+Subscript(S, e)/(2*π*τ*Subscript(C, e))*NIntegrate(BesselJ(0,y)*Exp(-((σ^2*Subscript(s, 1)^2)/4))*(-((a*Exp(-a*t))/((a+Subscript(p, 11))*(a+Subscript(p, 12))))+1/(Subscript(p, 11)-Subscript(p, 12))*((Subscript(p, 11)*Exp(Subscript(p, 11)*t))/(Subscript(p, 11)+a)-(Subscript(p, 12)*Exp(Subscript(p, 12)*t))/(Subscript(p, 12)+a)))*y/(σ*Subscript(w, 1))^2,{y,0,100})
Plot(Subscript(T, e),{t,0,1*10^-14})
Plot(Subscript(T, i),{t,0,1*10^-14})
```