10  Derivadas de funciones de varias variables

Ejercicio 10.1 Emparejar las siguientes funciones de dos variables con sus curvas de nivel.

  1. f(x,y)=cos(x)sen(y)

  2. g(x,y)=xy

  3. h(x,y)=2x+y

  4. i(x,y)=3x2+y2

La gráfica A corresponde a i(x,y).
La gráfica B corresponde a g(x,y).
La gráfica C corresponde a f(x,y).
La gráfica D corresponde a h(x,y).

Ejercicio 10.2 La ecuación de los gases perfectos relaciona la presión, el volumen y la temperatura de un gas perfecto. Esta ecuación se suele escribir de la forma

PV=nRT

donde P es la presión, V es el volumne, T la temperatura, R es la constante universal de los gases perfectos y n es el número de moles del gas, que también es constante para cada gas.

  1. ¿Cómo varía la presión de un gas perfecto cuando se aumenta la temperatura, manteniendo constante el volumen?

  2. ¿Cómo varía la presión de un gas perfecto cuando se aumenta el volumen, manteniendo constante la temperatura?

Podemos expresar la presión en función del volumen y la temperatura del gas perfecto mediante la función P(V,T)=nRTV.

  1. La variación de la presión con respecto a la temperatura, manteniendo constante el volumen es la derivada parcial de la presión con respecto a la temperatura, que vale

    PT=nRV.

    Como el volumen siempre es positivo, esto quiere decir que al aumentar la temperatura, manteniendo constante el volumen, la presión aumenta.

  2. La variación de la presión con respecto al volumen, manteniendo constante la temperatura es la derivada parcial de la presión con respecto al volumen, que vale

    PV=nRTV2.

    Como el volumen y la temperatura son siempre positivos, esto quiere decir que al aumentar el volumen, manteniendo constante la temperatura, la presión disminuye.

Ejercicio 10.3 La asimilación de CO2 de una planta depende de la temperatura ambiente (t) y de la intensidad de la luz (l), según la función

f(t,l)=ctl2,

donde c es una constante.

Estudiar cómo evoluciona la asimilación de CO2 para distintas intensidades de luz, cuando se mantiene la temperatura constante. Estudiar también cómo evoluciona para distintas temperaturas cuando se mantiene la intensidad de la luz constante.

La variación de la asimilación de CO2 con respecto a la intensidad de la luz viene dada por la derivada parcial con respeto a l, que vale

fl(t,l)=2ctl,

y la variación de la asimilación de CO2 con respecto a la temperatura viene dada por la derivada parcial con respecto a t, que vale

ft(t,l)=cl2.

Ejercicio 10.4 La función de producción de Cobb-Douglas

p(t,k)=atαk1α

se utiliza en Econometría para modelizar, de manera simplificada, la producción económica de un país, es decir, el valor monetario de los bienes que se producen en un año, en función de la cantidad de trabajo t y el capital invertido k, y donde a es una constante conocida como factor total de productividad y α es otra constante que mide la elasticidad del producto del trabajo y el capital.

  1. Tomando a=1 y α=0.5, ¿cómo son las curvas de nivel de esta función de producción?

  2. ¿Cómo varía la producción cuando aumenta la cantidad de trabajo y se mantiene constante el capital?

  3. ¿Cómo varía la producción cuando disminuye el capital?

  1. Para a=1 y α=0.5 se tiene la función p(t,k)=1t0.5k0.5=tk. Los puntos de la curva de nivel c cumplirán que tk=c, de donde se deduce que t=c2k y, por tanto, se trata de una función inversa. En la siguiente gráfica aparecen algunas curvas de nivel de esta función de producción.

  2. La variación de la producción con respecto a al trabajo, manteniendo constante el capital, lo da la derivada parcial con respecto al trabajo, que vale

    pt=aαtα1k1α.

    En el caso particular de a=1 y α=0.5 se tiene

    pt=12t1/2k1/2=12kt.

    Esto quiere decir que al aumentar el trabajo la producción económica aumenta a razón de esta cantidad. Se observa que la cantidad de trabajo aparece en el denominador, por lo que cuando la cantidad del trabajo es grande, el incremento de producción es pequeño.

  3. La variación de la producción con respecto a al capital, manteniendo constante el trabajo, lo da la derivada parcial con respecto al capital, que vale

    pk=atα(1α)k1α1=a(1α)tαkα.

    En el caso particular de a=1 y α=0.5 se tiene

    pt=12t1/2k1/2=12tk.

    Esto quiere decir que al aumentar el capital la producción económica aumenta a razón de esta cantidad. Al igual que antes, se observa que el capital aparece en el denominador, por lo que cuando el capital es grande, el incremento de producción es pequeño.

Ejercicio 10.5 Se tiene un cilindro de radio 5 cm y altura 3 cm. ¿Es el volumen del cilindro más sensible a una pequeña variación de su radio o de su altura? ¿Cuál debería ser la altura del cilindro para que su volumen fuese igual se sensible a una pequeña variación de su radio que a una pequeña variación de su altura?

El volumen de un cilindro depende de su radio y de su altura según la función v(r,h)=πr2h.

La tasa de variación del volumen con respecto al radio es la derivada parcial del volumen con respecto al radio, es decir,

vr=2πrh,

que para un cilindro de radio 5 cm y altura 3 cm es vr=2π53=30π.

Del mismo modo, la tasa de variación del volumen con respecto a la altura es la derivada parcial del volumen con respecto a la altura, es decir,

vh=πr2,

que para un cilindro de radio 5 cm y altura 3 cm es vh=π52=25π.

Así pues, el volumen de este cilindro es más sensible a una pequeña variación del radio que de la altura.

Para que fuese igualmente sensible a la variación del radio y de la altura, ambas derivadas parciales deberían ser iguales, es decir,

vr=vh2πrh=πr2h=r/2.

Ejercicio 10.6 Una empresa fabrica helados de tres sabores. El coste total de producción viene dada por la función

c(x,y,z)=110xyz3+x20+y40+z30+100,

donde x es el número de helados de chocolate, y el número de helados de fresa y z el número de helados de vainilla. Calcular el coste marginal para cada uno de los tipos de helados.

El coste marginal es la derivada parcial del coste total con respecto al número de helados producidos de cada tipo. Así, el coste marginal de los helados de chocolate es

cx=130yz(xyz)2/3+120.

El coste marginal de los helados de fresa es

cy=130xz(xyz)2/3+140.

Y el coste marginal de los helados de vainilla es

cz=130xy(xyz)2/3+130.

Ejercicio 10.7 La gráfica de una función f(x,y) se muestra a continuación.

¿Qué signo tienen las derivadas parciales en el punto (1,1). ¿Y en el punto (0,0)?

En el punto (1,1) derivada parcial con respecto a x es negativa, y la derivada parcial con respecto a y es positiva.

En el punto (0,0) ambas derivadas parciales se anulan.

