Questo script contiene una funzione normalcdf(a,mu,sigma) che calcola P(X≤a) nel caso di una distribuzione normale con media mu e con deviazione standard sigma. La funzione normalcdf2(a,b,mu,sigma) calcola 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)