vecGauss = -> theta = alea()*pi r = 2*racine(-2*ln(alea())) [2.2+r*cos(theta),1.7+r*sin(theta)] partfrac = (v) -> [v[0]-troncature(v[0]), v[1]-troncature(v[1])] histogramme ((partfrac vecGauss())[1] for n in [1..100000]), 0, 1, 80, 10000