Analiza matematyczna 2/Ćwiczenia 8: Ekstrema funkcji wielu zmiennych

Ekstrema funkcji wielu zmiennych. Ćwiczenia

Ćwiczenie 8.1.

a) Wyznaczyć wielomian Taylora rzędu drugiego funkcji $\displaystyle f(x,y)={\frac {\cos x}{\cos y}}$ w punkcie $\displaystyle (0,0)$.

b) Wyznaczyć wielomian Taylora rzędu drugiego funkcji $\displaystyle f(x,y)=\mathrm {arctg} \,({\frac {x-y}{x+y}})$ w punkcie $\displaystyle (1,1)$.

c) Wyznaczyć wielomian Taylora rzędu drugiego funkcji $\displaystyle f(x,y)={\frac {xy}{x^{2}+y^{2}}}$ w punkcie $\displaystyle (1,1)$.

d) Rozwinąć w szereg Taylora funkcję $\displaystyle f(x,y,z)=x^{3}+y^{3}+z^{3}-3xyz$ w punkcie $\displaystyle (1,1,1)$.

Wskazówka

Jak wyraża się wielomian Taylora za pomocą pochodnych cząstkowych?

d) Ile pochodnych cząstkowych niezerowych ma funkcja $\displaystyle f$?

Rozwiązanie

a) Obliczmy wartości pochodnych cząstkowych funkcji $\displaystyle f$ w punkcie $\displaystyle (0,0)$. Mamy

$\displaystyle {\begin{aligned}&{\frac {\partial f}{\partial x}}=-{\frac {\sin x}{\cos y}},&&\qquad {\frac {\partial f}{\partial x}}(0,0)=0;\\&{\frac {\partial f}{\partial y}}={\frac {\cos x\sin y}{\cos ^{2}y}},&&\qquad {\frac {\partial f}{\partial y}}(0,0)=0;\\&{\frac {\partial ^{2}f}{\partial x^{2}}}=-{\frac {\cos x}{\cos y}},&&\qquad {\frac {\partial ^{2}f}{\partial x^{2}}}(0,0)=-1;\\&{\frac {\partial ^{2}f}{\partial y^{2}}}={\frac {\cos x\cos y(1+\sin ^{2}y)}{\cos ^{3}y}},&&\qquad {\frac {\partial ^{2}f}{\partial y^{2}}}(0,0)=1;\\&{\frac {\partial ^{2}f}{\partial x\partial y}}=-{\frac {\sin x\sin y}{\cos ^{2}y}},&&\qquad {\frac {\partial ^{2}f}{\partial x\partial y}}(0,0)=0.\\\end{aligned}}$

Tak więc wielomian Taylora rzędu drugiego funkcji $\displaystyle f$ w punkcie $\displaystyle (0,0)$ ma postać

$\displaystyle T_{(0,0)}^{2}f(h_{1},h_{2})=1+{\frac {1}{2}}(-h_{1}^{2}+h_{2}^{2}).$

b) Obliczmy wartości pochodnych cząstkowych funkcji $\displaystyle f$ w punkcie $\displaystyle (1,1)$. Mamy

$\displaystyle {\begin{aligned}&{\frac {\partial f}{\partial x}}={\frac {y}{x^{2}+y^{2}}},&&\qquad {\frac {\partial f}{\partial x}}(1,1)={\frac {1}{2}};\\&{\frac {\partial f}{\partial y}}={\frac {-x}{x^{2}+y^{2}}},&&\qquad {\frac {\partial f}{\partial y}}(1,1)=-{\frac {1}{2}},\\&{\frac {\partial ^{2}f}{\partial x^{2}}}=-{\frac {2xy}{(x^{2}+y^{2})^{2}}},&&\qquad {\frac {\partial ^{2}f}{\partial x^{2}}}(1,1)=-1;\\&{\frac {\partial ^{2}f}{\partial y^{2}}}={\frac {2xy}{(x^{2}+y^{2})^{2}}},&&\qquad {\frac {\partial ^{2}f}{\partial y^{2}}}(1,1)=1,\\&{\frac {\partial ^{2}f}{\partial x\partial y}}={\frac {x^{2}-y^{2}}{(x^{2}+y^{2})^{2}}},&&\qquad {\frac {\partial ^{2}f}{\partial x\partial y}}(1,1)=0.\\\end{aligned}}$

Tak więc wielomian Taylora rzędu drugiego funkcji $\displaystyle f$ w punkcie $\displaystyle (1,1)$ ma postać

$\displaystyle T_{(1,1)}^{2}f(h_{1},h_{2})={\frac {1}{2}}h_{1}-{\frac {1}{2}}h_{2}+{\frac {1}{2}}\left(-h_{1}^{2}+h_{2}^{2}\right).$

c) Obliczmy wartości pochodnych cząstkowych funkcji $\displaystyle f$ w punkcie $\displaystyle (1,1)$. Mamy

$\displaystyle {\begin{aligned}&{\frac {\partial f}{\partial x}}={\frac {y^{3}-x^{2}y}{(x^{2}+y^{2})^{2}}},&&\qquad {\frac {\partial f}{\partial x}}(1,1)=0;\\&{\frac {\partial f}{\partial y}}={\frac {x^{3}-xy^{2}}{(x^{2}+y^{2})^{2}}},&&\qquad {\frac {\partial f}{\partial y}}(1,1)=0,\\&{\frac {\partial ^{2}f}{\partial x^{2}}}={\frac {2x^{3}y-6xy^{3}}{(x^{2}+y^{2})^{3}}},&&\qquad {\frac {\partial ^{2}f}{\partial x^{2}}}(1,1)=-{\frac {1}{2}};\\&{\frac {\partial ^{2}f}{\partial y^{2}}}={\frac {2xy^{3}-6x^{3}y}{(x^{2}+y^{2})^{3}}},&&\qquad {\frac {\partial ^{2}f}{\partial y^{2}}}(1,1)=-{\frac {1}{2}},\\&{\frac {\partial ^{2}f}{\partial x\partial y}}={\frac {-x^{4}-y^{4}+6x^{2}y^{2}}{(x^{2}+y^{2})^{3}}},&&\qquad {\frac {\partial ^{2}f}{\partial x\partial y}}(1,1)={\frac {1}{2}}.\\\end{aligned}}$

Tak więc wielomian Taylora rzędu drugiego funkcji $\displaystyle f$ w punkcie $\displaystyle (1,1)$ ma postać

$\displaystyle T_{(1,1)}^{2}f(h_{1},h_{2})={\frac {1}{2}}+{\frac {1}{2}}\left(-{\frac {1}{2}}h_{1}^{2}-{\frac {1}{2}}h_{2}^{2}+h_{1}h_{2}\right).$

d) Obliczmy wartości pochodnych cząstkowych funkcji $\displaystyle f$ w punkcie $\displaystyle (1,1,1)$. Mamy

