derivadas parciales ejercicios

FUNDAMENTOS MATEMÁTICOS DE LA INGENIERÍA Ingeriería Técnica Industrial. Especialidad en Mecánica. Boletin 6. Funciones d

Views 161 Downloads 1 File size 206KB

Report DMCA / Copyright

DOWNLOAD FILE

Recommend stories

Citation preview

FUNDAMENTOS MATEMÁTICOS DE LA INGENIERÍA Ingeriería Técnica Industrial. Especialidad en Mecánica. Boletin 6. Funciones de Varias Variables Curso 2003-2004

EJERCICIOS RESUELTOS 1. En cada apartado, calcular el límite y discutir la continuidad de la función de la que se calcula el límite ¡ ¢ x+y x √ (a) lim (c) lim x + 3y 2 (b) lim x+y (x,y)→(2,1) (x,y)→(2,4) x − y (x,y)→(1,1) ³ ´ arcsen xy y sen xy (f) lim (d) lim exy (e) lim 1 + xy (x,y)→(0,1) (x,y)→(0,0) (x,y)→( π ,2 ) 4 Solución: Todas las funciones que aparecen en cada uno de los apartados de este ejercicio son continuas en sus respectivos dominios de definición por ser composición, suma, producto y/o cociente de funciones continuas. Así, el cálculo de todos los límites se realiza por simple sustitución, ya que el punto donde deben ser calculado dichos límites siempre pertenece al dominio de la función. (a)

lim

(x,y)→(2,1)

x + 3y 2 = 2 + 3 · 12 = 5

Dominio: R2 .

2+4 x+y = = −3 Dominio: D = {(x, y) ∈ R2 : x 6= y}. 2−4 (x,y)→(2,4) x − y x 1 √ (c) lim Dominio: D = {(x, y) ∈ R2 : x + y > 0}. =√ x+y (x,y)→(1,1) 2 arcsen(x/y) (d) lim =0 Dominio: D = {(x, y) ∈ R2 : y 6= 0, |x/y| ≤ 1, 1 + xy 6= 0}. (x,y)→(0,1) 1 + xy (b)

(e) (f)

lim

lim

(x,y)→(π/4,2)

lim

(x,y)→(0,0)

y sen xy = 2 sen π/2 = 2

exy = 1

Dominio: R2 .

Dominio: R2 .

2. Calcular las derivadas parciales primeras y segundas de las siguientes funciones: (a) z = tg(2x − y) xy (d) w = x+y+z Solución:

x

(b) z = xe y

(c) z = x ln(xy) 1 (f) w = ln(xyz 2 ) (e)w = p 1 − x2 − y 2 − z 2

(a) z = tg(2x − y). Las derivadas parciales de z después de simplificar y factorizar quedan de la siguiente forma. Se ha utilizado que (sec x)0 = sec x tg x. Derivadas parciales de primer orden: ∂z ∂z = 2 sec2 (2x − y) = − sec2 (2x − y) ∂x ∂y Derivadas parciales de segundo orden: ∂2z ∂2z 2 = 8 sec (2x − y) tg(2x − y) = 2 sec2 (2x − y) tg(2x − y) ∂x2 ∂y 2 ∂2z ∂2z = = −4 sec2 (2x − y) tg(2x − y) ∂x∂y ∂y∂x x

(b) z = xe y Derivadas parciales de primer orden: ¶ µ x x 1 x ∂z x y y y 1+ = 1 · e + xe =e ∂x y y

x ∂z = xe y ∂y

1

µ

−x y2



=−

x2 xy e y2

Boletín 6. Ejercicios Resueltos. Fundamentos Matemáticos de la Ingeniería Mecánica. Curso 2003-04

2

Derivadas parciales de segundo orden: µ ¶ ¶ µ ¶ ¶ µ µ 2¶ µ x 1 x 1 x x ∂2z x 2 ∂2z −x x2 xy x 2x2 xy x x y y y y 1 + +e e = 2 + = e = e + e = e + − ∂x2 y y y y y2 ∂y 2 y3 y2 y2 y3 y ¶ µ ¶ ¶ µ µ 2 2 x x ∂ z x xy x x ∂ z x −x = = − 2e + 1+ ey = − 2 ey 2 + ∂x∂y ∂y∂x y y y2 y y

(c) z = x ln(xy) = x (ln x + ln y) = x ln x + x ln y Derivadas parciales de primer orden: ∂z 1 ∂z x = ln x + x + ln y = 1 + ln x + ln y = ∂x x ∂y y Derivadas parciales de segundo orden: 1 ∂2z ∂2z ∂2z 1 ∂2z −x = = 2 = = 2 2 ∂x x ∂y y ∂x∂y ∂y∂x y xy (d) w = . Las derivadas parciales de w después de simplificar y factorizar quedan de la siguiente x+y+z forma. Derivadas parciales de primer orden: y+z x+z y ∂w ∂w ∂w =y =x = −x 2 2 ∂x ∂y ∂z (x + y + z) (x + y + z) (x + y + z)2 Derivadas parciales de segundo orden: ∂2w ∂2w y+z ∂2w x+z y = −2y = = −2x = 2x 3 3 3 2 2 ∂x2 ∂y ∂z (x + y + z) (x + y + z) (x + y + z) ∂2w 2xy + xz + yz + z 2 ∂2w x−y−z ∂2w ∂2w = = = =y 3 3 ∂x∂y ∂y∂x ∂x∂z ∂z∂x (x + y + z) (x + y + z) ∂2w x−y+z ∂2w = = −x 3 ∂y∂z ∂z∂y (x + y + z) ¢−1/2 ¡ 1 (e) w = p = 1 − x2 − y 2 − z 2 . 1 − x2 − y 2 − z 2 Derivadas parciales de primer orden: ¢−3/2 ¢−3/2 ¡ x ∂w −1 ¡ 1 − x2 − y 2 − z 2 (−2x) = x 1 − x2 − y 2 − z 2 =q = ∂x 2 3 (1 − x2 − y 2 − z 2 ) ¢−3/2 ¢−3/2 ¡ ∂w y −1 ¡ 1 − x2 − y 2 − z 2 (−2y) = y 1 − x2 − y 2 − z 2 =q = ∂y 2 (1 − x2 − y 2 − z 2 )3 ¢−3/2 ¢−3/2 ¡ z ∂w −1 ¡ 1 − x2 − y 2 − z 2 (−2z) = z 1 − x2 − y 2 − z 2 =q = ∂z 2 (1 − x2 − y 2 − z 2 )3 Derivadas parciales de segundo orden: ¶ µ ¢ ¢ 2x2 + 1 − y 2 − z 2 ∂2w ¡ −3 ¡ 2 2 2 −3/2 2 2 2 −5/2 q 1 − x = 1 − x + x (−2x) = − y − z − y − z ∂x2 2 5 (1 − x2 − y 2 − z 2 ) ¶ µ ¢ ¢ ¡ ∂2w 2y 2 + 1 − x2 − z 2 −3 ¡ 2 2 2 −3/2 2 2 2 −5/2 q 1 − x = 1 − x + y (−2y) = − y − z − y − z ∂y 2 2 (1 − x2 − y 2 − z 2 )5 ¶ µ ¢ ¢ 2z 2 + 1 − x2 − y 2 ∂2w ¡ −3 ¡ 2 2 2 −3/2 2 2 2 −5/2 q 1 − x = 1 − x + z (−2z) = − y − z − y − z ∂z 2 2 (1 − x2 − y 2 − z 2 )5 ∂2w 3xy ∂2w = =q ∂x∂y ∂y∂x 5 (1 − x2 − y 2 − z 2 )

