TAREA DE CALCULO III

Leidy Katherin Narvaez (2140037-3753)

Universidad del Valle Cali, Colombia

E-mail:
[email protected]

Fecha: 11/06/2024

Use una computadora para graficar la superficie, el plano tangente y la recta

normal en la misma pantalla. Elija cuidadosamente el dominio para que evite
planos verticales extraños. Seleccione el punto de vista de tal manera que obtenga
una vista óptima de los tres objetos.

47. x y+ y z + zx=3 , ( 1 ,1 , 1 )

SOLUCIÓN:

Para el plano tangente necesitamos dos cosas un vector normal y el punto, un

vector perpendicular es el gradiente, que son las derivadas parciales de cada
componente.

∇ F= ( x , y , z )= ⟨ y + z , x + z , y + x ⟩ (1)

Se evalúa el gradiente en el punto ( 1 ,1 , 1 )

η⃗ =∇ F (1 , 1 ,1 ) =⟨ 2 , 2 ,2 ⟩

a ( x−x 0 ) + b ( y− y 0 ) + c ( z−z 0 ) =0 (2)

Reemplazamos los valores en la ecuación del plano tangente

2 ( x−1 ) +2 ( y−1 ) +2 ( z−1 )=0

Distribuimos el dos en los tres paréntesis

 2 x−2+2 y−2+2 z−2=0

2 x+ 2 y +2 z−6=0

x + y + z−3=0

x + y + z=3 Plano tangente

Se utilizo la ecuación de la recta normal

L :P+ t ⃗d (3)

( 1 ,1 , 1 ) +t ( 2 , 2 ,2 )

( 1 ,1 , 1 ) + ( 2 t , 2t , 2t )

⟨ 1+ 2t ,1+2 t , 1+2 t ⟩ Ecuación

paramétrica

Las ecuaciones paramétricas se plantean así:

