Задача условной оптимизации:
![$\[\left\{ \begin{gathered}
  f = 2\sqrt {xyz}  \to \max  \hfill \\
  \sqrt {x - {x^2}}  + \sqrt {y - {y^2}}  + \sqrt {z - {z^2}}  = 2\sqrt {xyz}  \hfill \\ 
\end{gathered}  \right.\]$ $\[\left\{ \begin{gathered}
  f = 2\sqrt {xyz}  \to \max  \hfill \\
  \sqrt {x - {x^2}}  + \sqrt {y - {y^2}}  + \sqrt {z - {z^2}}  = 2\sqrt {xyz}  \hfill \\ 
\end{gathered}  \right.\]$](https://dxdy-04.korotkov.co.uk/f/f/4/7/f4790514e534154f438629b85cc4ffe482.png)
![$\[\begin{gathered}
  L = 2\sqrt {xyz}  - \lambda \left( {\sqrt {x - {x^2}}  + \sqrt {y - {y^2}}  + \sqrt {z - {z^2}}  - 2\sqrt {xyz} } \right) \hfill \\
  \frac{{\partial L}}
{{\partial x}} = 0 \Leftrightarrow \frac{{\sqrt {yz} }}
{{\sqrt x }} - \lambda \frac{{1 - 2x}}
{{2\sqrt {x - {x^2}} }} + 2\lambda \frac{{\sqrt {yz} }}
{{\sqrt x }} = 0 \Leftrightarrow \frac{{1 - 2x}}
{{2\sqrt {1 - x} }}\sqrt x  = \frac{{f + 2\lambda f}}
{{2\lambda }} \hfill \\ 
\end{gathered} \]$ $\[\begin{gathered}
  L = 2\sqrt {xyz}  - \lambda \left( {\sqrt {x - {x^2}}  + \sqrt {y - {y^2}}  + \sqrt {z - {z^2}}  - 2\sqrt {xyz} } \right) \hfill \\
  \frac{{\partial L}}
{{\partial x}} = 0 \Leftrightarrow \frac{{\sqrt {yz} }}
{{\sqrt x }} - \lambda \frac{{1 - 2x}}
{{2\sqrt {x - {x^2}} }} + 2\lambda \frac{{\sqrt {yz} }}
{{\sqrt x }} = 0 \Leftrightarrow \frac{{1 - 2x}}
{{2\sqrt {1 - x} }}\sqrt x  = \frac{{f + 2\lambda f}}
{{2\lambda }} \hfill \\ 
\end{gathered} \]$](https://dxdy-01.korotkov.co.uk/f/8/7/5/8751ff9bb89ea9c5cc7c57182cb0721182.png)
Аналогично для производных по другим переменным. Далее я замечаю, что функция 
![$\[h\left( w \right) = \frac{{1 - 2w}}
{{2\sqrt {1 - w} }}\sqrt w \]$ $\[h\left( w \right) = \frac{{1 - 2w}}
{{2\sqrt {1 - w} }}\sqrt w \]$](https://dxdy-04.korotkov.co.uk/f/7/c/9/7c9bb8009ce7a75432915da0d952da3d82.png)
 монотонна при 
![$\[w \geqslant \frac{1}{2}\]$ $\[w \geqslant \frac{1}{2}\]$](https://dxdy-03.korotkov.co.uk/f/e/4/e/e4e06e41729bc549501ad61d7f63a38082.png)
. Следовательно, 
![$\[\frac{{1 - 2x}}
{{2\sqrt {1 - x} }}\sqrt x  = \frac{{1 - 2y}}
{{2\sqrt {1 - y} }}\sqrt y  = \frac{{1 - 2z}}
{{2\sqrt {1 - z} }}\sqrt z  \Leftrightarrow x = y = z\]$ $\[\frac{{1 - 2x}}
{{2\sqrt {1 - x} }}\sqrt x  = \frac{{1 - 2y}}
{{2\sqrt {1 - y} }}\sqrt y  = \frac{{1 - 2z}}
{{2\sqrt {1 - z} }}\sqrt z  \Leftrightarrow x = y = z\]$](https://dxdy-01.korotkov.co.uk/f/8/8/c/88cbc7ffe5e142a83f076612da1ae9aa82.png)
.
Далее, пользуясь условием-ограничением 
![$\[{\sqrt {x - {x^2}}  + \sqrt {y - {y^2}}  + \sqrt {z - {z^2}}  - 2\sqrt {xyz} }\]$ $\[{\sqrt {x - {x^2}}  + \sqrt {y - {y^2}}  + \sqrt {z - {z^2}}  - 2\sqrt {xyz} }\]$](https://dxdy-03.korotkov.co.uk/f/e/f/9/ef952a02fc58a6916bcadefa8ffec0bb82.png)
 получаю, что на множестве 
![$\[W = \left\{ {\left( {x,y,z} \right)|x,y,z \in \left[ {\frac{1}
{2},1} \right)} \right\}\]$ $\[W = \left\{ {\left( {x,y,z} \right)|x,y,z \in \left[ {\frac{1}
{2},1} \right)} \right\}\]$](https://dxdy-03.korotkov.co.uk/f/e/7/d/e7d0ef360ad489f725dc8c548f576d2a82.png)
 имеем стационарную точку Лагранжа 

.
На этом строгость заканчивается. Дальше я хотел сослаться на верхнюю выпуклость функции 

, что означало бы, что найденная стационарная точка есть точка максимума, более того, это точка глобального максимума, так что другие случае на координаты рассматривать не приходится. Но никакого док-ва выпуклости у меня к сожалению нет... 