∂2w 3yz ∂2w = =q ∂y∂z ∂z∂y (1 − x2 − y 2 − z 2 )5 ¢ ¡ (f) w = ln xyz 2 = ln x + ln y + 2 ln z

∂2w 3xz ∂2w = =q ∂x∂z ∂z∂x 5 (1 − x2 − y 2 − z 2 )

Boletín 6. Ejercicios Resueltos. Fundamentos Matemáticos de la Ingeniería Mecánica. Curso 2003-04

3

Derivadas parciales de primer orden: ∂w ∂w ∂w 1 1 2 = = = ∂x x ∂y y ∂z z Derivadas parciales de segundo orden: ∂2w 1 1 1 ∂2w ∂2w = − = − =− 2 2 2 2 2 2 ∂x x ∂y y ∂z z ∂2w ∂2w ∂2w ∂2w ∂2w ∂2w = = = = = =0 ∂x∂y ∂y∂x ∂x∂z ∂z∂x ∂y∂z ∂z∂y 1 1 + + xy, calcular todos los valores de x e y para los que se verifica que x y fx (x, y) = 0 y fy (x, y) = 0 simultáneamente.

3. Dada la función f (x, y) = Solución:

En primer lugar, obtenemos las derivadas parciales de primer orden de f : fx (x, y) = −

1 + y, x2

fy (x, y) = −

1 + x. y2

 1   − 2 + y = 0, x De la primera de las ecuaciones Ahora, resolvemos el sistema de ecuaciones (no lineales)   − 1 + x = 0. y2 1 obtenemos y = 2 y la sustitución en la segunda ecuación nos conduce a −x4 + x = 0 ⇐⇒ x(−x3 + 1) = x 0 ⇐⇒ x = 0, x = 1. Pero, nótese que los puntos de la recta x = 0 no pertenecen al dominio de f, por lo 1 que los valores pedidos son x = 1 e y = 2 = 1. 1 4. a) Comprobar que la función z = ex sen y satisface la ecuación de Laplace

∂2z ∂2z + = 0. ∂x2 ∂y 2

x b) Comprobar que la función z = e−t sen , donde c es una constante no nula, satisface la ecuación del c ∂2z ∂z − c2 2 = 0. calor ∂t ∂x Solución: (a) z = ex sen y. Calculemos las derivadas parciales segundas ∂2z ∂z = = z, 2 ∂x ∂x

∂z ∂z = ex sen y = z, = ex cos y, ∂x ∂y ∂2z ∂2z Así, + 2 = z − z = 0. 2 ∂x ∂y

∂2z ∂2z y . 2 ∂x ∂y 2 2 ∂ z = −ex sen y = −z. ∂y 2

∂z ∂ 2 z . y ∂t ∂x2 ∂z e−t z ∂z e−t ∂2z = − 2 sen (x/c) = − 2 . = −e−t sen (x/c) = −z, = cos (x/c) , 2 ∂t ∂x c ∂x c c ∂z ∂2z z De esta forma, − c2 2 = −z + c2 2 = −z + z = 0. ∂t ∂x c

(b) z = e−t sen (x/c) . Calculemos las derivadas parciales

v2 5. La aceleración centrípeta de una partícula que se mueve por un círculo es a = , donde v es la velocidad r y r es el radio del círculo. Estimar el máximo error porcentual en la aceleración si las cotas de error en v y r son del 2 por 100 y del 1 por 100, respectivamente. Solución: A partir de los datos del problema tenemos ∆a ' da =

dv ∆v dr ∆r = = ±0.02, = = ±0.01. Por consiguiente, v v r r

v2 ∂a 2v ∂a dv + dr = dv − 2 dr, ∂v ∂r r r

Boletín 6. Ejercicios Resueltos. Fundamentos Matemáticos de la Ingeniería Mecánica. Curso 2003-04

de donde teniendo en cuenta de nuevo que a = da ∆a ' = a a

2v r v2 r

v2 r2 v2 r

dv −

dr = 2

4

v2 , se tiene r dv dr − = 2 (±0.02) − (±0.01) = ±0.05. v r

Es decir, el máximo error porcentual en la medida de la aceleración centrípeta es del 5%. 6. Un leñador puede cometer un error del 5% al medir el diámetro y un 3% al medir la altura de un tronco cilíndrico, recto y circular. Estimar el máximo error porcentual en el volumen del tronco. Solución: El volumen del tronco es V =

π 2 D H, donde D es diámetro del tronco y H es la altura. Así, 4

∆V ' dV =

∂V ∂V π π dD + dH = DHdD + D2 dH ∂D ∂H 2 4

y por tanto, dV ∆V ' = V V

π 2 DH dD π 2 4D H

+

π 2 4D dH π 2 4D H

=2

dD dH + = 2 (±0.05) + (±0.03) = ±0.13. D H

Luego, el error porcentual al medir el volumen del tronco puede ser hasta del 1.3%. dz , siendo z = x sen y y x = e−t e y = π − t. dt dw b) Usando la regla de la cadena, hallar , siendo w = xy cos z, x = t, y = t2 y z = arcsen t. dt d2 w c) Usando la regla de la cadena, hallar 2 , siendo w = arccotan(2xy), x = cos t, y = sen t. Calcular, dt d2 w cuando t = 0. además, dt2 Solución:

7. a) Usando la regla de la cadena, hallar

(a) z = x sen y, x = e−t , y = π − t. dz ∂z dx ∂z dy = · + · = [sen y · (−e−t ) + x cos y · (−1)]|x=e−t ,y=π−t = dt ∂x dt ∂y dt = −e−t sen (π − t) − e−t cos (π − t) = e−t (cos t − sen t) .