$\displaystyle {\begin{aligned}&{\frac {\partial f}{\partial x}}=3x^{2}-3yz,&&\qquad {\frac {\partial f}{\partial x}}(1,1,1)=0;\\&{\frac {\partial f}{\partial y}}=3y^{2}-3xz,&&\qquad {\frac {\partial f}{\partial y}}(1,1,1)=0,\\&{\frac {\partial f}{\partial z}}=3z^{2}-3xy,&&\qquad {\frac {\partial f}{\partial z}}(1,1,1)=0;\\&{\frac {\partial ^{2}f}{\partial x^{2}}}=6x,&&\qquad {\frac {\partial ^{2}f}{\partial x^{2}}}(1,1,1)=6;\\&{\frac {\partial ^{2}f}{\partial y^{2}}}=6y,&&\qquad {\frac {\partial ^{2}f}{\partial y^{2}}}(1,1,1)=6;\\&{\frac {\partial ^{2}f}{\partial z^{2}}}=6z,&&\qquad {\frac {\partial ^{2}f}{\partial z^{2}}}(1,1,1)=6;\\&{\frac {\partial ^{2}f}{\partial x\partial y}}=-3z,&&\qquad {\frac {\partial ^{2}f}{\partial x\partial y}}(1,1,1)=-3;\\&{\frac {\partial ^{2}f}{\partial x\partial z}}=-3y,&&\qquad {\frac {\partial ^{2}f}{\partial x\partial z}}(1,1,1)=-3;\\&{\frac {\partial ^{2}f}{\partial y\partial z}}=-3x,&&\qquad {\frac {\partial ^{2}f}{\partial y\partial z}}(1,1,1)=-3;\\&{\frac {\partial ^{3}f}{\partial x^{3}}}=6,&&\qquad {\frac {\partial ^{3}f}{\partial x^{3}}}(1,1,1)=6;\\&{\frac {\partial ^{3}f}{\partial y^{3}}}=6,&&\qquad {\frac {\partial ^{3}f}{\partial y^{3}}}(1,1,1)=6;\\&{\frac {\partial ^{3}f}{\partial z^{3}}}=6,&&\qquad {\frac {\partial ^{3}f}{\partial z^{3}}}(1,1,1)=6.\end{aligned}}$

Pozostałe pochodne cząstkowe są równe zero. Tak więc rozwinięcie funkcji $\displaystyle f$ w szereg Taylora w punkcie $\displaystyle (0,0,0)$ ma postać

$\displaystyle {\begin{aligned}&f(x,y,z)=\\&{\frac {1}{2}}\left(6(x-1)^{2}+6(y-1)^{2}+6(z-1)^{2}-6(x-1)(y-1)-6(x-1)(z-1)-6(y-1)(z-1)\right)\\&+{\frac {1}{6}}\left(6(x-1)^{3}+6(y-1)^{3}+6(z-1)^{3}-18(x-1)(y-1)(z-1)\right).\end{aligned}}$

Ćwiczenie 8.2.

Wyznaczyć ekstrema lokalne funkcji

a) $\displaystyle f(x,y)=x^{4}+y^{4}-8x^{2}-2y^{2}+2006$,

b) $\displaystyle g(x,y)=x^{2}+8y^{3}-6xy+1,$

c) $\displaystyle \displaystyle h(x,y)=2xy+{\frac {1}{x}}+{\frac {2}{y}}$.

Wskazówka

Należy poszukać punktów krytycznych i zastosować kryterium Sylvestera do macierzy drugiej różniczki.

a) Jeśli $\displaystyle x_{1}$, $\displaystyle x_{2}$ i $\displaystyle x_{3}$ są pierwiastkami równania Parser nie mógł rozpoznać (SVG (MathML może zostać włączone przez wtyczkę w przeglądarce): Nieprawidłowa odpowiedź („Math extension cannot connect to Restbase.”) z serwera „https://wazniak.mimuw.edu.pl/api/rest_v1/”:): {\displaystyle \displaystyle p(x)=0} z jedną niewiadomą i Parser nie mógł rozpoznać (SVG (MathML może zostać włączone przez wtyczkę w przeglądarce): Nieprawidłowa odpowiedź („Math extension cannot connect to Restbase.”) z serwera „https://wazniak.mimuw.edu.pl/api/rest_v1/”:): {\displaystyle \displaystyle y_1} , Parser nie mógł rozpoznać (SVG (MathML może zostać włączone przez wtyczkę w przeglądarce): Nieprawidłowa odpowiedź („Math extension cannot connect to Restbase.”) z serwera „https://wazniak.mimuw.edu.pl/api/rest_v1/”:): {\displaystyle \displaystyle y_2} , Parser nie mógł rozpoznać (SVG (MathML może zostać włączone przez wtyczkę w przeglądarce): Nieprawidłowa odpowiedź („Math extension cannot connect to Restbase.”) z serwera „https://wazniak.mimuw.edu.pl/api/rest_v1/”:): {\displaystyle \displaystyle y_3} są pierwiastkami równania Parser nie mógł rozpoznać (SVG (MathML może zostać włączone przez wtyczkę w przeglądarce): Nieprawidłowa odpowiedź („Math extension cannot connect to Restbase.”) z serwera „https://wazniak.mimuw.edu.pl/api/rest_v1/”:): {\displaystyle \displaystyle q(y)=0} z jedną niewiadomą, to jakie rozwiązania ma układ dwóch równań (z dwoma niewiadomymi) Parser nie mógł rozpoznać (SVG (MathML może zostać włączone przez wtyczkę w przeglądarce): Nieprawidłowa odpowiedź („Math extension cannot connect to Restbase.”) z serwera „https://wazniak.mimuw.edu.pl/api/rest_v1/”:): {\displaystyle \displaystyle p(x)=0} i Parser nie mógł rozpoznać (SVG (MathML może zostać włączone przez wtyczkę w przeglądarce): Nieprawidłowa odpowiedź („Math extension cannot connect to Restbase.”) z serwera „https://wazniak.mimuw.edu.pl/api/rest_v1/”:): {\displaystyle \displaystyle q(y) = 0} ?

Rozwiązanie

a) Z warunku koniecznego istnienia ekstremum otrzymujemy układ dwóch niezależnych równań Parser nie mógł rozpoznać (SVG (MathML może zostać włączone przez wtyczkę w przeglądarce): Nieprawidłowa odpowiedź („Math extension cannot connect to Restbase.”) z serwera „https://wazniak.mimuw.edu.pl/api/rest_v1/”:): {\displaystyle \displaystyle 4x^3-16x=0} i Parser nie mógł rozpoznać (SVG (MathML może zostać włączone przez wtyczkę w przeglądarce): Nieprawidłowa odpowiedź („Math extension cannot connect to Restbase.”) z serwera „https://wazniak.mimuw.edu.pl/api/rest_v1/”:): {\displaystyle \displaystyle 4y^3-4y=0} . Pierwsze z nich ma rozwiązania Parser nie mógł rozpoznać (SVG (MathML może zostać włączone przez wtyczkę w przeglądarce): Nieprawidłowa odpowiedź („Math extension cannot connect to Restbase.”) z serwera „https://wazniak.mimuw.edu.pl/api/rest_v1/”:): {\displaystyle \displaystyle -2, 0, 2} , drugie Parser nie mógł rozpoznać (SVG (MathML może zostać włączone przez wtyczkę w przeglądarce): Nieprawidłowa odpowiedź („Math extension cannot connect to Restbase.”) z serwera „https://wazniak.mimuw.edu.pl/api/rest_v1/”:): {\displaystyle \displaystyle -1, 0, 1} . Punktami krytycznymi są więc pary $\displaystyle (0,0),(0,-1),(0,1),(2,0),(-2,0),(-2,-1),(-2,1),(2,-1),(2,1)$. Liczymy pochodne cząstkowe drugiego rzędu i budujemy macierz drugiej różniczki $\displaystyle d_{(x,y)}^{2}f$

$\displaystyle \left[{\begin{array}{cc}12x^{2}-16&0\\0&12y^{2}-4\end{array}}\right].$

Ponieważ ta macierz w punkcie $\displaystyle (0,0)$ ma postać $\displaystyle \left[{\begin{array}{cc}-16&0\\0&-4\end{array}}\right]$, w $\displaystyle (0,\pm 1)$ postać $\displaystyle \left[{\begin{array}{cc}-16&0\\0&8\end{array}}\right]$, w $\displaystyle (\pm 2,0)$ postać $\displaystyle \left[{\begin{array}{cc}32&0\\0&-4\end{array}}\right]$, wreszcie w $\displaystyle (\pm 2,1)$ i $\displaystyle (\pm 2,-1)$ postać $\displaystyle \left[{\begin{array}{cc}32&0\\0&8\end{array}}\right]$, więc funkcja $\displaystyle f$ nie ma ekstremów w punktach $\displaystyle (0,-1),(0,1),(2,0),(-2,0)$, ale ma maksimum w punkcie $\displaystyle (0,0)$ i ma minima w punktach $\displaystyle (-2,-1),\displaystyle (-2,1),(2,-1),(2,1)$.

b) Łatwo wyliczamy punkty krytyczne $\displaystyle (0,0)$ i $\displaystyle \left({\frac {9}{4}},{\frac {3}{4}}\right)$. Macierz drugiej różniczki $\displaystyle d_{(x,y)}^{2}g$ ma postać $\displaystyle \displaystyle \left[{\begin{array}{cc}2&-6\\-6&48y\end{array}}\right]$. Funkcja $\displaystyle g$ ma tylko jedno ekstremum -- minimum w punkcie $\displaystyle \left({\frac {9}{4}},{\frac {3}{4}}\right)$.

