# -*- coding: utf-8 -*-
"""
Created on Mon Jan  3 11:53:34 2022

@author: Sim
"""
import numpy as np, scipy.stats as st
a = [31, 29, 27, 25, 31, 30, 27, 35, 31]

lbd, ubd = st.t.interval(0.95, len(a)-1, loc=np.mean(a), 
                         scale=st.sem(a))
print(lbd, ubd)