(b) w = xy cos z, x = t, y = t2 , z = arcsen t. ¸¯ · ¯ ∂w dx ∂w dy ∂w dz 1 dw ¯ = = · + · + · = y cos z · 1 + x cos z · (2t) − xy sen z · √ 2 dt ∂x dt ∂y dt ∂z dt 1 − t ¯x=t,y=t2 ,z=arcsen t √ t4 3t2 − 4t4 t4 = 3t2 1 − t2 − √ = √ . = t2 cos (arcsen t) + 2t2 cos (arcsen t) − √ 2 2 1−t 1−t 1 − t2 (c) w = arccotg(2xy), x = cos t, y = sen t. d2 w En este caso, vamos a obtener la expresión de 2 aplicando la regla de la cadena de forma sucesiva dt y sin utilizar las expresiones explicítitas de w, x e y. Al final sustituiremos dichas expresiones. ∂w dx ∂w dy dw = · + · , luego Aplicando la regla de la cadena sabemos que dt ∂x dt ∂y dt µ ¶ µ ¶ d dw d ∂w dx ∂w dy d2 w = = · + · dt2 dt dt dt ∂x dt ∂y dt y aplicando la derivada para la suma y el producto se obtiene µ ¶ µ ¶ µ ¶ d ∂w d ∂w dx ∂w d2 x dy ∂w d2 y d ∂w dx ∂w dy . · + · = · + · 2 + · + · dt ∂x dt ∂y dt dt ∂x dt ∂x dt dt ∂y dt ∂y dt2

Boletín 6. Ejercicios Resueltos. Fundamentos Matemáticos de la Ingeniería Mecánica. Curso 2003-04

Ahora, utilizando la regla de la cadena para las funciones d dt

µ

∂w ∂x



=

∂ ∂x

µ

d dt

µ

∂w ∂x



·

dx ∂ + dt ∂y

µ

∂w ∂y



∂ = ∂x

µ

∂w ∂y



dx ∂ · + dt ∂y

µ

∂w ∂w y se tiene ∂x ∂y

∂w ∂x



·

dy ∂ 2 w dx ∂ 2 w dy = + · · 2 dt ∂x dt ∂y∂x dt

∂w ∂y



·

dy ∂ 2 w dx ∂ 2 w dy = · + . · dt ∂x∂y dt ∂y 2 dt

Por consiguiente, imponiendo la igualdad para las derivadas parciales cruzadas conseguimos ¶ ¶ µ 2 µ 2 dx ∂w d2 x dy ∂w d2 y ∂ w dx ∂ w dx ∂ 2 w dy d2 w ∂ 2 w dy = + = · · + · · + · · + · + · dt2 ∂x2 dt ∂y∂x dt dt ∂x dt2 ∂x∂y dt ∂x2 dt dt ∂y dt2 ¶ ¶ µ µ 2 2 ∂ 2 w dy dx ∂w d2 x ∂ 2 w dx dy ∂ 2 w ∂w d2 y ∂2w dx dy + + + = = · · + · · · + · · · ∂x2 dt ∂y∂x dt dt ∂x dt2 ∂x∂y dt dt ∂x2 dt ∂y dt2 ¶ ¶ µ µ 2 2 ∂ 2 w dy dx ∂w d2 x ∂ 2 w ∂w d2 y ∂2w dx dy + 2 · + + . = · · + · · · · ∂x2 dt ∂y∂x dt dt ∂x dt2 ∂y 2 dt ∂y dt2 Si ahora tenemos en cuenta que w = arccotg(2xy), x = cos t, y = sin t sen t, se obtiene ∂w 2y , =− ∂x 1 + 4x2 y 2 ∂2w xy 3 = 16 2, ∂x2 (1 + 4x2 y 2 ) dx = − sen t, dt

∂w 2x , =− ∂y 1 + 4x2 y 2

∂2w x3 y = 16 2, ∂y 2 (1 + 4x2 y 2 ) dy d2 x = − cos t, = cos t, dt2 dt

∂2w −1 + 4x2 y 2 =2 2 ∂y∂x (1 + 4x2 y 2 ) d2 y = − sen t, dt2

y por consiguiente, d2 w dt2

"

xy 3 −1 + 4x2 y 2 16 (− sen t)2 + 2 · 2 · cos t · (− sen t) + 2 2 y2 ) 2 y 2 )2 (1 + 4x (1 + 4x ¶ µ x3 y 2y + − · (− cos t) + +16 · (cos t)2 + 2 2 2 y 2 )2 1 + 4x y (1 + 4x ¶ µ ¸¯ ¯ 2x ¯ + − = · (− sen t) ¯ 1 + 4x2 y 2 x=cos t,y=sen t sen5 t −1 + 4 cos2 t sen 2 t = 16 cos t −4 2 2 cos t sen t+ (1 + 4 cos2 t sen 2 t) (1 + 4 cos2 t sen 2 t) ¢ ¡ sen t sen t cos t + 16 cos5 t +4 2 = 2 2 1 + 4 cos t sen t (1 + 4 cos2 t sen 2 t) 2 sen 4 t + 1 + 2 cos4 t = 8 cos t sen t (1 + 4 cos2 t sen 2 t)2 =

Ahora, para t = 0, tenemos

d2 w (0) = 0. dt2

∂z ∂z y usando la regla de la cadena, siendo z = x2 − y 2 y x = s cos t e y = t sen s. ∂t ∂s Solución: ¯ ∂z ∂z ∂x ∂z ∂y ¯¯ = 2x(−s sen t) − 2y sen s|x=s cos t,y=t sen s = = · + · ∂t ∂x ∂t ∂y ∂t ¯x=s cos t,y=t sen s

8. Hallar

= 2s cos s(−s sen t) − 2t sen s sen s = −2s2 cos t sen t − 2t sen2 s = −s2 sen(2t) − 2t sen2 s. ¯ ∂z ∂z ∂x ∂z ∂y ¯¯ = 2x cos t − 2yt cos s|x=s cos t,y=t sen s = = · + · ∂s ∂x ∂s ∂y ∂s ¯ 2

2

x=s cos t,y=t sen s 2 2

= 2s cos t − 2t cos s sen s = 2s cos t − t sen(2s).

5

Boletín 6. Ejercicios Resueltos. Fundamentos Matemáticos de la Ingeniería Mecánica. Curso 2003-04

6

9. a) Calcular las derivadas parciales de z por derivación implícita en los siguientes apartados: (a.1) x ln y + y 2 z + z 2 = 8.

(a.2) z = ex sen(y + z).

b) Calcular las derivadas parciales de w por derivación implícita en los siguientes apartados: √ √ (b.1) xyz + xzw − yzw + w2 = 5. (b.2) w − x − y − y − w + z = 0.

Solución:

(a) (a.1) Derivamos con respecto a x la identidad x ln y + y 2 z + z 2 = 8. ln y ∂z ∂z ∂z ln y + y 2 + 2z = 0, luego =− 2 ∂x ∂x ∂x y + 2z Si derivamos respecto de y se tiene x x ∂z ∂z x + 2y 2 z ∂z y + 2yz + 2yz + y 2 + 2z = 0 y por tanto, =− 2 =− y ∂y ∂y ∂y y + 2z y (y 2 + 2z) x (a.2) Repetimos el proceso para z = e sen(y + z). ∂z ∂z ∂z Derivando respecto de x se obtiene = ex sen(y + z) + ex cos(y + z) = z + ex cos(y + z) . Por ∂x ∂x ∂x ∂z z consiguiente, = . x ∂x 1 − e cos(y + z) ¶ µ ∂z ex cos(y + z) ∂z ∂z x = e cos(y + z) 1 + , luego = . Al derivar con respecto a y obtenemos ∂y ∂x ∂y 1 − ex cos(y + z)

