integrate x^a * (1 - x)^(b - 1) dx from 0 to 1