integrate x^(-1/2)*e^(-b*x^2)*sin(c*x) from 0 to infinity, real(b)>0