from scipy.stats import t
import numpy
dizi = [26, 33, 42, 45, 33, 53]
aralik = 0.95
diziOrt = numpy.mean(dizi)
stdSapma = numpy.std(dizi)
serbestlik_derecesi = len(dizi)-1
tkritigi = numpy.abs(t.ppf((1-aralik)/2, serbestlik_derecesi))
print((diziOrt-stdSapma*tkritigi/numpy.sqrt(len(dizi)), diziOrt+stdSapma*tkritigi/numpy.sqrt(len(dizi))) )