Finding the maximum value of the following product. 
Let $x,y,z$ be three reals which satisfies $x^2+y^2+z^2=1$. Find the maximum value of $$P = (x^2-yz)(y^2-zx)(z^2-xy).$$

I think expanding would be a bad idea. I tried to apply the Weirstrass product inequality but it may give the minimum not maximum value.

Lagrange multiplier is the correct tool for solving this type of constrained optimization problem. I don't know your backgroud, but the following is a completely elementary approach (which I only found out with the help of Lagrange).
First it's not hard to show that the maximum (which exists as the sphere is compact) is positive. If the maximum is achieved when two of the three factors are negative, say $x^2-yz<0$ and $y^2-xz<0$, then $yz>x^2\ge 0, xz>y^2\ge 0$, therefore $x,y,z$ all have the same sign. But then we may flip the sign of $z$, so that $(x^2+yz)(y^2+xz)(z^2-xy)\ge (x^2-yz)(y^2-xz)(z^2-xy)$. Therefore the max can only be achieved when all three factors are positive.
Use the geometric mean $\le$ the arithmetic mean, we have $f(x,y,z)^{1/3}\le\frac{x^2+y^2+z^2-xy-yz-xz}{3}=\frac{1-xy-yz-xz}{3}$. And equality holds only when $x^2-yz=y^2-xz=z^2-xy \Leftrightarrow (x-y)(x+y+z)=0, (y-z)(x+y+z)=0$. As $x,y,z$ cannot be all equal, we must have $x+y+z=0$. If we can show further that $\frac{1-xy-yz-xz}{3}$ achieves a maximum when $x+y+z=0$, then we can conclude maximum can be achieved for $f(x,y,z)$ for the same triple.
We are left to find when $\frac{xy+yz+xz}{3}$ achieves minimum on the sphere. This is equivalent to find the minimum of $(x+y+z)^2=(x^2+y^2+z^2) + 2(xy+yz+xz)=1+2(xy+yz+xz)$ on the sphere, and clearly the minimum of  $(x+y+z)^2$ is achieved when $x+y+z=0$.
Put all the work together and reverse the reasoning, we conclude that $f(x)$ achieves a maximum at $(x,y,z)\in \mathbb S^2$ if and only if $x^2-yz, y^2-xz, z^2-xy$ are all positive and $x+y+z=0$. (In fact, e.g. $x^2-yz = y^2+z^2+yz\ge 0$ on $x+y+z=0$, hence the condition of the three factors being positive is redundant. And the intersection of the sphere with the plane are exactly all the solutions to the optimization problem.)
Now we may use any point on the intersection of the plane and the sphere to calculate the maximum. E.g. $x=\frac{1}{\sqrt 2}, y = -\frac{1}{\sqrt 2}, z = 0$, then the maximum is just $\frac{1}{8}$ as pointed out in the comment.
