So far this is what I have done:
$$F(x,y) = \sqrt{xy} + λ(x+y-1) =0$$
$$F_x = \frac12(xy)^\left(\frac{-1}{2}\right).y + λ=0$$
$$F_y = \frac12(xy)^\left(\frac{-1}{2}\right).x + λ=0$$
$$Fλ = x+y-1=0$$
I found that $x=y$ and $x=y=\frac12$
I'm not sure if this makes sense and where to go from there.
This $3$D problem can be brought back to a $2$D issue:
From $x+y=1$, one has $y=1-x$. Then we have to find the maximum of $z=\sqrt{x(1-x)}$ which occurs at $x=\frac{1}{2}$ with value $z=\frac{1}{2}$ $\left(\text{using derivatives or the fact that the maximum value for the product of quantities whose sum is constant, here}\; x\; \text{and}\; 1-x, \text{is reached when they are equal}\; x=1-x \Rightarrow x=\frac12\right)$
Remark: In geometrical terms, we have sliced the surface with equation $z=\sqrt{xy}$ with vertical plane $x+y=1$, and have had a reasoning on the intersection curve (a circle).