Question: How do I speed up numerical integration of Bessel functions?

I am trying to evaluate the following integral to a float for various values of Zeta and rho. Evaluation of the integral seems to be very slow. Any suggestions on speeding it up? (I am a Maple newbie). "int(exp(-Zeta*eta)*BesselJ(0,rho*eta)*cos(eta),eta = 0 .. infinity)-int(exp(-Zeta*eta)*BesselJ(0,rho*eta)*sin(eta)/eta,eta = 0 .. infinity)+Zeta*int(eta*exp(-Zeta*eta)*BesselJ(0,rho*eta)*cos(eta),eta = 0 .. infinity)-Zeta*int(exp(-Zeta*eta)*BesselJ(0,rho*eta)*sin(eta),eta = 0 .. infinity)" Rgds, Saka.
Please Wait...