evalf(int(x*e^(-x^2)*J_0(2*x),x=0..infty))