c) Dla funkcji $\displaystyle h$ należy zrobić założenie $\displaystyle x\neq 0$ i $\displaystyle y\neq 0$. Łatwo wyliczyć, że jedynym punkt krytycznym jest $\displaystyle ({\frac {\sqrt[{3}]{2}}{2}},{\sqrt[{3}]{2}})$ i że w tym punkcie funkcja $\displaystyle h$ ma minimum (macierz drugiej różniczki $\displaystyle d_{(x,y)}^{2}h$ ma postać $\displaystyle \displaystyle \left[{\begin{array}{cc}{\frac {2}{x^{3}}}&2\\2&{\frac {4}{y^{3}}}\end{array}}\right]$).

Ćwiczenie 8.3.

Wyznaczyć ekstrema lokalne funkcji

a) $\displaystyle f(x,y)=e^{2x}(x+y^{2}+2y)$,

b) $\displaystyle g(x,y)=e^{x^{2}-y}(5-2x+y)$,

c) $\displaystyle h(x,y)=\ln |x+y|-x^{2}-y^{2}$,

d) $\displaystyle \displaystyle \phi (x,y)=x-2y+\ln {\sqrt {x^{2}+y^{2}}}+3\mathrm {arctg} \,{\frac {y}{x}}$.

Wskazówka

Należy poszukać punktów krytycznych i zastosować kryterium Sylvestera do macierzy drugiej różniczki.

a) Przy pochodnej cząstkowej $\displaystyle \displaystyle {\frac {\partial ^{2}f}{\partial x^{2}}}$ warto zauważyć, że jeden z jej składników jest równy $\displaystyle \displaystyle 2{\frac {\partial f}{\partial x}}$, zatem zeruje się w punktach krytycznych.

b) Skorzystać ze wskazówki do podpunktu a) przy wszystkich pochodnych cząstkowych drugiego rzędu.

d) Warto zapisać naszą funkcję w postaci $\displaystyle \phi (x,y)=x-2y+{\frac {1}{2}}\ln(x^{2}+y^{2})+3\mathrm {arctg} \,{\frac {y}{x}},$ łatwiej jest ją wtedy różniczkować.

Rozwiązanie

a) Warunek konieczny istnienia ekstremum sprowadza się do układu równań

$\displaystyle \left\{{\begin{array}{l}e^{2x}(2(x+y^{2}+2y)+1)=0\\e^{2x}(2y+2)=0\end{array}}\right.,$

którego rozwiązaniem jest tylko punkt $\displaystyle ({\frac {1}{2}},-1)$. Macierz drugiej różniczki $\displaystyle d_{(x,y)}^{2}f$ ma postać

$\displaystyle \left[{\begin{array}{cc}2e^{2x}(2(x+y^{2}+2y)+1)+2e^{2x}&2e^{2x}(2y+2)\\2e^{2x}(2y+2)&2e^{2x}\end{array}}\right].$

W naszym punkcie jest to macierz $\displaystyle \displaystyle \left[{\begin{array}{cc}2e&0\\0&2e\end{array}}\right]$, zatem funkcja $\displaystyle f$ ma w tym punkcie minimum.

b) Przekształcamy układ równań otrzymany z warunku koniecznego

$\displaystyle \left\{{\begin{array}{l}e^{x^{2}-y}[2x(5-2x+y)-2]=0\\e^{x^{2}-y}[-(5-2x+y)+1]=0\end{array}}\right.,$

zauważając, że z drugiego równania wynika, że $\displaystyle 5-2x+y=1$. Stąd z pierwszego równania $\displaystyle x=1$. Otrzymujemy jedyny punkt $\displaystyle (1,-2)$. Liczymy pochodne cząstkowe drugiego rzędu

$\displaystyle {\begin{aligned}{\frac {\partial ^{2}g}{\partial x^{2}}}=2xe^{x^{2}-y}[2x(5-2x+y)-2]+e^{x^{2}-y}[2(5-2x+y)-4x],\\{\frac {\partial ^{2}g}{\partial x\partial y}}=2xe^{x^{2}-y}[-(5-2x+y)+1]+2e^{x^{2}-y},\\{\frac {\partial ^{2}g}{\partial ^{2}}}=-e^{x^{2}-y}[-(5-2x+y)+1]-e^{x^{2}-y}.\end{aligned}}$

Pierwsze składniki każdej z nich zerują się w naszym punkcie krytycznym, zatem łatwo jest policzyć, że macierz drugiej różniczki $\displaystyle d_{(1,-2)}^{2}g$ ma postać $\displaystyle \left[{\begin{array}{cc}-2e^{3}&2e^{3}\\2e^{3}&-e^{3}\end{array}}\right]$. Stąd wnioskujemy, że $\displaystyle g$ nie ma ekstremum.

c) Dziedziną funkcji $\displaystyle h$ jest zbiór $\displaystyle \mathbb {R} ^{2}\setminus \{(x,y):y=-x\}$. Z warunku koniecznego istnienia ekstremum otrzymujemy układ równań

