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 (left tailed test)
print(stats.t.cdf(test_stat, n-1))
0.8000961247260316