Continuity estimates for $n$-harmonic equations Tadeusz IwaniecJani Onninen 30C6035J1535J70$p$-harmonic equationmaximal functionsOrlicz-Sobolev imbedding We investigate the nonhomogeneous $n$-harmonic equation $$\mbox{div}\, |\nabla u|^{n-2}\nabla u =f$$ for $u$ in the Sobolev space $\mathscr{W}^{1,n}(\Omega)$, where $f$ is a given function in the Zygmund class $\mathscr{L}\log^\alpha \mathscr{L}(\Omega)$. In dimension $n=2$ the solutions are continuous whenever $f$ lies in the Hardy space $\mathscr{H}^1(\Omega)$, so in particular, if $f\in \mathscr{L}\log \mathscr{L}(\Omega)$. We show in higher dimensions that within the Zygmund classes the condition $\alpha > n-1$ is both necessary and sufficient for the solutions to be continuous. We also investigate continuity of the map $f \rightarrow \nabla u$, from $\mathscr{L}\log^\alpha \mathscr{L}(\Omega)$ into $\mathscr{L}^n\log^\beta \mathscr{L}(\Omega)$, for \[ -1 < \beta < \frac{n\alpha}{n-1}-1. \] These and other results of the present paper, though anticipated by simple examples, are in fact far from routine. Certainly, they are central in the $p$-harmonic theory. Indiana University Mathematics Journal 2007 text pdf 10.1512/iumj.2007.56.2987 10.1512/iumj.2007.56.2987 en Indiana Univ. Math. J. 56 (2007) 805 - 824 state-of-the-art mathematics http://iumj.org/access/