He aquí un método menos ingenioso (pero que siempre funciona):
Maximicemos $F(x,y)=x^2+4y^2$ bajo las restricciones $x,y\ge 0$ y $x^3+y^3-x+y=0$ utilizando Multiplicadores de Lagrange. Sea $f(x,y,\lambda)=F(x,y)+\lambda(x^3+y^3-x+y),$ entonces $$\dfrac{\partial f}{\partial x}=2x+3\lambda x^2-\lambda,\,\,\,\,\,\,\dfrac{\partial f}{\partial y}=8x+3\lambda y^2+\lambda ,\,\,\,\,\,\,\dfrac{\partial f}{\partial\lambda}=x^3+y^3-x+y$$ y obtenemos un sistema de ecuaciones altamente no lineal. Con la ayuda de wolframalpha podemos tener soluciones aproximadas muy exactas (las soluciones exactas son un poco complicadas). Ahora para $$x\approx 0.9725956862081514773 ,\,\,\,\,\,\, y\approx 0.0524320766715842$$ obtenemos $$F(x,y)\approx 0.95693885948708459008$$ que es una muy buena aproximación del máximo de $F.$