First, because of the symmetric nature of this equation, we can factor every term. I'm going to go ahead and factor out the numerator of the first term.
(x5-x2) = x2(x3-1) = x2(x-1)(x2+x+1)
Then the first term will become this:
x2(x-1)(x2+x+1)/(x5+y2+z2)
If x is positive, then x2(x2+x+1) is also positive. (Also notice that these terms shrink quickly as you approach 0.) The problem reduces to solving whether
f(x,y,z) = (x-1)/(x5+y2+z2) + (y-1)/(x2+y5+z2) + (z-1)/(x2+y2+z5)
is positive. What makes you think that f(x,y,z) is positive in the interval 0≤{x,y,z}≤1?
Take the limit of the above as x, y, and z approach 0. You can make the transformation x=y=z in the limit since they all approach the same value, so
limx.y,z→0 f(x,y,z)
= limx→0 3(x-1)/(x5+x2+x2)
= limx→0 3(x-1)/[x2(x3+2)]
= limx→0 3/[x(x3+2)]-3/[x2(x3+2)]
(A polynomial of degree 5 grows faster than a polynomial of degree 4.)
= -∞
You're trying to prove something that's not true. Of course you're going to get stuck. Also notice that f(1,1,1) = 0, so you know that something is going to happen at (1,1,1). What happens is the following:
If (1,1,1) is a local minimum, then there is the possibility that the function is completely positive for (x,y,z)≥0. If not, then there are some points around (1,1,1) that are greater than zero and some points that are less than zero. I'll leave you to prove this for your general case, since you seem like you're the curious type. After proving this you will know why the conclusion you drew in the first post was erroneous.