Este script contém uma função normalcdf(a,mu,sigma)que calcula P(X≤a) no caso de uma distribuição normal de média mu e de desvio-padrão sigma. A função normalcdf2(a,b,mu,sigma) calcula P(a≤X≤b).
from math import * def normalcdf(a,mu,sigma): return (0.5)*(1+erf((a-mu)/(sqrt(2)*sigma))) def normalcdf2(a,b,mu,sigma): return normalcdf(b,mu,sigma)-normalcdf(a,mu,sigma)