Ejercicio 10.8 Obtener la ecuación del plano tangente y de la recta normal las siguientes superficies en los puntos indicados.

  1. xyz=8 en el punto (4,2,1).

  2. f(x,y)=yln(xy) en el punto (1/2,2,0).

  1. Tomando la función f(x,y,z)=xyz8, se trata de calcular la ecuación de la recta normal y el plano tangente a la superficie de nivel f(x,y,z)=0 en el punto P=(4,2,1). Por este teorema sabemos que el vector gradiente f(4,2,1) es perpendicular a la superficie de nivel en el punto (4,2,1) por lo que basta tomar este vector como vector director de la recta normal.

    f(x,y,z)=(fx,fy,fz)=(yz,xz,xy),

    y en el punto (4,2,1) vale

    f(4,2,1)=((2)(1),4(1),4(2))=(2,4,8).

    Así pues, la ecuación vectorial de la recta normal es

    P+tf(P)=(4,2,1)+t(2,4,8)=(4+2t,24t,18t).

    Y la ecuación del plano tangente

    ((x,y,z)P)f(P)=((x,y,z)(4,2,1))(2,4,8)=2x4y8z+24=0.

  2. La ecuación de la recta tangente a la superficie de f(x,y) en el punto (a,b,f(a,b)) es

    (a+tfx(a,b),b+tfy(a,b),f(a,b)t).

    Las derivadas parciales de f son

    fx=y1xyy=yx,fy=ln(xy)+y1xyx=ln(xy)+1,

    y en el punto (1/2,2) valen

    fx(1/2,2)=21/2=4,fy(1/2,2)=ln(122)+1=1.

    Así pues, sustituyendo en la ecuación de la recta tangente a f(x,y) en (1/2,2,0) se tiene

    (12+4t,2+t,t).

    Y la ecuación del plano normal es

    z=fx(1/2,2)(x12)+fy(1/2,2)(y2)+f(1/2,2)=4(x12)+(y2)+0=4x+y4=0.

Ejercicio 10.9 Calcular el gradiente de la función

f(x,y,z)=logxyz+arcsen(xz).

f(x,y,z)=(z1x2z2+12x,1y,x1x2z21z)

Ejercicio 10.10 Una nave espacial está en problemas cerca del sol. Se encuentra en la posición (1,1,1) y la temperatura de la nave cuando está en la posición (x,y,z) viene dada por T(x,y,z)=ex22y23z2 donde x,y,z se miden en metros. ¿En qué dirección debe moverse la nave para que la temperatura decrezca lo más rápidamente posible?

Para que la temperatura decrezca lo más rápidamente posible la nave debe moverse en la dirección opuesta al vector gradiente.

Como el gradiente de f vale

f(x,y,z)=(2xex22y23z2,4yex22y23z2,6zex22y23z2),

en el punto (1,1,1) se tiene f(1,1,1)=e6(2,4,6).

Ejercicio 10.11 Supongamos que la cantidad de agua almacenada en un pantano al final del año hidrológico, A en hectómetros cúbicos, viene dada por la función

A(p,t,c)=p3t1c2ecpt

donde p es la precipitación en litros/m2 caí­da durante el año hidrológico, t es la temperatura media del año hidrológico en ºC y c el consumo debido a abastecimiento de poblaciones cercanas y riego, en hectómetros cúbicos.

  1. Calcular el gradiente de la cantidad de agua almacenada.

  2. Suponiendo que hubiese algún año en el que el consumo fuese nulo, ¿qué condición tendría­ que cumplir la temperatura para que la derivada del agua almacenada con respecto a la temperatura fuese igual a la derivada con respecto a la precipitación?

  1. Las derivadas parciales de la función A son

Ap(p,t,c)=12(p3t1c2ecpt)1/23p2t1c2ecptctAt(p,t,c)=12(p3t1c2ecpt)1/2p3(t1)2c2ecptcpAc(p,t,c)=12(p3t1c2ecpt)1/22cecptc2ecptpt

Así pues, el gradiente de A es

A(c,p,t)=12(p3t1c2ecpt)1/2(3p2t1c2ecptct,p3(t1)2c2ecptcp,ecpt(2c+c2pt))

  1. Suponiendo c=0 las derivadas parciales con respecto a la precipitación y la temperatura valen

Ap(p,t,0)=12(p3t1)1/23p2t1At(p,t,0)=12(p3t1)1/2p3(t1)2

Así pues,

Ap(p,t,0)=At(p,t,0)3p2t1=p3(t1)23=pt1t=p3+1.

Ejercicio 10.12 La variable aleatoria bidimensional (X,Y) con función de densidad

f(x,y)=12πσxσye12((xμx)2σx2+(yμy)2σy2)

se conoce como normal bidimensional con X e Y independientes, de parámetros μ=(μx,μy) y σ=(σx,σy).

Calcular el gradiente de f e interpretarlo. ¿En qué punto se anula el gradiente? ¿Qué conclusiones sacas? ¿Cuál es la tasa de crecimiento de f cuando x?

f(x,y)=12πσxσye12((xμx)2σx2+(yμy)2σy2)(xμxσx2,yμyσy2).

El gradiente f(a,b) es un vector que indica la dirección en la que hay que moverse desde un punto (a,b) para que el valor de la función crezca lo más rápidamente posible. Como la función exponencial es siempre positiva y las desviaciones típicas también, la dirección del vector gradiente será la misma que la del vector (μxxσx2,μyyσy2), lo que indica que habría que moverse hacia el punto de medias.

Por otro lado, este vector se anula cuando (μxx,μyy)=(0,0), es decir, en el punto de medias (x=μx,y=μy).

Finalmente, si fijamos y como constante, como la exponencial tiene exponente negativo, se tiene

limx|f(x,y)|=((xμx)2σx4+(yμy)2σy4)2πσxσye12((xμx)2σx2+(yμy)2σy2).

Como en el numerador tenemos un polinomio en x y en el denominador una función exponencial, cuando x la función exponencial crece mucho más rápidamente que cualquier polinomio, y por tanto, limx|f(x,y)|=0.

Ejercicio 10.13 Tenemos dos objetos de masas m1 y m2 unidas por una cuerda que pasa a través de una polea como la de la siguiente figura.

Polea

Si m1m2, la aceleración del objeto de masa m1 viene dada por la ecuación

a=m1m2m1+m2g,

siendo g la aceleración de la gravedad.

Demostrar que se cumple la ecuación

m1am1+m2am2=0.

¿Qué ecuación cumplirán las fuerzas que se aplican sobre estos objetos?

Calculamos las derivadas parciales de la aceleración con respecto a las masas

am1=2gm2(m1+m2)2,am2=2gm1(m1+m2)2.

Sustituyendo ahora en la ecuación es inmediato ver que se cumple,

m1am1+m2am2=2gm1m2(m1+m2)2+2gm1m2(m1+m2)2=0.

Por el funcionamiento de la polea, si a es la aceleración que experimenta el primer objeto, la aceleración del segundo objeto será de la misma magnitud pero opuesta, es decir a. Por tanto, las fuerzas que deben actuar sobre los objetos para producir estas aceleraciones son f1=m1a y f2=m2a, respectivamente. Si calculamos la derivadas parciales de estas fuerzas con respecto a m1 y m2 respectivamente, se tiene

f1m1=a+m1am1,f2m2=am2am2.

Por tanto la ecuación que relaciona las fuerzas es

f1m1f2m1=a+m1am1+a+m2am2=2a.

Ejercicio 10.14 Utilizar la regla de la cadena para funciones de varias variables para demostrar las reglas de derivación de funciones de una variable real

(u+v)=u+v(uv)=uv(uv)=uv+uv(u/v)=uvuvv2

donde u y v son funciones de una variable real.

Sean u(t) y v(t) dos funciones de una variable. Si tomamos la función de dos variables f(x+y)=x+y y la función vectorial g(t)=(u(t),v(t)), entonces, aplicando la regla de la cadena se tiene

(u+v)=f(g(t))=f(g(t))g(t)=(1,1)(u,v)=u+v.

Aplicando el mismo procedimiento, si se toma f(x,y)=xy se tiene

(uv)=f(g(t))=f(g(t))g(t)=(1,1)(u,v)=uv.

Del igual modo, si se toma f(x,y)=xy se tiene

(uv)=f(g(t))=f(g(t))g(t)=(v,u)(u,v)=uv+uv.

Y finalmente, si se toma f(x,y)=xy se tiene

