Schervish's Method

This method uses a locally adaptive numerical integration algorithm, applied to the original problem. The details are given in the paper by Schervish (1984). One significant feature of this method is that it provides error bounds for the estimates of $P({\bf b})$ that it computes. The other methods only provide error estimates.