integrate ln(1-b^2 sin^2(x)) dx from 0 to pi/2