integrate exp(-b*P^2/2)*exp(i*P*x)*P dP from 0 to + inf