u=xyz-P(x, y, z)(x2+y2+z2-R2)
u=-[Pxx(x2+y2+z2-R2)+2xPx+2xPx+2P+Pyy(x2+y2+z2-R2)+2yPy+2yPy+2P+ Pzz(x2+y2+z2-R2)+2zPz+2zPz+2P]=0
Rearranging it we can get:
P(x2+y2+z2-R2)+4(xPx+yPy+zPz)+6P=0
As P is rotational symmetric, we first consider for P(x, x, x)=mx+n
Then 3Pxx(3x2-R2)+12Pxx+6P=0
Plug in P=mx+n, Px=m, Pxx=0, we can get m=n=0
So P(x, y, z)=0
Then u=xyz