$\displaystyle \left\{{\begin{array}{l}{\frac {1}{x+y}}-2x=0\\{\frac {1}{x+y}}-2y=0\end{array}}\right..$

W szczególności $\displaystyle x=y$, co po podstawieniu do pierwszego równania daje nam punkty $\displaystyle ({\frac {1}{2}},{\frac {1}{2}})$ i $\displaystyle (-{\frac {1}{2}},-{\frac {1}{2}})$. Macierz drugiej różniczki $\displaystyle d_{(x,y)}^{2}h$ ma postać

$\displaystyle \left[{\begin{array}{cc}-{\frac {1}{(x+y)^{2}}}-2&-{\frac {1}{(x+y)^{2}}}\\-{\frac {1}{(x+y)^{2}}}&-{\frac {1}{(x+y)^{2}}}-2\end{array}}\right].$

Stąd widać, że w obu punktach nasza funkcja ma maksimum.

d) Funkcja $\displaystyle \phi$ jest zdefiniowana poza prostą $\displaystyle x=0$. Warunek konieczny daje nam układ równań

$\displaystyle \left\{{\begin{array}{l}1+{\frac {x-3y}{x^{2}+y^{2}}}=0\\-2+{\frac {y+3x}{x^{2}+y^{2}}}=0\end{array}}\right..$

Redukując wyrażenie $\displaystyle x^{2}+y^{2}$, otrzymujemy $\displaystyle y+3x=-2(x-3y)$, czyli $\displaystyle y=x$. Wracając pierwszego równania otrzymujemy jedno rozwiązanie $\displaystyle (1,1)$. Macierzą drugiej różniczki $\displaystyle d_{(x,y)}^{2}\phi$ jest

$\displaystyle \left[{\begin{array}{cc}{\frac {y^{2}-x^{2}+6xy}{(x^{2}+y^{2})^{2}}}&{\frac {3y^{2}-3x^{2}-2xy}{(x^{2}+y^{2})^{2}}}\\{\frac {3y^{2}-3x^{2}-2xy}{(x^{2}+y^{2})^{2}}}&{\frac {x^{2}-y^{2}-6xy}{(x^{2}+y^{2})^{2}}}\end{array}}\right].$

W naszym punkcie macierz ta przyjmuje postać $\displaystyle \displaystyle \left[{\begin{array}{cc}{\frac {3}{2}}&-{\frac {1}{2}}\\-{\frac {1}{2}}&-{\frac {3}{2}}\end{array}}\right]$, zatem $\displaystyle \phi$ nie ma ekstremum.

Ćwiczenie 8.4.

Wyznaczyć ekstrema lokalne funkcji

a) $\displaystyle f(x,y)=\sin {x}\sin {y}\sin(x+y)$,

b) $\displaystyle h(x,y)=\sin {x}+\cos {y}+\cos(x-y)$
w zbiorze $\displaystyle (0,\pi )^{2}$.

Wskazówka

a) Warto pamiętać, że $\displaystyle \sin {\alpha }\cos \beta +\cos \alpha \sin \beta =\sin(\alpha +\beta )$.

Rozwiązanie

a) Mamy do rozwiązania układ równań

$\displaystyle \left\{{\begin{array}{l}0=\sin {y}(\cos {x}\sin(x+y)+\sin {x}\cos(x+y))=\sin {y}\sin(2x+y)\\0=\sin {x}(\cos {y}\sin(x+y)+\sin {y}\cos(x+y))=\sin {x}\sin(2y+x)\end{array}}\right..$

Ponieważ $\displaystyle x,y\in (0,\pi )$, zatem $\displaystyle \sin {y}\neq 0$ oraz $\displaystyle \sin(2x+y)=0$ wtedy i tylko wtedy, gdy $\displaystyle 2x+y=\pi$ lub $\displaystyle 2x+y=2\pi$. Wyliczamy stąd $\displaystyle y$ i wstawiamy do drugiego równania, w którym również zerować może się tylko drugi czynnik. Jeśli $\displaystyle y=\pi -2x$, to $\displaystyle 0=\sin(2\pi -3x)=-\sin {3x}$, czyli $\displaystyle x=\pi /3$ lub $\displaystyle x=2\pi /3$. Jeśli $\displaystyle y=2\pi -2x$, otrzymujemy te same punkty. Wobec założenia o dziedzinie punktami krytycznymi są $\displaystyle \left({\frac {\pi }{3}},{\frac {\pi }{3}}\right)$ i $\displaystyle \left({\frac {2\pi }{3}},{\frac {2\pi }{3}}\right)$. Macierzą drugiej różniczki $\displaystyle d_{(x,y)}^{2}f$ jest

$\displaystyle \left[{\begin{array}{cc}2\sin {y}\cos(2x+y)&\sin(2x+2y)\\\sin(2x+2y)&2\sin {x}\cos(x+2y)\end{array}}\right].$

Łatwo sprawdzić, że $\displaystyle f$ ma w $\displaystyle \left({\frac {\pi }{3}},{\frac {\pi }{3}}\right)$ maksimum i w $\displaystyle \left({\frac {2\pi }{3}},{\frac {2\pi }{3}}\right)$ minimum.

b) Tym razem należy rozwiązać układ

$\displaystyle \left\{{\begin{array}{l}0=\cos {x}-\sin(x-y)\\0=-\sin {y}+\sin(x-y)\end{array}}\right..$

Wynika stąd, że $\displaystyle \sin {y}=\cos {x}=\sin({\frac {\pi }{2}}-x)$. Ponieważ $\displaystyle x,y\in (0,\pi )$, więc $\displaystyle y={\frac {\pi }{2}}-x$ lub $\displaystyle y=\pi -({\frac {\pi }{2}}-x)={\frac {\pi }{2}}+x$. Otrzymujemy stąd jeden punkt krytyczny $\displaystyle \left({\frac {\pi }{3}},{\frac {\pi }{6}}\right)$, w którym funkcja $\displaystyle h$ osiąga maksimum.

Ćwiczenie 8.5.

Wyznaczyć ekstrema lokalne funkcji

a) $\displaystyle f(x,y)=1-{\sqrt {x^{2}+y^{2}}}$,

b) $\displaystyle g(x,y)={\sqrt[{5}]{x^{4}+y^{4}}}$,

c) $\displaystyle h(x,y)=x^{5}+y^{5}$.
Czy otrzymane ekstrema są też globalne?

Wskazówka

Należy poszukać punktów krytycznych.

a) b) Jak wyglądają poziomice danych funkcji? c) Jak wygląda zbiór zer danej funkcji i jak dzieli on dziedzinę na obszary, na których funkcja ta jest dodatnia lub ujemna?

Rozwiązanie

a) Dla naszej funkcji nie istnieją pochodne cząstkowe pierwszego rzędu w środku układu współrzędnych, a tam gdzie istnieją, nie zerują się. Zatem jedynym kandydatem na ekstremum jest punkt $\displaystyle (0,0)$. Zauważmy, że $\displaystyle f(0,0)=1$ i $\displaystyle {\sqrt {x^{2}+y^{2}}}>0$ dla dowolnego punktu $\displaystyle (x,y)$ na płaszczyźnie różnego od środka układu współrzędnych. W szczególności dla dowolnego $\displaystyle (x,y)\neq (0,0)$ mamy nierówność $\displaystyle f(x,y)<1$, co oznacza, że $\displaystyle f$ ma maksimum globalne w $\displaystyle (0,0)$. Warto także zauważyć, że wykres funkcji $\displaystyle f$ -- powierzchnia stożkowa -- powstaje przez obrót wykresu funkcji $\displaystyle z=\phi (x)=1-|x|=1-{\sqrt {x^{2}}}$ dookoła osi $\displaystyle 0z$.
wykres

