Add calculation of pivotal monte carlo confidence bounds to weibull.m