Internal function - put this function into uniroot to find the magnitude of the largest test statistic that will give a p-value of alpha.
GBJ_bound1_root(x, alpha, d, sig_vec)