plotting – Code that produces plot in V5 doesn’t work in later versions


I have problem in plotting Integral function.

I can compute/plot the graph of this integration below in Mathematica 5.0, but it is not possible to plot it in higher Mathematica versions.

My code is:

w = 100
b = 438095.2 (2.5*10^6 + 547826*y^2)
d = Sqrt(-8.8*10^21 *y^2 + 1.92*10^11 *(2.52*10^6 + 547826*y^2)^2)
Subscript(p, 1) = (-b + d)/2
Subscript(p, 2) = (-b - d)/2
T = 10^18 *Sqrt((Pi)/2)*NIntegrate(BesselJ(0, y)/Sqrt((y^2 + w^2)^3)*1/
d*(1 - Erf((Subscript(p, 2)*10^-15 - 1)/Sqrt(2)))*
Exp((Subscript(p, 2)*10^-15)^2/2 - Subscript(p, 2)*10^-15*(1 - (Tau)))*y*w, {y, 0, (Infinity)}, 
Method -> Oscillatory, MaxRecursion -> 12)
Plot(T, {(Tau), 0, 6*10^3}) 

code

The plot is generated by Mathematica 5.0:

plot

Error in later versions is NIntegrate::nconv

Can anybody please help? Any suggestions appreciated.