import scipy.stats as stats print(2*(1-stats.t.cdf(0.855, 29)))