(b) (b.1) Ahora obtenemos las derivadas parciales de w respecto de x, y, z siendo xyz+xzw−yzw+w2 = 5. ∂w ∂w ∂w Derivamos con respecto a x y obtenemos yz + zw + xz − yz + 2w = 0, y por tanto, ∂x ∂x ∂x yz + zw ∂w =− . ∂x xz − yz + 2w Derivando con respeto a y conseguimos xz + xz

∂w ∂w ∂w − zw − yz + 2w = 0, luego ∂y ∂y ∂y

∂w −xz + zw = . ∂y xz − yz + 2w Derivando con respeto a z se tiene xy + xw + xz

∂w ∂w ∂w − yw − yz + 2w = 0 y de aquí ∂z ∂z ∂z

xy + xw − yw ∂w =− . ∂z xz − yz + 2w

√ √ x − y − y − w + z = 0, obtenemos: µ ¶ ∂w 1 1 ∂w = 0, y por tanto, Derivando con respecto a x se consigue − √ − √ − ∂x 2 x−y 2 y−w ∂x √ ∂w y−w √ =√ . ∂x x − y (2 y − w + 1)

(b.2) De forma análoga para w −

µ ¶ ∂w 1 1 ∂w 1− = 0, luego + √ − √ ∂y 2 x−y 2 y−w ∂y √ √ ∂w − y−w+ x−y √ =√ . ∂y x − y (2 y − w + 1)

Derivando respecto de y obtenemos

µ ¶ 1 ∂w ∂w + 1 = 0 y por consiguiente, − √ − Ahora, al derivar respecto de z se obtiene ∂z 2 y−w ∂z √ ∂w y−w = −2 √ . ∂z 2 y−w+1

Boletín 6. Ejercicios Resueltos. Fundamentos Matemáticos de la Ingeniería Mecánica. Curso 2003-04

7

10. Determinar la derivada direccional de la función f en el punto dado y en la dirección indicada: (a) f (x, y) = xy,

P (1, 1), 2

(b) f (x, y, z) = ln(1 + xyz ),

u = i + j. P (2, 1, 0),

Solución:

u = −i + j + k.

µ ¶ 1 1 u = √ ,√ y conseguimos (a) Hacemos el vector u = i + j = (1, 1) unitario, tomando v = kuk 2 2 µ ¶ √ 1 1 2 Dv f (1, 1) = ∇f (1, 1) · v = (1, 1) · √ , √ = √ = 2. 2 2 2 µ ¶ −1 1 u 1 = √ ,√ ,√ (b) Tomamos v = . El gradiente de f es kuk 3 3 3 ¶ µ xz 2 2xyz yz 2 , , , ∇f (x, y, z) = (fx (x, y, z), fy (x, y, z), fz (x, y, z)) = 1 + xyz 2 1 + xyz 2 1 + xyz 2 luego Dv f (2, 1, 0) = ∇f (2, 1, 0) · v = (0, 0, 0) · v = 0.

π 11. a) Hallar la derivada direccional de f (x, y) = cos(x + y) en el punto P (0, π) en la dirección de Q( , 0). 2 b) Hallar la derivada direccional de f (x, y, z) = ln(x + y + z) en el punto P (1, 0, 0) en la dirección de Q(4, 3, 1). Solución: − → (a) Obsérvese que ∇f (x, y) = (− sen(x + y), − sen(x + y)) , luego ∇f (0, π) = 0 = (0, 0) y por tanto, − → Du f (0, π) = ∇f (0, π) · u = 0 · u = 0 para todo vector u unitario de R2 . −→ PQ −→ (b) Calculamos el vector unitario u en la dirección de P Q = (4, 3, 2) − (1, 0, 0) = (3, 3, 2), u = −→ = kP Qk ¶ µ 1 3 3 2 √ ,√ ,√ . Teniendo en cuenta que ∇f (x, y, z) = (1, 1, 1) , se tiene Du f (1, 0, 0) = x+y+z 22 22 22 ¶ µ 3 3 2 3 2 8 3 =√ +√ +√ =√ . ∇f (1, 0, 0) · u = (1, 1, 1) · √ , √ , √ 22 22 22 22 22 22 22 12. La temperatura en los puntos (x, y, z) de una pieza metálica está dada por T (x, y, z) = e2x+y+3z . ¿En qué dirección desde el origen se incrementa más rapidamente la temperatura? Solución: Como sabemos, viene dada por el vector gradiente, luego, ya que ¢ ¡ la dirección de máximo crecimiento ∇T (x, y, z) = 2e2x+y+3z , e2x+y+3z , 3e2x+y+3z , tenemos que el vector que marca el mayor crecimiento es ∇T (0, 0, 0) = (2, 1, 3). 2

13. a) Calcular el gradiente de g(x, y) = ye−x y el valor máximo de la derivada direccional en el punto P (0, 5). b) Calcular el gradiente de g(x, y, z) = xeyz y el valor máximo de la derivada direccional en el punto P (2, 0, −4). Solución:

´ ³ 2 2 (a) Puesto que ∇g(x, y) = −2xye−x , e−x , obtenemos que el valor máximo de la derivada direccional en el punto P (0, 5) viene dado por k∇g(0, 5)k = k(0, 1)k = 1.

(b) Análogamente, y utilizando que ∇g(x, y, z) = (eyz , xzeyz , xyeyz ) , el√valor máximo de la derivada direccional en el punto P (2, 0, −4) es k∇g(2, 0, −4)k = k(1, −8, 0)k = 65. 14. Hallar un vector normal y la recta tangente a la curva de nivel de f (x, y) = xy que pasa por el punto P (−1, 3). Solución:

Sabemos que el vector normal viene dado por el gradiente. Así, el vector normal es n = ∇f (−1, 3) = (3, −1) y la recta tangente viene dada, en forma implícita por, 3(x + 1) + (−1)(y − 3) = 0.

Boletín 6. Ejercicios Resueltos. Fundamentos Matemáticos de la Ingeniería Mecánica. Curso 2003-04

8

15. a) Calcular la ecuación del plano tangente y la ecuación de la recta normal a la superficie de ecuación z = x3 − 3xy + y 3 en el punto P (1, 2, 3). b) Calcular la ecuación del plano tangente y la ecuación de la recta normal a la superficie de ecuación xy 2 + 3x − z 2 = 4 en el punto P (2, 1, −2). Solución:

