import numpy
import scipy.stats
import math

def sample_stdev(list):
    sos = 0
    avg = numpy.mean(list)

    for i in list:
        sos += (i - avg)**2

    s = math.sqrt(1.00 / (len(list) - 1.00) * sos)
    return s

B = [16.34,16.88,15.10,14.30,17.31,15.56,13.84,16.14,16.10,14.81,13.39,12.29,19.28]
A = [19.06,13.94,11.52,16.33,12.61,12.20,14.60,14.39,17.70,11.85,13.00,18.49,15.00]


sa = sample_stdev(A)
sb = sample_stdev(B)

sp = math.sqrt(((12.00) * sb**2.00 + (12) * sa**2) / 24.00)
print(sp)
print math.sqrt(2.00/13.00)

ans = (15.4876923077 - 14.6684615385) + (2.064 * sp * math.sqrt(2.00/13.00))

# Upper bound
print(ans)

# Lower bound
ans = (15.4876923077 - 14.6684615385) - (2.064 * sp * math.sqrt(2.00/13.00))
print(ans)