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 and print the test statistic
print((x_bar - mu_null)/(s/math.sqrt(n)))
0.854544851976857