(a) Utilizamos que la ecuación del plano tangente a la superficie dada por z = f (x, y) en el punto (x0 , y0 , f (x0 , y0 )) tiene por ecuación z=

∂f ∂f (x0 , y0 )(x − x0 ) + (x0 , y0 )(y − y0 ) + f (x0 , y0 ). ∂x ∂y

∂f ∂f (1, 2) = 3 · 12 − 3 · 2 = −3, (1, 2) = ∂x ∂y −3 · 1 + 3 · 22 = 9, z0 = f (1, 2) = 3, la ecuación del plano tangente es z = −3(x − 1) + 9(y − 2) + 3, es decir, −3x + 9y − z = 12.   x = 1 − 3t, y = 2 + 9t, con t ∈ R. La recta normal en forma paramétrica viene dada por  z = 3 − t, Puesto que z = f (x, y) = x3 − 3xy + y 3 , se tiene

(b) En este apratado la superficie esta dada en forma implícita mediante xy 2 + 3x − z 2 = 4, luego tomando F (x, y, z) = xy 2 + 3x − z 2 − 4, la ecuación del plano tangente en el punto (2, 1, −2) viene dada por ∂F ∂F ∂F (2, 1, −2)(x − 2) + (2, 1, −2)(y − 1) + (2, 1, −2)(z + 2) = 0. ∂x ∂y ∂z ∂F ∂F ∂F (2, 1, −2) = 12 + 3 = 4, (2, 1, −2) = 2 · 2 · 1 = 4, (2, 1, −2) = −2 · (−2) = 4, ∂x ∂y ∂z el plano tangente tiene por ecuación 4(x − 2) + 4(y − 1) + 4(z + 2) = 0, i.e., x + y + z = 1.   x = 2 + t, La recta normal tiene por ecuaciones paramétricas y = 1 + t, con t ∈ R.  z = −2 + t, Así, usando que

16. Encontrar y clasificar los puntos críticos de las siguientes funciones: a) f (x, y) = x2 + y 2 + 2x − 6y + 6.

b) f (x, y) = 2x2 + 2xy + y 2 + 2x − 3. c) z = xy.

−4x . + y2 + 1 Solución:

d) z =

x2

Para la resolución de cada  ∂funo de los apartados de este ejercicio debemos resolver, en primer lugar, el  (x, y) = 0,  ∂x sistema de ecuaciones que nos proporciona los puntos críticos de f.   ∂f (x, y) = 0, ∂y  ∂f  (x, y) = 2x + 2 = 0  ∂f ∂f ∂x (a) Como (x, y) = 2x + 2, (x, y) = 2y − 6, la única solución del sistema  ∂x ∂y  ∂f (x, y) = 2y − 6 = 0 ∂y es (x, y) = (−1, 3). Es decir, el único punto crítico de f es el punto (−1, 3). ∂2f ∂2f ∂2f (x, y) = 2, (x, y) = 2, (x, y) = 0. Puesto que, Ahora calculamos las derivadas parciales ∂x2 ∂y 2 ∂x∂y µ 2 ¶2 ∂2f ∂ f ∂2f ∂2f (−1, 3) = 2 > 0, (−1, 3) · (−1, 3) − = 4 > 0, (−1, 3) ∂x2 ∂x2 ∂y 2 ∂x∂y

Boletín 6. Ejercicios Resueltos. Fundamentos Matemáticos de la Ingeniería Mecánica. Curso 2003-04

9

resulta que el punto crítico (−1, 3) es un mínimo relativo y la función toma en este punto el valor f (−1, 3) = (−1)2 + 32 + 2(−1) − 6 · 3 + 6 = −4. Nótese que f (x, y) = (x + 1)2 + (y − 3)2 − 4, luego en el punto (−1, 3) la función f toma su mínimo absoluto.  ∂f  (x, y) = 4x + 2y + 2 = 0  ∂x y obtenemos (x, y) = (−5/2, 4). Puesto que (b) Resolvemos el sistema   ∂f (x, y) = 2x + 2y − 3 = 0 ∂y ∂2f (−5/2, 4) = 4 > 0, ∂x2

∂2f ∂2f (−5/2, 4) · 2 (−5/2, 4) − 2 ∂x ∂y

µ

¶2 ∂2f (−5/2, 4) = 4 · 2 − 22 = 4 > 0, ∂x∂y

se deduce que (−5/2, 4) es un mínimo relativo. (c) Es fácil ver que el único punto crítico de f (x, y) = xy es el origen. Como ∂2f ∂2f (0, 0) · (0, 0) − ∂x2 ∂y 2

µ

∂2f (0, 0) ∂x∂y

¶2

= 0 · 0 − 12 = −1 < 0,

el origen es un punto de silla. (d) Calculamos las derivadas parciales de f (x, y) = ∂f x2 − y 2 − 1 (x, y) = 4 2, ∂x (x2 + y 2 + 1)

x2

−4x , + y2 + 1 ∂f x (x, y) = 8 2y 2 ∂y (x + y 2 + 1)

 x2 − y 2 − 1    4 2 2 = 0, (x + y 2 + 1) y resolvemos el sistema obteniendo los dos puntos críticos (x, y) = (1, 0), (x, y) = x   8 y = 0,  (x2 + y 2 + 1)2 (−1, 0). Ahora, para establecer la naturaleza de los dos puntos críticos, calculamos las derivadas parciales de segundo orden: ∂2f x2 − 3y 2 − 3 (x, y) = −8x 3, ∂x2 (x2 + y 2 + 1)

∂2f −3y 2 + x2 + 1 (x, y) = 8x 3, ∂y 2 (x2 + y 2 + 1)

∂2f 3x2 − y 2 − 1 (x, y) = −8y 3 ∂x∂y (x2 + y 2 + 1)

A partir de aquí, tenemos 1−3 ∂2f (1, 0) = −8 = 2, ∂x2 8

2 ∂2f (1, 0) = 8 = 2, ∂y 2 8

∂2f (1, 0) = 0, ∂x∂y

por lo que es inmediato observar que el punto (1, 0) es un mínimo relativo. Análogamente, 1−3 ∂2f (−1, 0) = 8 = −2, 2 ∂x 8

2 ∂2f (−1, 0) = −8 = −2, 2 ∂y 8

∂2f (−1, 0) = 0, ∂x∂y

y el punto (−1, 0) es un máximo relativo. 17. a) Calcular los extremos absolutos de f (x, y) = 2x − 2xy + y 2 en la región del plano acotada por las gráficas de las funciones y = x2 e y = 1. b) Calcular los extremos absolutos de f (x, y) = x2 + 2xy + y 2 en la región del plano acotada R = © ª 2 2 (x, y) : x + y ≤ 8 . Solución:

Boletín 6. Ejercicios Resueltos. Fundamentos Matemáticos de la Ingeniería Mecánica. Curso 2003-04

10

