Analyticity of the percolation density $theta$ in all dimensions

We prove that for Bernoulli bond percolation on $mathbb{Z}^d$, $dgeq 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.