Chapter 4: Partial Differentiation
Section 4.10: Optimization on Closed Domains
Example 4.10.4
Find the extreme values of the function fx,y=x2−x⁢y+y2+4⁢y+10/10 on the domain R consisting of the interior and boundary of the disk whose radius is 4, and whose center is at the origin.
Solution
Mathematical Solution
The unconstrained optimization problem requires solving the equations fx=fy=0, that is, the equations
15⁢x−110⁢y=0,−110⁢x+15⁢y+25=0
for x,y=−4/3,−8/3, a point interior to the disk. At this point, f=7/15. See Figure 4.10.4(a) where the points −4/3,−8/3,0 and −4/3,−8/3,7/15 are shown as black dots.
For the constrained problem, describe the boundary in polar coordinates xt=4 cost, yt=4 sint, t∈0,2 π, and find the extrema of
Ft=fxt,yt=13−8 cos⁡t⁢sin⁡t+8⁢sin⁡t/5
The distinct critical numbers are t=0,2 π/3,4 π/3, as corroborated by the graph of Ft in Figure 4.10.4(a).
use plots, plottools in module() local p1,p2,p3,p4,f,d; f:=(x^2-x*y+y^2+4*y+10)*(1/10); d:=sqrt(16-x^2); p1:=plot3d(f,x=-4..4,y=-d..d,shading=xy); p2:=plot3d(0,x=-4..4,y=-d..d,color=green); p3:=pointplot3d([[-4/3,-8/3,0],[-4/3,-8/3,7/15]],symbol=solidsphere,symbolsize=20,color=black); p4:=display(p1,p2,p3,labels=[x,y,z],tickmarks=[4,4,5],style=surfacecontour,lightmodel=none,orientation=[-55,75,0],axes=frame); print(p4); end module: end use:
Figure 4.10.4(a) Surface z=f over disk
plot((13-8*cos(t)*sin(t)+8*sin(t))/5,t=0..2*Pi,0..5);
Figure 4.10.4(b) Ft=fxt,yt
Table 4.10.4(a) lists the boundary extrema in Cartesian coordinates.
f4,0=13/5
f−2,23=13+63/5≐4.6785
f−2,−23=13+63/5≐0.5215
Table 4.10.4(a) Boundary extrema
Therefore, the absolute minimum is f−4/3,−8/3=7/15≐0.47, and the absolute maximum is f−2,23 = 13+63/5 ≐ 4.6785.
Maple Solution - Interactive
Initialize
Tools≻Load Package: Student Calculus 1
Loading Student:-Calculus1
Context Panel: Assign Function
fx,y=x2−x⁢y+y2+4⁢y+10/10→assign as functionf
Numerically find any unconstrained extrema
Context Panel: Optimization≻Minimize (local)
fx,y
110⁢x2−110⁢x⁢y+110⁢y2+25⁢y+1
→minimize
0.466666666666667,x=−1.33333333333333,y=−2.66666666666667
Analytically find any unconstrained extrema
Calculus palette: Partial-derivative operator Press the Enter key.
Context Panel: Solve≻Solve
Context Panel: Evaluate and Display Inline
∂∂ x fx,y=0,∂∂ y fx,y=0
x5−y10=0,−x10+y5+25=0
→solve
x=−43,y=−83
f−4/3,−8/3 = 715
Optimization Assistant: Optimize numerically subject to the boundary constraint
The Optimization Assistant is no longer listed in Tools≻Assistants. It is now found in Tools≻Tutors≻Optimization
Form the sequence fx,y,x2+y2=16 and press the Enter key. Context Panel: Optimization≻Optimization Assistant Use the to obtain extrema on the boundary.
Figure 4.10.4(c) shows the assistant finding a maximum on the circumference of the circle.
Figure 4.10.4(d) shows a graph of the surface determined by f, with the circular boundary of R lifted to this surface.
Select "Minimize" and/or add additional constraints on x and y to find other extrema on the boundary.
fx,y,x2+y2=16
110⁢x2−110⁢x⁢y+110⁢y2+25⁢y+1,x2+y2=16
Figure 4.10.4(c) Constrained maximum
Figure 4.10.4(d) Surface z=f and constraint
Apply the boundary constraint analytically
Write fxt,yt, where xt=4 cost,yt=4 sint describe the boundary of R; press the Enter key.
Context Panel: Simplify≻Simplify
Context Panel: Differentiate≻With Respect To≻t
Context Panel: Student Calculus1≻Solve≻Find Roots See Figure 4.10.4(e)
Context Panel: Evaluate and Display Inline (Evaluate fxt,yt at each distinct root.)
Figure 4.10.4(e) Roots dialog
f4 cost,4 sint
8⁢cos⁡t25−8⁢sin⁡t⁢cos⁡t5+8⁢sin⁡t25+8⁢sin⁡t5+1
= simplify
135+−8⁢cos⁡t+8⁢sin⁡t5
→differentiate w.r.t. t
8⁢sin⁡t25+−8⁢cos⁡t+8⁢cos⁡t5
→roots
0,2⁢π3,4⁢π3,2⁢π
f4,0 = 135→at 10 digits2.600000000
f4 cos2 π/3,4 sin2 π/3 = 135+6⁢35→at 10 digits4.678460970
f4 cos4 π/3,4 sin4 π/3 = 135−6⁢35→at 10 digits0.521539030
Maple Solution - Coded
Define the function f.
f≔x,y→x2−x⁢y+y2+4⁢y+10/10:
Find any unconstrained extrema
Obtain the equations fx=fy=0 via the diff command, and their solution via the solve command.
solvedifffx,y,x=0,difffx,y,y=0
Find extrema on the circumference of the disk
Representing the circumference of the disk in polar coordinates, use the diff command to obtain the derivative with respect to t.
Set that derivative equal to zero and, for the resulting equation, use the Roots command in the Student Calculus1 package to find all solutions in the interval t∈0,2 π.
Student:-Calculus1:-Rootsdifff4 cost,4 sint,t=0,t=0..2 π
Evaluate f at each distinct solution.
f4,0 = 135
f4 cos2 π/3,4 sin2 π/3 = 135+6⁢35
f4 cos4 π/3,4 sin4 π/3 = 135−6⁢35
<< Previous Example Section 4.10 Next Example >>
© Maplesoft, a division of Waterloo Maple Inc., 2024. All rights reserved. This product is protected by copyright and distributed under licenses restricting its use, copying, distribution, and decompilation.
For more information on Maplesoft products and services, visit www.maplesoft.com
Download Help Document