desvio_padrao_ex = 30
media_ex = 720
lim_inf_1 = 650
lim_sup_1 = 750
lim_2 = 800
lim_3 = 700
Z_inferior_1 = (lim_inf_1 - media_ex) / desvio_padrao_ex
Z_superior_1 = (lim_sup_1 - media_ex) / desvio_padrao_ex
Z_2 = (lim_2 - media_ex) / desvio_padrao_ex
Z_3 = (lim_3 - media_ex) / desvio_padrao_ex
probabilidade_1 = round((norm.cdf(Z_superior_1) - norm.cdf(Z_inferior_1)) * 100, 2)
probabilidade_2 = round((norm.cdf(-Z_2)) * 100, 2)
probabilidade_3 = round((norm.cdf(Z_3)) * 100, 2)
print(f"1) {probabilidade_1}% 2) {probabilidade_2}% 3) {probabilidade_3}% " )