import scipy.stats as st
a = range(10,14)
st.t.interval(0.95, len(a) - 1, loc = np.mean(a), 
              scale = st.sem(a))

(9.446, 13.554)