Tentei resolver o exercício da maneira mais concisa possível, consultando somente a documentação do scipy.
Talvez a maioria dos colegas já tenha se dado conta, mas como o loc
da função norm.interval
faz referência à media, portanto o parâmetro loc = 0
determina que a função calcule o intervalo de confiança para uma distribuição normal.
Logo norm.interval(alpha = 0.95, loc = 0, scale = sigma)
resolve o problema.