(uv)=f(g(t))=f(g(t))g(t)=(1v,uv2)(u,v)=uvuvv2.

Ejercicio 10.15 En un equilibrio químico, la concentración de una sustancia z depende de las concentraciones de otras dos sustancias x e y según la ecuación zxyln(x+yy)=0. Usar el diferencial para predecir la variación que experimentará la concentración de z si la concentración de x aumenta 0.1 mg/mm3 y la de y disminuye 0.2 mg/mm3 en el instante en que las concentraciones de x e y son 1 mg/mm3.

Despejando z de la ecuación se tiene la función de dos variables

z=f(x,y)=ln(x+yy)(xy)1/2=(ln(x+y)ln(y))(xy)1/2.

Las derivadas parciales de f son

fx=1x+y(xy)1/2+ln(x+yy)12(xy)1/2y,fy=(1x+y1y)(xy)1/2+ln(x+yy)12(xy)1/2x.

Estas funciones son continuas en un entorno del punto (1,1), por lo que la función f es diferenciable en (1,1) y podemos utilizar el diferencial total para aproximar la variación de z,

dz=f(1,1)(dx,dy)=(ln(2)+12,ln(2)12)(dx,dy)=(ln(2)+1)dx2+(ln(2)1)dy2.

Así pues, para variaciones de las concentraciones Δx=0.1 y Δy=0.2 mg/mm3, respectivamente, se tiene que la variación de z es

Δz(ln(2)+1)0.12+(ln(2)1)(0.2)20.1153 mg/mm3.

Ejercicio 10.16 Un cilindro metálico de radio 30cm y altura 50cm se dilata de manera que la tasa de variación de su radio es 1 cm/s y la tasa de variación de su altura es 2 cm/s. ¿Cuál es la tasa de variación de su volumen?

La fórmula del volumen de un cilindro de radio r y altura h es la función de dos variables V(r,h)=πr2h, pero como tanto el radio como la altura dependen del tiempo debido al proceso de dilatación, en realidad se tiene que el volumen depende del tiempo según la fórmula V(t)=πr(t)2h(t). Así pues, podemos calcular la derivada del volumen aplicando la regla de la cadena.

dVdt=Vrdrdt+Vhdhdt=2πrhdrdt+πr2dhdt.

Y en el instante en que r=30 y h=50 y empezamos a cambiar el radio con una tasa de variación drdt=1 cm/s y la altura con una tasa de variación dhdt=2 cm/s, la tasa de variación del volumen es

dVdt=2π30501+π3022=4800π cm3/s.

Ejercicio 10.17 Dos planetas se mueven en un mismo plano describiendo órbitas dadas por las funciones vectoriales f(t)=(3cos(t),sen(t)) y g(t)=(2sen(t/2),3cos(t/2)). Usar la regla de la cadena para funciones de varias variables para calcular la derivada de la distancia entre los planetas en el instante t=π. ¿Los planetas se están acercando o alejando?

La distancia euclídea entre dos puntos (x1,y1) y (x2,y2) del plano euclídeo es D(x,y)=(x2x1)2+(y2y1)2. Suponiendo que estos puntos son las posiciones de los planteas, que dependen del tiempo según las funciones vectoriales f(t) y g(t), aplicando la regla de la cadena, la derivada de la distancia entre los planetas es

dDdt=Dx1dx1dt+Dx2dx2dt+Dy1dy1dt+Dy2dy2dt=12((x2x1)2+(y2y1)2)1/2(2)(x2x1)dx1dt+12((x2x1)2+(y2y1)2)1/22(x2x1)dx2dt+12((x2x1)2+(y2y1)2)1/2(2)(y2y1)dy1dt+12((x2x1)2+(y2y1)2)1/22(y2y1)dy2dt=(x2+x1)dx1dt+(x2x1)dx2dt+(y2+y1)dy1dt+(y2y1)dy2dt(x2x1)2+(y2y1)2.

Como en el instante t=π tenemos

f(π)=(3cos(π),sen(π))=(3,0),g(π)=(2sen(π/2),3cos(π/2))=(2,0),f(π)=(3sen(π),cos(π))=(0,1),g(π)=(22cos(π/2),32sen(π/2))=(0,3/2),

sustituyendo en la expresión anterior se tiene

dDdt=(23)0+(2+3)0+(0+0)1+(00)(3/2)(2+3)2+(00)2=0.

Como la derivada es nula, en el instante t=π los planetas ni se están acercando ni alejando. De hecho, en ese instante la distancia entre los planetas es máxima, tal y como puede apreciarse en la siguiente gráfica.

Orbitas de los planetas.

Ejercicio 10.18 Para calcular el caudal de una fuente se ha llenado un depósito de base cuadrada 30±0.2 cm y altura 50±0.3 cm en un tiempo de 20±0.5 s. Expresar el caudal de la fuente incluyendo error en la medida.

El caudal de la fuente es el volumen dividido por el tiempo, pero el volumen del depósito depende de sus dimensiones, la función que da el caudal de la fuente es c(x,y,t)=x2yt, donde x es el lado de la base, y es la altura y t es el tiempo, y su gradiente vale

c(x,y,t)=(2xyt,x2t,x2yt2).

Para calcular el error en la medición del caudal de manera aproximada usando el diferencial del caudal,

dc=c(30,50,20)(dx,dy,dt)=(150,45,112.5)(dx,dy,dz)=150dx+45dy112.5dz,

que para unos errores dx=±0.2, dy=±0.3 y dt=±0.5 se tiene

dc=150(±0.2)+45(±0.3)112.5(±0.5)=±99.75.

Así pues, la expresión del caudal, incluyendo el error en la medida, es

c(30,50,20)±dc=2250±99.75 cm3/s.

Ejercicio 10.19 La potencia eléctrica P en un circuito depende de el voltaje E y la resistencia R según la fórmula P=E2/R. Si los errores cometidos al medir el voltaje y la resistencia son de un 1% y un 2% respectivamente, ¿cuál será el error relativo en el cálculo de la potencia eléctrica?

Usando el diferencial como una aproximación del error en la medición de las magnitudes, se tiene que el error relativo en la medición del potencial es

dPP=P(dE,dR)E2/R=(2ER,E2R2)(dE,dR)E2/R=R(2ERdEE2R2dR)E2=2dEEdRR.

Y como dEE y dRR son los errores relativos en las mediciones del voltaje y la resistencia, respectivamente, se tiene que el error relativo en la medición del potencial es

dPP=2(±0.01)(±0.02)=±0.04,

es decir, un 4%.

Ejercicio 10.20 La Quimiotaxis es el movimiento de los organismos dirigido por un gradiente de concentración, es decir, en la dirección en la que la concentración aumenta con más rapidez. El moho del cieno Dictyoselium discoideum muestra este comportamiento. En esta caso, las amebas unicelulares de esta especie se mueven según el gradiente de concentración de una sustancia química denominada adenosina monofosfato (AMP cíclico). Si suponemos que la expresión que da la concentración de AMP cíclico en un punto de coordenadas (x,y,z) es:

C(x,y,z)=4x2+y2+z4+1

y se sitúa una ameba de moho del cieno en el punto (1,0,1), ¿en qué dirección se moverá la ameba?

La dirección en la que más rápidamente aumenta la concentración de AMP es la dirección del gradiente de C, que vale

C(x,y,z)=(412(x2+y2+z4+1)3/22x,412(x2+y2+z4+1)3/22y,412(x2+y2+z4+1)3/24z3)=4(x2+y2+z4+1)3/2(x,y,2z3),

y en el punto (1,0,1) vale (4/27,0,8/27).