(a) Nótese que el problema posee solución puesto que f (x, y) = 2x − 2xy + y 2 es continua en la región dada, que a su vez es cerrada y acotada. Calculemos, en primer lugar, los puntos críticos de f (x, y) = 2x − 2xy + y 2 en la región abierta que acota las gráficas y = x2 e y = 1.  ∂f  (x, y) = 2 − 2y = 0  ∂x es (x, y) = (1, 1), que no pertenece a La única solución del sistema   ∂f (x, y) = −2x + 2y = 0 ∂y dicha región abierta. Ahora calculamos los puntos críticos sobre la frontera. Puesto que la frontera está dada por la unión {(x, y) ∈ R2 : −1 ≤ x ≤ 1, y = 1} ∪ {(x, y) ∈ R2 : −1 ≤ x ≤ 1, y = x2 }, debemos analizar la restricción de f a los conjuntos {(x, y) ∈ R2 : −1 ≤ x ≤ 1, y = 1}, {(x, y) ∈ R2 : −1 ≤ x ≤ 1, y = x2 }. Tomando y = 1, tenemos f (x, 1) = 2 · x − 2 · x · 1 + 12 = 1. Si y = x2 , entonces f (x, x2 ) = 2x − 2x3 + x4 , por lo que debemos analizar lo extremos absolutos de la función g(x) = 2x − 2x3 + x4 , con x ∈ [−1, 1]. Ahora, calculamos¡ los puntos críticos de g en el ¢ intervalo abierto (−1, 1). Para ello, resolvemos la ecuación g 0 (x) = 2 1 − 3x2 + 2x3 = 0. La raíces del anterior polinomio son: −1/2, 1 (doble) y sólo −1/2 pertenece al intervalo abierto (−1, 1). Ahora evaluamos g en x = −1/2 y en los extremos del intervalo y obtenemos g(−1/2) = −11/16, g(−1) = g(1) = 1, por lo que el máximo absoluto de g se alcanza en x = 1 y su mínimo absoluto en x = 0. Luego el máximo absoluto de f restringida a {(x, y) ∈ R2 : −1 ≤ x ≤ 1, y = x2 } se alcanza en el punto (1, 1) y el mínimo absoluto de f restringida a {(x, y) ∈ R2 : −1 ≤ x ≤ 1, y = x2 } se toma en el punto (−1/2, 1/4). En resumen, el mínimo absoluto de f en la región acotada por las gráficas de las funciones y = x2 e y = 1 es −11/16 y se alcanza en el punto (−1/2, 1/4) y su máximo valor es uno, alcanzándose en cualquiera de los puntos del segmento {(x, y) ∈ R2 : −1 ≤ x ≤ 1, y = 1}.

2 2 2 (b) Nótese que el problema posee solución © ª puesto que f (x, y) = x + 2xy + y es continua en todo R y 2 2 el círculo R = (x, y) : x + y ≤ 8 es cerrado y acotado. Téganse en cuenta que f puede expresarse en la forma f (x, y) = (x + y)2 , luego f (x, y) ≥ 0 y se anula en los puntosª(x, y) ∈ R2 tales que y = −x. Así, el menor valor de f restringida a R = © (x, y) : x2 + y 2 ≤ 8 es cero y se alcanza en los puntos del segmento y = −x restringido al circulo R. Es decir, f sólo es nula sobre el segmento S = {(x, y) ∈ R2 : y = −x, x2 + y 2 ≤ 8} = {(x, y) ∈ R2 : y = −x, 2x2 ≤ 8} = {(x, y) ∈ R2 : y = −x, |x| ≤ 2}. Por otra parte, los críticos de f restringida a R son todos los puntos del segmento S. En  puntos ∂f  (x, y) = 2x + 2y = 0  ∂x efecto, el sistema posee infinitas soluciones, todas de la forma y = −x.   ∂f (x, y) = 2x + 2y = 0 ∂y Ahora vamos estudiar ª los extremos absolutos de la función f restringida a la circunferencia C = © (x, y) : x2 + y 2 = 8 . Realizaremos este estudio de dos formas distintas:

el comportamiento de la función las semicircunferencias de Primera forma: Analizamos √ √ √ f sobre √ ecuaciones y = 8 − x2 e y = − 8 − x2 , con x ∈ [−2 2, 2 2]. Así, debemos estudiar las √ √ √ ¢ ¡ ¢2 ¢ ¡ ¡ √ funciones g(x) = f x, 8 − x2 = x + 8 − x2 = 2x 8 − x2 + 8 y h(x) = f x, − 8 − x2 = √ √ √ √ ¡ ¢2 x − 8 − x2 = −2x 8 − x2 + 8, con x ∈ [−2 2, 2 2]. √ √ √ ¢ ¡ −2x = 0 nos lleva a 2 8 − x2 − Para x ∈ (−2 2, 2 2), la ecuación g 0 (x) = 2 8 − x2 + 2x √ 2 2 ¡ 8√ −x ¢ ¡ √ ¢ 2x2 = 0, de donde x = ±2. Como g(2) = 16, g(−2) = 0, g 2 2 = g −2 2 = 8, el máximo absoluto de g es 16, alcanzándose en x = 2, y el mínimo √ absoluto es cero, tomándose en √ x = −2. Los corespondientes valores de y son: x = 2 =⇒ y = 8 − x2 = 2, x = 2 =⇒ y = − 8 − x2 = −2. Un análisis similar puede hacerse sobre la función h para deducir que el máximo absoluto de h es 16, alcanzándose en x = −2, y el mínimo absoluto es cero, tomándose en x√= 2. Los √ corespondientes valores de y son: x = 2 =⇒ y = − 8 − x2 = −2, x = −2 =⇒ y = − 8 − x2 = −2. Nótese que en los puntos donde g y h se anulan (x =√−2 y x = 2), los correspondiente valores √ para f (x = −2 =⇒ y = 8 − x2 = 2, x = 2, y = − 8 − x2 = −2) son (−2, 2) y (2, −2) que

Boletín 6. Ejercicios Resueltos. Fundamentos Matemáticos de la Ingeniería Mecánica. Curso 2003-04

11

