import scipy.stats as stats import math # Specify the number of occurences (x), the sample size (n), and the proportion claimed in the null-hypothesis (p) x = 10 n = 100 p = 0.5 # Calculate the sample proportion p_hat = x/n # Calculate the test statistic test_stat = (p_hat-p)/(math.sqrt((p*(1-p))/(n))) # Output the p-value of the test statistic (two-tailed test) print(2*stats.norm.cdf(test_stat))