Ejercicio 10.21 Se sabe que un gas perfecto a una presión de 1 atmósfera y una temperatura de 273.1 K ocupa un volumen de 22.4 litros.

  1. ¿En qué dirección deben cambiarse la presión y la temperatura para conseguir el mayor aumento del volumen?

  2. ¿Cómo varía el volumen si comenzamos a incrementar la temperatura a razón del doble del incremento de la presión?

  3. ¿Cuánto tendría que aumentar la temperatura por cada atmósfera de presión que se aumente para que el volumen permanezca constante?

La ecuación de los gases perfectos es PV=cT, donde c es una constante que depende del gas. En el caso particular del gas del ejercicio se tiene

c=PVT=12.4273.10.082 atm l/K.

Si de la ecuación de los gases perfectos despejamos el volumen, se tiene la función de dos variables V(P,T)=cTP.

  1. La dirección en la que hay que cambiar la temperatura y la presión para que el volumen aumente lo más rápidamente posible, es la dirección el vector gradiente,

    V(P,T)=(cTP2,cP),

    que para la temperatura y la presión dadas vale

    V(1,273.1)=(0.082273.112,0.0821)=(22.4,0.082).

  2. Para incrementar la temperatura a razón del doble del incremento de la presión, debemos cambiar la presión y la temperatura en la dirección del vector u=(1,2), y la variación que experimenta el volumen en esta dirección viene dado por la derivada direccional del volumen en la dirección de u, que vale

    Vu(1,273.1)=V(1,273.1)u|u|=(22.4,0.082)(1,2)5=22.41+0.082259.9442 l.

  3. Tomemos un vector unitario cualquiera u=(ux,uy). Como u es unitario, se cumple que ux2+uy2=1, por lo que, uy=1x2, y el vector tiene componentes u=(ux,1ux2). Para que el volumen permanezca constante, al cambiar la presión y la temperatura en la dirección de u, la derivada direccional del volumen en esa dirección debe ser nula, es decir,

    Vu(1,273.1)=V(1,273.1)u=(22.4,0.082)(ux,1ux2)=22.4ux+0.0821ux2=022.4ux=0.0821ux2273.1ux=1ux2273.12ux2=1ux2273.12ux2+ux2=1ux2=1273.12+1ux=±1273.12+1±0.0037.

    Por tanto, habrá que cambiar la presión y la temperatura en la dirección del vector u=(0.0037,10.00372)=(0.0037,0.9999).

    Podríamos haber resuelto el problema de manera más sencillas usando la propiedad de que el vector gradiente es ortogonal a la curva de nivel, y por tanto, para obtener la dirección en la que el volumen permanece constante bastaría con tomar un vector ortogonal al vector gradiente, como por ejemplo u=(0.082,22.4).

Ejercicio 10.22 Se han diseñado unas cápsulas con forma piramidal con base un rectángulo de lados a=3 cm, b=4 cm, y altura h=6 cm.

  1. ¿Cómo deben cambiar las dimensiones de la cápsula para que el volumen aumente lo más rápidamente posible? ¿Cuál sería la tasa de variación del volumen si cambian las dimensiones de la cápsula en la proporciones anteriores?

  2. Si se empiezan a cambiar las dimensiones de la cápsula de manera que el lado mayor del rectángulo disminuye la mitad de lo que aumenta el lado menor, y la altura aumenta el doble de lo que aumenta el lado menor, ¿cuál sería la tasa de variación del volumen de la cápsula en las condiciones anteriores?

El volumen de la cápsula viene dado por la función V(a,b,c)=13abc.

  1. La dirección de máximo crecimiento del volumen es la dirección del vector gradiente,

    V(a,b,c)=13(bc,ac,ab).

    que en el punto (3,4,6) vale V(3,4,6)=(8,6,4). La tasa de variación del volumen si se cambian sus dimensiones en la dirección del vector gradiente es

    |V(3,4,6)|=|(8,6,4)|=82+62+42=10.7703 cm3/s.

  2. La derivada direccional de V en (3,4,6) siguiendo la dirección del vector u=(1,1/2,2) es

Vu(3,4,6)=V(3,4,6)u|u|=(8,6,4)(1,1/2,2)12+(1/2)2+22=816(1/2)+425.255.6737 cm3/s.

Ejercicio 10.23 ¿En qué direcciones se anulará la derivada direccional de la función

f(x,y)=x2y2x2+y2

en el punto P=(1,1)?

Tomemos un vector unitario cualquiera u=(ux,uy). Como u es unitario, se cumple que ux2+uy2=1, por lo que, uy=1x2, y el vector tiene componentes u=(ux,1ux2). Para calcular la derivada direccional de f en la dirección del vector u, primero necesitamos el gradiente de f,

f(x,y)=(2x(x2+y2)(x2y2)2x(x2+y2)2,2y(x2+y2)(x2y2)2y(x2+y2)2)=(4xy2(x2+y2)2,4x2y(x2+y2)2),

que en el punto (1,1) vale

f(1,1)=(4112(12+12)2,4121(12+12)2)=(1,1).

Así pues, la derivada direccional de f en la dirección de u es

fu(1,1)=f(1,1)u=(1,1)(ux,1ux2)=0ux1ux2=0ux=1ux2ux2=1ux2ux2=1/2ux=±12.

Por tanto, la derivada direccional será nula en la dirección del vector (1/2,1/2).

Ejercicio 10.24 ¿Existe alguna dirección en la que la derivada direccional en el punto P=(1,2) de la función f(x,y)=x23xy+4y2 valga 14?

El gradiente de f es

f(x,y)=(2x3y,3x+8y),

que en el punto (1,2) vale

f(1,2)=(2132,31+82)=(4,13)

Como el máximo crecimiento de una función se da en la dirección del vector gradiente, se tiene que el máximo crecimiento de f es

|f(1,2)|=|(4,13)|=(4)2+13213.6015.

Por tanto, no existe ninguna dirección en la que la derivada direccional de f sea 14.

Ejercicio 10.25 La derivada direccional de una función f en un punto P es máxima en la dirección del vector (1,1,1) y su valor es 23. ¿Cuánto vale la derivada direccional de f en P en la dirección del vector (1,1,0)?

La derivada direccional de una función es máxima en la dirección del vector gradiente, por lo que el vector (1,1,1) tiene la misma dirección del vector gradiente, es decir, f(P)=k(1,1,1). Como además, su módulo debe ser 23 se tiene

|f(P)|=|k(1,1,1)|=k12+12+(1)2=k3=23,

de donde se deduce que k=2, y por tanto, el vector gradiente es f(P)=(2,2,2).

Así pues, la derivada direccional de f en la dirección del vector u=(1,1,0) es

fu(P)=f(P)u|u|=(2,2,2)(1,1,0)2=2+22=22.

Ejercicio 10.26 Dado el campo escalar

f(x,y,z)=x2y2+xyz3zx,

  1. Calcular la derivada direccional de f en (1,2,3) a lo largo del vector unitario u=12(1,1,0).

  2. ¿En qué dirección es máxima la derivada direccional de f en el punto anterior? Obtener el valor de dicha derivada direccional.

  1. El gradiente de f es

    f(x,y,z)=(2x+yz3z,2y+xz3,3xyz2x),

    que en el punto (1,2,3) vale

    f(1,2,3)=(21+2333,22+133,312321)=(53,23,53).

    Por tanto, la derivada direccional de f en la dirección del vector u es

    fu(1,2,3)=f(1,2,3)u=(53,23,53)(1,1,0)2=53232=302=152.

  2. La derivada direccional es máxima en la dirección del gradiente (53,23,53) y vale |f(1,2,3)|=6147.

