Na resolução da atividade consta, nos parâmetros do norm.interval() "scale = desvio_padrao / np.sqrt(n)". Acredito que esteja faltando a múltiplicação do Z referente ao nível de confiança de 90% (1.645), assim, o intervalo ficaria: (27.330435569981386, 28.669564430018614). É isso mesmo ou estou equivocado?
Olá José, a formula com o parâmetro scale, não necessita do valor do erro inferencial, no qual você precisa multiplicar o Z com o sigma (desvio_padrao / np.sqrt(n)) .
Para a formula só é necessário o valor do sigma.
norm.interval(alpha = 0.90,
loc = media_amostral,
scale = desvio_padrao / np.sqrt(n))
Espero ter ajudado, abraços.