{
x=1+ 2t
y=1+2 t
z=1+2 t

También se puede ver así:

x−1
x=1+2 t → t=
2

y−1
y=1+2 t → t=
2

z−1
z=1+2 t → t=
2

Siendoigual para x , y , z → t=x−1= y −1=z−1 Ecuación cartesiana

Ahora se procede a graficar los tres objetos en wólfram, como se observa en la
figura 1, en el que la recta interseca a la superficie y al plano, con el siguiente
código de programación:

f[x_, y_, z_] := x*y + y*z + z*x - 3;

pt = {1, 1, 1};

tangentPlane =
With[{grad = Grad[f[x, y, z], {x, y, z}] /. Thread[{x, y, z} -> pt]},
grad.{x, y, z} == grad.pt];

normalLine =
With[{grad = Grad[f[x, y, z], {x, y, z}] /. Thread[{x, y, z} -> pt]},
pt + t*grad];

surfacePlot =
ContourPlot3D[f[x, y, z] == 0, {x, -10, 10}, {y, -10, 10}, {z, -10, 10},
ContourStyle -> Opacity[0.5], Mesh -> None];

tangentPlanePlot =
ContourPlot3D[Evaluate[tangentPlane], {x, -10, 10}, {y, -10, 10}, {z, -10, 10},
ContourStyle -> Directive[Blue, Opacity[0.5]], Mesh -> None];

normalLinePlot =
ParametricPlot3D[Evaluate[normalLine], {t, -10, 10},
PlotStyle -> Directive[Black, Thick]];

Show[surfacePlot, tangentPlanePlot, normalLinePlot,

PlotRange -> All, Axes -> True, AxesLabel -> {"x", "y", "z"},
Boxed -> False, ViewPoint -> {1.5, -2.5, 1.5},
ViewVertical -> {0, 0, 1}]
 Figura 1. La superficie, el plano tangente y la recta normal en la misma pantalla

59. Suponga que un científico tiene razones para creer que dos cantidades x & y
están relacionadas linealmente, es decir que y=mx+b , al menos
aproximadamente, para algunos valores de m y b. El científico realiza un
experimento y recolecta datos bajo la forma de puntos, ( x 1 , y 1 ) , ( x 2 , y 2 ) , … , ( x n , y n ),
tras de lo cual traza estos puntos. Los puntos no se sitúan exactamente en una
recta, así que el científico quiere determinar constantes m y b de tal forma que la
recta y=mx+b “se ajuste” lo más posible a los puntos (véase la figura).
Sea d i= yi −( mxi +b ) la desviación vertical del punto respecto a la recta. El método
que se utilizo fue el de mínimos cuadrados y el cual determina a m y b de tal modo
n
que se minimice ∑ d 2i , la suma de los cuadrados de estas desviaciones.
i=1

Demuestre que, de acuerdo con este método, la recta que mejor se ajuste se
obtiene cuando

n n
m ∑ x1 +bn=∑ y 1
i=1 i=1

(4)

n n n
m ∑ x +b ∑ x1=∑ xi y i
2
i
i=1 i=1 i=1

(5)

SOLUCIÓN:

Así, la recta se determina resolviendo en estas dos ecuaciones y las dos

incógnitas m y b. (Véase la sección 1.2 para un análisis adicional y aplicaciones
del método de mínimos cuadrados)

Para iniciar se observa el gráfico donde tenemos cierta cantidad de puntos de x i y

y i, que no están en la misma recta, entonces el método de los mínimos cuadrados
nos permite encontrar una recta que se aproxime mejor a los puntos como lo es la
ecuación y=mx+b .

Para el valor de x i se tiene dos puntos, es decir punto inicial ( x i , y i ) y otro punto
( x i , m x i +b ). Ahora se define a d i, como la distancia del punto (x ¿ ¿ i, y i )¿ al punto
n
( x i , m x i +b ). Ahora se quiere minimizar ∑ d 2i , como ya se conoce que
i=1
n
d i= yi −( m x i +b ), se reemplaza en ∑ ( y i−( m x i+ b ) )2, que depende de m y b, y así
i=1
n
encontrar el valor de m y b, tal que f ( m, b ) =∑ ( y i−( m x i +b ) ) tome el valor mínimo.
2

i=1

Para esto se calculan las derivadas parciales con respecto a m y b:

 n n
∂f
=∑ 2 ( y i−( m xi +b ) )∗(−x i )=¿−2 ∑ ( x i yi −m x 2i −b x i ) =0 ( 6 ) ¿
∂ m i=1 i=1

n n
∂f
=∑ 2 ( y i− ( m x i +b ) )∗(−1 )=¿−2 ∑ ( y i−m x i−b )=0 ¿
∂ b i=1 i=1

Por esto se obtiene un punto crítico (m, b).

Después se realizó el cálculo del determinante, para esto se determinan las

derivadas parciales de segundo orden.

n n
fmm=−2 ∑ −x2i =2 ∑ x 2i >0
i=1 i=1

n n
fbb=−2 ∑ (−1 )=2 ∑ 1=2 n
i=1 i=1

n n
fmb=−2 ∑ −x i=2 ∑ x i
i=1 i=1

( ) (∑ )
n n 2 n n 2
D=
fmm fbm
fmb fbb( i=1
)
=2 ∑ x i ∗ (2 n )− 2 ∑ x i =4 n ∑ x i −4
2

i=1 i=1
2

i=1
x i >0

Se identifica que para el determinante se puede utilizar la desigualdad de Cauchy-

Schwarz.

( ) ( )( )
n 2 n n

∑ ai bi ≤ ∑ a2i ∑ b2i
i=1 i=1 i=1

(7)
Para esto se definen los términos:

( ) ( )( ) ( )
n 2 n n n n 2 n

∑ 1∗x i ≤ ∑ 12 ∑ x 2i =n∗∑ x 2i → ∑ x i ≤ n ∑ x 2i
i=1 i=1 i=1 i=1 i=1 i =1

Finalmente se identifica que el parámetro (m, b) es un mínimo, por medio de la

resolución de la prueba de la segunda derivada, que indica que si D > 0 y fmm (m,
b) > 0, entonces f (m, b) es un mínimo local, en resumen, de lo mencionado
anteriormente:
fmm > 0, D > 0

Definición Una función de dos variables tiene un máximo local en (a, b) si f(x, y) ≤
f(a, b) cuando (x, y) está cerca de (a, b). [Esto significa que f (x, y) ≤ f (a, b) para
todos los puntos (x, y) en algún disco con centro (a, b)] El número f(a, b) se llama
valor máximo local. Si f(x, y)≥ f(a, b) cuando (x, y) está cerca de (a, b) f tiene un
mínimo local en (a, b) y f(a, b) es un valor mínimo local.

46. Grafique el sólido que se ubica entre las superficies:

z=e− x cos ( x2 + y 2 ) ∧z=2−x 2− y 2, para |x|≤1 , | y|≤ 1. Use un sistema algebraico
2

computacional para aproximar el volumen de este sólido con cuatro decimales.

SOLUCIÓN:

Se tienen las dos superficies y los límites para los parámetros de x & y

z=e− x cos ( x2 + y 2 )
2

2 2
z=2−x − y

|x|≤1 ⟶−1≤ x ≤ 1

| y|≤ 1 ⟶−1≤ y ≤ 1

Con estos criterios establecidos se calculó el volumen

1 1
v=∫ ∫ ( 2−x − y −e cos ( x + y ) ) dy dx
2
2 2 −x 2 2

−1 −1

Se separaron las dos integrales

1 1 1 1

∫ ∫ ( 2−x 2− y 2 ) dy dx−∫ ∫ ( e−x cos ( x 2+ y 2 ) ) dy dx

2

−1 −1 −1 −1

Primero, se realiza la primera integral dada por

1 1

∫ ∫ ( 2−x 2− y 2 ) dy dx
−1 −1
 )|