Ejercicio 10.27 Sea f(x,y,z)=0 la ecuación que define una superficie en el espacio real. Demostrar que el vector gradiente de f en cualquier punto (a,b,c) es normal a la superficie en ese punto.

Consideremos una función vectorial g(t) cuya trayectoria circule sobre la superficie f(x,y,z)=0 y tal que g(0)=(a,b,c). Entonces, como los puntos de g están en la superficie, cumplirán la ecuación f(g(t))=0, de manera que si calculamos su derivada se tiene

f(g(t))=f(g(t))g(t)=0.

En particular en t=0 se tiene

f(g(0))g(0)=f(a,b,c)g(0)=0,

de donde se deduce que los vectores f(a,b,c) y g(0) son normales. Pero como g(0) es el vector tangente a la trayectoria de g en el punto (a,b,c) y esto ocurre para cualquier trayectoria inscrita en la superficie que pase por (a,b,c) en t=0, se concluye que f(a,b,c) es normal a la superficie en (a,b,c).

Ejercicio 10.28 Un cuerpo se mueve en el plano a través de los puntos de coordenadas (x,y) relacionadas mediante la siguiente expresión 2exysen(x)+ycos(x)=2. Calcular las ecuaciones de las rectas tangente y normal a la trayectoria del cuerpo cuando x=0.

Veamos primero cuáles son las coordenadas del punto donde hay que calcular las rectas tangente y normal a la trayectoria. Sustituyendo x=0 en la ecuación anterior se tiene

2e0sen0+ycos0=2y=2,

y, por tanto, se trata del punto (0,2).

Sabemos que la ecuación de la recta tangente a la gráfica de la función f en el punto de coordenadas (a,f(a)) es:

ya=f(a)(xa)

En nuestro caso, no tenemos la expresión explícita de la función f, pero sí que tenemos la ecuación

F(x,y)=2exysen(x)+ycos(x)2=0,

que define a y como función implícita de x. Podemos afirmar que esta función existe en un entorno de x=0 ya que se cumplen las condiciones del teorema de la función implícita, es decir,

  1. F(0,2)=2e0sen(0)+2cos(0)2=0.
  2. Fx=2exyysen(x)+2exycos(x)ysen(x).
    Fy=2exyxsen(x)+cos(x).
    Son continuas en un entorno de (0,2).
  3. Fy(0,2)=2e020sen(0)+cos(0)=10.

Así pues, la derivada de y con respecto a x es

dydx(0,2)=Fx(0,2)Fy(0,2)=2e022sen(0)+2e02cos(0)ysen(0)2e020sen(0)+cos(0)=21=2.

Por lo tanto, la ecuación de la recta tangente es y2=2(x0)y=2x+2,

y la ecuación de la recta normal es

y2=12(x0)y=12x+2.

Podríamos haber llegado al mismo resultado sabiendo que F(x,y) es normal a la trayectoria F(x,y)=0. Como F(0,2)=(2,1), la ecuación vectorial de la recta normal es

(0,2)+t(2,1)=(2t,2+t),

y la ecuación vectorial de la recta tangente es

(0,2)+t(1,2)=(t,22t),

que son las mismas rectas calculadas antes.

Ejercicio 10.29 Hallar la ecuación de la recta tangente y normal a la curva x2+y2=3xy1 en los puntos en que x=1. Calcular también los extremos relativos y decir si son máximos o mínimos.

Consideremos y como función de x. Veamos primero los puntos de la curva en los que x=1:

12+y2=31y1y23y+2=0.

Resolviendo la ecuación obtenemos dos soluciones y=1 e y=2, de modo que existen dos puntos para los que x=1, que son el (1,1) y el (1,2).

Para calcular las ecuaciones de las rectas tangente y normal en estos puntos, necesitamos calcular la derivada dydx en dichos puntos, pero antes hay que asegurarse de que la ecuación

F(x,y)=x2+y23xy+1=0

define a y con función implícita de x en un entorno de los puntos (1,1) y (1,2). Veamos que se cumplen las condiciones del teorema de la función implícita.

  1. F(1,1)=12+12311+1=0 y F(1,2)=12+243121=0.
  2. Fx=2x3y y Fy=2y3x, que son continuas en todo su dominio, en particular en (1,1) y (1,2).
  3. Fy(1,1)=2131=10 y Fy(1,2)=2231=10.

Así pues, F(x,y)=0 define a y como función implícita de x en un entorno de los puntos (1,1) y (1,2), de modo que la derivada de y con respecto a x en estos puntos es

dydx(1,1)=Fx(1,1)Fy(1,1)=21312131=11=1,dydx(1,2)=Fx(1,2)Fy(1,2)=21322231=41=4.

Así pues, la ecuación de la recta tangente en el punto (1,1) es

y1=dydx(1,1)(x1)y=2x,

y la ecuación de la recta normal es

y1=1dy/dx(1,1)(x1)y=x,

mientras que la ecuación de la recta tangente en el punto (1,2) es

y2=dydx(1,2)(x1)y=4x2,

y la ecuación de la recta normal es

y2=1dy/dx(1,2)(x1)y=9x4.

Por otro lado, para calcular los extremos relativos, primero calculamos los puntos críticos, que son los que anulan la primera derivada.

dydx=3y2x2y3x=03y2x=0y=2x/3.

Pero además deben pertenecer a la curva de la función, y por tanto deben satisfacer la ecuación de la función.

x2+(2x/3)2=3x(2x/3)1x2+4x2/9=2x21x2=9/5x=±3/5.

Así pues, existen dos puntos críticos que son el (3/5,2/5) y (3/5,2/5). Para ver si son puntos de máximo o mínimo relativos, necesitamos calcular la segunda derivada en dichos puntos:

d2ydx2=ddx(3y2x2y3x)=(3dydx2)(2y3x)(2dydx3)(3y2x)(2y3x)2.

En el punto, (3/5,2/5) la segunda derivada vale

d2ydx2(3/5,2/5)=(302)(225335)(203)(325235)(225335)2=25,

que al ser positiva, indica que el punto (3/5,2/5) es un punto de mínimo relativo.

En el punto, (3/5,2/5) la segunda derivada vale

d2ydx2(3/5,2/5)=(302)(225335)(203)(325235)(225335)2=25,

que al ser negativa, indica que el punto (3/5,2/5) es un punto de máximo relativo.

Ejercicio 10.30 En una reacción química se cumple la ecuación

xln(y)+2ey2+zxxz2=1

donde x, y y z son las concentraciones de tres sustancias. ¿Cómo cambia la concentración de z si empezamos a disminuir x al mismo ritmo que se incrementa y en el instante en que x=2, y=1 y z=1?

Consideremos la función F(x,y,z)=xln(y)+2ey2+zxxz2+1 y veamos si se cumplen las condiciones del teorema de la función implícita en un entorno del punto (2,1,1).

  1. F(2,1,1)=2ln(1)+2e12122(1)2+1=0.

  2. Las derivadas parciales de F son

    Fx=ln(y)2ey2+zx21z2,Fy=xy4yey2+zx,Fz=2ey2+zx+2xz3,

    que son continuas en un entorno del punto (2,1,1).

  3. Fz(2,1,1)=2e1212+22(1)3=30.

Así pues, se cumplen la condiciones del teorema de la función implícita y F(x,y,z)=0 define a z como función implícita z=f(x,y) en un entorno de (2,1,1).

Las derivadas parciales de z son

fx(2,1,1)=Fx(2,1,1)Fz(2,1,1)=3/23=12,fy(2,1,1)=Fy(2,1,1)Fz(2,1,1)=43=43,

por lo que el vector gradiente de f es f(2,1)=(12,43).

Finalmente, la derivada parcial de f en la dirección del vector u=(1,1) es