b) Podobnie jak w poprzednim punkcie funkcja $\displaystyle g$ ma niezerowe pochodne cząstkowe pierwszego rzędu poza środkiem układu współrzędnych, gdzie te pochodne w ogóle nie istnieją. Tym razem $\displaystyle g(0,0)=0$, a dla $\displaystyle (x,y)\neq (0,0)$ wartość $\displaystyle g(x,y)$ jest dodatnia. Zatem w punkcie $\displaystyle (0,0)$ funkcja $\displaystyle g$ ma globalne minimum.
wykres

c) Pochodne cząstkowe pierwszego rzędu funkcji $\displaystyle h$ zerują się tylko w punkcie $\displaystyle (0,0)$, jednakże tym razem funkcja $\displaystyle h$ nie ma ekstremum w punkcie $\displaystyle (0,0)$. Mamy bowiem $\displaystyle h(0,0)=0$, $\displaystyle h(a,0)=a^{5}>0$ dla $\displaystyle a>0$ i $\displaystyle h(a,0)=a^{5}<0$ dla $\displaystyle a<0$, zatem dowolnie blisko środka układu współrzędnych funkcja przyjmuje i wartości dodatnie i ujemne, zatem i mniejsze i większe od wartości w tym punkcie.
wykres

Ćwiczenie 8.6.

a) Pokazać, że funkcja $\displaystyle f(x,y)=(1+e^{x})\cos {y}+xe^{x}$ ma nieskończenie wiele minimów, natomiast nie ma żadnego maksimum.

b) Pokazać, że funkcja $\displaystyle f(x,y)=3x^{4}-4x^{2}y+y^{2}$ nie ma minimum w punkcie $\displaystyle (0,0)$, ale jej zacieśnienie do dowolnej prostej przechodzącej przez początek układu współrzędnych ma silne minimum w tym punkcie.

Wskazówka

a) Należy poszukać punktów krytycznych i zastosować kryterium Sylvestera do macierzy drugiej różniczki. Warto pamiętać, że $\displaystyle F(x)=e^{x}$ przyjmuje tylko wartości dodatnie.

b) Należy zbadać znak funkcji na osiach układu współrzędnych i w punktach postaci $\displaystyle (a,2a^{2})$. Jak wygląda zacieśnienie funkcji do prostej przechodzącej przez początek układu współrzędnych?

Rozwiązanie

a) Warunek konieczny istnienia ekstremum sprowadza się do układu

$\displaystyle \left\{{\begin{array}{l}e^{x}(\cos {y}+1+x)=0\\-(1+e^{x})\sin {y}=0\end{array}}\right..$

Z drugiego równania wynika, że $\displaystyle y=k\pi$ dla pewnego $\displaystyle k\in \mathbb {Z}$. Jeśli $\displaystyle k$ jest parzyste, to z pierwszego równania $\displaystyle x=-2$, jeśli nieparzyste, to $\displaystyle x=0$. Tworzymy macierz drugiej różniczki $\displaystyle d_{(x,y)}^{2}f$

$\displaystyle \left[{\begin{array}{cc}e^{x}(\cos {y}+2+x)&-e^{x}\sin {y}\\-e^{x}\sin {y}&-(1+e^{x})\cos {y}\end{array}}\right].$

Niech $\displaystyle m$ będzie dowolną liczbą całkowitą. W punkcie $\displaystyle (0,(2m+1)\pi )$ rozważana powyżej macierz ma postać $\displaystyle \left[{\begin{array}{cc}1&0\\0&2\end{array}}\right],$ natomiast w punkcie $\displaystyle (-2,2m\pi )$ postać $\displaystyle \left[{\begin{array}{cc}e^{-2}&0\\0&-(1+e^{-2})\end{array}}\right],$ zatem funkcja $\displaystyle f$ ma minimum w każdym punkcie postaci $\displaystyle (0,(2m+1)\pi )$, a nie ma ekstremum w żadnym z punktów postaci $\displaystyle (-2,2m\pi )$.

b) Zauważmy, że $\displaystyle f(0,0)=0$ oraz $\displaystyle f(0,b)=b^{2}>0$ dla dowolnej niezerowej liczby $\displaystyle b$. Z drugiej strony $\displaystyle f(a,2a^{2})=3a^{4}-8a^{4}+4a^{4}=-a^{4}<0$ dla dowolnej niezerowej liczby $\displaystyle a$. Z tych dwóch faktów funkcja nie może mieć minimum w swoim miejscu zerowym $\displaystyle (0,0)$, bo dowolnie blisko tego miejsca przyjmuje zarówno wartości dodatnie, jak i ujemne. Widzimy, że zawężenie funkcji $\displaystyle f$ do prostej $\displaystyle x=0$, czyli funkcja $\displaystyle g(y)=f(0,y)=y^{2}$, ma globalne minimum w punkcie $\displaystyle 0$. Podobnie dla dowolnego $\displaystyle m\in \mathbb {R}$ zawężenie funkcji $\displaystyle f$ do prostej $\displaystyle y=mx$, czyli funkcja $\displaystyle h_{m}(x)=f(x,mx)=3x^{4}-4mx^{3}+m^{2}x^{2}$, ma minimum w punkcie $\displaystyle 0$ (zob. ćwiczenia z Analizy matematycznej I do modułu 10).

Ćwiczenie 8.7.

Wyznaczyć ekstrema lokalne funkcji

a) $\displaystyle f(x,y,z)=x^{4}-y^{3}+2z^{3}-2x^{2}+6y^{2}-3z^{2}$,

b) $\displaystyle g(x,y,z)=x^{3}+xy+y^{2}-2zx+2z^{2}+3y-1$,

c) $\displaystyle h(x,y,z)=xyz(4-x-y-z)$.

Wskazówka

Należy poszukać punktów krytycznych i zastosować kryterium Sylvestera do macierzy drugiej różniczki.

a) Warto skorzystać ze wskazówki ćwiczenia 8.1. a).

c) Zwróćmy uwagę, że funkcja $\displaystyle h$ zeruje się na czterech płaszczyznach: $\displaystyle x=0$, $\displaystyle y=0$, $\displaystyle z=0$ i $\displaystyle x+y+z=4$. Najpierw należy pokazać, że w żadnym z punktów trzech pierwszych płaszczyzn funkcja $\displaystyle h$ nie osiąga ekstremum, bo dowolnie blisko każdego takiego punktu przyjmuje zarówno wartości dodatnie lub ujemne. (Ten fakt jest też prawdziwy dla ostatniej płaszczyzny, ale sprawdzenie tego nie jest konieczne, co będzie widoczne w dalszym postępowaniu). Następnie szukamy punktów krytycznych pod założeniem $\displaystyle x\neq 0,y\neq 0,z\neq 0$.

Rozwiązanie

a) Każda z pochodnych cząstkowych rzędu pierwszego funkcji $\displaystyle f$ zależy tylko od tej zmiennej, względem której jest liczona. Z warunku koniecznego istnienia ekstremum otrzymujemy układ trzech niezależnych równań $\displaystyle 4x^{3}-4x=0$, $\displaystyle -3y^{2}+12y=0$ i $\displaystyle 6z^{2}-6z=0$. Punkty krytyczne zatem to $\displaystyle (0,0,0),\displaystyle (0,0,1),\displaystyle (1,0,0),\displaystyle (1,0,1),\displaystyle (-1,0,0),\displaystyle (-1,0,1),\displaystyle (0,4,0),\displaystyle (0,4,1),\displaystyle (1,4,0),\displaystyle (1,4,1)$, $\displaystyle (-1,4,1)$, $\displaystyle (-1,4,1)$. Macierz drugiej różniczki $\displaystyle d_{(x,y,z)}^{2}f$ ma postać