1

(
1 3
y
¿∫
2
2 y −x y− dx
−1 3 −1

(( ))
1
¿ ∫ 2−x 2−
−1
1
3 )(
− −2+ x 2+
1
3
dx

−1
2
(
1
3
2 1
¿ ∫ 2−x − + 2−x − dx
3 )
1

−1
(
¿ ∫ 4−2 x −
2 2
3
dx )
1
¿∫
−1
( 103 −2 x ) dx2

)|
1

(
3
10 2x
¿ x− dx
3 3 −1

¿ ( 103 − 23 )−( −103 + 32 )

¿ ( 83 )−(−83 )= 163
Se determino la integral por medio de wólfram ya que es una integral compleja
para desarrollarla a mano.

1 1

∫ ∫ ( e−x cos ( x 2+ y 2 ) ) dy dx
2

−1 −1

Este es el código en wólfram Cloud para encontrar el resultado de esta integral:

(* Definición de las funciones *)

f1[x_, y_] := Exp[-x^2] Cos[x^2 + y^2]
f2[x_, y_] := 2 - x^2 - y^2

(* Cálculo del volumen *)

volume = NIntegrate[f2[x, y] - f1[x, y], {x, -1, 1}, {y, -1, 1},
 Method -> {"GlobalAdaptive", "SymbolicProcessing" -> 0,
"SingularityHandler" -> None}, MaxRecursion -> 100,
PrecisionGoal -> 6]

(* Mostrar el resultado *)
volume

El resultado final de esta integral es el volumen:

1 1 1 1

∫ ∫ ( 2−x 2− y 2 ) dy dx−∫ ∫ ( e−x cos ( x 2+ y 2 ) ) dy dx=3.02707

2

−1 −1 −1 −1

Este es el código en Wólfram Cloud de la gráfica (Figura 2), Gráfica del sólido:
Plot3D se utiliza para graficar las dos superficies y rellenar el área entre ellas, con
diferentes niveles de opacidad para visualizar claramente el sólido.

Plot3D[{Exp[-x^2] Cos[x^2 + y^2], 2 - x^2 - y^2}, {x, -1, 1}, {y, -1, 1},
PlotStyle -> {Directive[Opacity[0.5], Blue], Directive[Opacity[0.5], Red]},
PlotLegends -> {"z = Exp[-x^2] Cos[x^2 + y^2]", "z = 2 - x^2 - y^2"},
AxesLabel -> {"x", "y", "z"},
PlotRange -> {{-1, 1}, {-1, 1}, Automatic}, ImageSize->Large ]

Figura 2. Solido entre superficies