fu(2,1)=f(2,1)u|u|=(12,43)(1,1)|(1,1)|=12+432=1162,

lo que quiere decir que por cada unidad que cambiemos las concentraciones de x e y en la dirección del vector (1,1), la concentración de z aumentará 1162 unidades.

Ejercicio 10.31 Una función f(x,y) se llama homogénea de grado n si satisface que f(tx,ty)=tnf(x,y). Demostrar que si f es homogénea de grado n se cumplen las ecuaciones

xfx+yfy=nf(x,y),x22fx2+2xy2fxy+y22fy2=n(n1)f(x,y).

En primer lugar, veamos que una función homogénea de grado n cumple que sus derivadas parciales son homogéneas de grado n1. Derivando ambos lados de la igualdad f(tx,ty)=tnf(x,y) con respecto a x se tiene

f(tx,ty)x=tnf(x,y)xf(tx,ty)txd(tx)dx=tnf(x,y)xf(tx,ty)xt=tnf(x,y)xf(tx,ty)x=tn1f(x,y)x,

por lo que fx es homogénea de grado n1, y del mismo modo se prueba que fy también lo es.

Veamos ahora que se cumple la primera ecuación. Consideremos la función vectorial g(t)=(tx,ty), entonces f(g(t))=f(tx,ty)=tnf(x,y) y derivando con respecto a t ambos lados de la ecuación se tiene

df(tx,ty)dt=d(tnf(x,y))dtf(tx,ty)txd(tx)dt+f(tx,ty)tyd(ty)dt=ntn1f(x,y)tn1f(x,y)txx+tn1f(x,y)tyy=ntn1f(x,y)f(x,y)txx+f(x,y)tyy=nf(x,y).

Como esto es cierto para cualquier t, en particular, tomando t=1 se concluye

f(x,y)xx+f(x,y)yy=nf(x,y).

La segunda ecuación se prueba de forma análoga.

Ejercicio 10.32 Demostrar que si f(x,y) es una función de varias variables tal que x e y dependen, a su vez, de otras dos variables, s y t, es decir, x=g(s,t) e y=h(s,t), entonces

fs=fxxs+fyysft=fxxt+fyyt

Dado un campo escalar f(x,y), donde (x,y) son las coordenadas cartesianas de un punto en el plano, usar el resultado anterior para calcular fr y fθ, donde r y θ son las coordenadas polares del punto.

Para calcular ft suponemos s constante, de manera que las funciones g(s,t) y h(s,t) son funciones de una variable t. Aplicando la regla de la cadena para una función de dos variables, se tiene

ft=fxdxdt+fydydt=fxxt+fyyt

Del mismo modo, para calcular fs suponemos t constante, de manera que las funciones g(s,t) y h(s,t) son funciones de una variable s, y aplicando de nuevo la para una función de dos variables, se tiene

fs=fxdxds+fydyds=fxxs+fyys

Supongamos ahora que f(x,y) es un campo escalar donde x e y son las coordenadas cartesianas de un punto en el plano euclídeo. Como el cambio de coordenadas polares a cartesianas viene dado por las funciones

x=g(r,θ)=rcos(θ)y=h(r,θ)=rsen(θ),

aplicando el resultado anterior se tiene

fr=fxxr+fyyr=fxcos(θ)+fysen(θ),fθ=fxxθ+fyyθ=fx(rsen(θ))+fyrcos(θ).

Ejercicio 10.33 La relación que modeliza el potencial eléctrico V de un punto del plano en función de su distancia, es V=logD, donde D=x2+y2.

  1. Calcular el gradiente de V.

  2. Hallar la dirección de máxima variación del potencial eléctrico en el punto (x,y)=(3,3).

  3. Calcular la matriz Hessiana y el Hessiano de V en el punto anterior.

  4. Aplicar el ejercicio anterior para calcular el gradiente y el Hessiano en coordenadas polares.

  5. Si nos movemos a lo largo de la curva y=x+1, cuál será el máximo potencial alcanzado?

  1. V(x,y)=(xx2+y2,yx2+y2).

  2. V(3,3)=3/6(1,1).

  3. 2V(x,y)=(y2x2y4+2x2y2+x42xyy4+2x2y2+x42xyy4+2x2y2+x4x2y2y4+2x2y2+x4)

    2V(3,3)=(01/61/60)

    y |2V(3,3)|=1/36.

  4. Tomando las funciones de cambio a coordenadas polares x=rcos(θ) e y=rsen(θ), por el ejercicio anterior se tiene

    Vr=Vxcos(θ)+Vysen(θ)=xx2+y2cos(θ)+yx2+y2sen(θ)=rcos(θ)2r2+rsen(θ)2r2=1rVθ=Vx(rsen(θ))+Vyrcos(θ)=xx2+y2(rsen(θ))+yx2+y2rcos(θ)=rsen(θ)cos(θ)r2+rsen(θ)cos(θ)r2=0

    Podríamos haber llegado a este mismo resultado haciendo el cambio a coordenadas polares directamente en la función del potencial ya que V(r,θ)=ln(r).

    A partir de aquí, obtener la matriz hessiana resulta más sencillo.

    2V(r,θ)=(1r2000)

  5. Sustituyendo y=x+1 en la función del potencial se tiene que

    V(x,y)=ln(x2+(x+1)2)=ln(2x2+2x+1),

    que en realidad es una función de una variable. Si calculamos su derivada se tiene

    dVdx=12x2+2x+14x+222x2+2x+1=2x+12x2+2x+1.

    por lo que el único punto crítico está en x=1/2. Estudiando el signo de la derivada a la izquierda y la derecha de este punto se observa fácilmente que se trata de un máximo relativo, por lo que el potencial máximo vale V(1/2)=log22.

Ejercicio 10.34 Una barra de metal de un metro de largo se calienta de manera irregular y de forma tal que a x metros de su extremo izquierdo y en el instante t minutos, su temperatura en grados centígrados esta dada por H(x,t)=100e0.1tsen(πxt) con 0x1.

  1. Calcular Hx(0.2,1) y Hx(0.8,1). ¿Cuál es la interpretación práctica (en términos de temperatura) de estas derivadas parciales? Explicar por qué cada una tiene el signo que tiene.

  2. Calcular la matriz hessiana de H.

  1. La derivada parcial de H con respecto a x es

    Hx(x,t)=100e0.1tcos(πxt)πt

    y en los puntos que nos piden vale

    Hx(0.2,1)=100e0.1cos(0.2π)π=229.9736Hx(0.8,1)=100e0.1cos(0.8π)π=229.9736

    La derivada parcial Hx(x0,t0) indica la variación instantánea que experimenta la temperatura con respecto a la variación de la distancia al extremo izquierdo en el punto. El signo de la derivada parcial indica si la variación de la temperatura es creciente (aumenta la temperatura) o decreciente (disminuye). Así en el punto (0.2,1) la temperatura aumentará a razón de 229.9736 grados centígrados por cada metro que nos alejemos del extremo izquierdo de la barra de metal, mientras que en el (0.8,1) la temperatura disminuirá a razón de 229.9736 grados centígrados por cada metro que nos alejemos del extremo izquierdo de la barra de metal.

  2. Para calcular la matriz Hessiana necesitamos las derivadas parciales de segundo orden:

    Ht(x,t)=100(xe0.1tsen(πxt)+e0.1txsen(πxt))==100(0.1e0.1tsen(πxt)+e0.1tcos(πxt)πx)==100e0.1t(0.1sen(πxt)+πxcos(πxt)),2Hx2(x,t)=x(100e0.1tπtcos(πxt))=100e0.1tπt(sen(πxt)πt)==100e0.1tπ2t2sen(πxt),2Htx(x,t)=t(100e0.1tπtcos(πxt))==100(te0.1tπtcos(πxt)+e0.1t(t(πt)cos(πxt)+πttcos(πxt)))==100(0.1e0.1tπtcos(πxt)+e0.1t(πcos(πxt)πtsen(πxt)πx))==100e0.1t(0.1πtcos(πxt)+πcos(πxt)π2xtsen(πxt))==100e0.1t((0.1πt+π)cos(πxt)π2xtsen(πxt)),

    2Hxt(x,t)=2Htx(x,t)(igualdad de las derivadas cruzadas por el teorema de Schwartz)2Ht2(x,t)=t(100e0.1t(0.1sen(πxt)+πxcos(πxt)))==100(te0.1t(0.1sen(πxt)+πxcos(πxt))++e0.1t(t(0.1sen(πxt))+t(πxcos(πxt))))==100(0.1e0.1t(0.1sen(πxt)+πxcos(πxt))++e0.1t(0.1cos(πxt)πxπxcos(πxt)πx))==100e0.1t(0.01sen(πxt)0.1πxcos(πxt)0.1πxcos(πxt)π2x2cos(πxt))==100e0.1t(0.01sen(πxt)(0.2+π2x2)cos(πxt)).

    Así pues, la matriz Hessiana es (100e0.1tπ2t2sen(πxt)100e0.1t((0.1πt+π)cos(πxt)π2xtsen(πxt))100e0.1t((0.1πt+π)cos(πxt)π2xtsen(πxt))100e0.1t(0.01sen(πxt)(0.2+π2x2)cos(πxt)))

