Here's a counterexample in $n\ge3$ dimensions. With $\lVert\cdot\rVert$ denoting the Euclidean norm, set $$ f(x)=\begin{cases} x_1x_2x_3\lVert x\rVert^{-d-4},&{\rm if\ }x\not=0\cr 0,&{\rm if\ }x=0. \end{cases} $$ This is harmonic on $x\not=0$, has first and second order derivatives at the origin, but $f$ along with all its derivatives are discontinuous at the origin.
Note that $f=\partial^3(\lVert x\rVert^{-d+2})/\partial x_1\partial x_2\partial x_3$ is a third order derivative of a harmonic function, so is harmonic (away from the origin). As $f$ vanishes whenever all but two coordinates are zero, it along with its first order derivatives all vanish on the coordinate axes. So, $f$ along with all its first and second order derivatives vanish at the origin.