Add pivotal monte carlo confidence bounds to weibull.m