integrate r e^(-ar^2+b cos t r) dr