pertenecen, como no podía ser de otra forma, al segmento S = {(x, y) ∈ R2 : y = −x, |x| ≤ 2}, donde la función f toma el valor cero. Uniendo los resultados anteriores, tenemos que el máximo absoluto de f es 16, tomándose en los puntos (2, 2) y (−2, 2), y el mínimo absoluto de f es cero y se alcanza en todos los puntos de segmento S. ª © Segunda forma: La circunferencia C = (x, y) : x2 + y 2 = 8 puede escribirse, en coordenadas √ √ polares, en la forma x = 2 2 cos θ, y = θ, con θ ∈ [0, 2π). Así, la función f restringida √ 2 2 sen √ a C puede expresarse en la forma f (2 2 cos θ, 2 2 sen θ) = 8(cos θ + sen θ)2 , con θ ∈ [0, 2π), y por tanto, debemos analizar la función α(θ) = 8(cos θ + sen θ)2 ,con θ ∈ [0, 2π). (No incluimos el punto θ = 2π, pues α(0) = α(2π)). No analizaremos la función α en los puntos donde ésta se anula (los puntos θ ∈ [0, 2π) tales que cos θ + sen θ = 0, es decir, θ = 3π/4 y θ = 7π/4), ya que sabemos que en los correspondientes puntos de la circunferencia, la función f es nula y toma su valor mínimo. Derivando la función α, obtenemos α0 (θ) = 16(cos θ + sen θ)(− sen θ + cos θ), y los valores θ ∈ [0, 2π) que no anulan a α, pero sí a su derivada, son aquellos que verifican − sen θ + cos θ = 0. Los únicos valores en el intervalo [0, 2π) que resuleven esta ecuación son θ = π/4 y θ = 5π/4. Así, puesto que α(0) = α(2π) = 8 y α(π/4) = α(5π/4) = 16, el valor máximo de α es 16 y se alcanza en θ = π/4 y θ = 5π/4. Por consiguiente, el valor máximo de√f restrigida a la circunferencia C es 16 y dicho máximo √ √ máximo se√alcanza en los valores x = 2 2 cos π/4 = 2, y = 2 2 sen π/4 = 2 y x = 2 2 cos 5π/4 = −2, y = 2 2 sen 5π/4 = −2. Uniendo los resultados anteriores, tenemos que el máximo absoluto de f en R es 16, tomándose en los puntos (2, 2) y (−2, 2), y el mínimo absoluto de f en R es cero y se alcanza en todos los puntos de segmento S. 18. Calcular la mínima distancia desde el punto P (5, 5, 0) al paraboloide de ecuación z = x2 + y 2 . Solución: La distancia de un punto cualquiera (x, y, z) al punto (5, 5, 0) viene dada por p d(x, y, z) = (x − 5)2 + (y − 5)2 + z 2 .

Ahora bien, si el punto (x, y, z) pertenece al paraboloide de ecuación z = x2 + y 2 , entonces la función distancia se reduce a la función de dos variables q 2 g(x, y) = (x − 5)2 + (y − 5)2 + (x2 + y 2 ) . Pero obsérvese que optimizar esta función equivale a la optimización de la función que resulta de suprimir la raíz cuadrada. Así, calcularemos los extremos de la función ¢2 ¡ f (x, y) = (x − 5)2 + (y − 5)2 + x2 + y 2

Comenzamos calculando los puntos críticos resolviendo el sistema de ecuaciones  ∂f  (x, y) = 2x − 10 + 4x3 + 4xy 2 = 0,  ∂x   ∂f (x, y) = 2y − 10 + 4yx2 + 4y 3 = 0. ∂y

¢ ¡ Restando ambas ecuaciones se obtiene 0 = 2x−2y+4x3 −4y 3 +4xy 2 −4yx2 = 2 (x − y) 2x2 + 1 + 2y 2 , de donde x = y, y sustituyendo en cualquiera de las dos ecuaciones del sistema debemos resolver la ecuación 2x − 10 + 4x3 + 4x3 = 0. Es decir, debemos encontrar las raíces del polinomio p(x) = 2x − 10 + 8x3 = 0. La única solución real de p es x = 1, por lo que el único punto crítico de f es x = 1, y = x = 1, esto es, el punto (1, 1). Ahora veremos que este punto crítico es mínimo relativo para la función f. Calculamos la derivadas parciales de segundo orden de f. ∂2f ∂2f (x, y) = 2 + 12x2 + 4y 2 =⇒ (1, 1) = 18 2 ∂x ∂x2

Boletín 6. Ejercicios Resueltos. Fundamentos Matemáticos de la Ingeniería Mecánica. Curso 2003-04

12

∂2f ∂2f (x, y) = 2 + 12y 2 + 4x2 =⇒ (1, 1) = 18 2 ∂y ∂y 2 ∂2f ∂2f (x, y) = 8yx =⇒ (1, 1) = 8 ∂x∂y ∂x∂y

µ 2 ¶2 ∂2f ∂2f ∂ f ∂2f (1, 1) = 18 > 0 y (1, 1) · (1, 1) − = 18 · 18 − 82 = 260 > 0 y por (1, 1) ∂x2 ∂x2 ∂y 2 ∂x∂y tanto, el punto (1, 1) es un mínimo relativo para la función f. Obsérvese que el correspondiente punto del paraboloide es el punto (1, 1, 2) (z = x2 + y 2 = 12 + 12 = 2) Luego

Representando el paraboloide y el punto (5, 5, 0) podemos asegurar que el problema posee una única solución y por tanto, la mínima distancia del punto (5, 5, 0) al paraboloide de ecuación z = x2 + y 2 podemos obtenerla calculando la distancia del punto (5, 5, 0) al punto (1, 1, 2) de dicho paraboloide, p d = (1 − 5)2 + (1 − 5)2 + 22 = 6.

19. Un tanque metálico rectangular sin tapa debe contener 256 metros cúbicos de líquido. ¿Cuáles son las dimensiones del tanque que requiere menos material para su construcción? Solución: Para la construcción del tanque debemos construir la base y sus cuatro caras. Si x e y son las longitudes de los lados de la base y z es la altura, entonces debemos minimizar la función g(x, y, z) = xy + 2xz + 2yz Puesto que el volumen del tanque es 256 m3 , se tiene que xyz = 256, luego z =

256 y la función a xy

minimizar es f (x, y) = xy + 2x

256 512 512 256 + 2y = xy + + xy xy y x

con

x, y > 0

En primer lugar, obtenemos los puntos críticos de la función f resolviendo el sistema  ∂f yx2 − 512   = 0, (x, y) =  ∂x x2 ∂f xy 2 − 512    = 0, (x, y) = ∂y y2

cuya única solución real es x = 8, y = 8. Para clasificar el punto crítico (8, 8) de f, calculamos sus derivadas parciales de segunco orden: 1024 ∂2f 1024 ∂2f (x, y) = 3 =⇒ (8, 8) = 3 = 2, 2 ∂x x ∂x2 8 1024 ∂2f 1024 ∂2f (x, y) = 3 =⇒ (8, 8) = 3 = 2, ∂y 2 y ∂y 2 8 ∂2f ∂2f (x, y) = 1 =⇒ (8, 8) = 1. ∂x∂y ∂x∂y µ 2 ¶2 ∂2f ∂2f ∂ f ∂2f (8, 8) = 2 > 0 y (8, 8) · 2 (8, 8)− (8, 8) = 2 ·2− 12 = 3 > 0 y el punto Por consiguiente, ∂x2 ∂x2 ∂y ∂x∂y crítico (8, 8) es un mínimo relativo para f. Así, las dimensiones de la caja que requieren menor material 256 son x = 8, y = 8, z = = 4. 8·8