Ejercicio 10.35 La ecuación diferencial parcial

2ux2+2uy2+2uz2=0,

se conoce como ecuación de Laplace y se aplica a multitud de fenómenos relacionadas con la conducción de calor, el flujo de fluidos o el potencial eléctrico.

  1. Comprobar que la función u(x,y,z)=1x2+y2+z2 satisface la ecuación de Laplace.

  2. ¿Existe algún punto en el que el crecimiento de la función sea nulo?

  3. Si fijamos z=1, calcular 4ux2y2.

  1. Para comprobar que u(x,y,z) satisface la ecuación de Laplace calculamos las tres derivadas parciales segundas que intervienen en la ecuación. Comenzando con las derivadas parciales con respecto a la variable x, obtenemos:

    ux=12(x2+y2+z2)3/22x=x(x2+y2+z2)3/22ux2=x(x(x2+y2+z2)3/2)=(x2+y2+z2)3/2+3x2(x2+y2+z2)5/2uy=y(x2+y2+z2)3/22uy2=(x2+y2+z2)3/2+3y2(x2+y2+z2)5/2uz=z(x2+y2+z2)3/22uz2=(x2+y2+z2)3/2+3z2(x2+y2+z2)5/2

    Por lo tanto,

    2ux2+2uy2+2uz2=3(x2+y2+z2)3/2+3(x2+y2+z2)(x2+y2+z2)5/2=3(x2+y2+z2)3/2+3(x2+y2+z2)3/2=0.

  2. Una condición necesaria para que el crecimiento de una función de varias variables en un punto sea nulo es que el gradiente en dicho punto se anule, y el gradiente se anula si se anulan sus tres componentes:

    u=0(ux,uy,uz)=(0,0,0)

    Por lo tanto, tenemos un sistema no lineal de tres ecuaciones con tres incógnitas:

    x(x2+y2+z2)3/2=0

    y(x2+y2+z2)3/2=0

    z(x2+y2+z2)3/2=0

    Y teniendo en cuenta que el término (x2+y2+z2), por tratarse de una suma de cuadrados, únicamente puede ser 0 si x=y=z=0; y a igual conclusión llegamos si suponemos que es distinto de 0, ya que entonces la primera ecuación implica que necesariamente x=0, la segunda implica que y=0, y la tercera implica que z=0. Por lo tanto, concluimos que el único punto en el que el crecimiento puede ser nulo es (x,y,z)=(0,0,0), pero dicho punto no pertenece al dominio de definición de la función (tendríamos un cero como denominador de una fracción), por lo que no hay ningún punto en el que la función presente un crecimiento nulo.

  3. Suponiendo z=1, la función resultante presenta únicamente dos variables:

    u(x,y,1)=1x2+y2+1=(x2+y2+1)1/2

    La derivada propuesta es:

    4ux2y2=x(x(y(uy)))

    Operando como ya hicimos en los cálculos previos de las derivadas segundas, obtenemos:

    uy=y(x2+y2+1)3/22uy2=y(uy)=(x2+y2+1)3/2+3y2(x2+y2+1)5/23uxy2=x(2uy2)=3x(x2+y2+1)5/215y2x(x2+y2+1)7/24ux2y2=x(3uxy2)=3(x2+y2+1)5/215(x2+y2)(x2+y2+1)7/2=+105x2y(x2+y2+1)9/2.

Ejercicio 10.36 Hallar los extremos relativos y los puntos de silla de la función f(x,y)=(x2+y2)22a2(x2y2) con a0.

No tiene máximos relativos.
Mínimos relativos en (a,0) y (a,0).
Punto de silla en (0,0).

Ejercicio 10.37 Determinar los extremos relativos y los puntos de silla del campo escalar h(x,y)=xy+xy222x2.

Máximo relativo en (1/8,1).
No tiene mínimos relativos.
Puntos de silla en (0,0) y (0,2).

Ejercicio 10.38 El rendimiento de una cosecha, R, depende de las concentraciones de nitrógeno, n, y fósforo, p, presentes en el suelo según la función R(n,p)=npe(n+p). ¿Cuáles deben ser las concentraciones de nitrógeno y fósforo para que el rendimiento de la cosecha sea máximo?

Se trata de un problema de optimización, así que, primero calculamos los puntos donde se anula el gradiente

Rn=e(n+p)p(1n)=0n=1 o p=0.Rp=e(n+p)n(1p)=0n=0 o p=1.

Así pues, existen dos puntos críticos, que son (0,0) y (1,1). Para ver en cuál de ellos hay un máximo relativo, calculamos el hessiano en cada uno de ellos. La matriz hessiana vale

2R(n,p)=(e(n+p)p(n2)e(n+p)(npnp+1)e(n+p)(npnp+1)e(n+p)n(p2))

El hessiano en el punto (0,0) vale

|2R(0,0)|=|0110|=1<0,

por lo que en el punto (0,0) hay un punto de silla.

El hessiano en el punto (1,1) vale

|2R(0,0)|=|e200e2|=e4>0,

por lo que en el punto (1,1) hay un extremo relativo, y como 2Rn2=e2<0 se trata de un máximo relativo. Por tanto, el rendimiento de la cosecha será máximo para n=p=1.

Ejercicio 10.39 Una empresa fabrica ordenadores portátiles y de sobremesa. El precio de los ordenadores portátiles es de 1200€ y el de los de sobremesa 800€. Si el coste de producir x ordenadores portátiles e y ordenadores de sobremesa semanales viene dado por la función C(x,y)=2x2+y2+10000ln(xy) en euros, ¿qué cantidad de ordenadores de cada tipo deberá producir semanalmente para maximizar el beneficio?

Los ingresos obtenidos por la venta de x ordenadores portátiles es 1200x y por la venta de y ordenadores de sobremesa es 800y, de modo que la función de beneficio es

B(x,y)=1200x+800y2x2y210000ln(xy).

Para obtener el máximo de esta función primero calculamos los puntos críticos.

B(x,y)=(12004x10000x,8002y10000y)=(0,0).

Resolviendo la primera ecuación se tiene x=8.58 o x=291.42, y resolviendo la segunda y=12.92 o y=387.08, por lo que hay 4 puntos críticos (8.58,12.92), (8.58,387.08), (291.42,12.92) y (291.42,387.08).

