Integrate exp( -x**2/2/a**2 - y**2/2/b**2 ) from -Infinity to +Infinity