Let $u(x,y)$ be a harmonic function on domain s.t all the partial derivatives of $u(x,y)$ vanish at the same point of , then $u(x,y)$ is constant.
Now the thing is if the harmonic conjugate of $u(x,y)$ exists say $v(x,y)$ then $f=u+iv$ is analytic and $f^m(z)$ vanishes for all $z \in D$ then $f(z)$ is const so is $u(x,y)$.
But the question is I know that for a star shaped domain the harmonic conjugates exists not for any domain.. So what will be the solution of this?
For each $z\in D$ there is a star-shaped region $B$ with $z\in B\subset D$. (For example, a small enough open ball centered at $z$.) By your harmonic conjugate reasoning $u$ is constant on $B$. This (and the connectedness of $D$) is enough to show that $u$ is constant on all of $D$.