Vtgyjfy

I'm asked to find the global (=absolute) extremas of the function $$f(x,y,z)=xy-yz+xz$$ on the domain $$x^2+y^2 leq zleq1$$ which is obviously a paraboloid.

First, I find the extremas in the inside of the domain $$nabla f=(y,x-z,x-y)=(0,0,0)$$ and I get the trivial solution $$x=y=z=0$$. Then I thought about using Lagrange multiplier for the border with constraint $$g(x,y,z)=x^2+y^2-z$$ where where $z leq 1$. So :
$$nabla f=(y,x-z,x-y)=lambda nabla g=lambda(2x,2y,-1)$$ and $$x^2+y^2=z$$ After a bit of solving, I get $y=frac{x}{2x-1}$ and $z=x^2(1+frac{1}{(2x-1)^2})$.

At that point, the problem is that expecially the $x$ equivalence of $z$ seems a bit tricky, so I'm not sure if a made any mistake, took the wrong path or if there is an easier way to proceed in this situation. So my question is if my approach is correct, and if so, how should I proceed next, and also if there is an easier way to proceed on that conic domain than with Lagrange multiplier ?

Thanks for your help !

We first find any critical point in the interior of the domain, which is $x^2 + y^2 < 1.$ So we have $f'(x, y, z) = (y + z, x - z, x - y)$ and $f'(x, y, z) = 0$ can only happen when $x = y = z$ and $y = -z,$ so $x = y = z = 0,$ but this is not in the interior of the domain. Hence, there are no critical points.

Consider now one part of the frontrier, the one defined by $x^2 + y^2 = z = 1.$ Hence we have to maximise $u(x, y) = f(x, y, 1) = x + xy - y$ with $x^2 + y^2 = 1.$ Set $g(x, y) = x^2 + y^2 - 1.$ Since $g'(x, y) = (2x, 2y),$ the derivative of $g$ is never zero and the Lagrange multipliers method apply. We have to solve $u'(x, y) = dfrac{lambda}{2} g'(x, y),$ which signifies $(1 + y, x - 1) = lambda (x, y).$ Observe that $x = 1$ implies $y = 0$ and so we can consider the point $(1, 0)$ as one of the critical points in the restriction. Putting it aside, we have (by division), $dfrac{1 + y}{x - 1} = dfrac{x}{y}$ or else $y + y^2 = x^2 - x = 1 -y^2 - x$ or $x = 1 - 2y^2 - y.$ From the relation $x^2 + y^2 = 1$ we get $1 + 4y^4 + y^2 - 4y^2-2y+y^2 = 1,$ that is, $4y^4-2y^2-2y=0$ which leads to $y = 0$ (and hence $x = 1$) or else $2y^3-y-1=0.$ By inspection, we can factor $(y - 1)(2y^2+2y+1) = 0$ and the only solution is $y = 1$ with $x = -3,$ which is not in the restriction. Hence, we got only one point so far $(1, 0, 1).$

We finally tackle the other part of the restriction, which is $x^2 + y^2 - z = 0.$ Here we solve $f'(x, y, z) = lambda (2x, 2y, -1)$ and this leads to (upon summing) $y + x = 2lambda(x + y)$ and $y - x = lambda.$ If $x = -y,$ then $z = 2x^2 leq 1$ and we are optimising $h(x) = f(x, -x, 2x^2) = 4x^3 - x^2$ for $|x| leq dfrac{1}{2}$ and the optimisers for $h$ are $x = -dfrac{1}{sqrt{2}}$ and $x = sqrt{2} - dfrac{1}{2},$ this allow deducing the critical points for $f$ to be $left(-dfrac{1}{sqrt{2}}, dfrac{1}{sqrt{2}}, 1right)$ and $left(dfrac{1}{sqrt{2}}, -dfrac{1}{sqrt{2}}, 1right).$ If $x neq -y,$ then $lambda = dfrac{1}{2}$ and $y = x + dfrac{1}{2}.$ Again, substituting $x^2 + y^2 = z leq 1$ leads to the point (details ommited) $left(-dfrac{1+sqrt{7}}{4}, dfrac{1 - sqrt{7}}{4}, 1right).$

To terminate the problem, recall the region is compact and evaluate in each of the four points.

Take my calculations with a grain of salt, I did them and checked them twice, but with calculus of several variables, I always commit arithmetic errors.

Substitute $x = u - v$ and $y = u + v$ (thus $x^2 + y^2 = 2u^2 + 2v^2$). We extremize $$xy - yz + xz = u^2 - v^2 - 2vz = u^2 - (v + z)^2 + z^2$$
subject to the constraints $2u^2 + 2v^2 leq z leq 1$. Denote this rewritten function $f(u, v, z)$.

Analysis of the minimum.

