# Irwin–Hall distribution

In probability and statistics, the Irwin–Hall distribution, named after Joseph Oscar Irwin and Philip Hall, is probability distribution for a random variable defined as sum of a number of independent random variables, each having a uniform distribution. For this reason it is also known as the uniform sum distribution.

The generation of pseudo-random numbers having an approximately normal distribution is sometimes accomplished by computing the sum of a number of pseudo-random numbers having a uniform distribution; usually for the sake of simplicity of programming. Rescaling the Irwin–Hall distribution provides the exact distribution of the random variates being generated.

This distribution is sometimes confused with the Bates distribution, which is the mean (not sum) of n independent random variables uniformly distributed from 0 to 1.

## Definition

The Irwin–Hall distribution is the continuous probability distribution for the sum of n independent and identically distributed U(0, 1) random variables:

$X=\sum _{k=1}^{n}U_{k}.$ The probability density function (pdf) is given by

$f_{X}(x;n)={\frac {1}{2\left(n-1\right)!}}\sum _{k=0}^{n}\left(-1\right)^{k}{n \choose k}\left(x-k\right)^{n-1}\operatorname {sgn} (x-k)$ where sgn(x − k) denotes the sign function:

$\operatorname {sgn} \left(x-k\right)={\begin{cases}-1&xk.\end{cases}}$ Thus the pdf is a spline (piecewise polynomial function) of degree n − 1 over the knots 0, 1, ..., n. In fact, for x between the knots located at k and k + 1, the pdf is equal to

$f_{X}(x;n)={\frac {1}{\left(n-1\right)!}}\sum _{j=0}^{n-1}a_{j}(k,n)x^{j}$ where the coefficients aj(k,n) may be found from a recurrence relation over k

$a_{j}(k,n)={\begin{cases}1&k=0,j=n-1\\0&k=0,j0\end{cases}}$ The coefficients are also A188816 in OEIS. The coefficients for the cumulative distribution is A188668.

The mean and variance are n/2 and n/12, respectively.

## Special cases

$f_{X}(x)={\begin{cases}1&0\leq x\leq 1\\0&{\text{otherwise}}\end{cases}}$ $f_{X}(x)={\begin{cases}x&0\leq x\leq 1\\2-x&1\leq x\leq 2\end{cases}}$ • For n = 3,
$f_{X}(x)={\begin{cases}{\frac {1}{2}}x^{2}&0\leq x\leq 1\\{\frac {1}{2}}\left(-2x^{2}+6x-3\right)&1\leq x\leq 2\\{\frac {1}{2}}\left(x^{2}-6x+9\right)&2\leq x\leq 3\end{cases}}$ • For n = 4,
$f_{X}(x)={\begin{cases}{\frac {1}{6}}x^{3}&0\leq x\leq 1\\{\frac {1}{6}}\left(-3x^{3}+12x^{2}-12x+4\right)&1\leq x\leq 2\\{\frac {1}{6}}\left(3x^{3}-24x^{2}+60x-44\right)&2\leq x\leq 3\\{\frac {1}{6}}\left(-x^{3}+12x^{2}-48x+64\right)&3\leq x\leq 4\end{cases}}$ • For n = 5,
$f_{X}(x)={\begin{cases}{\frac {1}{24}}x^{4}&0\leq x\leq 1\\{\frac {1}{24}}\left(-4x^{4}+20x^{3}-30x^{2}+20x-5\right)&1\leq x\leq 2\\{\frac {1}{24}}\left(6x^{4}-60x^{3}+210x^{2}-300x+155\right)&2\leq x\leq 3\\{\frac {1}{24}}\left(-4x^{4}+60x^{3}-330x^{2}+780x-655\right)&3\leq x\leq 4\\{\frac {1}{24}}\left(x^{4}-20x^{3}+150x^{2}-500x+625\right)&4\leq x\leq 5\end{cases}}$ 