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.
 
            