$\displaystyle \left[{\begin{array}{ccc}12x^{2}-4&0&0\\0&-6y+12&0\\0&0&12z-6\end{array}}\right].$

Wobec tego w punkcie $\displaystyle (0,0,0)$ macierzą tą jest $\displaystyle \left[{\begin{array}{ccc}-4&0&0\\0&12&0\\0&0&-6\end{array}}\right]$, w $\displaystyle (0,0,1)$ - $\displaystyle \left[{\begin{array}{ccc}-4&0&0\\0&12&0\\0&0&6\end{array}}\right]$, w $\displaystyle (\pm 1,0,0)$ - $\displaystyle \left[{\begin{array}{ccc}8&0&0\\0&12&0\\0&0&-6\end{array}}\right]$, w $\displaystyle (0,4,0)$ - $\displaystyle \left[{\begin{array}{ccc}-4&0&0\\0&-12&0\\0&0&-6\end{array}}\right]$, w $\displaystyle (0,4,1)$ - $\displaystyle \left[{\begin{array}{ccc}-4&0&0\\0&-12&0\\0&0&6\end{array}}\right]$, w $\displaystyle (\pm 1,4,0)$ - $\displaystyle \left[{\begin{array}{ccc}8&0&0\\0&-12&0\\0&0&-6\end{array}}\right]$, wreszcie w $\displaystyle (\pm 1,4,1)$ - $\displaystyle \left[{\begin{array}{ccc}8&0&0\\0&-12&0\\0&0&6\end{array}}\right]$. Stąd widać na mocy kryterium Sylvestera, że funkcja $\displaystyle f$ ma minima w punktach $\displaystyle (1,0,1)$ i $\displaystyle (-1,0,1)$ i maksimum w punkcie $\displaystyle (0,4,0)$ oraz, że są to jedyne ekstrema tej funkcji.

b) Warunek konieczny istnienia ekstremum prowadzi do układu

$\displaystyle \left\{{\begin{array}{l}3x^{2}+y-2z=0\\x+2y+3=0\\-2x+4z=0\end{array}}\right.,$

którego rozwiązaniami są dwie trójki liczb $\displaystyle (-{\frac {1}{2}},-{\frac {5}{4}},-{\frac {1}{4}})$ i $\displaystyle (1,-2,{\frac {1}{2}})$. Macierz drugiej różniczki $\displaystyle d_{(x,y,z)}^{2}g$ ma postać

$\displaystyle \left[{\begin{array}{ccc}6x&1&-2\\1&2&0\\-2&0&4\end{array}}\right].$

Ponieważ

$\displaystyle {\rm {det}}\left[{\begin{array}{ccc}-3&1&-2\\1&2&0\\-2&0&4\end{array}}\right]=-36\quad {\rm {i}}\quad {\rm {det}}\left[{\begin{array}{cc}-3&1\\1&2\end{array}}\right]=-7,$

funkcja $\displaystyle g$ nie ma ekstremum w punkcie $\displaystyle (-{\frac {1}{2}},-{\frac {5}{4}},-{\frac {1}{4}})$, natomiast wobec

$\displaystyle {\rm {det}}\left[{\begin{array}{ccc}6&1&-2\\1&2&0\\-2&0&4\end{array}}\right]=36\quad {\rm {i}}\quad {\rm {det}}\left[{\begin{array}{cc}6&1\\1&2\end{array}}\right]=11$

funkcja $\displaystyle g$ ma minimum w punkcie $\displaystyle (1,-2,{\frac {1}{2}})$.

c) Funkcja $\displaystyle h$ zeruje się na czterech płaszczyznach: $\displaystyle x=0$, $\displaystyle y=0$, $\displaystyle z=0$ i $\displaystyle x+y+z=4$. Pokażmy najpierw, że w żadnym punkcie pierwszych trzech z nich nie ma ekstremum. Weźmy punkt $\displaystyle (x_{0},y_{0},z_{0})$ leżący na płaszczyźnie $\displaystyle x=0$ oraz zdefiniujmy funkcję $\displaystyle s(x,y,z)=yz(4-x-y-z)$. Mamy $\displaystyle x_{0}=0$. Ponieważ częścią wspólną każdych dwóch z naszych płaszczyzn jest tylko prosta, więc dowolnie blisko punktu $\displaystyle (x_{0},y_{0},z_{0})$ możemy znaleźć taki punkt $\displaystyle (x_{1},y_{1},z_{1})$, że $\displaystyle x_{1}=0$ oraz $\displaystyle s(x_{1},y_{1},z_{1})\neq 0$. Niech np. $\displaystyle s(x_{1},y_{1},z_{1})>0$ (drugi przypadek jest symetryczny). Z ciągłości funkcji $\displaystyle s$ dla dostatecznie małej liczby dodatniej $\displaystyle \delta$ zachodzi $\displaystyle s(\delta ,y_{1},z_{1})>0$ oraz $\displaystyle s(-\delta ,y_{1},z_{1})>0$ (bo $\displaystyle x_{1}=0$). Ale wtedy $\displaystyle h(\delta ,y_{1},z_{1})=\delta s(\delta ,y_{1},z_{1})>0$ oraz $\displaystyle h(-\delta ,y_{1},z_{1})=-\delta s(-\delta ,y_{1},z_{1})<0$, zatem funkcja $\displaystyle h$ nie ma minimum w punkcie $\displaystyle (x_{0},y_{0},z_{0})$ (bo jest to miejsce zerowe, a dowolnie blisko tego miejsca funkcja $\displaystyle h$ przyjmuje zarówno wartości dodatnie, jak i ujemne). Analogicznie postępujemy z punktami z płaszczyzn $\displaystyle y=0$ i $\displaystyle z=0$.

Wobec tego wystarczy poszukać punktów krytycznych pod założeniem $\displaystyle x\neq 0,y\neq 0,z\neq 0$. Wtedy warunek konieczny istnienia ekstremum prowadzi do układu Cramera

$\displaystyle \left\{{\begin{array}{l}2x+y+z=4\\x+2y+z=4\\x+y+2z=4\end{array}}\right.,$

którego rozwiązaniem jest jedna trójka liczb $\displaystyle (1,1,1)$. Macierz drugiej różniczki $\displaystyle d_{(x,y,z)}^{2}h$ ma postać

$\displaystyle \left[{\begin{array}{ccc}-2yz&z(4-2x-2y-z)&y(4-2x-y-2z)\\z(4-2x-2y-z)&-2xz&x(4-x-2y-2z)\\y(4-2x-y-2z)&x(4-x-2y-2z)&-2xy\end{array}}\right].$

Ponieważ

$\displaystyle {\rm {det}}\left[{\begin{array}{ccc}-2&-1&-1\\-1&-2&-1\\-1&-1&-2\end{array}}\right]=-4\quad {\rm {i}}\quad {\rm {det}}\left[{\begin{array}{cc}-2&-1\\-1&-2\end{array}}\right]=3$

funkcja $\displaystyle h$ ma maksimum w punkcie $\displaystyle (1,1,1)$. Jest to jedyne ekstremum tej funkcji.

Ćwiczenie 8.8.

a) Wyznaczyć ekstrema lokalne funkcji

