  We prove that for Bernoulli bond percolation on $\mathbb{Z}^d$, $d\geq 2$ the
percolation density is an analytic function of the parameter in the
supercritical interval $(p_c,1]$. This answers a question of Kesten from 1981.