A continuación calculamos el hessiano en cada uno de ellos.

|2B(x,y)|=|4x2+10000x2002y2+10000y2|=(4x2+10000)(2y2+10000)x2y2

En el punto (8.58,12.92) se tiene |2B(8.58,12.92)|=7634.37>0 y 2Bx2(8.58,12.92)=131.8393>0, por lo que en se trata de un mínimo relativo.

En el punto (8.58,387.08) se tiene |2B(8.58,387.08)|=254.88<0, por lo que se trata de un punto de silla.

En el punto (291.42,12.92) se tiene |2B(291.42,12.92)|=224.81<0, por lo que se trata de otro punto de silla.

En el punto (291.42,387.08) se tiene |2B(291.42,387.08)|=7.51>0 y 2Bx2(291.42,387.08)=3.88<0, por lo que en se trata de un máximo relativo.

Así pues, para maximizar el beneficio deben fabricarse semanalmente 291.42 ordenadores portátiles y 387 ordenadores de sobremesa.

Ejercicio 10.40 En el ajuste de regresión de una recta y=a+bx, se suele utilizar la técnica de mínimos cuadrados que consisten en buscar los valores de a y b que hacen mínima la función

f(a,b)=i=1n(yiabxi)2,

donde el sumatorio abarca a todos los pares de la muestra (xi,yi) para i=1,,n, siendo n el tamaño de la muestra.

Demostrar que esta función alcanza el mínimo en el punto

a=y¯bx¯ y b=sxysx2.

Ejercicio 10.41 Una empresa que fabrica cajas de cartón quiere construir una caja rectangular sin tapa con un volumen de 25 litros. ¿Qué dimensiones debe tener la caja para emplear la menor cantidad de cartón posible?

El volumen de una caja de base rectangular se obtiene multiplicando sus dimensiones V(x,y,z)=xyz, donde x e y son las dimensiones del rectángulo de la base y z la altura.

Por otro lado, el area de su superficie es S(x,y,z)=xy+2xz+2yz. Así pues, en el punto donde se maximize esta función debe cumplirse, por el método de los multiplicadores de Lagrange, las siguientes condiciones:

V(x,y,z)=25000 cm3S(x,y,z)=kV(x,y,z)

Si calculamos los gradientes de V y S se tiene

S(x,y,z)=(y+2z,x+2z,2x+2y)V(x,y,z)=(yz,xz,xy),

de manera que se deben cumplir las ecuaciones

y+2z=kyzx+2z=kxz2x+2y=kxy

Para resolver el sistema, multiplicamos la primera ecuación por x, la segunda por y y la tercera por z, para que las tres ecuaciones tengan el mismo lado derecho, de donde se deduce que

xy+2xz=xy+2yz=2xz+2yz.

De la primera ecuación se deduce fácilmente que x=y, y sustituyendo en la segunda, se tiene que x2=2xz, que tiene solución x=0 o z=x/2. Como la primera opción daría un volumen nulo, no puede ser, así que solo tendría sentido la segunda. Finalmente como el volumen debe ser 25000 cm3, sustituyendo las soluciones en la fórmula del volumen se tiene

V(x,y,z)=25000xxx2=25000x3=50000x36.8403,

por lo que las dimensiones de la caja son x=y=36.8403 cm y z=36.8403/2=18.4201 cm.

Ejercicio 10.42 Un cometa sigue la órbita dada por la función 3x2+y2=9, donde el origen de coordenadas es la posición del sol. Si una nave espacial se mantiene fija en la posición (4,1), ¿en qué punto de su trayectoria el cometa estará más cerca del nave?

La distancia del cualquier punto (x,y) de la trayectoria del cometa a la posición de la nave (4,2) es f(x,y)=(x4)2+(y2)2, y si llamamos g(x,y)=3x+y2, por el método de los multiplicadores de Lagrange, se tienen que cumplir las ecuaciones

g(x,y)=9f(x,y)=kg(x,y)

Calculando los gradientes de f y g,

f(x,y)=((x4)(x4)2+(y2)2,(y2)(x4)2+(y2)2)g(x,y)=(3,2y),

se tienen que cumplir las ecuaciones

(x4)(x4)2+(y2)2=k3(y2)(x4)2+(y2)2=k2y,

de donde se deduce que

(x4)3(x4)2+(y2)2=(y2)2y(x4)2+(y2)2(x4)3=(y2)2yx=1123y.

Como además tiene que cumplirse la ecuación de la trayectoria del cometa, se tiene que

g(1123y,y)=93329y+y2=9y3+15y18=0y1.047,

y, por tanto, x11231.047=2.6347.

Ejercicio 10.43 La función de producción de un determinado producto es f(k,l)=kl3 donde k son las unidades de capital y l las unidades de mano de obra. Si el coste por unidad de capital es 7€ y el coste por unidad de mano de obra es 6€, determinar la cantidad de capital y mano de obra que minimiza el coste si deben producirse 290 unidades del producto.

La función de coste, para k unidades de capital y l unidades de mano de obra es c(k,l)=7k+6l, así que se trata de minimizar esta función con la restricción f(k,l)=290. Aplicando el método de los multiplicadores de Lagrange, debe cumplirse

f(k,l)=290c(k,l)=λf(k,l)

Si calculamos los gradientes de f y g se tiene

f(k,l)=(l3,3kl2)c(k,l)=(7,6),

de manera que deben cumplirse las ecuaciones

l3=7λ3kl2=6λ,

de donde se deduce que

l37=3kl26k=2l7.

Finalmente, sustituyendo en la restricción se tiene

f(2l7,l)=2902l7l3=290l4=1015l=101545.6444,

y k25.64447=1.6127.

Ejercicio 10.44 Calcular los polinomios de Maclaurin de segundo grado para las siguientes funciones.

  1. f(x,y)=ln(x+y+1).
  2. g(x,y)=ex+y.
  3. h(x,y)=sen(x)sen(y).
  1. pf,(0,0)2(x,y)=x+yx22xyy22.

  2. pg,(0,0)2(x,y)=1+x+y+x22+xy+y22.

  3. ph,(0,0)2(x,y)=xy.

Ejercicio 10.45 Un modelo ecológico explica el número de individuos de una población mediante la función

f(x,t)=etx,

donde t es el tiempo y x el número de predadores en la región.

Calcular el número aproximado de individuos en la población para t=0.1 y x=0.9 utilizando el polinomio de Taylor de segundo grado de la función en el punto (1,0).

La fórmula del polinomio de Taylor de grado 2 para la función f(x,t) en el punto (1,0) es

Pf,(1,0)2(x,t)=f(1,0)+f(1,0)(x1,y)+12(x1,y)2f(1,0)(x1y)

Calculando las primeras derivadas parciales obtenemos fácilmente el gradiente

f(x,t)=(etx2,etx),

que en el punto (1,0) vale

f(1,0)=(e0,e0)=(1,1).

Calculando las segundas derivadas parciales obtenemos la matriz hessiana

2f(x,t)=(2etx3etx2etx2etx)

que en el punto (1,0) vale

2f(1,0)=(2e0e0e0e0)=(2111).

Sustituyendo en la fórmula del polinomio de Taylor de segundo grado se tiene

Pf,(1,0)2(x,t)=1(x1)+t+(x1)2(x1)t+12t2=x23x+2txt+12t2+3.

Finalmente, sustituyendo x=0.9 y t=0.1 se tiene que el número aproximado de individuos para estos valores es

Pf,(1,0)2(0.9,0.1)=0.9230.9+20.10.90.1+120.12+3=1.2245.