$\displaystyle f(x,y,z)=4-x^{2}-{\frac {y}{x}}-{\frac {z^{2}}{y}}-{\frac {1}{z}}.$

b) Wyznaczyć ekstrema lokalne funkcji

$\displaystyle \Phi (x,y,z)=\sin(x+y+z)-\sin {x}-\sin {y}-\sin {z}$

w zbiorze

$\displaystyle (0,\pi )^{2}$.

Wskazówka

Należy poszukać punktów krytycznych i zastosować kryterium Sylvestera do macierzy drugiej różniczki.

b) W wyliczaniu punktów krytycznych przyda się tożsamość trygonometryczna - wzór na różnicę cosinusów.

Rozwiązanie

a) Zakładamy, że $\displaystyle x,y,z\neq 0$. Otrzymany z warunku koniecznego układ równań

$\displaystyle \left\{{\begin{array}{l}-2x+{\frac {y}{x^{2}}}=0\\-{\frac {1}{x}}+{\frac {z^{2}}{y^{2}}}=0\\-2{\frac {z}{y}}+{\frac {1}{z^{2}}}=0\end{array}}\right.$

ma jedyne rozwiązanie - punkt $\displaystyle \left({\frac {\sqrt[{3}]{2}}{2}},{\frac {1}{2}},{\frac {\sqrt[{3}]{2}}{2}}\right)$. Liczymy pochodne cząstkowe drugiego rzędu

$\displaystyle {\begin{aligned}{\frac {\partial ^{2}f}{\partial x^{2}}}=-2-2{\frac {y}{x^{3}}},\quad {\frac {\partial ^{2}f}{\partial y^{2}}}=-2{\frac {z^{2}}{y^{3}}},\quad {\frac {\partial ^{2}f}{\partial z^{2}}}=-2{\frac {1}{y}}-2{\frac {1}{z^{3}}},\\{\frac {\partial ^{2}f}{\partial x\partial y}}={\frac {1}{x^{2}}},\quad {\frac {\partial ^{2}f}{\partial y\partial z}}=2{\frac {z}{y^{2}}},\quad {\frac {\partial ^{2}f}{\partial z\partial x}}=0\end{aligned}}$

i budujemy macierz drugiej różniczki $\displaystyle d_{\left({\frac {\sqrt[{3}]{2}}{2}},{\frac {1}{2}},{\frac {\sqrt[{3}]{2}}{2}}\right)}^{2}f$ w tym punkcie, która ma postać

$\displaystyle A=\left[{\begin{array}{ccc}-6&2{\sqrt[{3}]{2}}&0\\2{\sqrt[{3}]{2}}&-4{\sqrt[{3}]{4}}&4{\sqrt[{3}]{2}}\\0&4{\sqrt[{3}]{2}}&-12\end{array}}\right].$

Mamy det$\displaystyle A=-144{\sqrt[{3}]{4}}<0,$ det$\displaystyle \left[{\begin{array}{cc}-6&2{\sqrt[{3}]{2}}\\2{\sqrt[{3}]{2}}&-4{\sqrt[{3}]{4}}\end{array}}\right]=20{\sqrt[{3}]{4}}>0$ oraz $\displaystyle -6<0$. Zatem z kryterium Sylvestera funkcja $\displaystyle f$ ma maksimum w punkcie $\displaystyle \left({\frac {\sqrt[{3}]{2}}{2}},{\frac {1}{2}},{\frac {\sqrt[{3}]{2}}{2}}\right)$.

b) Otrzymujemy układ równań

$\displaystyle \left\{{\begin{array}{l}\cos(x+y+z)=\cos {x}\\\cos(x+y+z)=\cos {y}\\\cos(x+y+z)=\cos {z}\end{array}}\right..$

W szczególności $\displaystyle \cos {x}=\cos y=\cos z$, czyli $\displaystyle x=y=z$, ponieważ $\displaystyle x,y,z\in (0,\pi )$. Zatem $\displaystyle 0=\cos {3x}-\cos {x}=-2\sin {2x}\sin {x}$, a stąd $\displaystyle x={\frac {\pi }{2}}$. Macierz drugiej różniczki $\displaystyle d_{(x,y,z)}^{2}\Phi$ ma postać

$\displaystyle \left[{\begin{array}{ccc}\sin {x}-\sin(x+y+z)&-\sin(x+y+z)&-\sin(x+y+z)\\-\sin(x+y+z)&\sin {y}-\sin(x+y+z)&-\sin(x+y+z)\\-\sin(x+y+z)&-\sin(x+y+z)&\sin {z}-\sin(x+y+z)\end{array}}\right].$

Ponieważ

$\displaystyle {\rm {det}}\left[{\begin{array}{ccc}2&1&1\\1&2&1\\1&1&2\end{array}}\right]=4\quad {\rm {i}}\quad {\rm {det}}\left[{\begin{array}{cc}2&1\\1&2\end{array}}\right]=3,$

funkcja $\displaystyle \Phi$ ma w punkcie $\displaystyle ({\frac {\pi }{2}},{\frac {\pi }{2}},{\frac {\pi }{2}})$ minimum.

Ćwiczenie 8.9.

(Zadanie Huygensa) Pomiędzy liczby dodatnie $\displaystyle a$ i $\displaystyle b$ ($\displaystyle a\leq b$) wstawić liczby dodatnie $\displaystyle x_{1},...,x_{n}$ tak, aby ułamek

$\displaystyle f(x_{1},...,x_{n})={\frac {x_{1}x_{2}...x_{n}}{(a+x_{1})(x_{1}+x_{2})...(x_{n-1}+x_{n})(x_{n}+b)}}$

miał największą wartość.

Wskazówka

Przyrównując pochodne cząstkowe pierwszego rzędu do zera należy pamiętać, że dziedziną naszej funkcji jest przedział $\displaystyle (0,+\infty )$. Wychodząc z układu równań otrzymanego z warunku koniecznego istnienia ekstremum, proszę wyrazić zależność między kolejnymi punktami $\displaystyle a,x_{1},x_{2},...,x_{n},b$ za pomocą liczby $\displaystyle q={\frac {x_{1}}{a}}$. Jakiego rodzaju jest to zależność?

By pokazać, że funkcja osiąga maksimum globalne w punkcie krytycznym, rozważamy najpierw prosty przypadek $\displaystyle n=1$ (mamy wtedy funkcję jednej zmiennej rzeczywistej). Następnie, jeśli $\displaystyle n>1$ ustalamy dowolne $\displaystyle n-1$ zmiennych i rozważamy zacieśnienie naszej funkcji do półprostej -- mamy wtedy znowu funkcję jednej zmiennej. Co o niej możemy powiedzieć? Jak z tego wywnioskować naszą tezę?

Rozwiązanie

Zauważmy, że $\displaystyle f$ jest dobrze zdefiniowaną funkcją $\displaystyle n$ zmiennych dodatnich. Jeśli przyjmiemy oznaczenia

$\displaystyle x'=(x_{2},...,x_{n})\quad {\rm {i}}\quad p(x')=(x_{2}+x_{3})...(x_{n-1}+x_{n})(x_{n}+b),$

