Olá, essa é minha solução.
media = 300
desvio_padrao = 50
Z_inferior1 = (250 - media) / desvio_padrao
Z_superior1 = (350 - media) / desvio_padrao
Z_inferior2 = (400 - media) / desvio_padrao
Z_superior2 = (500 - media) / desvio_padrao
probabilidade1 = norm.cdf(Z_superior1) - norm.cdf(Z_inferior1)
probabilidade2 = norm.cdf(Z_superior2) - norm.cdf(Z_inferior2)
print(f'Probabilidade entre 250 e 350 = {round(probabilidade1 * 100, 2)}% | Probabilidade entre 400 e 500 = {round(probabilidade2 * 100, 2)}%')