integrate x/ (a * sqrt(2 * pi)) exp(-(x -b)**2/(2*a**2))dx