If $v$ and $z$ are fixed, then $f$ is minimized at $u=0$. Minimizing $f$ for a fixed $z$ (but not $v$) is thus equivalent to maximizing $(v + z)^2$ over $0 leq v leq sqrt{z/2}$. This happens when $v$ is as large as possible, giving a minimum for fixed $z$ of $f(0, sqrt{z/2}, z) = - sqrt{2} z^{3/2} - frac{z}{2}$. The minimum over the whole region is thus $-frac{1}{2} - sqrt{2}$, attained at $(u, v, z) = (0, sqrt{1/2}, 1)$, viz. $(x, y, z) = (-sqrt{1/2}, sqrt{1/2}, 1)$.

Analysis of the maximum.

The maximum of $f(u, v, z)$ for fixed $z$ is on the boundary $u^2 + v^2 = z/2$; otherwise, we could increase $f$ by making $|u|$ larger. Thus we can substitute $u^2 = z/2 - v^2$ in
begin{align*}
f(u, v, z) &= u^2 - v^2 - 2vz \
&= frac{z}{2} - 2v^2 - 2vz \
&= -2 left(v + frac{z}{2} right)^2 + frac{z^2 + z}{2}.
end{align*}
For a fixed $z$, the maximum of $f$ is $frac{z^2 + z}{2}$, given at $v = -z/2$. (This point always falls within the constraints: $v^2 = z^2/4 < z/2$ for all $z in [0, 2]$.) The maximum over the entire region, therefore, is $1$, given by $(u, v, z) = (1/2, -1/2, 1)$, viz. $(x, y, z) = (1, 0, 1)$.

Interior Critical Points

Let
$$
f(x,y,z)=xy-yz+xztag1
$$
To find an interior critical point, we need
$$
begin{align}
0
&=delta f\
&=(y+z),delta x+(x-z),delta y+(x-y),delta ztag2
end{align}
$$
for all $delta x$, $delta y$, and $delta z$. Thus, the only interior critical point is $x=y=z=0$ and
$$
bbox[5px,border:2px dashed #C0A000]{f(0,0,0)=0}tag3
$$

Surface of the Paraboloid

On the paraboloid
$$
x^2+y^2-z=0tag4
$$
the allowable variations ($delta x$, $delta y$, and $delta z$) satisfy
$$
2x,delta x+2y,delta y-delta z=0tag5
$$
Orthogonality says that at the points for which the variations that satisfy $(5)$ also satisfy $(3)$, we have
$$
frac{y+z}{2x}=frac{x-z}{2y}=y-xtag6
$$
Applying $x^2+y^2=z$ and $frac{y+z}{2x}=y-x$, we get
$$
3x^2-2xy+y^2+y=0tag7
$$
Applying $x^2+y^2=z$ and $frac{x-z}{2y}=y-x$, we get
$$
x^2-2xy+3y^2-x=0tag8
$$
Subtracting $(8)$ from $(7)$ yields
$$
(x+y)(2x-2y+1)=0tag9
$$
Plugging $y=x+frac12$ into $(7)$ gives
$$
2x^2+x+tfrac34=0tag{10}
$$
which has no real solutions.

Plugging $y=-x$ into $(7)$ gives
$$
6x^2-x=0tag{11}
$$
which gives the critical points $(0,0,0)$, which was considered in $(3)$, and $left(frac16,-frac16,frac1{18}right)$ where
$$
bbox[5px,border:2px dashed #C0A000]{f!left(tfrac16,-tfrac16,tfrac1{18}right)=-tfrac1{108}}tag{12}
$$

Interior of the Base

On the base, where $z=1$ and $delta z=0$, $(2)$ says that the only interior critical point is $(1,-1,1)$, which is outside the base.

Edge of the Base

Using $(2)$ and $(5)$ on the edge of the base, orthogonality says that the critical points are where $frac{y+1}{2x}=frac{x-1}{2y}$; that is,
$$
y(y+1)=x(x-1)tag{13}
$$
Combining $(4)$ and $(13)$ yields
$$
x(x-1)left(4x^2-2right)tag{14}
$$
Combining $(4)$, $(13)$, and $(14)$ says that the critical points are in the set
$$
left{(0,-1,1),(1,0,1),left(tfrac1{sqrt2},-tfrac1{sqrt2},1right),left(-tfrac1{sqrt2},tfrac1{sqrt2},1right)right}tag{15}
$$
where
$$
begin{align}
&bbox[5px,border:2px solid #C0A000]{f!(0,-1,1)=1}\
&bbox[5px,border:2px solid #C0A000]{f!(1,0,1)=1}\
&bbox[5px,border:2px dashed #C0A000]{f!left(tfrac1{sqrt2},-tfrac1{sqrt2},1right)=sqrt2-tfrac12}\
&bbox[5px,border:2px solid #C0A000]{f!left(-tfrac1{sqrt2},tfrac1{sqrt2},1right)=-sqrt2-tfrac12}
end{align}tag{16}
$$
The extreme points are given in the solid boxed equations in $(16)$:
$$
-sqrt2-tfrac12le f(x,y,z)le1tag{17}
$$