import scipy.stats as stats
import math
# Specify the sample mean (x_bar), the sample standard deviation (s), the mean claimed in the null-hypothesis (mu_null), and the sample size (n)
x_bar = 62.1
s = 13.46
mu_null = 60
n = 30
# Calculate the test statistic
test_stat = (x_bar - mu_null)/(s/math.sqrt(n))
# Output the p-value of the test statistic (two tailed test)
print(2*(1-stats.t.cdf(test_stat, n-1)))
0.3998077505479367