20. Una caja rectangular, cuyos lados son paralelos a los ejes de coordenadas, se inscribe en el elipsoide de ecuación 96x2 + 4y 2 + 4z 2 = 36. ¿Cuál es el mayor volumen posible para tal caja? Solución:

Boletín 6. Ejercicios Resueltos. Fundamentos Matemáticos de la Ingeniería Mecánica. Curso 2003-04

13

El volumen de la caja puede expresearse en la forma V = 8xyz, donde x, y, z > 0 ypel punto (x, y, z) está en la superficie del elipsoide de ecuación 96x2 + 4y 2 + 4z 2 = 36. Por tanto, z = −24x2 − y 2 + 9 y la función a maximizar es p V (x, y) = 8xy −24x2 − y 2 + 9 con − 24x2 − y 2 + 9 ≥ 0, x > 0, y > 0. Luego, vamos a resolver© un problema de máximo absoluto paraªuna función continua V en un conjunto cerrado y acotado R = (x, y) ∈ R2 : 24x2 + y 2 ≤ 9, x ≥ 0, y ≥ 0 ..

Comenzamos calculando los puntos críticos de V en el interior del conjunto R. La única solución del sistema  ∂V 48x2 + y 2 − 9   p = 0, (x, y) = −8y √   ∂x √ −24x2 − y 2 + 9 2 , y = 3. Puesto que V = 0 en las fronteras de R, es x = 2 2  4 ∂V 24x + 2y − 9   = 0, (x, y) = −8x p  2 2 ∂y −24x − y + 9 ! Ã√ √ 2 x = 0, y = 0 y 24x2 + y 2 − 9 = 0, el máximo absoluto de V se alcanza en el punto , 3 y su máximo 4 v ! Ã√ à √ !2 u √ u √ √ √ ¡√ ¢2 2 2 2 t 3 + 9 = 6 2. Además, las dimensiones de la , 3 =8· · 3 · −24 valor es V − 4 4 4 v à √ !2 u √ √ u p √ √ ¡√ ¢2 2 2 2 = , 2y = 2 3, 2z = 2 −24x2 − y 2 + 9 = 2t−24 caja son: 2x = 2 − 3 + 9 = 2 3. 4 2 4

21. Una caja de R3 tiene la forma del cubo 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, 0 ≤ z ≤ 1. Se coloca dentro de la caja la placa que es la porción del plano x + y + z = 1 que cabe dentro de ella. Se calienta la placa de tal forma que la temperatura en el punto (x, y, z) es T (x, y, z) = 4 − 2x2 − y 2 − z 2 . Hallar el punto más caliente y más frío de la placa. Solución: El problema en cuestión se trata de un problema de extremos absolutos en región cerrada y acotada. Puesto que x + y + z = 1, tenemos z = 1 − x − y y la función a optimizar es T (x, y) = 4 − 2x2 − y 2 − (1 − x − y)2 © ª restringida a la región R = (x, y) ∈ R2 : x ≥ 0, y ≥ 0, x + y ≤ 1 . Puesto que T es continua en R y R es cerrada y acotada, T alcanza su máximo absoluto y su mínimo absoluto en puntos de R. Ahora calculamos los puntos críticos de T en el interior de R. Para ello, resolvemos el sistema  ∂T   (x, y) = −6x + 2 − 2y = 0, ∂x ∂T   (x, y) = −4y + 2 − 2x = 0. ∂y

1 2 , y = . Y este punto se encuentra en el interior de R, pues 5 5 ¶ µ ¶2 µ ¶2 µ 1 2 1 2 3 1 2 1 2 = 4−2 x = > 0, y = > 0, x + y = + = < 1. Obsérvese que T , − − 5 5 5 5 5 5 5 5 5 µ ¶2 1 2 18 1− − = . 5 5 5 ª © Ahora analizamos la restricción de T al conjunto (x, y) ∈ R2 : 0 ≤ x ≤ 1, y = 0 . Es decir, estudiamos la función f (x) = T (x, 0) = 4 − 2x2 − (1 − x)2 = 3 − 3x2 + 2x, con x ∈ [0, 1]. En primer lugar, calculamos los puntos críticos de f en el intervalo abierto (0, 1), es decir, resolvemos la ecuación f 0 (x) = −6x + 2 = 0 en el intervalo (0, 1). La única solución de esta ecuación es x = 1/3 ∈ (0, 1). Evaluamos f en el punto crítico y en los extremos del intervalo [0, 1]. Como f (0) = 3, f (1/3) = 3 − 3 · (1/3)2 + 2 · 1/3 = 10 3 y f (1) = 2, el máximo absoluto de f es 10/3 y se alcanza en x = 1/3 y el mínimo absoluto de f es 2, tomándose en x = 1.

La única solución de este sistema es x =

14

Boletín 6. Ejercicios Resueltos. Fundamentos Matemáticos de la Ingeniería Mecánica. Curso 2003-04

© ª Estudiamos a continuación la restricción de T al conjunto (x, y) ∈ R2 : 0 ≤ y ≤ 1, x = 0 . Tenemos que analizar la función g(y) = T (0, y) = 4 − y 2 − (1 − y)2 = 3 − 2y 2 + 2y, para 0 ≤ y ≤ 1. La ecuación g 0 (y) = −4y + 2 = 0 posee una única solución y = 1/2 ∈ (0, 1). Evaluando g en los extremos del intervalo 7 2 (0, 1) y en el punto crítico y = 1/2, se obtiene g(0) = 3, g(1/2) = 3 − 2 · (1/2) + 2 · 1/2 = , g(1) = 3. 2 Por consiguiente, el máximo absoluto de g es 7/2 y se alcanza en y = 1/2 y el mínimo absoluto se toma en y = 0 e y = 1 y su valor es 3. Por último, estudiamos el comportamiento de la función T sobre el conjunto © ª (x, y) ∈ R2 : 0 ≤ x ≤ 1, x + y = 1 .

2

Tomando y = 1−x, debemos analizar la función h(x) = T (x, 1−x) = 4−2x2 −(1−x)2 −(1 − x − (1 − x)) = 3 − 3x2 + 2x, con x ∈ [0, 1]. Obsérvese que las funciones f y h coinciden, luego el máximo absoluto de h es 10/3 y se alcanza en x = 1/3 y el mínimo absoluto de h es 2, tomándose en x = 1. Teniendo en consideración los cálculos anteriores, podemos establecer que la función T toma su máximo ¶ µ 18 1 2 y su valor es absoluto en el punto , y su mínimo absoluto es 2, tomándose en los puntos (1/3, 0) 5 5 5 y (1/3, 2/3). ¶ µ 1 2 2 1 2 2 (pues z = 1 − − = ) y los puntos Luego, el punto más caliente de la placa es el punto , , 5 5 5 5 5 5 (1/3, 0, 2/3) y (1/3, 2/3, 0) son los puntos de menor temperatura.