from scipy import mean, std, sqrt, absolute
import scipy.stats as st
minta=[1193, 1198, 1203, 1191, 1195, 1196, 1199, 1191,\
1201, 1196, 1193, 1198, 1204, 1196, 1198, 1200]
n=len(minta)
matlag=mean(minta)
mu0=1200
sig=3
alfa=0.05 # 1-alfa=0.95 a feladatban

print("u:")
u=sqrt(n)*(matlag-1200)/sig
if absolute(u)<st.norm.ppf(1-0.5*alfa): print("accept")
else: print("reject")



print("t:")
s=std(minta)
t=sqrt(n)*(matlag-1200)/s
if absolute(t)<st.t.ppf(1-0.5*alfa,df=n-1): print("accept")
else: print("reject")

# ketoldali kritikus ertekek a szokasos szinteken:
# for alfa in [0.1,0.05,0.01]:
#   print(st.norm.ppf(1-0.5*alfa))