This function computes the coverage of the interval
probet(a, b, c, alpha)
numerical value between 0 and 1 corresponding to the coverage
lower bound of the prior 95%~confidence interval
upper bound of the prior 95%~confidence interval
mean parameter of the prior distribution
scale parameter of the prior distribution
solbeta