Note: $(\frac{1}{\sqrt 2} , -\frac{1}{\sqrt 2} , 0)$ and $(\frac{1}{\sqrt 6} , \frac{1}{\sqrt 6} ,- \frac{2}{\sqrt 6})$ satisfy the constraint, and both are orthogonal to one annother.
We could parameterize the constraints.
$x = \frac {1}{\sqrt{2}} \cos \theta + \frac{1}{\sqrt 6} \sin\theta\\
y = -\frac {1}{\sqrt{2}} \cos \theta + \frac{1}{\sqrt 6} \sin\theta\\
z = - \frac{2}{\sqrt 6} \sin\theta$
then we have a problem in one variable.
maximize $(\frac {1}{\sqrt{2}} \cos \theta + \frac{1}{\sqrt 6} \sin\theta)(-\frac {1}{\sqrt{2}} \cos \theta + \frac{1}{\sqrt 6} \sin\theta)(- \frac{2}{\sqrt 6} \sin\theta)$
$(-\frac 12 \cos^2\theta + \frac 16\sin^2\theta)(-\frac{2}{\sqrt 6}\sin\theta)\\
\frac 1{\sqrt 6} \cos^2\theta\sin\theta - \frac 1{3\sqrt 6}\sin^3\theta\\
\frac {\sin 3\theta}{3\sqrt 6}$
This is maximized when $\sin 3\theta = 1,$ and minimized when it equals $-1$
$\pm\frac {1}{3\sqrt 6}$
Calculus is not actually required.
But you wanted to use Lagrange Multipliers.
$F(x,y,z,\lambda,\mu) = xyz - \lambda(x^2+y^2 + z^2 - 1) - \mu(x+y+z)\\
\frac{\partial F}{\partial x} = yz - 2\lambda x -\mu = 0\\
\frac{\partial F}{\partial y} = xz - 2\lambda y -\mu = 0\\
\frac{\partial F}{\partial z} = xy - 2\lambda z -\mu = 0\\
\frac{\partial F}{\partial \lambda} = x^2 + y^2 + z^2 - 1 = 0\\
\frac{\partial F}{\partial \lambda} = x+y+z = 0$
Setting line 1 equal to line 2.
$yz - 2\lambda x = xz - 2\lambda y\\
(x-y) z = 2\lambda (y-x)$
$y = x$ or $-2\lambda = {z}$
first case: $x = y$
$z = -2x\\
x^2 +y^2 + z^2 = x^2 + x^2 + (2x)^2 = 6x^2 = 1\\
(x,y,z) = (\pm \frac{1}{\sqrt{6}}, \pm\frac{1}{\sqrt{6}}, \mp \frac {2}{\sqrt 6})\\
xyz = \pm \frac{1}{3\sqrt 6}$
else $-2\lambda = z$
$yz + zx = \mu\\
xy + z^2 = \mu\\
z^2 -yz - zx +xy = 0\\
(z-x)(z-y) = 0$
And we get a similar result.