Reorganize some code, fix beta PDF for a,x = 1,0 and b,x = 1,1
Added Brent-Dekker floating-point root finder
Made `real-dist-hpd-interval' return type more precise; documented
Added HPD interval estimates and Monte Carlo variance
Defer building tables; drops 600ms from (require math) on my machine