to licznik pochodnej cząstkowej funkcji $\displaystyle f$ po $\displaystyle x_{1}$ wyraża się wzorem

$\displaystyle {\begin{aligned}x_{2}...x_{n}(a+x_{1})(x_{1}+x_{2})p(x')-x_{1}x_{2}...x_{n}[(x_{1}+x_{2})p(x')+(a+x_{1})p(x')]=\\=x_{2}...x_{n}p(x')[(a+x_{1})(x_{1}+x_{2})-x_{1}(x_{1}+x_{2}+a+x_{1})]=x_{2}...x_{n}p(x')(ax_{2}-x_{1}^{2}),\end{aligned}}$

zatem zeruje się wtedy i tylko wtedy, gdy $\displaystyle ax_{2}-x_{1}^{2}=0$. Postępując analogicznie dla pozostałych zmiennych, uzyskamy układ równań

$\displaystyle \left\{{\begin{array}{l}ax_{2}-x_{1}^{2}=0\\x_{1}x_{3}-x_{2}^{2}=0\\\vdots \\x_{n-2}x_{n}-x_{n-1}^{2}=0\\x_{n-1}b-x_{n}^{2}=0\end{array}}\right..$

Przekształcając te równania i porównując je, otrzymamy zależność

$\displaystyle {\frac {x_{1}}{a}}={\frac {x_{2}}{x_{1}}}={\frac {x_{3}}{x_{2}}}=...={\frac {x_{n-1}}{x_{n-2}}}={\frac {x_{n}}{x_{n-1}}}={\frac {b}{x_{n}}},$

co oznacza, że ciąg $\displaystyle a,x_{1},x_{2},...,x_{n},b$ jest geometryczny o ilorazie $\displaystyle \displaystyle q={\frac {x_{1}}{a}}$. W konsekwencji $\displaystyle b=aq^{n+1}$ i stąd $\displaystyle q={\sqrt[{n+1}]{\frac {b}{a}}}$. Punktem krytycznym jest zatem

$\displaystyle M({\sqrt[{n+1}]{a^{n}b}},{\sqrt[{n+1}]{a^{n-1}b^{2}}},{\sqrt[{n+1}]{a^{n-2}b^{3}}},...,{\sqrt[{n+1}]{a^{2}b^{n-1}}},{\sqrt[{n+1}]{ab^{n}}})$

oraz

$\displaystyle {\begin{aligned}f(M)={\frac {\sqrt[{n+1}]{a^{1+2+...+n}b^{1+2+...+n}}}{{\sqrt[{n+1}]{a^{n}}}({\sqrt[{n+1}]{a}}+{\sqrt[{n+1}]{b}}){\sqrt[{n+1}]{a^{n-1}b}}({\sqrt[{n+1}]{a}}+{\sqrt[{n+1}]{ab}})...{\sqrt[{n+1}]{b^{n}}}({\sqrt[{n+1}]{a}}+{\sqrt[{n+1}]{b}})}}=\\={\frac {1}{({\sqrt[{n+1}]{a}}+{\sqrt[{n+1}]{b}})^{n+1}}}.\end{aligned}}$

Wystarczy jeszcze pokazać, że dla dowolnego punktu $\displaystyle P(x_{1},...,x_{n})\in (0,+\infty )^{n}$ różnego od punktu $\displaystyle M$ zachodzi $\displaystyle f(M)>f(P)$. Zanim to udowodnimy, zauważmy, że (co wynika jeszcze z układu równań uzyskanego z warunku koniecznego ekstremum) punkt

$\displaystyle P(x_{1},...,x_{k})=M$ wtedy i tylko wtedy, gdy $\displaystyle x_{1}={\sqrt {ax_{2}}},x_{2}={\sqrt {x_{1}x_{3}}},...,x_{n-1}={\sqrt {x_{n-2}x_{n}}},x_{n}={\sqrt {x_{n-1}b}}.$

Rozważmy teraz przypadek $\displaystyle n=1$. Szukamy wtedy maximum funkcji $\displaystyle \displaystyle h(x)={\frac {x}{(a+x)(x+b)}}$ jednej zmiennej dodatniej $\displaystyle x$. Z powyższego rozumowania wiemy, że punktem krytycznym jest punkt $\displaystyle {\sqrt {ab}}$. Chcemy teraz pokazać, że $\displaystyle h(x)<h({\sqrt {ab}})$ dla dowolnego $\displaystyle x\in (0,{\sqrt {ab}})\cup ({\sqrt {ab}},+\infty )$. Można to udowodnić, wykorzystując rachunek różniczkowy jednej zmiennej rzeczywistej (zachęcamy do tego ćwiczenia jako przypomnienia z Analizy matematycznej I), ale można też zrobić to bardziej elementarnie. Mamy mianowicie $\displaystyle (x-{\sqrt {ab}})^{2}\geq 0$ dla dowolnej liczby rzeczywistej $\displaystyle x$, a równość zachodzi dokładnie wtedy, gdy $\displaystyle x={\sqrt {ab}}$. Przekształcając tę nierówność, otrzymujemy kolejno $\displaystyle x^{2}+ab\geq 2{\sqrt {ab}}x$, a stąd $\displaystyle ax+x^{2}+ab+bx\geq (a+2{\sqrt {ab}}+b)x$, czyli $\displaystyle (a+x)(x+b)\geq ({\sqrt {a}}+{\sqrt {b}})^{2}x$, co jest równoważne nierówności $\displaystyle h(x)\leq h({\sqrt {ab}})$ i równość zachodzi dokładnie wtedy, gdy $\displaystyle x={\sqrt {ab}}$, co dowodzi naszej tezy w tym przypadku.

Teraz, jeśli $\displaystyle n$ jest dowolną liczbą naturalną większą od 1, ustalmy dowolną liczbę $\displaystyle k\in \{1,2,...,n\}$ oraz $\displaystyle n-1$ dowolnie wybranych liczb dodatnich $\displaystyle x_{1},...,x_{k-1},x_{k+1},...,x_{n}$ i rozważmy funkcję

$\displaystyle g(x)=f(x_{1},x_{2},...,x_{k-1},x,x_{k+1},...,x_{n}).$

Zauważmy, że jest to funkcja $\displaystyle h$ z poprzedniego przypadku, dla przedziału $\displaystyle (x_{k-1},x_{k+1})$ (lub $\displaystyle (x_{k+1},x_{k-1})$, jeśli liczba $\displaystyle x_{k+1}$ jest mniejsza od $\displaystyle x_{k-1}$), pomnożona przez stałą dodatnią. Zatem z poprzedniego rozumowania funkcja $\displaystyle g$ osiąga silne maksimum w punkcie $\displaystyle x_{k}={\sqrt {x_{k-1}x_{k+1}}}$. Zatem ogólnie funkcja $\displaystyle f$ osiąga silne maksimum w punkcie $\displaystyle P(x_{1},...,x_{k})$, dokładnie wtedy, gdy zachodzą związki

$\displaystyle x_{1}={\sqrt {ax_{2}}},x_{2}={\sqrt {x_{1}x_{3}}},...,x_{n-1}={\sqrt {x_{n-2}x_{n}}},x_{n}={\sqrt {x_{n-1}b}},$

czyli dokładnie wtedy, gdy $\displaystyle P=M$.