The best way to approach this question is to use lagrange multipliers
Here is a link if you would like to read an example: http://tutorial.math.lamar.edu/Classes/CalcIII/LagrangeMultipliers.aspx
To begin,
$$f(x, y, z) = x^3 + y^3 + z^3 - 3xyz $$
$$g(x, y, z) = x^2 + y^2 + z^2 $$
The constraint is that
$$g(x, y, z) = x^2 + y^2 + z^2 = 1 $$
Now we create our 4 equation, using the constraint and by taking the partial derivatives of the function.
$$f_x = \lambda g_x, 3x^2 = 2\lambda x ....(1)$$
$$f_y = \lambda g_y, 3y^2 = 2\lambda y ....(2)$$
$$f_z = \lambda g_z, 3z^2 = 2\lambda z ....(3)$$
$$x^2 + y^2 + z^2 = 1 ....(4)$$
Substituting equations (1), (2), and (3) into (4):
$$\frac{2\lambda x}{3} + \frac{2\lambda y}{3} + \frac{2\lambda z}{3} = 1$$
Also, notice that $$3x = 2\lambda$$
$$3y = 2\lambda$$
$$3z = 2\lambda$$
Meaning
$$x = y = z$$
So our maximum will occur when x and y and z are all equal.
Plugging this into our constraint:
$$3x^2 = 1$$ $$x = y = z = +/-\frac{\sqrt{3}}{3}$$
Notice that maximum value will only occur when x, y, and z are all positive.
Therefore, maximum value is:
$$3x^3 - 3x^3 = 0$$