Calculo Iii - M. Chacon PDF

Material complementario
para el curso
C
alculo III
con uso de software
MA1003
Escuela de Matem
atica
Universidad de Costa Rica
Lic. Mynor Chac

on Daz
2012
0.0
c
2012
Mynor Chac
on Daz
MA1003
Indice general
1. Introducci
on General
2. Operaciones con vectores
7
2
2.1. Representacion de un punto en R y R
. . . . . . . . . . . . . . .
2.2. Representacion de un vector en R2 y R3 . . . . . . . . . . . . . . .
10
2.3. Adicion, sustraccion y magnitud(o norma) de vectores en la TI-89
12
2.4. Vector unitario . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
15
2.5. Producto Punto o Escalar . . . . . . . . . . . . . . . . . . . . . . .
16
2.6. Producto Vectorial o Producto Cruz . . . . . . . . . . . . . . . . .
17
2.7. Uso de la tecla STO> . . . . . . . . . . . . . . . . . . . . . . . . . .
18
3. Operaciones con matrices
21
3.1. Representacion, adicion y producto de matrices . . . . . . . . . . .
21
3.2. Transpuesta de una matriz
. . . . . . . . . . . . . . . . . . . . . .
23
3.3. Determinante de una matriz . . . . . . . . . . . . . . . . . . . . . .
24
3.4. Metodo de Gauss-Jordan
. . . . . . . . . . . . . . . . . . . . . . .
26
3.5. Uso del comando soluc() en el caso de infinitas soluciones . . . .
28
3.6. Valores Propios y Vectores Propios de una matriz . . . . . . . . . .
31
4. Superficies en R3
35
4.1. Representacion de superficies generadas por curvas planas, con

GenGCF . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.2. Representacion de superficies cuadricas, con GenGCF . . . . . . .
x2 y 2 z 2
4.2.1. Elipsoide: 2 + 2 + 2 = 1 . . . . . . . . . . . . . . . . . .
a
b
c
x2 y 2 z 2
4.2.2. Hiperboloide de una hoja: 2 + 2 2 = 1 . . . . . . . .
a
b
c
1
35
37
37
39
0.0
MA1003
x2 y 2 z 2
2 2
a2
b
c
2
x
y2
Paraboloide Elptico: z = 2 + 2 . . .
a
b
x2 y 2
Paraboloide Hiperbolico: z = 2 2 .
a
b
x2 y 2
2
Cono Elptico: z = 2 + 2 . . . . . .
a
b
Cilindros . . . . . . . . . . . . . . . . .
4.2.3. Hiperboloide de dos hojas:
=1 . . . . . . . .
40
4.2.4.
. . . . . . . . . . .
43
. . . . . . . . . . .
44
. . . . . . . . . . .
46
. . . . . . . . . . .
48
4.2.8. Cilindro con Directriz en forma parametrica . . . . . . . . .
51
4.2.9. Conos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
54
4.2.10. Superficie de Revolucion . . . . . . . . . . . . . . . . . . . .
58
4.2.11. Un caso particular . . . . . . . . . . . . . . . . . . . . . . .
61
4.2.5.
4.2.6.
4.2.7.
5. Funciones Vectoriales
63
5.1. Definicion y propiedades . . . . . . . . . . . . . . . . . . . . . . . .
63
5.2. Derivacion e integracion . . . . . . . . . . . . . . . . . . . . . . . .
63
5.3. Curvas en R2 , R3 . . . . . . . . . . . . . . . . . . . . . . . . . . . .
68
5.4. Parametrizacion
. . . . . . . . . . . . . . . . . . . . . . . . . . . .
69
5.4.1. Parametrizacion de y = f (x) . . . . . . . . . . . . . . . . .
70
5.4.2. Parametrizacion de F (x, y) = 0 . . . . . . . . . . . . . . . .
71
5.4.3. Parametrizacion de curvas en
R3
. . . . . . . . . . . . . . .
72
5.5. Dibujo de curvas con Winplot y GenGCF . . . . . . . . . . . . . .
75
5.6. El Triedro Intrnseco o Marco de Frenet . . . . . . . . . . . . . . .
78
5.7. Curvatura, Torsion y Longitud de arco . . . . . . . . . . . . . . . .
82
6. Funciones de Varias Variables o Campos Escalares
89
6.1. Graficacion de funciones de dos variables . . . . . . . . . . . . . . .
89
6.2. Dominio de una funcion de dos variables . . . . . . . . . . . . . . .
91
6.3. Lmite de una funcion de dos variables . . . . . . . . . . . . . . . .
93
6.4. Derivadas Parciales . . . . . . . . . . . . . . . . . . . . . . . . . . .
97
6.5. Derivadas Parciales de orden superior . . . . . . . . . . . . . . . . 100

6.6. Diferenciabilidad de z = f (x, y) . . . . . . . . . . . . . . . . . . . . 101
6.6.1. Aproximacion lineal de z = f (x, y) . . . . . . . . . . . . . . 104
6.6.2. Diferencial Total de z = f (x, y) . . . . . . . . . . . . . . . . 106
6.7. Regla de la cadena . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
6.7.1. Regla de la Cadena involucrando funciones no explcitas . . 111
6.7.2. Derivacion Implcita . . . . . . . . . . . . . . . . . . . . . . 117
6.8. Derivada direccional, interpretacion geometrica. . . . . . . . . . . . 118
6.9. Plano Tangente y vector normal a una superficie . . . . . . . . . . 122
6.9.1. Valores extremos de la derivada direccional . . . . . . . . . 125
c
2012
Mynor Chac
on Daz
MA1003
0.0
7. Extremos de Campos Escalares

129
7.1. Extremos Relativos . . . . . . . . . . . . . . . . . . . . . . . . . . . 129
7.2. Extremos Absolutos . . . . . . . . . . . . . . . . . . . . . . . . . . 144
7.3. Extremos Condicionados . . . . . . . . . . . . . . . . . . . . . . . . 148
8. Integraci
on M
ultiple
8.1. Integral doble . . . . . . . . . . . . . . . . . . . .
8.1.1. Integral doble en region rectangular . . .
8.1.2. Integral doble en una region general R . .
8.1.3. Integral doble en coordenadas polares . .
8.2. Integral triple . . . . . . . . . . . . . . . . . . . .
8.2.1. Coordenadas cilndricas . . . . . . . . . .
8.2.2. Coordenadas esfericas . . . . . . . . . . .
8.3. Cambio de Variables mas generales . . . . . . . .
8.3.1. Trasformaciones con Winplot . . . . . . .
8.4. Formulas para momentos, masa y centro de masa
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
9. Integral de Lnea e Integral de Superficie

9.1. Integral de Lnea . . . . . . . . . . . . . . . . . . . . . . . .
9.1.1. Integral de lnea de un campo escalar . . . . . . . .
9.1.2. Aplicaciones fsicas . . . . . . . . . . . . . . . . . . .
9.1.3. Integral de lnea de un campo vectorial . . . . . . .
9.1.4. Campos conservativos. . . . . . . . . . . . . . . . . .
9.2. Teorema de Green . . . . . . . . . . . . . . . . . . . . . . .
9.2.1. Area
de una region plana, usando Teorema de Green
9.2.2. Teorema de Green para una region con agujero . . .
9.3. Integral de superficie . . . . . . . . . . . . . . . . . . . . . .
9.3.1. Consideraciones basicas y Producto Fundamental . .
9.3.2. Ejemplos de parametrizaciones de algunas superficies
9.3.3. Calculo de areas de superficies . . . . . . . . . . . .
9.3.4. Calculo de Integrales de superficie . . . . . . . . . .
9.3.5. Superficies orientables . . . . . . . . . . . . . . . . .
9.3.6. Calculo de Integrales de Flujo . . . . . . . . . . . . .
9.4. Teoremas de Stokes y de Gauss . . . . . . . . . . . . . . . .
9.4.1. Teorema de Stokes . . . . . . . . . . . . . . . . . . .
9.4.2. Teorema de Gauss o de la divergencia . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
159
159
159
163
178
188
193
202
206
209
219
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
227
227
227
229
236
248
255
261
263
271
271
272
278
285
290
291
294
294
305
c
2012
Mynor Chacon Daz
0.0
c
2012
Mynor Chac
on Daz
MA1003
CAPITULO
Introduccion General
Este material representa un esfuerzo por incorporar programas computacionales

de uso libre a la ense
nanza del Calculo en varias variales. Es un intento mas por
hacer mas comprensible y por mejorar el aprendizaje de las tecnicas de Calculo.
El software que se utiliza esta compuesto por:
El graficador Winplot
Un Emulador de la calculadora cientfica Texas Instrument TI-89
Una version light de MuPAD
El graficador GenGFC
El uso de tales programas esta en funcion de los temas que cubre el curso
MA1003-Calculo III. La idea de este trabajo es introducir el uso del software de
acuerdo a la necesidad del tema que se este abarcando en determinado momento.
En los captulos 2 y 3 se hace un repaso de conceptos de Algebra Lineal, haciendo
uso de la calculadora cientfica Texas Instrument TI-89 y MuPAD para la solucion
de sistemas de ecuaciones, y las operaciones con matrices. El graficador Winplot
es de gran utilidad en la graficacion de vectores.
El captulos 4 abarca el estudio de las superficies cuadricas, cilindros, conos y superficies de revolucion, en esta parte se hace imprescindible el uso del graficador
GenGCF, ademas del Winplot.
En el captulo 5, se hace el estudio de las funciones vectoriales. Para la derivacion
e integracion de tales funciones es de gran utilidad la TI-89. Seguidamente se da
enfasis a la parametrizacion de curvas planas y en el espacio. La graficacion de
curvas se hace con la ayuda del GenGCF y Winplot. Luego se trata el Triedro
5
1.0
MA1003
Intrnseco o Marco de Frenet, la curvatura, longitud de arco y torsion de una curva. Conociendo la parametrizacion de una curva, Winplot brinda la posibilidad
de obtener, de forma automatica, el valor de la longitud de arco y la curvatura.
En el captulo 6, se estudian los campos escalares. De manera breve se tratan
los temas de dominio maximo, continuidad y lmites de funciones de dos variables z = f (x, y), en esta parte los graficadores son de gran ayuda para visualizar regiones que representan dominios y para la graficacion de funciones de dos
variables. La TI-89 y Mupad brindan herramientas valiosas para el calculo, comprension y revision de derivadas parciales de orden superior, regla de la cadena
y derivadas de funciones implcitas. En el tema de optimizacion de funciones sin
y con condiciones (metodo de Lagrange), el uso integrado de todos los programas mencionados permite enfocar los problemas de una manera diferente (m
as
geometrica)a la tradicional, permite tambien el uso de metodos que generalmente
no se utilizan, debido a lo largo de los procedimientos.
El estudio de la integracion de campos escalares se realiza en el captulo 7, incluye
integrales dobles, integrales triples y sus aplicaciones; ademas se estudian los cambios de variables en integrales m
ultiples. En este captulo se hace un uso discreto
de la TI-89 y Mupad con respecto al calculo de integrales m
ultiples. Winplot brinda una valiosa herramienta para la transformacion automatica de regiones dado
un cambio de varibales.
El octavo y u
ltimo captulo, comprende el estudio de Campos vectoriales, incluye
integrales de lnea, independencia de la trayectoria y el teorema de Green, as como integrales de superficie, teorema de la divergencia y el teorema de Stokes.
En esta parte se puede utilizar todas las herramientas tecnologicas que se hayan
aprendido durante los captulos anteriores.
Se inciste en que, lo expuesto en este material no es un tutorial del software utilizado, consiste mas bien en dar las pautas y reglas para utilizar las herramientas.
Como descargar el software y todo lo referente a su uso inicial y basico se hara en
las lecciones presenciales en el laboratorio.
Ademas del contenido matematico, el estudiante solamente tendra que aprender
a utilizar los programas, no necesita de conocimientos en programacion.
c
2012
Mynor Chac
on Daz
CAPITULO
Operaciones con vectores
En este captulo se revisaran las operaciones basicas con vectores. Se utilizara para
ello el emulador de la TI-89 y el graficador Winplot.
2.1.
Representaci
on de un punto en R2 y R3
Para representar el punto A = (2, 4) en el graficador Winplot, siga la

ruta Ventana >2-dim .
Luego aparece la ventana interior conteniendo el sistema de ejes cartesianos,

desde la barra del men
u siga la ruta Ecua >Punto >(x,y)...
7
2.1
MA1003
Aparece la ventana donde se digitan las coordendas del punto, se elige su

tama
no y color.
El punto debe verse como sigue.
Analogamente para representar el punto A = (1, 2, 3), se sigue la ruta

Ventana >3-dim y luego Ecua >Punto >Cartesiana ....
c
2012
Mynor Chac
on Daz
MA1003
2.2
Se eligen las caractersticas del punto.
El punto debe verse como sigue.
c
2012
Mynor Chacon Daz
2.2
2.2.
MA1003
Representaci
on de un vector en R2 y R3
Sean A = (2, 4) y B = (3, 1) dos puntos, el vector v =AB= B A =

(1, 3) puede ser representado geometricamente utilizando Winplot. En la
siguiente ventana apareceran representados los puntos y el vector v =AB,

desde A hasta B, ademas el mismo vector anclado en el origen cuya representacion es v = (1, 3). Una vez representados los puntos, se sigue la ruta
Ecua >Segmento >(x,y)...
En la ventana segmento, las coordenadas de los puntos (x1 , y1 ) y (x2 , y2 )

son las de los puntos A y B en cualquier orden , en este caso (x1 , y1 ) = (2, 4)
y (x2 , y2 ) = (3, 1). Por otro lado el ancho del l
apiz se elige como 2, para
aumentar el grosor de la lnea del segmento y que se aprecie mejor. Puesto
que el vector v = (1, 3) tiene punto final B = (3, 1), entonces en la opcion
flechas: se elige la alternativa p2, esto hace que el segmento lleve la punta
de flecha en B.
c
2012
Mynor Chac
on Daz
10
MA1003
2.2
Se obtiene as la imagen siguiente:
Para graficar el vector anclado se repite el proceso anterior, dibujando el

segmento que va de (0,0) hasta (1,-3).
La imagen final es la siguiente:
11
c
2012
Mynor Chacon Daz
2.3
MA1003
Actividad 1
Sean A = (2, 2, 2) y B = (3, 3, 3). Represente los puntos A,B, el vector
v =AB y el vector anclado en el origen, esto es v = (1, 1, 1). Haga girar el

sistema de ejes de manera que se vea como en la siguiente figura:
2.3.
Adici
on, sustracci
on y magnitud(o norma) de vectores en la TI-89
En los siguientes ejmplos se utilizaran vectores en R2 y R3 indistintamente. El vector (a,b,c) en la TI-89 se representa como [a,b,c]. La siguiente
imagen representa la suma (-1,2,3)+(4,0,-2)=(3,2,1)
5
El producto del vector (-1,2,3) por el escalar 5/2, es decir (1, 2, 3) =
2

5
15
, 2,
, se aprecia como sigue:
2
2
c
2012
Mynor Chac
on Daz
12
MA1003
2.3
Seguidamente se va a calcular la magnitud

de unvector, kvk. Como ejemplo

15
5
. Para ello se debe seguir
se va a calcular la norma del vector , 2,
2
2
la ruta 2nd >Math >Matrices >Normas >norm(
Luego de presionar Enter, se obtiene:
Para no tener que digitar de nuevo el vector use la flecha hacia arriba para
marcarlo
13
c
2012
Mynor Chacon Daz
2.3
MA1003
Una vez marcado el vector al que se le va a calcular la magnitud.
Se presiona Enter, y se cierra el parentesis redondo que se abrio automaticamente con el comando norm(.
c
2012
Mynor Chac
on Daz
14
MA1003
2.4

5
15
5 14

Seguidamente se obtiene la norma , 2,
=
2
2
2
2.4.
Vector unitario
Dado el vector v Rn , un vector con la misma direccion pero de magniv

tud uno es el vector u =
. Para obtener el vector unitario con la direccion
kvk
del vector (1,1,1), siga la ruta
2nd >Math >Matrices >Operac Vectores >vectUnid(
digte el vector, cierre el parentesis redondo y luego de presionar Enter se
obtiene, el vector unitario
v
u=
=
kvk
1 1 1
, ,
3 3 3
15
c
2012
Mynor Chacon Daz
2.5
MA1003
2.5.
Producto Punto o Escalar
Dado los vectores v Rn , a = (a1 , a2 , ..., an ) y b = (b1 , b2 , ..., bn ) se define

el producto punto entre ellos como
ab=
n
Y
ai b i
i=1
Para efectuar a b, siendo a = (2, 3, 5) y b = (0, 2, 1), siga la ruta

2nd >Math >Matrices >Operac Vectores >producEsc(
Luego digite los vectores separados por una coma, producEsc([2,3,-5],[0,2,1]),

se presiona Enter y se obtiene el resultado -1.
c
2012
Mynor Chac
on Daz
16
MA1003
2.6.
2.6
Producto Vectorial o Producto Cruz
Dado los vectores u = (1, 2, 3), v = (3, 1, 7) R3 , para calcular el

producto vectorial entre ellos se sigue la ruta
2nd >Math >Matrices >Operac Vectores >producVect(
Luego se digitan los vectores separados por una coma, producVect([1,2,3],[3,-1,7]),

se presiona Enter y se obtiene el resultado (1, 2, 3) (3, 1, 7) = (17, 2, 7).
17
c
2012
Mynor Chacon Daz
2.7
2.7.
MA1003
Uso de la tecla STO>
A veces es necesario realizar distintos calculos utilizando algunos vectores fijos. Para no digitarlos cada vez que se utilicen, se van a almacenar
con un nombre determinado. Por ejemplo se va a almacenar los vectores
(1, 2, 3), (3, 1, 7) a las letras (nombres) u, v respectivamente. Para esto dig, seguido de la tecla
y
ite el vector [1,2,3] luego presione la tecla
la letra u, ubicada en la esquina superior derecha de la tecla
. Luego de
presionar Enter, se debe ver la siguiente pantalla.
De igual manera se almacena el vector (3, 1, 7) en la letra v.

c
2012
Mynor Chac
on Daz
18
MA1003
2.7
Si se desea el valor 3u5v+2(uv), simplemente se digita 3*u-5*v+2*prodVect(u,v),

luego Enter y se obtiene el resultado deseado (22, 15, 40).
Se sugiere al estudiante, borrar la asignacion de nombres que se acaba de

hacer, pues estas quedan en la memoria y pueden causar alg
un problema al
querer usar los nombres para otra tarea. Para esto presione 2nd >F1, y elija
la primera opcion Borrar a-z.
19
c
2012
Mynor Chacon Daz
2.7
MA1003
Actividad 2
Sean u = (1, a, 1) y v = (a, 1, a) vectores en R3 . Use la TI-89 para calcular
ku vk.
c
2012
Mynor Chac
on Daz
20
CAPITULO
Operaciones con matrices
En este captulo se revisa el concepto de matriz, las operaciones elementales

con matrices, el determinante de una matriz, la definicion de matriz reducida
escalonada. Ademas el metodo de Gauss-Jordan para resolver sistemas de
ecuaciones lineales y el calculo de valores propios de una matriz.
3.1.
Representaci
on, adici
on y producto de matrices
Para representar la matriz A =
!
2 3 5
en la TI-89 se digita [[2,3,5][-7,0,8]]
7 0 8
y luego se presiona Enter
Para efectuar operaciones aritmeticas resulta comodo almacenar las matrices

y luego realizar los calculos. Por ejemplo si se desea calcular 3 A 5 B,
21
3.1
MA1003
!
1 0 3
siendo A la matriz mencionada anteriormente y B =
2 1 4
Primero se almacenan las dos matrices con sus respectivos nombres: a y b.
Luego se digita 3*a-5*b y se presiona Enter.
!
1
9 0
As se obtiene el resultado 3A 5B =
31 51 4
Sea
1 2
C = 1 5
3 0
Suponga que se desea realizar AC, pero no se quiere almacenar la matriz C.
Para efectuar esta operacion, digite a*[[1,2][-1,5][3,0]], presione Enter
y se obtiene el resultado deseado:
AC =
c
2012
Mynor Chac
on Daz
!
1 2
2 3 5
1 5 =
7 0 8
3 0
22
!
14
19
17 14
MA1003
3.2.
3.3
Transpuesta de una matriz
Si
A=
!
2 3 5
7 0 8
la Trasnpuesta de A es
2 7
AT = 3
0
5 8
Para realizar esta operacion con la TI-89, digite [[2,3,5][-7,0,8]],luego
siga la ruta
2nd >Math >Matrices >T
Y presione Enter.
23
c
2012
Mynor Chacon Daz
3.3
MA1003
3.3.
Determinante de una matriz
Suponga que se necesita calcular
5
1
A=
0
2
el determinante de la matriz
1 0 3
4 2 2
1 0 5
6 1 1
Para realizar este calculo siga la ruta

2nd >Math >Matrices >det(
En seguida ingrese la matriz respectiva, cierre el parentesis ), que se abrio automaticamente con el comando det( y finalice presionando Enter.
Se obtiene de esta forma el resultado det(A) = 194.
c
2012
Mynor Chac
on Daz
24
MA1003
3.3
Actividad 1
Determine los valores de x para los cuales

x + 4 1
1

det A = 1 x 5
1 =0

0
1
x + 3
Obtenga el determinante con la TI-89, luego para resolver la ecuacion use
el comando soluc, el cual puede ser accesado presionando la tecla F2 y la
opcion 1:soluc(.
El formato de dicho comando es soluc(ecuaci
on, variable)
Las soluciones deberan apreciarse como en las siguientes imagenes:
25
c
2012
Mynor Chacon Daz
3.4
3.4.
MA1003
M
etodo de Gauss-Jordan
Sistema con soluci

on u
nica.
El Metodo de Gauss-Jordan, se utiliza para resolver un sistema de ecuciones
lineales. Por ejemplo, para resolver el sistema
x + 2y + 3z = 9
4x + 5y + 6z = 24
3x + y 2z = 4
Ingrese a la calculadora TI-89 la matriz aumentada del sistema.
ahora siga la ruta

2nd >Math >Matrices >escalred(
Se marca la matriz
c
2012
Mynor Chac
on Daz
26
MA1003
3.4
Se presiona Enter, se cierra el parentesis redondo y se obtiene la matriz

Escalonada Reducida correspondiente.
De la matriz escalonada reducida
1 0 0 | 4
0 1 0 | 2
0 0 1 | 3
se obtiene se obtiene la solucion, x = 4, y = 2, z = 3
Sistema con infinitas soluciones
Considerese ahora el sistema
(
x+ y+ z =2
x 2y + 3z = 0
Al calcular la matriz escalonada reducida correspondiente a la matriz aumentada
27
c
2012
Mynor Chacon Daz
3.5
MA1003
1 1 1 |
1 2 3 |
!
2
0
Se obtiene
De donde
4
5
x+ z =
3
3
2
2
y z=
3
3
Si se denomina z = t entonces las infinitas soluciones se representan de
manera parametrica como sigue:
x = 45t
y = 2(1+t)
3
z=t
3.5.
tR
Uso del comando soluc() en el caso de infinitas

soluciones
Para ejemplificar el comando soluc, se va a resolver de nuevo el sistema

(
x+ y+ z =2
x 2y + 3z = 0
Para utilizar el comando soluc presione F2 y elija la primer opcion

c
2012
Mynor Chac
on Daz
28
MA1003
3.5
Ahora digite solucion(x+y+z=2 and x-2y+3z=0,{x,y,z}). Para digitar la

palabra and presione la tecla CATALOG
Luego de presionar Enter se obtienen los siguientes resultados
29
c
2012
Mynor Chacon Daz
3.5
MA1003
Si se denomina @1=t, se obtienen las mismas soluciones:
x = 45t
tR
y = 2(1+t)
3
z=t
Note que coinciden con las obtenidas con el metodo de Gauss-Jordan.
Actividad 2
Encuentre la recta intersecci
on de los planos
(
2x + 3y 2z = 4
x y + 2z = 3
Actividad 3
El objetivo en esta actividad es resolver sistemas de ecuaciones utilizando
MuPAD
1. Para resolver el sistema
(
x+y =5
x y = 1
Digte solve([x + y = 5, x y = 1], {x, y})

2. Ahora considere el sistema de la secci
on 3.4 , 3.5:
c
2012
Mynor Chac
on Daz
30
MA1003
3.6
x+y+z =2
x 2y + 3z = 0
intente resolverlo utilizando el mismo comando del inciso anterior y

compruebe que no funciona, es decir, el comando utilizado no resuelve
sistemas con soluciones infinitas.
3. Para resolver tal sistema digite en MuPAD:
equations:= {x + y + z = 2, x 2 y + 3 z = 0}
solve(equations)
3.6.
Valores Propios y Vectores Propios de una matriz
Recuerde que si A es una matriz cuadrada, 6= 0 es un valor propio

(eigenvalor) de A si v 6= 0, tal que A I = 0, en tal caso v es llamado
vector propio (eigenvector) de A asociado al valor propio .
En este curso solo interesa determinar los valores propios y no los vectores
propios. Por ejemplo ingrese la matriz
12 4 4
A = 4 12 4
4 4 12
en la TI-89, luego de seguir la ruta
2nd >Math >Matrices >valPr(
debera verse como en la siguiente imagen

31
c
2012
Mynor Chacon Daz
3.6
MA1003
De aqu se obtiene los tres valores propios: 16,4,16. Debido a que el valor
16 aparece dos veces y el 4 una vez, se dira que el valor propio = 16
tiene multiplicidad algebraica 2 y el valor propio = 4 tiene multiplicidad
algebraica 1.
Ahora se resolvera el mismo problema utilizando el software MuPAD Light,
para ello haga clic en el cono
Seguidamente digite A:=matrix([[12, 4, 4], [4, 12, 4], [4, 4, 12]]) para
definir la matriz A y luego digite linalg::eigenvectors(A)
c
2012
Mynor Chac
on Daz
32
MA1003
3.6
Los vectores columna que aparecen a la derecha de = 16, (1, 1, 0), (1, 0, 1)
son los vectores propios correspondientes, analogamente (1, 1, 1) es el vector
propio de = 4. Si se desean solamente los valores propios y no los vectores
propios se digita numeric::eigenvalues(A)
Actividad 4
Determine los valores propios de
2 1
A = 1
2
2
1
1
6
1. Utilizando la TI-89.
2. Utilizando la instruccion numeric::eigenvalues(A) en MuPAD.
3. Utilizando la instruccion linalg::eigenvectors(A)en MuPAD.
4. Noto alguna situacion extra
na?, a que se deber
a?
33
c
2012
Mynor Chacon Daz
3.6
c
2012
Mynor Chac
on Daz
MA1003
34
CAPITULO
Superficies en R3
Este captulo comprende el trazo de superficies simples generadas por curvas

planas, superficies cuadricas y otras mas generales. Para esto se hace el uso
de GenGCF y del Winplot.
4.1.
Representaci
on de superficies generadas por curvas planas, con GenGCF
Para iniciar el estudio de superficies en el espacio R3 . Considere una

curva del tipo y = f (x), el conjunto de todos los puntos del espacio del tipo
(x, f (x), z) forma una superficie cuya grafica se obtiene al tomar la curva
y = f (x) y moverla paralelamente al eje z. La ecuacion de tal superficie se
escribe y = f (x), z R. Analogamente, si la curva tiene el formato F (x, y) =
0, la superficie se escribe F (x, y) = 0, z R. Sin embargo en ambos casos es
costumbre no presentar la variacion de z, por lo que el lector debe deducir si
se refiere a la curva plana o a la curva espacial, esto normalmente depende
del contexto en el que se escriba la ecuacion.
Ejemplo con Winplot

La grafica de y = x3 generada por Winplot es la siguiente,
35
4.1
MA1003
Debido al nivel de dificultad, el procedimiento para dibujar la superficie

y = x3 , z R en Winplot se explicara en clase, la gràfica de la superficie es
la siguiente:
Ejemplo con GenGCF

La grafica de z 2 y 2 = 1 es la siguiente:
c
2012
Mynor Chac
on Daz
36
MA1003
4.2
Para obtener la grafica de la superficie z 2 y 2 = 1, x R en el GenGCF,

haga clic en Generador de archivos para Graphic Calculator cuyo cono es
. Luego se digita la ecuacion de manera que se vea como en la imagen
siguiente:
La superficie generada por el GenGCF es la siguiente:
4.2.
4.2.1.
Representaci
on de superficies cu
adricas, con GenGCF
Elipsoide:
x2 y 2 z 2
+ 2 + 2 =1
a2
b
c
Abra el Generador de archivos para Graphic Calculator cuyo cono es

x2 y 2
z2
. Para graficar
+
+
= 1, digite la ecuacion de manera que se
5
8
16
vea como en la imagen siguiente:
37
c
2012
Mynor Chacon Daz
4.2
MA1003
La grafica generada debe verse as:
Note que la interseccion del elipsoide con el plano z = 2 genera una elipse.
A esta elipse se le denomina Traza del elipsoide con el plano z = 2.
La traza con el plano y = 1 tambien genera una elipse, observe que la grafica
del plano tiene un opacidad del 0.5.
c
2012
Mynor Chac
on Daz
38
MA1003
4.2.2.
4.2
Hiperboloide de una hoja:
x2 y 2 z 2
+ 2 2 =1
a2
b
c
x2 y 2
z2
La grafica de
+
= 1, es
5
8
16
La traza con el plano x = 1 genera una hperbola.

39
c
2012
Mynor Chacon Daz
4.2
MA1003
La traza con el plano z = 1 genera una elipse.
4.2.3.
Hiperboloide de dos hojas:
x2
La grafica de
y 2 z 2 = 1, es
6
c
2012
Mynor Chac
on Daz
40
x2 y 2 z 2
2 2 =1
a2
b
c
MA1003
4.2
La traza con el plano z = 2 genera una hperbola.
La traza con el plano y = 1 genera una hiperbola.

41
c
2012
Mynor Chacon Daz
4.2
c
2012
Mynor Chac
on Daz
MA1003
42
MA1003
4.2.4.
4.2
Paraboloide Elptico: z =
x2 y 2
+ 2
a2
b
La grafica de 5z = 3x2 + 5y 2 , es
La traza con el plano x =
1
genera una parabola.
2
43
c
2012
Mynor Chacon Daz
4.2
MA1003
4.2.5.
x2 y 2
Paraboloide Hiperb
olico: z = 2 2
a
b
La grafica de 2z = x2 3y 2 , es
c
2012
Mynor Chac
on Daz
44
MA1003
4.2
1
La traza con el plano z = genera una hiperbola con eje de simetra paralelo
2
al eje X.
1
La traza con el plano z = genera una hiperbola con eje de simetra
2
paralelo al eje Y.
45
c
2012
Mynor Chacon Daz
4.2
MA1003
La traza con el plano x = 1/2 genera una parabola.
4.2.6.
x2 y 2
Cono Elptico: z = 2 + 2
a
b
2
La grafica de z 2 = x2 + y 2 , es
c
2012
Mynor Chac
on Daz
46
MA1003
4.2
La traza con el plano x = 1 genera una hiperbola.
La traza con el plano y = 2x genera dos rectas.
47
c
2012
Mynor Chacon Daz
4.2
MA1003
La traza con el plano z =
4.2.7.
1 1
+ x genera una elipse
2 2
Cilindros
Un cilindro es la superficie formada por una familia de rectas paralelas llamadas generatrices que pasan por los puntos de una curva espacial C llamada curva directriz.
Sea
P0 un punto sobre la curva C.
P un punto generico sobre una recta generatriz L.
v = (a, b, c) un vector director de cualquier generatriz.

Ecuaciones Parametricas de L
c
2012
Mynor Chac
on Daz
48
MA1003
4.2
x = x0 + a t
y = y0 + b t
z =z +ct
0
tR
Supongase que C esta determinada por la interseccion de dos superficies:

(
F (x, y, z) = 0
G(x, y, z) = 0
Para determinar la ecuacion H(x, y, z) = 0 del cilindro se debe eliminar las
variables x0 , y0 , z0 , t de las ecuaciones:
x = x0 + a t
y = y0 + b t
z = z0 + c t
F (x0 , y0 , z0 ) = 0
G(x0 , y0 , z0 ) = 0
Ejemplo
Hallar la ecuacion del cilindro elptico cuya directriz es la elipse de ecuaciones
x = cos , y = sen , z = sen + cos , y cuyas generatrices son perpendiculares al plano que contiene dicha elipse.
Soluci
on
Observe que cualquier punto de la curva satisface las ecuaciones z = x + y y
x2 + y 2 = 1, de manera que dicha curva es la interseccion de las superficies:
(
F (x, y, z) = x + y z = 0
G(x, y, z) = x2 + y 2 1 = 0
Puesto que las generatrices son perpendiculares al plano x + y z = 0,
entonces deben ser paralelas al vector (1, 1, 1). De de forma que ya es
posible construir el sistema que va a generar la ecuacion del cilindro, a
saber:
x = x0 + t
y = y0 + t
z = z0 t
x0 + y0 z0 = 0
2
x0 + y02 1 = 0
49
c
2012
Mynor Chacon Daz
4.2
MA1003
De las tres primeras ecuaciones se tiene x0 = x t, y0 = y t, z0 = z t,

x+yz
al sustitur en la ecuacion del plano se obtiene t =
, luego, de la
3
u
ltima ecuacion del sistema se obtiene:
(x t)2 + (y t)2 =

x+yz
x
3
2
x+yz
+ y
3
2
=
(2x y + z)2 + (2y x + z)2 =
5x2 + 5y 2 + 2z 2 8xy + 2xz + 2yz =
El desarrollo de los cuadrados puede hacerse utilizando MuPAD, para ello

haga clic en el cono
.
Seguidamente digite expand((2x y + z)2 + (2y x + z)2 )
Ahora utilizaremos GenGCF para ver las superficies z = x + y, x2 + y 2 = 1

que dan origen a la curva directriz
La superficie cilndrica generada junto con las dos superficies anteriores

se aprecian en las siguiente imagen:
c
2012
Mynor Chac
on Daz
50
MA1003
4.2
El grafico aislado del cilindro se aprecia as:
4.2.8.
Cilindro con Directriz en forma param

etrica
Suponga que la curva Directriz es:
x = x(s)
y = y(s)
z = z(s)
51
sR
c
2012
Mynor Chacon Daz
4.2
MA1003
Siendo v un vector director de una recta generatriz, la ecuacion parametrica

del cilindro es R(s, t) = r(s) + tv.
x = x(s) + t v1
R(s, t) =
y = y(s) + t v2
z = z(s) + t v
3
s, t R
Eliminando los parametros s y t, se obtiene la ecuacion cartesiana correspondiente.
Ejemplo
Resolver el ejemplo de la seccion anterior pero ahora considerando solamente
la representacion parametrica de la elipse directriz:
x = cos , y = sen , z = sen + cos
y cuyas generatrices son paralelas al vector director v = (1, 1, 1).
c
2012
Mynor Chac
on Daz
52
MA1003
4.2
Soluci
on
La ecuacion parametrica del cilindro es
R(, t) = r() + tv = (cos , sen , sen + cos ) + t (1, 1, 1)
el sistema que se debe construir queda de la siguiente manera
x = cos + t
s, t R
R(, t) :=
y = sen + t
z = sen + cos t
Despejando se obtiene x t = cos , y t = sen , esto implica que (x t)2 +
(y t)2 = 1, luego al sustituir sen y cos en z = sen + cos t se obtiene
z =yt+xttt=
x+yz
3
De nuevo
(x t)2 + (y t)2 =

x+yz
x
3
2
x+yz
+ y
3
2
=
(2x y + z)2 + (2y x + z)2 =
5x2 + 5y 2 + 2z 2 8xy + 2xz + 2yz =
53
c
2012
Mynor Chacon Daz
4.2
4.2.9.
MA1003
Conos
Un cono es una superficie formada por una familia de rectas no paralelas

llamadas generatrices, que pasan por un vertice(a, b, c) y por los puntos de
una curva C llamada curva directriz.
Ecuaciones Parametricas de L
x = (1 t) a + t x0
y = (1 t) b + t y0
z = (1 t) c + t z
0
tR

(
F (x, y, z) = 0
G(x, y, z) = 0
Para determinar la ecuacion H(x, y, z) = 0 del cono se debe eliminar las
variables x0 , y0 , z0 , t de las ecuaciones:
x = (1 t) a + t x0
y = (1 t) b + t y0
z = (1 t) c + t z0
F (x0 , y0 , z0 ) = 0
G(x0 , y0 , z0 ) = 0
Ejemplo
Halle la ecuacion de la superficie conica que tiene por vertice el punto (1, 2, 3)
y cuya directriz es la interseccion de las superficies x2 + y 2 = 4, z = 0
Soluci
on
Sea (x0 , y0 , z0 ) un punto de la directriz, esto significa que se debe cumplir
x20 + y02 = 4, z0 = 0. As el sistema a trabajar es el siguiente:
c
2012
Mynor Chac
on Daz
54
MA1003
4.2
x = (1 t) + t x0
y = 2(1 t) + t y0
z = 3(1 t) + t z0
x20 + y02 4 = 0
z0 = 0
Despejando se obtiene:
x1+t
x0 = t
y0 = y2+2t
t
z = z3+3t
0
t
Luego
x20
y02
(x 1 + t)2 (y 2 + 2t)2
=4
+
=4
t2
t2
Por otro lado, dado que z0 = 0, entonces z 3 + 3t = 0, lo cual implica

z
que t = 1 .
3
En resumen se tienen las siguientes ecuaciones:
(
(x 1 + t)2 + (y 2 + 2t)2 = 4t2
t = 1 z3
De manera que

x1+1

z 2
z 2
+ y2+2 1
=
3
3
(3x z)2 + (3y 2z)2 =

z 2
4 1
3
4 (3 z)2
Luego de desarrollar y simplificar se obtiene la ecuacion del cono:

9x2 + 9y 2 + z 2 6xz 12yz + 24z 36 = 0
El gràfico generado por tal GenGCF correspondiente a tal ecuacion es:
55
c
2012
Mynor Chacon Daz
4.2
MA1003
Ejemplo
Determine la ecuacion del cono que aparece en la imagen siguiente.
Soluci
on
Puesto que el vertice es (1, 1, 2) y las superficies que determinan la curva
son z = x2 + y 2 y x + y + z = 2, entonces el sistema correspondiente es el
siguiente:
x = (1 t) + tx0
y = (1 t) + ty0
z = (1 t)(2) + tz0
x20 + y02 = z0
x0 + y0 + z0 = 2
de aqu x0 =
x1+t
y1+t
z 2 + 2t
, y0 =
, z0 =
.
t
t
t
Por otro lado
x0 + y0 + z0 = 2
Sustituyendo en esta ecuacion, se obtiene:
c
2012
Mynor Chac
on Daz
56
MA1003
4.2
x 1 + t y 1 + t z 2 + 2t
+
+
=2
t
t
t
x + y + z 4 + 4t = 2t
x + y + z 4 = 2t
t=
x y z + 4
2
Analogamente al sustituir en
x20 + y02 = z0
se obtiene:
z 2 + 2t
=
t
x1+t
t
2

+
y1+t
t
2
(z 2 + 2t)t = (x 1 + t)2 + (y 1 + t)2

(z 2 x y z + 4)t =
(2 x y)t =
x y z + 4
x1+
2
2
x y z + 4
+ y1+
2
(2x 2 x y z + 4)2 (2y 2 x y z + 4)2

+
4
4
4(2 x y)t = (x y z + 2)2 + (y x z + 2)2

4(2 x y)
x y z + 4
2
= (x y z + 2)2 + (y x z + 2)2
2(2 x y)(4 x y z) = (x y z + 2)2 + (y x z + 2)2
57
c
2012
Mynor Chacon Daz
2
4.2
4.2.10.
MA1003
Superficie de Revoluci
on
Es una superficie formada por una familia de crculos generatrices, centrados

en los puntos de un eje E, que pasan por los puntos de una curva directriz,
C.
Note que cualquier crculo perpendicular al eje E puede ser considerado como
la interseccion de un plano perpendicular a E, v1 x0 +v2 y0 +v3 z0 = t con una
esfera centrada en el punto (a, b, c) de E, (x0 a)2 +(y0 b)2 +(z0 c)2 = R2 .

(
F (x, y, z) = 0
G(x, y, z) = 0
Si (x0 , y0 , z0 ) es un punto en la superficie, debe cumplirse el sistema
(x0 a)2 + (y0 b)2 + (z0 c)2 = R2
v1 x0 + v2 y0 + v3 z0 = t
F (x0 , y0 , z0 ) = 0
G(x0 , y0 , z0 ) = 0
c
2012
Mynor Chac
on Daz
58
MA1003
4.2
De estas ecuaciones se eliminan (x0 , y0 , z0 ) y se obtiene la ecuacion h(R2 , t) =

0, as la ecuacion de la superfice de revolucion es
h((x a)2 + (y b)2 + (z c)2 , v1 x + v2 y + v3 z) = 0
Ejemplo
Determine la ecuacion de la superficie de revolucion que se obtiene al girar
la recta
x z = 1, x y + z = 0
alrededor del eje x = y = z.
Soluci
on
y0
z0
x0
=
=
es decir es la
El eje x = y = z representa la recta
1
1
1
recta que pasa por (0, 0, 0) y tiene vector director v = (1, 1, 1), de manera
que las ecuaciones de los crculos generatrices son
(
(x0 0)2 + (y0 0)2 + (z0 0)2 = R2

1 x0 + 1 y0 + 1 z0 = t
Las superficies F (x, y, z), G(x, y, z) son respectivamente xz = 1, xy+z =

0, de manera que el sistema que se debe analizar es el siguiente:
x20 + y02 + z02 = R2
x +y +z = t
0
0
0
1
=0
0
0
x0 y0 + z0 = 0
(1)
(2)
(3)
(4)
De la ecuacion (3) se obtiene z0 = x0 1 y al sustituirla en (4) se obtiene la

nueva ecuacion x0 y0 + x0 1 = 0, lo cual implica que y0 = 2x0 1.
Hasta aqu se tiene z0 = x0 1 y y0 = 2x0 1, al reemplazar en las ecuaciones
(1) y (2) del sistema se obtiene
(
x20 + (2x0 1)2 + (x0 1)2 = R2

x0 + (2x0 1) + (x0 1) = t
Desarrollando y simplificando se obtiene:

(
6x20 6x0 + 2 = R2
4x0 2 = t
59
c
2012
Mynor Chacon Daz
4.2
MA1003
t+2
De este sistema se obtiene x0 =
y al sustitur en la primera de estas
4
ecuaciones se llega a la identidad:
3 2 1
t + = R2
8
2
Una vez que se obtiene una ecuacion en terminos de t y R y dado que
todos los puntos (x, y, z) de la superficie de revolucion deben satisfacer las
ecuaciones (1) y (2) del sistema original:
x + y + z = t, x2 + y 2 + z 2 = R2
se sustituyen t y R2 para llegar a la ecuacion de la superficie de revolucion:
3
1
(x + y + z)2 + = x2 + y 2 + z 2
8
2
Desarrollando y simplificando se obtiene la ecuacion final:
5x2 + 5y 2 + 5z 2 6xy 6xz 6yz 4 = 0
La superficie generada por esta ecuacion es la siguiente:
La recta de color rojo es el eje de revolucion y la de color morada es la recta

directriz.
La superficie sin la generatriz ni el eje es la siguiente:
c
2012
Mynor Chac
on Daz
60
MA1003
4.2.11.
4.2
Un caso particular
Si C := f (y, z) = 0 se hace girar sobre el eje Y , la Superficie de Revolu

cion tiene como ecuacion f (y, x2 + z 2 ) = 0.
En resumen si:
f (x, y)
p = 0 gira alrededor del eje X, la ecuacion de la superficie es
f (x, y 2 + z 2 ) = 0.
f (x, y) = 0 gira alrededor del eje Y , la ecuacion de la superficie es
f ( x2 + z 2 , y) = 0.
f (y, z) = 0 gira alrededor del eje Y , la ecuacion de la superficie es
f (y, x2 + z 2 ) = 0.
f (y,
pz) = 0 gira alrededor del eje Z, la ecuacion de la superficie es
f ( x2 + y 2 , z) = 0.
f (x, z)
p = 0 gira alrededor del eje X, la ecuacion de la superficie es
f (x, y 2 + z 2 ) = 0.
f (x,
pz) = 0 gira alrededor del eje Z, la ecuacion de la superficie es
f ( x2 + y 2 , z) = 0.
61
c
2012
Mynor Chacon Daz
4.2
MA1003
Ejemplo
2
Si z = ex se hace girar alrededor del eje Z, en este caso f (x, z) =

2
z ex = 0 y la ecuacion de la superficie de revolucion es
2
p
2
2
x2 +y 2
2
2
= 0 z = e(x +y )
f ( x + y , z) = 0 z e
2
El grafico de z = ex generado por Winplot es:
El grafico de z = e(x
2 +y 2
) generado por Winplot es:
El grafico de z = e(x
2 +y 2
) generado por GenGCF es:
c
2012
Mynor Chac
on Daz
62
CAPITULO
Funciones Vectoriales
En este captulo se hace un estudio breve de las funciones vectoriales. Dada la

parametrizacionde una fucion f : Rn R, se utiliza Winplot y el GenGCF
para el trazado de la curva respectiva.
5.1.
Definici
on y propiedades
Una funcion f : R Rn se llama funcion vectorial, se denota

f (t) = (x1 (t), . . . , xn (t))
cada componente xi : R R, i = 1 . . . n. En el caso particular n = 3, se
acostumbra a escribir f (t) = (x(t), y(t), z(t)).
5.2.
Derivaci
on e integraci
on
Dada f (t) = (x1 (t), . . . , xn (t)), se define

lm f (t) = lm x1 (t), . . . , lm xn (t)
t7a
t7a
t7a
Se deduce entonces que

f 0 (t) = (x01 (t), . . . , x0n (t))
siempre que cada xi sea diferenciable, analogamente
Z
Z
f (t) dt =
Z
x1 (t) dt, . . . ,
63

xn (t) dt
5.2
MA1003
siempre que cada xi sea integrable.

Seguidamente se muestra como derivar e integrar una funcion vectorial utilizando el emulador de la TI-89, como ejemplo se va a usar la funcion

f (t) =
1
t , sen t,
1 + t2
3
Para definir la funcion siga la ruta F4 >Definir
Para ingresar la funcion en la TI-89,se digita

Definir f (t) = [[t 3, sen t, 1/(1 + t 2)]]
luego se digita ENTER

c
2012
Mynor Chac
on Daz
64
MA1003
5.2
.
De esta manera es posible utilizar la funcion para poder evaluarla, derivarla
o integrarla sin tener que digitarla de nuevo, para invocarla digite f (t).
Por ejemplo, para calcular f (0) simplemente se digita f(0) y se obtiene el
resultado (0, 0, 1), es decir f (0) = (0, 0, 1).
.
Ahora para determinar f 0 (x), siga la ruta F3 >d(derivar
65
c
2012
Mynor Chacon Daz
5.2
MA1003
.
Digite d(f(t),t) y luego ENTER de esta manera se obtiene:
0
f (t) =
2t
3t , cos t,
1 + t2
2
.
R
Analogamente para determinar f (x)dx, siga la ruta
R
F3 > (integrar
de manera que se vea como sigue
obtiene:
Z
f (t) dt =
c
2012
Mynor Chac
on Daz
(f (t), t) y luego se presiona ENTER as se

t4
, cos t, arctan t
4
66
MA1003
5.2
f (t) dt. En el pre-
La instruccion (f(t), t, a, b) da como resultado

a
sente ejemplo se tiene

Z

f (t) dt =
, 1 cos(1),
4
4
.
Seguidamente se exponen algunas propiedades de las derivacion del producto
por un escalar, del producto punto y del producto cruz de funciones vectoriales.
Propiedades
Sean r(t), u(t) dos funciones vectoriales derivables, R y f (t) una funcion
real derivable, entonces:
1.
d
[ r(t)] = r0 (t)
dt
2.
d
[r(t) u(t)] = r0 (t) u0 (t)
dt
3.
d
[f (t) u(t)] = f 0 (t) u(t) + f (t) u0 (t)
dt
4.
d
[r(t) u(t)] = r0 (t) u(t) + r(t) u0 (t)
dt
5.
d
[r(t) u(t)] = r0 (t) u(t) + r(t) u0 (t)
dt
67
c
2012
Mynor Chacon Daz
5.3
5.3.
MA1003
Curvas en R2 , R3
Al grafico de una funcion r : R Rn , se le denomina curva C, y se escribe

C := r(t), si r : R R2 , su grafica es una curva plana y si r : R R3 , su
grafica es una curva en el espacio. El vector r(t) es el vector posicion.
Si un objeto se mueve a lo largo de una curva C := r(t), entonces r(t0 )

indica la posicion del objeto en en instante t0 .
Considere el siguiente grafico de la curva C, los vectores r(t), r(t + h) y
r(t + h) r(t).
Geometricamente es claro que si h tiende a cero el vector r(t + h) r(t)

tiende a ser tangente a la curva en el punto r(t) y sigue la direccion de la
curva, esto puede expresarse como sigue:
lm
h 70
r(t + h) r(t)
= r0 (t)
h
donde r0 (t) un vector tangente a la curva con la direccion de la misma.

c
2012
Mynor Chac
on Daz
68
MA1003
5.4
En la siguiente imagen se aprecian vectores tangentes en distintos puntos
En general si existe r0 (t) y no es nula, la recta que pasa por r(t) y es paralela a r0 (t) se llama Recta Tangente a C en r(t). Al vector v(t) = r0 (t)
se le conoce como vector velocidad y al vector a(t) = r00 (t) vector aceleracion.
5.4.
Parametrizaci
on
Considere la curva C definida por la funcion vectorial r(t) = (x(t), y(t), z(t))
el sistema de ecuaciones siguiente:
x = x(t)
y = y(t)
z = z(t)
69
c
2012
Mynor Chacon Daz
5.4
MA1003
representa la parametrizacion de C y t es el parametro.
5.4.1.
Parametrizaci
on de y = f (x)
Una parametrizacion inmediata de la curva determinada por el grafico de

la funcion y = f (x), desde (a, f (a)) hasta (b, f (b)), se obtiene al utilizar la
variable independiente x como parametro t, es decir la parametrizacion es
la siguiente
x=t
axb
y = f (t)
Es decir, expresado de otra forma, la funcion real univariable y = f (x) y
la funcion vectorial r(t) = (t, f (t)) representan la misma curva. Se advierte
que existen otras posibles parametrizaciones de la misma curva, la opcion
anterior es simplemente la mas simple.
Ejemplo
La funcion lineal y = 3x + 2, 1 x 1, se parametriza como
x=t
1 t 1
y = 3t + 2
Sin embargo, desde el punto de vista del algebra lineal, el segmento de recta
que va del punto (1, 1) al punto (1, 5), tiene vector director (2, 6) y pasa
por el punto (1, 1), por lo que las ecuaciones parametricas de tal segmento
de recta son:
x = 1 + 2t
0t1
y = 1 + 6t
En resumen se obtuvieron dos parametrizaciones distintas para el mismo
segmento de recta.
c
2012
Mynor Chac
on Daz
70
MA1003
5.4.2.
5.4
Parametrizaci
on de F (x, y) = 0
Por ejemplo, considere la elipse con ecuacion x2 + xy + y 2 1 = 0, para

obtener una parametrizacion complete cuadrados utilizando los dos primeros
teminos, y se procede como sigue:
1 = x2 + xy + y 2
y2 y2
= x2 + xy +
+ y2
4
4

y 2 3 2
= x+
+ y
2
4
!2

y 2
3
= x+
y
+
2
2
!2
y 2
3
La ecuacion x +
+
y
= 1, sugiere escribir
2
2
x + 2 = cos t
0 t 2
3 y = sen t
2

x = cos t
y=
2
3
y
2
0 t 2
sen t
x = cos t
y=
2
3
1
3
sen t
0 t 2
sen t
Note que la parametrizacion anterior genera la grafica de la funcion vectorial

r(t) =

1
2
cos t sen t, sen t
3
3
en particular

1 2
r(0) = (1, 0),
r
= ,
2
3 3
Al observar estos dos puntos en la curva, se deduce que la trayectoria de la
curva va en direccion contraria a la de las agujas del reloj. Esta

se conoce
como la direccion positiva

71
c
2012
Mynor Chacon Daz
5.4
5.4.3.
MA1003
Parametrizaci
on de curvas en R3
A continuacion se consideran varios ejemplos que cubren lo necesario y basico

de este tema.
Ejemplo
Recuerde que la ecuacion x2 + y 2 = 9 representa un crculo de radio 3, en el
plano XY, en coordenadas polares las ecuaciones correspondientes son
x = 3 cos t
0 t 2
y = 3 sen t
que pasa si se a
nade z = 2 a las ecuaciones anteriores?
La imaginacion sugiere que la parametrizacion
x = 3 cos t
y = 3 sen t
z=2
0 t 2
representa la misma curva elevada a una altura de 2 unidades en el eje Z,

como se aprecia en la siguiente imagen:
c
2012
Mynor Chac
on Daz
72
MA1003
5.4
La funcion vectorial correspondiente es r(t) = (3 cos t, 3 sen t, 2).

Ejemplo
Cual es la parametrizacion de la elipse que se observa en el grafico siguiente?
La curva se encuentra situada en el plano y = 1, luego la elipse tiene semieje

mayor de longitud 4(paralelo al eje Z), semieje menor de longitud 2 (paralelo
al eje X). Esto significa que la ecuacion de la elipse es
x2 z 2
+
=1
4
16
x 2
z 2
=1
2
4
Esto sugiere la parametrizacion
x = 2 cos t
y=1
0 t 2
z = 4 sen t
73
c
2012
Mynor Chacon Daz
5.4
MA1003
Como funcion vectorial se escribe r(t) = (2 cos t, 1, 4 sen t).

Ejemplo
Cual es la parametrizacion de la curva interseccion de la superfice esferica
x2 + y 2 + z 2 = 25 con el plano 2x + y z = 0?
Se necesita entonces parametrizar la curva

(
x2 + y 2 + z 2 = 25
C:
2x + y z = 0
Despejando z de la segunda ecuacion y reemplazando en la primera se obtiene:
25 = x2 + y 2 + (2x + y)2
= 5x2 + 4xy + 2y 2
4
2
5 = x2 + xy + y 2
5
5
4
4
2
4
= x2 + xy + y 2 y 2 + y 2
5
25
25
5

2
2
6
=
x + y + y2
5
25
!2

2
2
6
=
x+ y +
y
5
5
c
2012
Mynor Chac
on Daz
74
MA1003
5.5

la ecuacion
2
x+ y
5
2
6
y
5
!2
= 5, sugiere:
x + 52 y = 5 cos t
6
y = 5 sen t
5
x=
y=
5 cos t 25 y
sen t
5 5
6
x = 5 cos t
y=
5 5
6
2 5
6
sen t
0 t 2
sen t
Puesto que z = 2x + y, se tiene:
z = 2x + y
!
2 5
5 5
5 cos t sen t + sen t
= 2
6
6
r
5
sen t
= 2 5 cos t +
6
De manera que la parametrizacion de la curva C es:
x
=
5
cos
t
2
sen t
q
y = 5 56 sen t
z = 25 cos t + 5 sen t
0 t 2
5.5.
Dibujo de curvas con Winplot y GenGCF
Para trazar una curva en R3 , utilizando Winplot, siga la ruta

Ventana >3-dim
En seguida desde la barra de menu de Winplot siga el camino
Ecua >Curva
75
c
2012
Mynor Chacon Daz
5.5
MA1003
En seguida aparece la ventana siguiente, en la cual , aparece la curva en

forma parametrica correspondiente a
r(t) = (cos 7t, sen 7t, t) , 0 t 2
c
2012
Mynor Chac
on Daz
76
MA1003
5.5
La curva aparece sin ejes, para agregarlos siga la ruta Ver >Ejes >Ejes
De esta manera se obtiene su grafica:
Para trazar una curva en R3 , utilizando GenGCF, desde la barra de menu,

siga la ruta
Dib 3D >Curva param
etrica
77
c
2012
Mynor Chacon Daz
5.6
MA1003
Por default aparece la curva en forma parametrica correspondiente a

r(t) = (2 sen t, 3 cos t, t)
Cambiando el intervalo de variacion de t, desde 0 hasta 8 25,132, se
obtiene
Cuyo grafico se aprecia en la imagen siguiente:
5.6.
El Tri
edro Intrnseco o Marco de Frenet
Suponga que r(t) = (x(t), y(t), z(t)) representa la trayectoria de un objeto

en el espacio tridimensional. El sistema de ejes formado por los vectores
unitarios i, j, k sirve para determinar la posicion que ocupa el objeto en el
instante t0 , esto es r(t0 ) = (x(t0 ), y(t0 ), z(t0 )), sin embargo en cada punto
de la curva existe otro sistema de ejes ortogonales que brinda informacion
valiosa. Este sistema es llamado el Marco de Frenet o Triedro Intrnseco,
esta constituido por los vectores:
T (t) =
r0 (t)
T 0 (t)
,
N
(t)
=
, B(t) = T (t) N (t)
kr0 (t)k
kT 0 (t)k
Concretamente, el vector Tangente unitario T (t) indica el movimiento hacia

adelante del objeto en el instante t, el vector Normal unitario N (t) indica la
c
2012
Mynor Chac
on Daz
78
MA1003
5.6
direccion en que se dobla la curva, el vector Binormal unitario B(t) indica

la tendencia de su movimiento a girarfuera del plano creado por T (t) y
N (t) en la direccion ortogonal a dicho plano. El plano formado por T (t)
y N (t) es llamado Plano Osculador, el plano formado por T (t) y B(t) es
llamado Plano Rectificante, el plano formado por B(t) y N (t) es llamado
Plano Normal.
La imagen anterior corresponde al Marco de Frenet,

calculado en el punto

correspondiente a t = , para la curva r(t) =

2 cos t, 2 sen t, 2 sen t .
4
Winplot brinda un opcion en donde se puede apreciar el Marco de Frenet
de forma automatica, sin tener que hacer calculos manuales. Para observar
esto, una vez dibujada la curva r(t) anterior, desde la barra de men
u siga la
ruta Una >Traza curvas...
79
c
2012
Mynor Chacon Daz
5.6
MA1003
y marque las
En la proxima ventana digite la curva, el valor de t =
4
ventanas correspondientes a Frenet y ancho de vectores, como se aprecia
en la imagen siguiente:
As se obtiene el Marco de Frenet en t =

c
2012
Mynor Chac
on Daz
80
:
4
MA1003
5.7
Al mover el cursor en la barra de desplazamiento, el triedro se desplaza sobre

la curva:
La siguiente es una imagen del Marco de Frenet en un punto t 6=
81
:
4
c
2012
Mynor Chacon Daz
5.7
5.7.
MA1003
Curvatura, Torsi
on y Longitud de arco
Dada una curva r(t) = (x(t), y(t)), recuerde que, T (t) indica el movimiento
hacia adelante del objeto en el instante t, su longitud es constante, solo
cambia su direccion cuando el objeto se mueve a lo largo de la curva. Se
llama curvatura a la razon con que T (t) girapor unidad de longitud, a lo
largo de la curva, por esta razon se define

dT
curvatura = =
ds

dT
En el caso de una curva tridimensional = , puede ser considerada
ds
como la razon(por unidad de longitud) con que el plano normal gira cuando
el objeto se mueve a lo largo de la curva.
Analogamente, para una curva tridimensional, B(t) indica la tendencia de su
movimiento a girarfuera del plano creado por T (t) y N (t)(plano osculador)
en la direccion ortogonal a dicho plano, se define torsion como la razon(por
unidad de longitud) con que el plano osculador gira alrededor de T (t), trata
de medir cuanto se tuerce la curva, se definie como
dB
N
ds
Por otro lado, retomando la definicion del vector unitario
torsion = =
T (t) =
v(t)
r0 (t)
=
0
kr (t)k
kv(t)k
Se puede escribir v = kvk T , ahora al derivar a ambos lados se obtiene:
a = kvk0 T + kvk T 0
= kvk0 T + kvk T 0 1
kT 0 k
= kvk0 T + kvk T 0
kT 0 k
0
T
= kvk0 T + kvk
kT 0 k
kT 0 k
= kvk0 T + kvk N kT 0 k
= kvk0 T + kvk kT 0 k N
|{z}
| {z }
aT
aN
Es decir a = aT T + aN N , a los valores aT , aN se les llama componente

tangencial y componente normal de la aceleracion. En terminos algebraicos
c
2012
Mynor Chac
on Daz
82
MA1003
5.7
esta ecuacion dice que el vector aceleracion a(t) = r00 (t) es combinacion
lineal de los vectores T (t) y N (t), es decir a(t) esta contenido en el plano
osculador.
Seguidamente se dan otras formulas mas comodas y u
tiles para obtener:
T, N, B, , , aT y aN .
Se sugiere utilizar las siguientes opciones:
r0 (t)
kr0 (t)k
r0 (t) r00 (t)
B(t) =
kr0 (t) r00 (t)k
(r0 (t) r00 (t)) r0 (t)
N (t) =
k(r0 (t) r00 (t)) r0 (t)k
T (t) =
Para calcular la curvatura, se sugiere:
|y 00 |
para C dada por y=f(x)
[1 + (y 0 )2 ] 2
|x0 y 00 y 0 x00 |
kr0 (t)k3
kr0 (t) r00 (t)k

kr0 (t)k3
para C dada por r(t)=(x(t),y(t))
para C dada por r(t)=(x(t),y(t),z(t))
Para calcular la torsion, se sugiere:

=
k[r0 (t) r00 (t)] r000 (t)k

kr0 (t) r00 (t)k2
Para obtener las componentes tangencial y normal de la aceleracion, se sugiere:

va
kvk
kv ak
=
kvk
aT =
aN
O tambien se puede utilizar la ecuacion

kak2 = a2T + a2N
83
c
2012
Mynor Chacon Daz
5.7
MA1003
Ejemplo
En el caso de la helice circular, r(t) = (2 cos t, 2 sen t, 3t), es facil verificar
2
0,153846153. Winplot
que, en cualquier punto de la curva se tiene =
13
da una opcion para obtener el valor de la curvatura directamente; para esto
primero dibuje la curva, luego siga la ruta
Una >Traza curvas...
En la ventana de nombre traza se obtiene el valor de la curvatura.
curvatura: 0.15385
Note que el valor obtenido de Winplot y el obtenido del calculo manual

coinciden.(Aplicando redondeo a cuatro cifras decimales)
c
2012
Mynor Chac
on Daz
84
MA1003
5.7
Por otro lado, recuerde que la Longitud de Arco de la curva r(t), t [a, b],
se obtiene con la formula
Z
L=
kr0 (t)k dt
De nuevo, Winplot brinda la opcion para determinar la longitud de arco de

manera automatica. Para esto siga la ruta
Una >Longitud de arco...
Para seguir con el ejemplo de la helice circular; en la ventana longitud de
la curva y en la casilla correspondiente a l
mite inferior, digite 0. En la
casilla correspondiente a l
mite superior, digite pi. Por u
ltimo haga clic
en longitud=. El valor 11.32717 representa la longitud de la curva desde
t = 0 hasta t = .
Para verificar este valor note que r0 (t) = (2 sin t, 2 cos t, 3), as
kr0 (t)k =
(4 sin2 t + 4 cos2 t + 9) =
13
por lo tanto
Z
L=
13 = 13 11,3271734
Observe que el valor obtenido con Winplot es una excelente aproximacion a

la longitud de arco obtenida del calculo directo.
85
c
2012
Mynor Chacon Daz
5.7
MA1003
Ejemplo
Considere la figura siguiente:
1. Sea la parte de la curva C que pasa por el punto P (0, 2, 2 2).

Determine la parametrizacion de .
2. Determine el vector Binormal B(t) en el punto P .
3. Determine la ecuacion cartesiana del plano osculador en el punto P .
h i
4. Determine la longitud de arco de , en 0, .
2
5. Determine la curvatura de en t = .
2
Soluci
on
1. La curva a parametrizar es
(
x2 + y 2 = 4
z 2 = 4y
De x2 + y 2 = 4, se deduce x = 2 cos t, y = 2 sen t, por lo que z 2 =
8 sen t z = 2 2 sen t. La parte de la curva que pasa por el punto
P (0, 2, 2 2) es la que tiene como parametrizacion
x = 2 cos t
:=
y = 2 sen t
z = 22 sen t
c
2012
Mynor Chac
on Daz
86
MA1003
2.
5.7

r(t) = 2 cos t, 2 sen t, 2 2 sen t r( ) = 0, 2, 2 2

2

2 cos t
2 sen t, 2 cos t,
r0 ( ) = (2, 0, 0)
2
sen t
r0 (t) =
r00 (t) =
2 cos t, 2 sen t,
2 (sen2 t + 1)
3
2 sen 2 t

r00 ( ) = 0, 2, 2
2
Luego

r0 ( ) r00 ( ) = (2, 0, 0) 0, 2, 2 = 0, 2 2, 4
2
2
De manera que

2
r0 ( 2 ) r00 ( 2 )

=
= 0
r ( 2 ) r00 ( 2 )
!
3 6
0,
,
3 3
3. La ecuacion del Plano Osculador es

0x + 2 2y + 4z = 0, 2 2, 4 0, 2, 2 2
2 2y + 4z = 4 2
4. Es difcil realizar los calculos correpondientes a la longitud de arco, por
lo que se utiliza Winplot
87
c
2012
Mynor Chacon Daz
5.7
MA1003
5. Analogamente la curvatura es:
c
2012
Mynor Chac
on Daz
88
CAPITULO
Funciones de Varias Variables o Campos Escalares
6.1.
Gr
aficaci
on de funciones de dos variables
Una funci
on de varias variables o Campo Escalar es una funcion que
tiene por dominio un subconjunto D R n , n 2, esto es f : D R n R .
Por ejemplo, una funcion de dos variables es
f (x, y) =
sen(x2 + y 2 )
x2 + y 2
Se puede obtener su grafica utilizando Winplot y GenGCF. En Winplot,

haga clic en Ventana y luego en 3-dim
Luego siga la ruta Ecua >Expl

cita..., en la ventana emergente z=f(x,y)
se digita la funcion y se hace clic en ok
89
6.1
MA1003
(Note que la funcion f (x, y) =
sen(x2 + y 2 )
es la que aparece de oficio
x2 + y 2
al abrir la ventana)
La grafica de la funcion es la siguiente:
Para generar la grafica anterior en GenGCF, luego de abrir el software digite

la funcion tal como se aprecia en la imagen siguiente:
c
2012
Mynor Chac
on Daz
90
MA1003
6.2.
6.2
Dominio de una funci

on de dos variables
El dominio maximo es el conjunto de puntos del plano en donde la funcion

esta bien definida, a saber Df = {(x, y)tqf (x, y) R }
Ejemplo
Determine y grafique el dominio maximo de f : R 2 R definida por
1
f (x, y) = p
2
2
(x + y 1) (4 x2 y 2 )
Soluci
on
Determinar el dominio equivale en este caso a resolver la desigualdad

x2 + y 2 1 4 x2 y 2 > 0
Esto se cumple cuando ambos factores son positivos o ambos son negativos.
En el segundo caso esto significa x2 + y 2 < 1 (interior del disco de radio 1)
y 4 < x2 + y 2 < 1 (exterior del disco de radio 2), si embargo no existe un
punto que satisfaga ambas desigualdades a al vez.
91
c
2012
Mynor Chacon Daz
6.2
MA1003
En el primer caso se debe cumplir x2 + y 2 > 1 (exterior del disco de radio

1) y 4 > x2 + y 2 < 1 (interior del disco de radio 2), por lo que los puntos
que satisfagan ambas desigualdades son aquellos que se encuentren fuera del
disco de radio 1 e interior al disco de radio 2.
Al escribir la desigualdad en el GenGCF se obtiene directamente la region
correspondiente al dominio maximo de la funcion:
El grafico de la funcion generado por GenGCF es el siguiente:
c
2012
Mynor Chac
on Daz
92
MA1003
6.3.
6.3
Lmite de una funci

on de dos variables
Para analizar el lm f (x), solamente existen dos opciones para acercarnos al

xa
valor a. El acercamiento puede hacerse por la derecha o por la izquierda.
De hecho se tiene que

lm f (x) = L lm+ f (x) = L = lm f (x)
xa
xa
xa
En el caso de una funcion de dos variables, para estudiar el

lm
f (x, y)
(x,y)(a,b)
existen infinidad de caminos para acercarnos al punto (a, b). En la siguiente

imagen se muestra alguans caminos para acercarse a punto (a, b), cada trazo
con distinto color es una camino posible por el cual (x, y) puede acercarse a
(a, b). Obviamente existen infinitas opciones.
93
c
2012
Mynor Chacon Daz
6.3
MA1003
Este hecho hace muy difcil el calculo de lmites de funciones de dos variables,
ya que el lmite existe si y solo si tiene el mismo valor por cualquier camino
que pase por (a, b). A veces se puede probar que
lm
f (x, y) = L
(x,y)(a,b)
pero esto requiere de la definicion formal de lmite la cual se omite en este

curso. a menos que f(a,b)=c, entonces es claro que:
lm
f (x, y) = f (a, b) = c
(x,y)(a,b)
Si esto sucede se dice que f (x, y) es continua en (a, b), y si f (x, y) es continua
en todo punto de una region D R 2 , entonces f (x, y) es continua en D.
Ejemplo
Si f (x, y) =
x2 + y 2
5
, es claro que f (1, 2) = , de manera que
x+y
3
x2 + y 2
5
=
(x,y)(1,2) x + y
3
lm
5
De manera que f (x, y) es continua en , de hecho los puntos de discon3
tinuidad de la funcion se encuentran sobre la recta y = x, en la sguiente
imagen se aprecia el grafico de la fucnion y la lnea que contiene los puntos
de discontinuidad.
c
2012
Mynor Chac
on Daz
94
MA1003
6.3
Ejemplo
Algunas tecnicas conocidas tambien pueden utilizarse en el calculo de lmites,
como se puede apreciar en el siguiente lmite:
x4 y 4
(x2 y 2 ) (x2 + y 2 )
=
l
m
= lm x2 y 2 = 0
(x,y)(0,0) x2 + y 2
(x,y)(0,0)
(x,y)(0,0)
x2 + y 2
lm
x4 y 4
es discontinua en (0, 0), sin embargo, se puede
x2 + y 2
redefinir de la siguiente manera para que sea continua:
En este caso f (x, y) =
4
4
x y , si x 6= 0
f (x, y) = x2 + y 2
0,
si x = 0
95
c
2012
Mynor Chacon Daz
6.3
MA1003
Ejemplo
En algunos casos, para verificar que el lmite no existe, se calcula el lmite
por distintos caminos que se acerquen al punto (a, b), un ejemplo de esto lo
brinda el
x2 y
(x,y)(0,0) x4 + y 2
lm
Al acercanos al (0, 0) por cualquier recta que pase por el origen, es decir con
puntos (x, y) que se encuentran sobre la recta y = mx, m R , se tiene que
si x 0 entonces obviamente y 0, por lo que:
mx3
mx
x2 y
=
l
m
= lm 2
=0
4
2
4
2
2
x0 x + m x
x0 x + m2
(x,y)(0,0) x + y
lm
Analogamente al acercanos al (0, 0) con puntos (x, y) que se encuentran

sobre la parabola y = x2 , se tiene que:
x2 y
x4
1
1
=
l
m
= lm =
4
2
4
4
x0 x + x
x0 2
(x,y)(0,0) x + y
2
lm
Puesto que el lmite tiene resultados 0 y
1
, se concluye que
2
x2 y
(x,y)(0,0) x4 + y 2
lm
Otra forma de llegar a la misma conclucion, es considerar las curvas y =

kx2 , k R , sobre estos caminos se tiene:
x2 y
kx4
k
=
l
m
=
x0 x4 + k 2 x4
(x,y)(0,0) x4 + y 2
1 + k2
lm
Es decir sobre cada camino y = kx2 , el valor del lmite cambia de acuerdo
al valor de k, por lo que el lmite no existe.
c
2012
Mynor Chac
on Daz
96
MA1003
6.4
x2 y
es discontinua en (0, 0). En la
x4 + y 2
siguiente imagen, la discontinuidad se ve reflejada en la irregularidad de la
superficie z = f (x, y) cerca del origen de coordenadas.
En particular entonces f (x, y) =
6.4.
Derivadas Parciales
Sea f (x, y),la derivada parcial de f con respecto a x es aquella funcion,

f
o por D1 f , tal que su valor, en cualquier punto (x, y)
representada por
x
en el dominio de respectivo esta dado por
f (x + h, y) f (x, y)
h0
h
D1 f = lm
si este lmite existe. Analogamente, la derivada parcial de f con respecto a

f
y es aquella funcion, representada por
o por D2 f , tal que su valor, en
y
cualquier punto (x, y) en el dominio de f esta dado por
f (x, y + h) f (x, y)
h0
h
D2 f = lm
si este lmite existe.
97
c
2012
Mynor Chacon Daz
6.4
MA1003
En general si f (x1 , x2 , . . . , xn ) es un campo escalar, entonces

f (x1 , . . . , xi + h, . . . , xn) f (x1 , . . . , xi , . . . , xn )
h0
h
representa la derivada de f con respecto a su i-esima variable.
Ejemplo
Di f = lm
Sea f (x, y) = 3x2 y y 2 , se va a calular la derivada parcial con respecto a y:

f
y
f (x + h, y) f (x, y)
h0
h
2
[3x (y + h) (y + h)2 ] [3x2 y y 2 ])
lm
h0
h
2
2
3x y + 3x h (y + h)2 3x2 y + y 2
lm
h0
h
3x2 h y 2 2yh h2 + y 2
lm
h0
h
3x2 h 2yh h2
lm
h0
h
2
lm 3x 2y h
= lm
=
=
=
=
=
h0
= 3x2 2y
Note que, el resultado se puede obtener al derivar la expresion f (x, y) =

3x2 y y 2 con respecto a la variable y y considerando a x como si fuera una
constante. Analogamente, procediendo de esta manera se deduce que:
f
= 6xy
x
Ejemplo

Sea f (x, y, z) = ln x2 + y 2 + z 2 , entonces:
f
2z
= 2
z
x + y2 + z2
c
2012
Mynor Chac
on Daz
98
MA1003
6.4
Ejemplo
En este ejemplo se utiliza la TI-89 y MuPAD para obtener derivadas parciales
2
de f (x, y) = ex /y .
En la TI-89, siga la ruta F3>d( derivar, luego digite la funcion e indique la
variable con respecto a la cual se va a derivar
Para derivar con MuPAD, digite diff(exp( x 2 / y),x), de esta manera

se obtiene
En ambos casos se obtiene

x2
2xe y
f
=
x
y
99
c
2012
Mynor Chacon Daz
6.5
6.5.
MA1003
Derivadas Parciales de orden superior
Sea f (x, y),la derivada parcial de f con respecto a x o a y tambien son

funciones de dos variables. Y si las derivadas de estas funciones existen,
reciben el nombre de segundas derivadas parciales de f . La derivada de f
x
con respecto a y se denota de varias formas:
f
x

=
2f
= fxy = D12 f
y x
Si se deriva de nuevo con respecto de la variable x, se obtiene la siguiente

derivada parcial de orden 3:
2f
y x

=
3f
= fxyx = D121 f
xy x
Sea W = f (x, y, z), al derivar W con respecto a z, luego con respecto a y,

luego con respecto a x y por u
ltimo con respecto a x, se obtiene:
4W
= Wzyxx = D3211 W
x2 yz
Ejemplo
Sea W = f (x, y, z) = arctan(3xyz), se va a calcular
3f
3f
=
.
zx x
zx2
Para realizar este calculo con MuPAD, digite diff(arctan(3*x*y*z),x,x,z)

y obtenga el resultado inmediatamente:
c
2012
Mynor Chac
on Daz
100
MA1003
6.6
En la TI debe digitar d(d(d(tan1 (3*x*y*z),x),x),z), el resultado obtenido

es
3f
162z 2 (3z 2 x2 y 2 1) xy 3
=
zx2
(9z 2 x2 y 2 + 1)3
Tal como se aprecia en la siguientes imagenes:
El lector debe verificar que ambas expresiones algebraicas, la obtenida de la

TI y la de MuPAD, son iguales.
6.6.
Diferenciabilidad de z = f (x, y)
Recuerdese que si y = f (x) se define y = f (x0 + h) f (x0 ), en tal caso

x = (x0 + h) x0 = h. Ademas la definicion de derivada establece:
lm
x0
y
= f 0 (x0 )
x
Ese lmite puede rescribirse de la siguiente forma

y
f 0 (x0 )
x
O, de otra manera si 0, cuando x 0, entonces

101
c
2012
Mynor Chacon Daz
6.6
MA1003
y
= f 0 (x0 ) +
x
Esto significa que y = f 0 (x0 )y + x, siempre que f 0 (x0 ) exista.

La generalizacion de lo anterior a una funcion de dos variables da como
consecuencia la siguiente definicion.
Definici
on
El campo escalar z = f (x, y) es diferenciable en P0 = (x0 , y0 ), si D1 f (P0 )
y D2 f (P0 ) existen y se cumple
z = D1 f (P0 ) x + D2 f (P0 ) y + 1 x + 2 y
siempre que i 0, i = 1, 2, cuando x 0, y 0, y ademas

z = f (x, y) f (x0 , y0 ), y = x x0 , x = y y0
En tal caso se dice que z = f (x, y) es diferenciable en P0 = (x0 , y0 ).
Se puede probar que si D1 f (x, y) , D2 f (x, y), son funciones continuas en

una region abierta R, entonces f es diferenciable en R. Ademas si f es
diferenciable en R se puede garantizar que f es continua en R.
Se debe tener cuidado pues, la sola existencia de D1 f (P0 ) , D2 f (P0 ) NO

garantiza la diferenciabilidad de f en P0 , un ejemplo de esto lo brinda el
siguiente caso.
c
2012
Mynor Chac
on Daz
102
MA1003
6.6
Ejemplo
Sea
xy
, si (x, y) 6= (0, 0)
2
2
f (x, y) = x + y
0,
si (x, y) = (0, 0)
Si nos acercamos a (x, y) por el eje X, es decir, con puntos de la forma (x, 0),
se tiene
lm
f (x, y) =
(x,y)(0,0)
lm
x2
(x,y)(0,0)
xy
0
= lm 0 = 0
= lm 2
2
x0
x0 x + 0
+y
Mientras que acercandonos a (x, y) por la recta y = x se tiene
lm
(x,y)(0,0)
f (x, y) =
lm
(x,y)(0,0)
x2
1
1
xy
=
l
m
= lm =
2
2
2
x0 2x
x0 2
x +y
2
lm
(x,y)(0,0)
f (x, y) @
Esto permite garantizar que f no es continua en (0, 0), en consecuencia f

no es diferenciable en (0, 0).
Ahora analizaremos la existencia de las derivadas parciales en (0, 0). Seg

un
la definicion
f (x + h, y) f (x, y)
f
(0, 0) = lm
h0
x
h
f (0 + h, 0) f (0, 0)
= lm
h0
h
f (h, 0) 0
= lm
h0
h
f (h, 0)
= lm
h0
h
0
= lm
h0 h
= 0
103
c
2012
Mynor Chacon Daz
6.6
MA1003
f (x + h, y) f (x, y)
f
(0, 0) = lm
h0
y
h
f (0, 0 + h) f (0, 0)
= lm
h0
h
f (0, h) 0
= lm
h0
h
f (0, h)
= lm
h0
h
0
= lm
h0 h
= 0
En resumen, las derivadas parciales de primer orden
f
f
(0, 0),
(0, 0) exisx
y
ten, sin embargo f no es diferenciable en (0, 0).
De nuevo se hace incapie en que, ademas de la existencia de las derivadas

parciales D1 f (x, y), D2 f (x, y) se requiere que sean funciones continuas para
poder garantizar la diferenciabilidad. Resumiendo se tiene que:
f es diferenciable en R
6.6.1.
f f
,
existen y son continuas en R
x y
Aproximaci
on lineal de z = f (x, y)
Supongase que z = f (x, y) es diferenciable en P0 = (x0 , y0 ), esto significa

que
z = D1 f (P0 ) x + D2 f (P0 ) y + 1 x + 2 y
siempre que i 0, i = 1, 2, cuando x 0, y 0
Puesto que
z = f (x, y) f (x0 , y0 ), x = x x0 , y = y y0
c
2012
Mynor Chac
on Daz
104
MA1003
6.6
entonces
z = D1 f (P0 ) x + D2 f (P0 ) y + 1 x + 2 y
f (x, y) f (x0 , y0 ) = D1 f (P0 ) (x x0 )+D2 f (P0 ) (y y0 )+1 x +2 y
f (x, y) f (x0 , y0 ) + D1 f (P0 ) (x x0 ) + D2 f (P0 ) (y y0 )

{z
}
|
L(x,y)
la expresion L(x, y), representa una Aproximacion Lineal de f (x, y)

f (x, y) L(x, y)
.
Si retomamos la pen
ultima expresion, se tiene:
f (x, y) f (x0 , y0 ) + D1 f (P0 ) (x x0 ) + D2 f (P0 ) (y y0 )
z f (x0 , y0 ) + D1 f (P0 ) (x x0 ) + D2 f (P0 ) (y y0 )
0 f (x0 , y0 ) + D1 f (P0 ) (x x0 ) + D2 f (P0 ) (y y0 ) z
0 D1 f (P0 ) (x x0 ) + D2 f (P0 ) (y y0 ) z + f (x0 , y0 )

Se puede decir que para valores de (x, y) muy cercanos a (x0 , y0 ) y tomando
z0 = f (x0 , y0 ), se puede escribir:
D1 f (P0 ) (x x0 ) + D2 f (P0 ) (y y0 ) z + z0 = 0
Entonces cualquier punto (x, y, z) que se encuentre en un peque
no vecindario
de (x0 , y0 , z0 ), satisface la ecuacion del plano siguiente:
D1 f (P0 ) (x x0 ) + D2 f (P0 ) (y y0 ) (z z0 ) = 0
Este es el Plano Tangente a f (x, y) en el punto (x0 , y0 , z0 ). La ecuacion
del plano puede reescribirse as:
105
c
2012
Mynor Chacon Daz
6.6
MA1003
f
f
(P0 ) x +
(P0 ) y z =
x
y

f
f
(P0 ) ,
(P0 ) , 1 (x0 , y0 , z0 )
x
y

f
f
De lo anterior se deduce que
(P0 ) ,
(P0 ) , 1 es un vector normal
x
y
al plano tangente a f (x, y) en el punto (x0 , y0 , z0 ).
6.6.2.
Diferencial Total de z = f (x, y)
De nuevo recuerde que, si z = f (x, y) es diferenciable entonces se cumple

z = D1 f (P0 ) x + D2 f (P0 ) y + 1 x + 2 y
siempre que i 0, i = 1, 2, cuando x 0, y 0.
Definici
on
Si z = f (x, y) es diferenciable en (x, y), se define como la Diferencial
Total de f, a la funcion df = dz definida como sigue
dz = D1 f (x, y) x + D2 f (x, y) y
En particular si z = x, D1 f (x, y) = 1, D2 f (x, y) = 0 y la ecuacion anterior se escribe dz = dx = x. De forma analoga, si z = y, D1 f (x, y) =

0, D2 f (x, y) = 1 el diferencial se escribe dz = dy = y. en conclusion las
diferenciales de las varibles independientes se definen como:
dx = x,
dy = y
Entonces
La Diferencial Total de f (x, y) es la funcion
dz =
c
2012
Mynor Chac
on Daz
f
f
(x, y) dx +
(x, y) dy
x
y
106
MA1003
6.6
Debido a que dz suele ser mas facil de calcular que z, empleamos el hecho
de que z dz en algunas aproximaciones.
Ejemplo
Se desea construir un recipiente de metal, con tapa, en forma de cilindro
circular recto con una altura de 17 cm, radio de la base de 8 cm y un espesor
de 0.1 cm. Si el costo del metal que se va a utilizar es de 52 colones por cm3 .
Use diferenciales para determinar el costo del metal en la manufacturacion
del recipiente.
Soluci
on
El volumen exacto del metal del recipiente, es la diferencia entre los vol
umenes
de dos cilindros. El mas peque
no(el interno) tiene altura y0 = 17 cm y radio
de la base x0 = 8 cm, el cilindro mas grande tiene:
altura = y = y0 + espesor de la tapa + espesor de la base = 17+0,1+0,1 = 17,2
radio = x = x0 + espesor de la pared cilndrica = 8 + 0,1 = 8,1

De lo anterior se tiene que
dx = x x0 = 0,1
dy = y y0 = 0,2
El volumen del recipiente esta dado por la funcion, V (x, y) = x2 y y las
V
V
derivadas parciales estan dadas por
= 2xy,
= x2 , por lo tanto el
x
y
diferencial total de la funcion volumen es:
V
V
(x0 , y0 ) dx +
(x0 , y0 ) dy
x
y
= 2 x0 y0 dx + x20 dy
dV (x0 , y0 ) =
= 2 (17) (8) (0,1) + (17)2 (0,2)

= 85
107
c
2012
Mynor Chacon Daz
6.7
MA1003
Puesto que, V 85 cm3 , el costo total de la manufacturacion es aproximadamente 52(85) = 4420 colones.
Definici
on
Si z = f (x1 , . . . , xn ) y P = (e
x1 , . . . , x
en ) entonces el incremento de f en P
esta dado por
f (P ) = f (e
x1 + x1 , . . . , x
en + x1 ) f (P )
Si i , i = 1, 2...n y xi , i = 1, 2...n tales que i 0 cuando xi 0, y

tal que es posible escribir
z =
n
X
Di f (P ) xi +
i=1
n
X
i xi
i=1
se dice que f es diferenciable en P . en tal caso la Diferencial total de

z = f (x1 , . . . , xn ) es
dz =
6.7.
f
f
f
f
dx1 +
dx2 + . . . +
dxn1 +
dxn
x1
x2
xn1
xn
Regla de la cadena
Teorema
Sea z = f (x, y) una funcion diferenciable de x e y, y sean x = F (s, t) y
y = G(s, t). Si existen las derivadas parciales de primer orden de x y y con
respecto a s y t, entonces z es funcion de s y t y se cumple que
f x
f y
z
=
+
s
x s
y s
z
f x
f y
=
+
t
x t
y t
c
2012
Mynor Chac
on Daz
108
MA1003
6.7
En realidad el teorema anterior se puede generalizar a mas o menos varibales,

por ejemplo si w = f (x, y, z) es una funcion diferenciable de y x = F (s, t),
y = G(s, t) y z = G(s, t) y en caso de existir las derivadas parciales de
primer orden de x, y y z con respecto a s y t, entonces w es funcion de s y
t y se cumple que
w
f x
f y
f z
=
+
+
s
x s
y s
z s
w
f x
f y
f z
=
+
+
t
x t
y t
z t
Si w = f (x, y, z), pero x = F (t), y = G(t), z = G(t), entonces la regla de la

cadena se escribe de la siguiente manera
f dx
f dy
f dz
dw
=
+
+
dt
x dt
y dt
z dt
A continuacion, dos ejemplos representativos.
Ejemplo
s
z z
Sea z = f (x, y) = x2 + y, donde x = st, y = , obtenga
, .
t
s t
Soluci
on
Se tiene en este caso que
x
x
y
1 y
s f
f
= t,
= s,
= ,
= 2,
= 2x,
=1
s
t
s
t t
t
x
y
z
f x
f y
=
+
s
x s
y s
1
= 2xt +
t
1
= 2(st)t +
t
= 2st2 +
109
1
t
c
2012
Mynor Chacon Daz
6.7
MA1003
z
f x
f y
=
+
t
x t
y t
s
= 2xs 2
t
s
2
= 2s t 2
t
Se puede obtener el resultado sin recurrir a la regla de la cadena, sustituyendo las funciones x, y directamente en la funcion z
z = x2 + y = (st)2 +
s
s
= s2 t2 +
t
t
Derivando directamente se obtiene:

1
z
= 2st2 + ,
s
t
z
s
= 2s2 t 2
t
t
Ejemplo
Sea z = f (x, y) = x2 y 2 , donde x = cosh t, y = senh t, obtenga
Soluci
on
Se tiene que
dx
dy
f
f
= senh t,
= cosh t,
= 2x,
= 2y
dt
dt
x
y
dz
f dx
f dy
=
+
dt
x dt
y dt
= 2 cosh t senh t 2 senh t cosh t
= 0
De hecho z = x2 y 2 = cosh2 t senh2 t = 1

c
2012
Mynor Chac
on Daz
110
dz
=0
dt
z
.
t
MA1003
6.7
6.7.1.
Regla de la Cadena involucrando funciones no explcitas

y
R
R
diferenciable, se
Sea g : derivable y sea z = f (x, y) = xy + xg
x
f f
va a calcular
,
, luego se verificara que se cumple la igualdad
x y
z
z
+y
= xy + z
x
x
y
Note que f en primer instancia depende de dos variables pero g es una funcion univariable.
y
y
z
0 y
= y+g
+xg
x
x
x
x2
y y
y
= y+g
g0
x
x
x
Analogamente

z
0 y
= x+xg
y
x
y
= x + g0
x

1
x
De esta manera se tiene:
h
y y
y i
h
y i
z
z
+y
= x y+g
g0
+ y x + g0
x
y
x
x
x
x
y
y
y
= xy + xg
yg 0
+ xy + yg 0
x
x
x
y
= 2xy + xg
x
y
= xy + xy + xg
|
{z x }
z
= xy + z
111
c
2012
Mynor Chacon Daz
6.7
MA1003
MuPAD da la opcion de obtener derivadas parciales que involucran funciones

no explcitas,como se aprecia en la siguiente imagen:
Note que los

son los mismos obtenidos anteriormente, entendiendo
y resultados
0 y
=g
. Incluso se puede hacer el calculo algebraico de
que D(g)
x
x
x
z
z
+y
x
y
, usando MuPAD:
c
2012
Mynor Chac
on Daz
112
MA1003
6.7
Ahora, se va a considerar un ejemplo en donde w = f (u, v), diferenciable

con u u(x, y), v v(x, y)
Ejemplo

Sea w = f
1. Obtenga
x
, xy
y

diferenciable, con u =
x
, v = xy, entonces:
y
w w 2 w
,
,
, sin utilizar software.
x y xy
3w
2. Obtenga
, usando Mupad, debe escribir el resultado en el
x3
formato acostumbrado.

w
w
3. Determine w(2, 2) =
(2, 2),
(2, 2) , conociendo que D1 f (1, 4) =
x
y
4, D2 f (1, 4) = 1. Mas adelante se vera que al vector

w(2, 2) =

w
w
(2, 2),
(2, 2)
x
y
se le conoce como el gradiente y que el vector

w
w
(2, 2),
(2, 2), 1
x
y

representa un vector normal (perpendicular) a la superficie w = f
x
, xy
y
en el punto x = 2, y = 2.
Soluci
on

1. Puesto que w = f
x
, xy
y

diferenciable, con u =
x
, v = xy, entonces:
y
w
1
= fu + yfv
x
y
x
w
= 2 fu + xfv
y
y
113
c
2012
Mynor Chacon Daz
6.7
MA1003

2w
2w
1
1 x
x
=
= 2 fu +
fuu + xfuv + fv + y
fvu + xfvv
xy
yx
y
y y2
y2
x
x
x
1
fvu + xyfvv
fu 3 fuu + fuv + fv +
2
y
y
y
y
1
x
= 2 fu + fv 3 fuu + xyfvv
y
y
2. Usando Mupad se obtiene
3w
:
x3
En notacion usual esto significa:
"

yfuvv +
3w
1
=
y
y
yf
+
f
+
vvv
uvv
x3
y
y
fuuv
y
# "

yfuuv +
1
+ y yfuvv + fuuv +
y
y
3
1
= y 3 fvvv + 3yfuvv + fuuv + 3 fuuu
y
y
Observese que MuPAD, tambien hace el trabajo de acomodar los terminos, como se puede apreciar en la imagen siguiente:
c
2012
Mynor Chac
on Daz
114
fuuu
y
1
y
MA1003
6.7

x
3. Determine un vector perpendicular a la superficie w = f
, xy en
y
el punto x = 2, y = 2, conociendo que D1 f (1, 4) = 4, D2 f (1, 4) = 1
El vector es:

w(2, 2) =

w
w
(2, 2),
(2, 2)
x
y
Se sabe que:
w
1
= fu + yfv
x
y
w
x
= 2 fu + xfv
y
y
Ademas observe que si x = 2, y = 2 u = 1, v = 4, de manera que:
w
1
(2, 2) = fu (1, 4) + 2fv (1, 4)
x
2
1
= D1 f (1, 4) + 2D2 f (1, 4)
2
1
= (4) + 2(1) = 0
2
analogamente
w
2
(2, 2) = 2 fu (1, 4) + 2fv (1, 4)
y
2
115
c
2012
Mynor Chacon Daz
6.7
MA1003
1
= D1 f (1, 4) + 2D2 f (1, 4)
2
1
= (4) + 2(1) = 4
2
As entonces
w(2, 2) = (0, 4)

x
y el vector normal a la superficie w = f
, xy es
y

w
w
(2, 2),
(2, 2), 1 = (0, 4, 1)
x
y
Ejemplo: Regla de la cadena aplicada a un sistema de ecuaciones

Considere el siguiente sistema de ecuaciones, en donde se supone que u =
u(x, y), v = v(x, y) son diferenciables y donde x,y son varibles independientes.
(
x+y =u+v
xu + yv = 1
u v u v
, , , , para esto se derivan las dos ecuaciones del
x x y y
sistema con respecto a la variable independiente x, as se obtiene:
(
1 + 0 = ux + v x
Se va a calcular
u + xux + yvx = 0
Es decir se obtiene el sistema:
(
ux + vx = 1
xux + yvx = u
Resolviendo para ux y vx se obtiene:

ux =
y+u
x+u
, vx =
yx
xy
analogamente si se derivan las dos ecuaciones del sitema original con respecto
a la variable, y, entonces se obtiene:
c
2012
Mynor Chac
on Daz
116
MA1003
6.7
0 + 1 = uy + v y
xuy + v + yvy = 0
Es decir se obtiene el sistema:

(
uy + v y = 1
xuy + yvy = v
Resolviendo para ux y vx se obtiene:

uy =
6.7.2.
x+v
y+v
, vx =
yx
xy
Derivaci
on Implcita
Supongase que F (x, y, z) = 0 es diferenciable y define implcitamente la

funcion z = f (x, y), f diferenciable. Derivando con respecto a, x, aplicando
la regla de la cadena a la ecuacion F (x, y, f (x, y)) = 0, se obtiene:
Fx 1 + Fy 0 + Fz zx = 0 zx =
Fx
Fz
Analogamente derivando con respecto a , y, se obtiene:

Fx 0 + Fy 1 + Fz zy = 0 zy =
Fy
Fz
De otro modo:
f
Fy
=
y
Fz
f
Fx
= ,
x
Fz
Ejemplo
Sea z 3 xy + yz + y 3 2 = 0 una ecuacion que define implcitamente la
f
f
funcion z = f (x, y), f diferenciable., Determine
(1, 1),
(1, 1).
x
y
Soluci
on
En este caso, F (x, y, z) = z 3 xy + yz + y 3 2, de manera que Fx = y,
Fy = x + z + 3y 2 y Fz = 3z 2 + y, as
f
y
= 2
,
x
3z + y
f
x + z + 3y 2
=
y
3z 2 + y
Obviamente en el punto (1, 1), x = 1, y = 1, al sustituir en

117
c
2012
Mynor Chacon Daz
6.8
MA1003
z 3 xy + yz + y 3 2 = 0

se obtiene z 3 + z 2 = 0 (z 1) z 2 + z + 2 = 0, cuya u
nica solucion
es z = 1, por lo tanto los valores a sustituir en las derivadas obtenidas son:
x = 1, y = 1, z = 1, el resultado obtenido es:
f
1
= ,
x
4
f
3
=
y
4
Ejemplo

La ecuacion g x2 + y 2 , y 2 + z 2 + x + y + z = 0, g diferenciable, define
z z
.
implcitamente z = f (x, y). Determine
x y
Soluci
on
Se tiene g x2 + y 2 , y 2 + z 2 +x+y+z = 0, g diferenciable, define implcita| {z } | {z }

u
v

mente z = f (x, y). Denote g x2 + y 2 , y 2 + z 2 + x + y + z = 0, entonces:
{z
}
|
F (x,y,z)
Fx
gu ux + gv vx + 1
2xgu + 1
z
=
=
=
x
Fz
gu uz + gv vz + 1
2zgv + 1
Analogamente
z
Fy
2ygu + 2ygv + 1
gu uy + gv vy + 1
=
=
=
y
Fz
gu uz + gv vz + 1
2zgv + 1
6.8.
Derivada direccional, interpretaci

on geom
etrica.
Recuerdese que si P = (x0 , y0 ), y f es una funcion de dos variables entonces

f (x0 + h, y0 ) f (x0 , y0 )
f
(P ) = lm
h0
x
h
f
f (x0 , y0 + h) f (x0 , y0 )
(P ) = lm
h0
y
h
Note que estos dos lmites pueden reescribirse como en la siguiente forma:
c
2012
Mynor Chac
on Daz
118
MA1003
6.8
f
f (x0 + h, y0 ) f (x0 , y0 )
(P ) = lm
h0
x
h
f [(x0 , y0 ) + (h, 0)] f (x0 , y0 )
= lm
h0
h
f [(x0 , y0 ) + h (1, 0)] f (x0 , y0 )
= lm
h0
h
f (P + h i) f (P )
= lm
h0
h
Analogamente
f
f (x0 , y0 + h) f (x0 , y0 )
(P ) = lm
h0
y
h
f [(x0 , y0 ) + (0, h)] f (x0 , y0 )
= lm
h0
h
f [(x0 , y0 ) + h (0, 1)] f (x0 , y0 )
= lm
h0
h
f (P + h j) f (P )
= lm
h0
h
Al derivar f con respecto a X, se considera la funcion z = f (x, y), como
una funcion de una sola varible z = f (x, y) = F (x) (pues Y se toma como
constante, y = y0 ), cuya grafica genera una curva en la superficie, z =
f (x, y), paralela al plano XZ. De forma identica si se asume x = x0 al
f
calcular
, se genera una curva plana z = f (x, y) = G(y) en un plano
y
paralelo al plano YZ. En ambos casos, las curvas planas, vistas como lo que
son en realidad, curvas en el espacio son las siguientes:
s(t) = (t, y0 , F (t)) ,
r(t) = (x0 , t, G(t))
Puesto que s(t), se encuentra en un plano paralelo al XZ, su trayectoria se

hace en la direccion del eje X, o sea siguiendo la direccion del vector i, de
forma similar la trayectoria de r(t) se hace en la direccion del eje Y, en la
direccion del vector j.
f (P + h i) f (P )
f
, se obtiene la
(P ) = lm
h0
h
x
variacion de z cuando a partir de Q = (x0 , y0 , f (x0 , y0 )) el punto (x, y0 , f (x, y0 ))
se mueve en direccion del vector i.
En resumen, al calcular
119
c
2012
Mynor Chacon Daz
6.8
MA1003
f
f (P + h j) f (P )
(P ) = lm
, se obtiene la
h0
y
h
variacion de z cuando a partir de Q = (x0 , y0 , f (x0 , y0 )) el punto (x0 , y, f (x0 , y))
se mueve en direccion del vector j.
Analogamente,, al calcular
f
f
(P ) ,
(P ), pueden ser interprex
y
tadas como la derivada de f (x, y), en el punto P (x0 , y0 ) en la direccion del
vector i o j respectivamente, una forma de resaltar este hecho es utilizando
la siguiente notacion:
Por lo tanto, las derivadas parciales
Di f (P ) =
f (P + h i) f (P )
f
(P ) = lm
h0
x
h
Dj f (P ) =
f (P + h j) f (P )
f
(P ) = lm
h0
y
h
Si los vectores: i, j; son sustituidos por un vector unitario u, entonces se

obtiene:
f (P + h u) f (P )
h
y recibe el nombre de Derivada Direccional de z = f (x, y) en el punto
P (x0 , y0 ) y en la direcci
on del vector unitario u.
Du f (P ) = lm
h0
Se puede probar que, para efectos de calculo, se procede como en la siguiente

definicion:
f (P + h u) f (P )
= f (P ) u
h0
h
Du f (P ) = lm
Como ejemplo considere la funcion z= f (x, y)= 4 x2 y 2 , z 0, el punto

1 1
P = (1, 1) y el vector director u = , . Observe que u es el vector
2 2
unitario en la direccion i + j y que P = (1, 1) implica Q = (1, 1, f (1, 1)) =
(1, 1, 2).
c
2012
Mynor Chac
on Daz
120
MA1003
6.8
Si se fija x = x0 = 1, la funcion z = f (1, y) = G(y) = 4 12 y 2 = 3 y 2 , de

manera que r(t) = (1, t, G(t)) = (1, t, 3 t2 ), la siguiente imagen muestra la
curva r(t):
Si se fija y = y0 = 1, la funcion z = f (x, 1) = F (x) = 4 x2 12 = 3 x2 , de

manera que s(t) = (t, 1, F (t)) = (t, 1, 3 t2 ), la siguiente imagen muestra
ambas curvas: r(t), s(t):
Al movernos sobre la superficie en la direccion del vector i + j, pasando

por Q = (1, 1, 2) se genera la curva w(t) = (t, t, 4 t2 ). Luego las derivadas
direccionales en P , en la direccion de i, j y de i + j son respectivamente:
f
(P ) = 2,
x
f
(P ) = 2
y

1 1
Du f (P ) = f (P ) u = (2, 2) ,
= 2 2
2 2
121
c
2012
Mynor Chacon Daz
6.9
6.9.
MA1003
Plano Tangente y vector normal a una superficie
De nuevo se va a utilizar el ejemplo anterior para determinar la ecuacion

del plano tangente y un vector normal este. Recuerde las dos curvas cuyas
trayectorias son paralelas a los planos YZ y XZ, respectivamente:

r(t) = 1, t, 3 t2 ,
s(t) = t, 1, 3 t2
Los vectores tangentes repectivos son

r0 (1) = u = (0, 1, 2)
s0 (1) = v = (1, 0, 2)
La imagen siguiente muestra las curvas y sus vectores tangentes
Puesto que los vectores u, v son tangentes a las curvas, es natural que el
vector w = v u = (2, 2, 1) sea el vector normal al Plano Tangente a la
superficie en el punto Q
c
2012
Mynor Chac
on Daz
122
MA1003
6.9
Por lo tanto la ecuacion del plano es 2x + 2y + z = 6, la siguiente imagen

muestra parte del plano tangente:
Note que esto coincide con lo expuesto en la seccion 6.6.1, donde se dedujo
que el vector normal del plano tangente a z = f (x, y) en el punto (x0 , y0 , z0 )
esta dado por

f
f
(P0 ) ,
(P0 ) , 1
x
y
o su vector paralelo:

f
f
(P0 ) ,
(P0 ) , 1
x
y
123
c
2012
Mynor Chacon Daz
6.9
MA1003
En este ejemplo

f
f
(P0 ) ,
(P0 ) , 1 = ((2), (2), 1) = (2, 2, 1)
x
y
En el caso en que la superficie este dada por la ecuacion F (x, y, z) = 0,

siendo z = f (x, y) entonces, la derivacion implcita dice que el normal al
plano tangente es:

f
f
Fx
Fy
(P0 ) ,
(P0 ) , 1 = (P0 ), (P0 ), 1
x
y
Fz
Fz
Evidentemente un vector paralelo que se puede utilizar en lugar del anterior

es:
(Fx (P0 ), Fy (P0 ), Fz (P0 )) = F (P0 )
En resumen:

f
f
z = f (x, y) tiene como normal al plano tangente el vector , , 1
x y

f f
o
,
, 1 .
x y
F (x, y, z) = 0 tiene como normal al plano tangente el Gradiente de
F , a saber el vector F = (Fx , Fy , Fz ) o F = (Fx , Fy , Fz )
Ejemplo
Calcular la derivada direccional de f (x, y, z) = x2 + y 2 z 2 en el punto
(1,1,1), a lo largo de la curva de interseccion de las superficies:
(
x2 + y 2 + 2z 2 = 4
z = exy
Soluci
on
Las superficies involucradas que determinan la curva son F (x, y, z) = x2 +
y 2 + 2z 2 4 = 0 y G(x, y, z) = z exy = 0, por lo tanto los vectores
normales a sus respectivos planos tangentes son:
F = (2x, 2y, 4z) F (1, 1, 1) = (2, 2, 4)

c
2012
Mynor Chac
on Daz
124
MA1003
6.9
G = (exy , exy , 1) G(1, 1, 1) = (1, 1, 1)

Un vector director a lo largo de la curva puede obtenerse por
F (1, 1, 1) G(1, 1, 1) = (2, 2, 4) (1, 1, 1) = (2, 6, 4) = 2(1, 3, 2)

El vector u es
1
u=
(1, 3, 2) =
k(1, 3, 2)k
1
3
2
, ,
14 14 14
Por otro lado, se tiene f = (2x, 2y, 2z) f (1, 1, 1) = (2, 2, 2), de
manera que la derivada direccional es:

Du f (1, 1, 1) = f (1, 1, 1) u = (2, 2, 2)
6.9.1.
1
3
2
, ,
14 14 14
12
=
14
Valores extremos de la derivada direccional
Recuerde la definicion de la derivada direccional de f en el punto P en la

direccion del vector unitario u
Du f (P ) = f (P ) u
Gracias a propiedades del producto punto se puede escribir:
Du f (P ) = f (P ) u = kf (P )k kuk cos = kf (P )k cos

|{z}
1
Siendo el angulo entre los vectores kf (P )k y kuk 0 , como se

2
aprecia en la imagen siguiente
125
c
2012
Mynor Chacon Daz
6.9
MA1003
El valor maximo que puede alcanzar el coseno es 1, esto ocurre en el momento

en que el valor del angulo sea 0, es decir cuando los vectores kf (P )k y
kuk son tienen la misma direccion.
Por lo tanto el valor maximo de la derivada direccional es, kf (P )k, se

alcanza en la direccion del vector gradiente de f en P, f (P ) .
Analogamente el valor mnimo de la derivada direccional es , kf (P )k,
se alcanza en la direccion contraria al vector gradiente de f en P.
c
2012
Mynor Chac
on Daz
126
MA1003
6.9
Ejemplo
La temperatura en una caja es aproximada por
T (x, y, z) = xyz(1 x)(2 y)(3 z)
x [0, 1] , y [0, 2] , z [0, 3]

1 1 1
Si un mosquito se localiza en P1
, ,
en que direccion debe volar para
5 3 4
calentarse mas rapido, cual es la mayor razon de cambio en esa direccion?
(debe utilizar el software adecuado)
Soluci
on
Debido a que las derivadas de T son un poco incomodas de calcular, se
utilizara la TI-89, primero se define la funcion temperatura como se aprecia
en la imagen siguiente:
Luego se define la funcion f (x, y, z), de la siguiente forma:

127
c
2012
Mynor Chacon Daz
6.9
MA1003

d
d
d
g(x, y, z) =
(t(x, y, z), x) ,
(t(x, y, z), y) ,
(t(x, y, z), z)
dx
dx
dx

1 1 1
11 11 2
1 1 1
, ,
=g
, ,
=
, ,
De manera que f
5 3 4
5 3 4
48 75 9

11 11 2 7 32641

1 1 1

, ,
=
, ,
=
y f
5 3 4 48 75 9
3600

En resumen, el mosquito debe viajar en la direccion del vector
7 32641
y la mayor razon de cambio de la temperatura es
3600
c
2012
Mynor Chac
on Daz
128
11 11 2
, ,
48 75 9
CAPITULO
Extremos de Campos Escalares
7.1.
Extremos Relativos
Definici
on
Un campo escalar f tiene m
aximo relativo en a S Rn si f (x) f (a)
para todo x en un entorno abierto B(a) S, analogamente tiene un mnimo
relativo en a si f (x) f (a) para todo x en un entorno abierto B(a) S.
En ambos casos se dice que a es un extremo de f.
Definici
on
Sea f diferenciable en a, si f (a) = 0 se dice que a es un punto estacionario o crtico de f.
Definici
on
Sea f diferenciable en un punto estacionario a. Si todo entorno abierto
B(a) contiene x tal que f (x) < f (a) y otros para los que f (x) > f (a) se dice
que a es un punto de ensilladura de f.
Recuerde que si f es diferenciable, entonces
f (a + y) f (a) = f (a)y + kyk E(a, y)
donde
a = (x1 , ...xn ), y = (x1 , ..., xn ), lm E(a, y) = 0
y0
129
7.1
MA1003
Si a es un punto estacionario, entonces f (a) = 0, por lo que:

f (a + y) f (a) = 0 y + kyk E(a, y)
f (a + y) f (a) = kyk E(a, y)
En consecuencia, la naturaleza del punto estacionario queda determinada
por el signo del termino kyk E(a, y).
Si
E(a, y) > 0 kyk E(a, y) = f (a + y) f (a) > 0
entonces
f (a + y) > f (a)
En resumen, se tiene:
E(a, y) > 0 f (a + y) > f (a) f (a) es mnimo
E(a, y) < 0 f (a + y) < f (a) f (a) es maximo
F
ormula de Taylor de segundo orden
Si f es un campo escalar con derivadas parciales segundas Dij f continuas en
un entorno B(a), entonces para todo y Rn , tal que (a + y) B(a) tenemos
f (a + y) f (a) = f (a) y +
1
yH(a)y t + kyk2 E2 (a, y)
2!
en donde E2 (a, y) 0 cuando y 0.

Este teorema, garantiza que
kyk E(a, y) =
1
yH(a)y t + kyk2 E2 (a, y)
2!
Por lo tanto la naturaleza del punto estacionario queda determinada por el

signo del termino 2!1 yH(a)y t
La expresion yH(a)y t , representa una forma cuadr
atica, H(a) es una
matriz simetrica con entradas hij = Dij f (a), el siguiente teorema da una luz
acerca del signo de la forma cuadratica.
c
2012
Mynor Chac
on Daz
130
MA1003
7.1
Teorema
Dada la formula cuadratica yH(a)y t , se cumple:
yH(a)y t > 0, y 6= 0 si y solo s todos los valores propios de H(a)
son positivos.
yH(a)y t < 0, y 6= 0 si y solo s todos los valores propios de H(a)

son negativos.
Recuerde que el signo de E(a, y) coincide con el signo de yH(a)y t y que

este signo determina si a es maximo, mnimo o punto de silla de f, por lo
tanto se tiene el siguiente teorema.
Naturaleza del punto estacionario
Sea f un campo escalar con derivadas parciales segundas continuas Dij f en
una entorno B(a), y designemos con H(a), la matriz hessiana en un punto
estacionario a.
Si todos los autovalores de H(a) son positivos, f tiene un mnimo
relativo en a.
Si todos los autovalores de H(a) son negativos, f tiene un m

aximo
relativo en a.
Si H(a) tiene autovalores positivos y negativos, f tiene un punto

de ensilladura en a.
131
c
2012
Mynor Chacon Daz
7.1
MA1003
Para el caso en que f : R2 R, se tiene
H(a) =
2f
(a)
x2
2f
(a)
xy
2f
(a)
yx
2f
(a)
y 2
En este caso el teorema anterior se traduce en el siguiente Criterio de las

derivadas de segundo orden, donde det(H(a)) = :
Si > 0, fxx (a) > 0 f posee mnimo relativo en a
Si > 0, fxx (a) < 0 f posee m

aximo relativo en a
Si < 0 f posee punto de ensilladura en a
Si = 0 No hay criterio
Para el caso en que f : R3 R, se tiene
fxx (a) fxy (a) fxz (a)
fyx (a) fyy (a) fyz (a)
H(a) =
f
(a)
f
(a)
f
(a)
zx
zy
zz
Si 1 = fxx (a), 2 =
!
fxx (a) fxy (a)
, 3 = det(H(a)), entonces
fyx (a) fyy (a)
Si 1 > 0, 2 > 0, 3 > 0 f posee mnimo relativo en a
Si 1 < 0, 2 > 0, 3 < 0 f posee m

aximo relativo en a
En cualquier otro caso f no tiene extremo en a
c
2012
Mynor Chac
on Daz
132
MA1003
7.1
Ejemplo 1
Considere la imagen siguiente:
1. Determine el sistema que permite determinar los puntos P y Q.

2. Resuelva el sistema y deduzca las coordenadas de P y Q.
3. Determine la matriz Hessiana en el punto P y verifique que P es un
mnimo relativo de f(x,y)
Soluci
on
1. Sea f (x, y) =
2x + 2y + 1
x2 + y 2 + 1
2. Obteniendo las derivadas con MuPAD, se obtienen:
133
c
2012
Mynor Chacon Daz
7.1
MA1003
De manera que el sistema simplificado a resolver es:

(
x 2xy x2 + y 2 + 1 = 0
y 2xy + x2 y 2 + 1 = 0
3. De nuevo MuPAD nos da la solucion:

Los puntos crticos reales son (1, 1),
pondientes son:
c
2012
Mynor Chac
on Daz
134

1 1
,
, las imagenes corres2 2
MA1003
7.1

1 1
Es decir los candidatos a P y Q son (1, 1, 1),
, , 2 , de acuerdo
2 2
a la posici
en la grafica de f(x,y), se deduce que P = (1, 1, 1) y
on
1 1
, ,2
Q=
2 2
4. Las derivadas de segundo orden obtenidas con MuPad son:
135
c
2012
Mynor Chacon Daz
7.1
MA1003
Utilizando la TI-89 para la evaluacion de las derivadas en el punto

P = (1, 1, 1), se obtiene la siguiente matriz Hessiana:
H(1, 1, 1) =
2
3
2
3
2
4
Puesto que = det(H) = > 0 y fxx (P ) = > 0, se concluye que el
3
3
punto P representa un mnimo relativo de f.
c
2012
Mynor Chac
on Daz
136
MA1003
7.1
Ejemplo 2
Determine los extremos relativos de f (x, y) = x4 + y 4 2x2 + 4xy 2y 2 .
Soluci
on
1. Sea f (x, y) = x4 + y 4 2x2 + 4xy 2y 2 , entonces:
fx = 4x3 4x + 4y = 0
fy = 4y 3 + 4x 4y = 0
simplificando se obtiene:
x3 x + y = 0
y3 + x y = 0
como se aprecia en la imagen siguiente los soluciones son: (0, 0), ( 2, 2), ( 2, 2)
Usando MuPAD se obtiene:

137
c
2012
Mynor Chacon Daz
7.1
MA1003
En (0, 0), se tiene
H=
!
4 4
4 4
Puesto que det(H) = 0, se tiene que en (0, 0) no hay criterio.
En ( 2, 2), ( 2, 2), se tiene
H=
!
20 4
4 20
Puesto que det(H) = 384 > 0, fx = 20 > 0, se tiene que en ( 2, 2), ( 2, 2)

existen m`nimos relativos.
Por otro lado los valores propios de H correspondiente a los puntos
( 2, 2), ( 2, 2), son: 24, 16 como se aprecia en la siguiente imagen:

c
2012
Mynor Chac
on Daz
138
MA1003
7.1
Por tener valores propios positivos se deduce que en ( 2, 2), ( 2, 2)

existen m`nimos relativos.
Los valores propios de la matriz correspondiente en (0,0) son: -8, 0, por
lo tanto no hay criterio
139
c
2012
Mynor Chacon Daz
7.1
MA1003
Ejemplo 3
Determine los extremos relativos de
f (x, y, z) = x2 + y 2 + 3z 2 + yz + 2zx xy
Clasificar los puntos crticos utilizando los dos metodos expuestos en el material multimedia. Para obtener los valores propios de la matriz Hessiana
utilice MuPAD.
Soluci
on
Para determinar los extremos relativos de
f (x, y, z) = x2 + y 2 + 3z 2 + yz + 2zx xy
Luego de calcular las derivadas parciales de primer orden se plantea el sistema siguiente:
fx = 2x + 2z y = 0
fy = 2y + z x = 0
f = 6z + y + 2x = 0
z
Este es un sistema lineal con solucion u
nica (0,0,0).
Para el primer metodo se necesita la matriz Hessiana la cual se puede obtner
con MuPAD, la imagen siguiente muestra el resultado:
c
2012
Mynor Chac
on Daz
140
MA1003
7.1
Al evaluar en (0,0,0), obviamente los valores de la matriz no sufren ninguna

alteracion, entonces:
2 1 2
H(0, 0, 0) = 1 2 1
2
1 6
En tal caso
1 = fxx = 2,
2 = det
=3
1 2
2
3 = det(H) = 4
Puesto que
1 > 0, 2 > 0, 3 > 0
Se concluye que en (0,0,0), f alcanza un mnimo relativo.
Por otro lado, los valores propios de H(0,0,0) son todos positivos, por lo
tanto, el segundo de los metodos garantiza que f alcanza un mnimo relativo
en (0,0,0).
141
c
2012
Mynor Chacon Daz
7.1
MA1003
Ejemplo 4
Determine los extremos relativos de z = x2 2xy + y 2 2x + 2y. el sistema
a resolver para determinar puntos crticos es:
(
fx = 2x 2y 2 = 0
fy = 2x + 2y + 2 = 0
Evidentemente las soluciones de este sistema son infinitas, de hecho todos

los puntos de la recta x y 1 = 0 son puntos crticos. Por otro lado note
que fxy = 2, fxx = 2, fyy = 2, de manera que
H=
!
2 2
2
2
Puesto que = det(H) = 0, se tiene que no hay criterio para tales puntos
crticos. Los valores propios de H, son 0 y 4, de manera que tampoco hay
criterio:
c
2012
Mynor Chac
on Daz
142
MA1003
7.1
Que se puede hacer en tal caso?. Note que en este ejemplo se puede rescribir
z = x2 2xy + y 2 2x + 2y = (x y)2 2 (x y) = (x y) (x y 2)
De manera que sobre cualquier punto en la recta x y 1 = 0, se tiene
z = (x y) (x y 2) = (1)(1 2) = 1, ahora se calculara la imagen de
un punto arbitrario cercano a la recta, es decir, supongase que (x, y) cumple
x y = 1 y sea (x + h, y + k) un punto cercano, es decir h, k suficientemente
peque
nos. Entonces
f (x + h, y + k) (x + h y k)(x + h y k 2)
(x y + h k)(x y + h k 2)
(1 + h k)(1 + h k 2)
(h k + 1)(h k 1)
(h k)2 1
como (h k)2 0 (h k)2 1 + 1 (h k)2 1 1, entonces

f (x + h, y + k) 1 h, k peque
nos. Por lo tanto, se concluye que
f (x + h, y + k) f (x, y) x, y en la recta x y = 1
Por lo que se concluye que los puntos de la recta constituyen mnimos de la
funcion. En la siguiente imagen se aprecia que son mnimos absolutos.
143
c
2012
Mynor Chacon Daz
7.2
7.2.
MA1003
Extremos Absolutos
Recuerde que una funcion real f : [a, b] R alcanza un maximo absoluto

y un mnimo absoluto en [a, b], este resultado se puede extender a campos
escalares, f : R Rn R, siendo R una region del espacio n-dimensional
que contiene a los puntos de su frontera o borde y que ademas satisface que:
x R, M 0, tq |f (x)| M , a una region de este tipo se le conoce como
regi
on cerrada y acotada. el teorema correspondiente dice:
Teorema
Sea f : R Rn R, un campo escalar diferenciable, R una region cerrada y
acotada, entonces f alcanza extremos absolutos en un punto estacionario
o en un punto de la frontera de R.
Ejemplo
Determine los valores que representan el maximo y mnimo absoluto de la
funcion f (x, y) = x2 + y 2 xy + x + y en la region
c
2012
Mynor Chac
on Daz
144
MA1003
7.2
R = {(x, y, z)/ x 0, y 0, x + y 3}
Soluci
on
La region se observa en la imagen siguiente:
Del sistema
(
fx = 2x y + 1 = 0
fy = 2y x 4y + 1 = 0
se obtiene la solucion P = (1, 1), el cual representa el u

nico punto crtico.
Note que P R. Luego se debe analizar el comportamiento de f en cada una
de la curvas que constituyen la frontera o borde de la region R, en este caso
las curvas frontera se denotaran por: , , .
An
alisis sobre
(
:=
y=0
z = x2 + x ;
3 x 0

1
Luego z = 2x + 1 = 0 x =
en tal caso el punto crtico es , 0
2
(recuerde que sobre , y = 0), ahora considerando
los
extremos
del
intervalo

1
x = 3, x = 0, se obtienen tres puntos crticos , 0 ; (3, 0); (0, 0).
2
0
21 ,
An
alisis sobre
145
c
2012
Mynor Chacon Daz
7.2
MA1003
(
:=
x=0
z = y2 + y ;
3 y 0

1
Luego z = 2y + 1 = 0 y =
en tal caso el punto crtico es 0,
2
(recuerde que sobre , x = 0), ahora considerando
del intervalo

los extremos
1
; (0, 3); (0, 0).
y = 3, y = 0, se obtienen tres puntos crticos 0,
2
0
21 ,
An
alisis sobre
:=
y = 3 x
z = (x 3)2 + x2 x(x 3) 3
= 3x2 + 9x + 6 ; 3 x 0
Luego z 0 = 6x + 9= 0 x = 3
, puesto que sobre , y = 3
x,
2
3
3 3
entonces y = 3 ; de manera que el punto crtico es , ,
2
2 2
ahora considerando los extremos
del
intervalo
x
=
3,
x
=
0,
se
obtienen

3 3
tres puntos crticos ,
; (0, 3); (0, 3).
2 2
Por u
ltimo, se toman todos los candidatos obtenidos y se calcula las imagenes
respectivas bajo la funcion f (x, y) = x2 + y 2 xy + x + y.
c
2012
Mynor Chac
on Daz
146
MA1003
7.2
DE PUNTOS
CLASIFICACION
( x, y )
(1, 1)
( 21 , 0)
(3, 0)
(0, 0)
(0, 21 )
(0, 3)
( 32 , 23 )
f (x, y)
1
14
6
0
14
6
34
Condicion
Mnimo valor
Maximo valor
Maximo valor
-
En consecuencia, el valor maximo(Absoluto) de f (x, y) es 6 y se alcanza en

los puntos (3, 0, 6); (0, 3, 6); el valor mnimo(Absoluto) de f (x, y) es -1 y
se alcanza en el punto (1, 1, 1).
La superficie definida por la funcion f : R R2 R, es la siguiente:
La imagen siguiente muestra los puntos maximos y mnimos:

147
c
2012
Mynor Chacon Daz
7.3
MA1003
7.3.
Extremos Condicionados
Dada la funcion f : Rn R, estamos interesados en determinar los extremos

que satisfagan las m condiciones: g1 (x1 , x2 , ..., xn ) = 0, g2 (x1 , x2 , ..., xn ) =
0, ..., gm (x1 , x2 , ..., xn ) = 0. Para resolver este problema se construye la Funci
on de Lagrange:
L (x1 , x2 , ..., xn , 1 , 2 , ..., m ) = f (x1 , x2 , ..., xn ) +
m
X
i gi (x1 , x2 , ..., xn , )
i=1
Se puede demostrar que los puntos crticos condicionados de f (x1 , x2 , ..., xn )

coinciden con los puntos crticos de L (x1 , x2 , ..., xn , 1 , 2 , ..., m ); los valores i reciben el nombre de Multiplicadores de Lagrange. En resumen,
resolver el problema:
Optimizar: f (x1 , x2 , ..., xn , )
con las conciones: gi (x1 , x2 , ..., xn , ) = 0, i = 1, ...m
equivale a resolver:
Optimizar: L (x1 , x2 , ..., xn , 1 , 2 , ..., m )
Para clasificar los puntos crticos se utilizara uno de los dos metodos que el
Doctor Joseph C. Varilly1 expone en las notas elaboradas para uso de los estudiantes y profesores de la catedra de MA1003-Calculo III. A continuacion
1
Catedr
atico de la Escuela de Matem
atica, UCR
c
2012
Mynor Chac
on Daz
148
MA1003
7.3
se presenta un resumen del metodo 2

M
etodo 2
Sea f : Rn R, con las m condiciones:
g1 (x1 , x2 , ..., xn ) = 0, g2 (x1 , x2 , ..., xn ) = 0, ..., gm (x1 , x2 , ..., xn ) = 0
Estas
determinan una porcion V Rn de dimension n-m. Por ejemplo si
f : R3 R y las condiciones son g1 = x+y +z 3 = 0, g2 = x2 +y 2 z = 0,
entonces este sistema de ecuaciones determina una curva (interserccion de
ambas superficies) y los puntos de la curva constituyen una porcion de R3
de dimension n m = 3 2 = 1. Si fuese solo una condicion g1 = 0, los
puntos sobre el plano x + y + z 3 = 0 constituyen un subespacio de R3 de
dimension n m = 3 1 = 2.
En el supuesto caso anterior de las dos condiciones, recuerde que dos vectores
ortogonales a las dos superficies, el plano y el paraboloide, son sus vectores
gradientes g1 y g2 , de manera que los vectores tangentes(v) a la curva
de interseccion(porcion V) deben ser ortogonales a g1 y g2 , es decir se
debe cumplir v g1 = 0 y v g2 = 0.
En general, se tiene que en un punto a V , existen nm vectores tangentes
independientes v1 , v2 , ..., vnm que satisfacen:
vi g1 (a) = 0, vi g2 (a) = 0, ..., vi gn (a) = 0
i = 1, ..., n m
Analogamente a lo visto para extremos relativos, la matriz Hessiana de L
Lx1 x1 (a) Lx1 x2 (a) Lx1 xn (a)
A = D2 L(a) =
Lxn x1 (a) Lxn x2 (a) Lxn xn (a)

va a se
nalar la naturaleza del punto a V . Para eso, de la matriz A se
construye la matriz H = (hij ) de una forma cuadratica, definida por:
(hij ) = vi (A vj )
La matriz H es de tama
no (n m) (n m), sus menores principales
1 , 2 , ...nm se definen como:
149
c
2012
Mynor Chacon Daz
7.3
MA1003
Ak = matriz obtenida de H al eliminar las primeras k1 filas y columnas de H

k = det (Ak )
De manera que
1 = det (A1 ) = det (H) ; nm = det (Anm ) = hnm,nm
Luego se aplican los siguientes criterios:
1 > 0, 2 > 0, 3 > 0..., nm > 0
> 0, < 0, > 0...,

1
2
3
nm < 0, n m par
1 =
6 0
1 < 0, 2 > 0, 3 < 0..., nm < 0, n m impar
signos distintos no alternos y

nm < 0
Mnimo
M
aximo
M
aximo
Puntos de silla
1 = 0 No hay criterio
Ejemplo 1
En la siguiente imagen se aprecian los cuatro puntos: P, Q, R, S ubicados en
la interseccion del cilndro elptico 3x2 + 4xy + 6y 2 = 140 y del paraboloide
circular 8z = x2 + y 2 .
c
2012
Mynor Chac
on Daz
150
MA1003
7.3
1. Utilice Multiplicadores de Lagrange para plantear el problema que

le permita determinar los puntos P, Q, R, S.(NO LO RESUELVA)
2. Determine las coordenadas de los puntos P, Q, R, S.
3. Utilizar el metodo visto en clase para clasificar el punto P .
Soluci
on
1. Maximizar la funcion f (x, y) =
4xy + 6y 2 = 140

1 2
x + y 2 sujeto a la condicion 3x2 +
8
2. Para determinar las coordenadas de los puntos P, Q, R, S se resuelve

el sistema siguiente, basado en la funcion de Lagrange:
L(x, y, ) =

1 2
x + y 2 + 3x2 + 4xy + 6y 2 140
8
Lx = 4 x + 6x + 4y = 0
Ly = 41 y + 4x + 12y = 0
L = 3x2 + 4xy + 6y 2 140 = 0
Usando MuPAD se obtiene:
Es decir
1
, x = 2, y = 4
56
1
= , x = 2, y = 4
56
1
= , x = 2 14, y = 14
56
=
151
c
2012
Mynor Chacon Daz
7.3
MA1003
1
, x = 2 14, y = 14
56
Para obtener las coordenadas de P, Q, R, S se calcula las imagenes de

tales puntos:
f (2, 4) =
f (2, 4) =
5
= 2,5
2
5
= 2,5
2

35
f (2 14, 14) =
= 8,75
4
35
f (2 14, 14) =
= 8,75
4
De la imagen
Se deduce que los signos de las coordenadas de los puntos son P (+, ), Q(, ), R(, +),
por lo que
c
2012
Mynor Chac
on Daz
152
MA1003
7.3
35
2 14, 14,
4

35
R = 2 14, 14,
4

5
S = 2, 4,
2

5
Q = 2, 4,
2

P =

3. De la imagen siguiente se obtiene D2 L(x, y, ) =
6 +
4
1
4
4
12 +
!
1
4
y el gradiente g(x, y) = (6x + 4y, 4x + 12y)

En el punto , P =
35
1
2 14, 14,
, recuerde que se tiene = ,
4
16
por lo que
D2 L =
81 14
14 12
g = (8 14, 4 14)
= A,
Puesto que n m = 1, se necesita determinar un vector v tal que

v g = 0. En este caso sirve tomar v = (1, 2). La matriz H debe tener
tama
no 1 1, la u
nica entrada es:
h11
18 14
= (1 2)
14 12
153
!
25
1
= = 1
8
2
c
2012
Mynor Chacon Daz
7.3
MA1003
Entonces
1 = nm <0 y n m es impar, esto significa que el punto
35
es un maximo condicionado.
P = 2 14, 14,
4
Ejemplo 2
Determine todos los puntos crticos de f (x, y, z) = x4 + y 4 + z 4 , sujeta
a la condicion x2 + y 2 + z 2 = 1. Luego clasifique el punto estacionario

1
1 1
, , .
3
3 3
Soluci
on
La funcion de Lagrange correspondiente es

L(x, y, z, ) = x4 + y 4 + z 4 + x2 + y 2 + z 2 1
Luego:
Lx = 4x3 + 2x = 0 2x (x2 + ) = 0 x = 0 = 2x2
L = 4y 3 + 2y = 0 2y (y 2 + ) = 0 y = 0 = 2y 2
y
3
2
2
L
z = 4z + 2z = 0 2z (z + ) = 0 z = 0 = 2z
L = x 2 + y 2 + z 2 1 = 0
Observe que (0, 0, 0) no satisface la ecuacion L = 0, por lo que no es punto
crtico. En seguida se tomaran varios casos que se deducen del sistema:
Caso x = 0, y = 0, z = 0, pero el punto (0, 0, 0) no satisface la ecuacion
f = x2 + y 2 + z 2 1 = 0, entonces no es un punto crtico.
Caso x = 0, y = 0, z 6= 0, en tal caso 2z 2 + = 0, entonces z = 1
= 2, como consecuencia se tienen dos puntos crticos (0, 0, 1, 2).
Caso x = 0, y 6= 0, z = 0, en tal caso 2y 2 + = 0, entonces y = 1
Caso x 6= 0, y = 0, z = 0, en tal caso 2x2 + = 0, entonces x = 1
Caso x = 0, y 6= 0, z 6= 0, en tal caso
y 2 + z 2 = 1 2y 2 + = 0 2z 2 + = 0
De las u
ltimas dos ecuaciones se tiene y 2 = z 2 = , entonces
2
c
2012
Mynor Chac
on Daz
154
MA1003
7.3

1
y 2 + z 2 = 1 = 1 = 1 y = z =
2
2
2
1
1
como consecuencia se tienen cuatro puntos crticos (0, , , 1).
2
2
Caso x 6= 0, y = 0, z 6= 0, procediendo como en el caso anterior se
1
1
obtienen cuatro puntos crticos ( , 0, , 1).
2
2
Caso x 6= 0, y 6= 0, z = 0, analogamente se obtienen cuatro puntos
1
1
crticos ( , , 0, 1).
2
2
Caso x 6= 0, y 6= 0, z 6= 0, entonces se deduce que
x2 = y 2 = z 2 x = y, x = z
de manera que
1
x + y + z = 1 x = = 2x2 = 2
3
2

1
2
=
3
3
1
1
1
2
En consecuencia se tienen ocho puntos crticos ( , , , ).
3
3
3 3
En total se obtuvieron 26 puntos crticos:
(0, 0, 1, 2); (0, 1, 0, 2); (1, 0, 0, 2) 6 puntos crticos
1
1
1
1
1
1
(0, , , 1); ( , 0, , 1); ( , , 0, 1) 12 puntos crticos
2
2
2
2
2
2
1
1
1
2
( , , , ) 8 puntos crticos
3
3
3 3
155
c
2012
Mynor Chacon Daz
7.3
MA1003
Utilizando MuPAD, se obtienen los 26 puntos crticos:
c
2012
Mynor Chac
on Daz
156
MA1003
7.3
Ahora se va a analizar la naturaleza del punto

a=
1
1 1
, ,
3
3 3

, =
2
3
En este ejemplo se tiene n m = 3 1 = 2, por lo tanto se deben buscar 2

vectores independientes v1 , v2 tales que vi g(a) = 0, para esto observe que

g(x, y, z) = (2x, 2y, 2z) g(a) =
2 2
2
, ,
3
3 3
A simple vista se eligen los vectores v1 = (1, 1, 0) y v2 = (0, 1, 1), por otro
lado se tiene
12x2 + 2
0
0
D2 L(x, y, z; ) =
0
12y 2 + 2
0
2
0
0
12z + 2)

2
En el punto (a; ) = a;
se tiene:
3
157
c
2012
Mynor Chacon Daz
7.3
MA1003
0 0
A = D2 L(a; ) = 0 83 0
0 0 38
8
3
Luego se construye la matriz H:

h11 = v1 (Av1 ) = 1 1

h12 = v1 (Av2 ) = 1 1

h21 = v2 (Av1 ) = 0 1
h22 = v2 (Av2 ) = 0 1
H=
16
3
8
3
8
3
16
3
8
=
3
8
3
0 0
0
8
3
0 0
0
8
3
1 0
0
8
3
1 0
0

0 0
1
16
8
0 1 =
3
3
0
0 83

0 0
0
8
8
0 1 =
3
3
1
0 38

1
0 0
8
8
0 1 =
3
3
0 38
0

0 0
0
16
8
0 1 =
3
3
0 83
1
!
2 1
1 2
16
Como 1 = |H| = 8 > 0, nm = 2 =
> 0 se concluye que existe un
3

1
1 1
mnimo condicionado en a = , , .
3
3 3
c
2012
Mynor Chac
on Daz
158
CAPITULO
Integracion Multiple
Como es lo usual en un curso de calculo en una variable, toca ver ahora

el concepto de integral para funciones de dos y tres variables, esto nos lleva
a las integrales dobles y triples. Se estudiara los cambios de variables mas
usuales y las aplicaciones inmediatas a algunos conceptos propios de los
cursos de Fsica tales com: masa, centroide, densidad de un alambre, etc.
8.1.
8.1.1.
Integral doble
Integral doble en regi

on rectangular
En el curso de calculo en una variable dada una funcion continua

y posiZ
b
tiva definida en el intervalo [a, b], para poder definir la integral
f (x) dx,
a
como el area bajo la curva, se divide el intervalo en sub-intervalos y luego

se aproxima el area total con la suma de las areas de todos los rectangulos
generados por los sub-intervalos.
159
8.1
MA1003
f (x) dx
a
n
X
f (xi ) xi = lm
n7
i=1
n
X
f (xi ) xi
i=1
Analogamente en el caso de un campo positivo y continuo, f : R R2 R,

siendo R el rectangulo [a, b] [c, d], se subdivide este en n sub-rectangulos
de area xi yi , los rectangulos de la imagen anterior se convierten ahora
en prismas rectangulares de volumen f (x , y ) xi yi y el area bajo la
curva se convierte en el volumen entre el plano XY y la superficie
z = f (x, y). Siendo (x , y ) un punto en el sub-rectangulo correspondiente.
c
2012
Mynor Chac
on Daz
160
MA1003
8.1
ZZ
f (x) dx dy
X X
i
f (x , y ) xi yj
Ejemplo 1
Considere la funcion f (x, y) = 2x2 y 2 , en el intervalo R = [1, 1][1, 1]
se tiene que f es positiva. Utilice Winplot para calcular el volumen limitado
por el plano XY y la superficie z = 2 x2 y 2 dado por
ZZ
V =
f (x) dx dy
R
Utilizando una particion de R en 1002 = 10000 sub-rectangulos.
161
c
2012
Mynor Chacon Daz
8.1
MA1003
Soluci
on
Una vez graficada la funcion, desde la barra de men
u de Winplot se sigue
la ruta Una >Integrar f(x,y), seguidamente se completa la ventana tal
como aparece en la imagen siguiente:
Es decir se tiene:
ZZ
V =

2 x2 y 2 dx dy 5,33333
ZZ
Nota: Cuando la funcion no sea positiva, la integral
f (x) dx dy se define
R
de igual manera pero en tal caso el resultado no coincide con el volumen bajo
la superficie.
Teorema de Fubinni (1)
Sea f : R R2 R, un campo escalar diferenciable, R la region definida
por [a, b] [c, d]. Entonces
Z bZ
ZZ
f (x) dx dy =
R
f (x, y) dy dx
a
Aplicando el teorema anterior, para la funcion f (x, y) = 2 x2 y 2 en el

intervalo R = [1, 1] [1, 1], se tiene:
c
2012
Mynor Chac
on Daz
162
MA1003
ZZ
8.1
2x y

2 x2 y 2 dy dx
dx dy =
1

2x

2x
=
1
y3
y
3
1
dx
1

1
1
2
2x +
dx
3
3
10
2x2
3

dx
16
5,33333
3
Note que el resultado coincide con lo obtenido con Winplot y con lo que se
obtiene tambien con la TI-89, como se aprecia en la siguiente imagen:
8.1.2.
Integral doble en una regi

on general R
Teorema de Fubinni (2)

Sea f : R R2 R, un campo escalar diferenciable, R la region definida
por a x b, h(x) y g(x).
163
c
2012
Mynor Chacon Daz
8.1
MA1003
ZZ
h(x)
g(x)
f (x, y) dy dx
|
{z
}
primero repecto a y
f (x) dx dy =
Entonces
Analogamente si R es la region definida por h(y) x g(y), c y d.
ZZ
Entonces
Z
f (x) dx dy =
c
g(y)
f (x, y) dx dy
|
{z
}
primero repecto a x
h(y)
Ejemplo 2
ZZ
Sea R = [0, ] [0, ]. Calcule
|cos(x + y)| dx dy.

R
Primero se analizara el valor absoluto:
|cos(x + y)| =
c
2012
Mynor Chac
on Daz
cos(x + y),
si 0 x + y
cos(x + y),
si
cos(x + y), si
164
3
x + y 2
2
< x + y < 3
2
2
MA1003
8.1
Graficando las regiones que corresponden a las tres desigualdades, se obtiene
Si nos restringimos al cuadrado R = [0, ] [0, ] se observa tres subregiones

de R a saber R1 , R2 , R3 ,, las cuales representan el dominio de integracion de
tres integrales dobles.
ZZ
ZZ
|cos(x + y)| dx dy =
cos(x + y) dx dy
R1
ZZ
ZZ
cos(x + y) dx dy +
+
R2
cos(x + y) dx dy
R3
Al observar la imagen de la superficie f (x, y) = |cos(x + y)|, tambien se

puede ver la necesidad de partir la region en tres partes.
165
c
2012
Mynor Chacon Daz
8.1
MA1003
De manera que
ZZ
x
2
cos(x + y) dx dy =
cos(x + y) dy dx
0
R1
[sen(x + y)]02
dx
sen
1
2
Analogamente
c
2012
Mynor Chac
on Daz
166
i
sen x dx
[1 sen x] dx

MA1003
8.1
ZZ
cos(x + y) dx dy =
cos(x + y) dy dx
R3
3 2 x
1
2
Como se aprecia en la siguiente imagen:
Por u
ltimo para calcular la integral sobre R2 se debe partir esta en dos
pedazos:
De manera que se obtiene:

167
c
2012
Mynor Chacon Daz
8.1
MA1003
ZZ
ZZ
ZZ
cos(x + y) dx dy
cos(x + y) dx dy +
cos(x + y) dx dy =
R2
/2
cos(x + y) dy dx +
x
2
/2
3 2 x

+1 +
+1
2
2
=+2
Por lo tanto se tiene:
ZZ
|cos(x + y)| dx dy =

1 + ( + 2) +
1 = 2
2
2
Ejemplo 3
Z
Considere la integral
0
2xx2
84x
f (x, y) dydx
1. Dibuje la region de integacion.
2. Reescriba la integral con el orden de integracion invertido.
Soluci
on
La region de integracion es la siguiente:
c
2012
Mynor Chac
on Daz
168
cos(x + y) dy dx
MA1003
8.1
Para cambiar el orden de integracion se deben distinguie tres regiones como

se ve en la figura siguiente:
169
c
2012
Mynor Chacon Daz
8.1
MA1003
Note que
p
y = 2x x2 x2 + y 2 = 2x (x 1)2 = 1 y 2 x = 1 1 y 2
8 y2
y = 8 4x y 2 = 8 4x x =
4
Por lo tanto la integral en el nuevo orden de integracion se convierte en:
2 2
8y 2
4
Z
f (x, y)dxdy+
1y 2
Z
f (x, y)dxdy+
8y 2
4
1+
f (x, y)dxdy
1y 2
Ejemplo 4
Z
2ex 1
f (x, y) dy dx +
Dada I =
ln(2)
f (x, y) dy dx. Dibuje

0
2x3
x3
la region R de integracion y cambie el orden.

Soluci
on
La region R de integracion es la siguiente:
Cada sub-region corresponde a una de las dos integrales.

Para cambiar el orden de integracion se deben ver dos sub-regiones las cuales
se aprecian en la imagen siguiente:
c
2012
Mynor Chac
on Daz
170
MA1003
8.1
Las integrales correspondientes son:

Z
1
2
3y3
y2
f (x, y) dx dy +
I=
0
ln
y+1
2
f (x, y) dx dy
1
2
ln
y+1
2
Ejemplo 5
Cambiar el orden de integracion en
10
Z
0
20xx2
Z
f (x, y)dydx+
x(20x)
20
18
10
4x
x6
Z
f (x, y)dydx+
x(20x)
20
20
18
20xx2
f (x, y)dydx
x(20x)
20
Debe presentar el grafico de la region generado con GenGCF o Winplot.

Debe presentar todos los calculos que fueron necesarios para determinar
los lmites en el nuevo orden de integracion.
Soluci
on:
x(20 x)
, es una parabola que corta el eje x = 0
20
y en x = 20 y tiene vertice cuando x = 10 en el punto (10, 5).
Observe que la curva y =
Luego si y = 20x x2 x2 + y 2 = 20x (x 10)2 + y 2 = 100 lo cual

representa un crculo de radio r = 10.
La primera impresion del dibujo es la siguiente
171
c
2012
Mynor Chacon Daz
8.1
MA1003
La interseccion entre la parabola y la funcion f (x) =

resolver el sistema
4x
se obtiene al
x6
y = 20x x2
4x
y=
x6
Igualando ambas ecuaciones y elevando al cuadrado se obtiene 16x2 = (x
6)(20x x2 ), resolviendo tal ecuacion se obtienen los valores x = 0, 4, 10, 18,
4(4)
= 8 y por el
para x = 4 no existe solucion pues por un lado y =
46
p
otro y = 20(4) 42 = 8, de manera que los puntos de interseccion son
(0, 0), (10, 10), (18, 6). En este caso tales puntos pueden ser obtenidos con
claridad de la grafica generada por GenGCF.
De nuevo utilizando las intersecciones y el GenGCF se obtiene:
Para verificar que la regiòn es la correcta note que las tres integrales originales producen la siguiente particion de la region entre las curvas:
c
2012
Mynor Chac
on Daz
172
MA1003
8.1
Para cambiar la integracion al orden dxdy, se debe cambiar la forma de

apreciar el dibujo , ahora se debe de ver de la siguiente manera:
173
c
2012
Mynor Chacon Daz
8.1
MA1003
De aqu se nota que la variacion de y sera en los intervalos de 0 a 5, de 5 a 6 y

de 6 a 10. Para los lmites de integracion con respecto a x , se debe despejar
la variable x en funcion de la variable y, para las tres funciones involucradas.
yp
= 20x x2 (x 10)2 + y 2 = 100 (x 10)2 = 100 y 2 x =
10 100 y 2
x(20 x)
20y = 20x x2 20y = x2 20x + 100 100
20
p
x = 10 100 20y
y=
4x
6y
xy 6y = 4x xy 4x = 6y x =
x6
y4
y=
Luego de un breve analisis se concluye que hay que partir la region en cuatro
sub-regiones como se aprecia en las siguientes imagenes:
De manera que las tres integrales originales se convierten en cuatro que son:
10 10020y
10
f (x, y) dxdy +
100y 2
10+
10
10+
100y 2
10+ 10020y
100y 2
10
f (x, y) dxdy+
6y
y4
f (x, y) dxdy +
100y 2
c
2012
Mynor Chac
on Daz
174
10
f (x, y) dxdy
100y 2
MA1003
8.1
Ejemplo 6
Dada la expresion:
Z
I=
0
1x2
1+4xx2
f (x, y) dydx +
1
x1
1+4xx2
Z
+
3
2+ 5
f (x, y) dydx
1+4xx2
1+4xx2
f (x, y) dydx
Elabore una grafica con Winplot o el GenGCF que muestre la region de

integracion sombreada. Escriba I seg
un el orden de integracion dxdy.
Soluci
on:
De acuerdo con los lmites de integracion se debe considerar las siguientes
funciones para la graficacion de la region de integracion. De la ecuacion
y = 1 + 4x x2 se obtiene la igualdad (x2)2 +y 2 = 5, de y = 1 x2

se obtiene la igualdad x2 +y 2 = 1, por u
ltimo de se tiene la funcion y = x1.
Para determinar los puntos de interseccion de las dos circunferencias se sustituye x2 + y 2 = 1 en la ecuacion x2 + y 2 4x = 1 se obtiene 1 4x = 1 o
sea x = 0, esto significa que y 2 = 1 y = 1; por lo tanto los puntos de
interseccion son (0, 1). Analogamente los puntos de interseccion entre la
recta y = x 1 y la circunferencia (x 2)2 + y 2 = 5 son los (3, 2).
En la siguiente imagen aparece la region dividida en tres subregiones coloreadas de distintos colores, cada una de ellas corresponde a cada una de las
tres integrales iniciales del enunciado.
175
c
2012
Mynor Chacon Daz
8.1
MA1003
Para invertir el orden debemos ver la region como se aprecia en la siguiente

imagen:
Notepque al despejar x en la ecuacion (x 2)2 + y 2 = 5 p

se obtiene x =
2
2
2
2 5 y , analogamente de x + y = 1 se obtiene x = 1 y 2 , por lo
tanto con el cambio de orden quedan las siguientes integrales:
Z
I=
2+
5y 2
f (x, y) dxdy +
5y 2
2+
5y 2
f (x, y) dxdy+
1y 2
2+
5y 2
f (x, y) dxdy
0
c
2012
Mynor Chac
on Daz
176
y+1
MA1003
8.1
Ejemplo 7
Encuentre el volumen del solido V que se encuentra sobre el plano XY, bajo
el plano z = x e interior al cilindro x2 + y 2 = 4.
Soluci
on:
La siguiente figura muestra la mitad del solido
De manera que la mitad del volumen es esta dada por:
V
=
2
4x2
x dy dx
0
=
0
1
=
2
x
2
4x2
dx
0
4 x2 dx

2
1
x3
=
4x
2
3 0

2
1
8
=
8
2
3 0
8
3
El volumen total del solido es V =
16
.
3
177
c
2012
Mynor Chacon Daz
8.1
MA1003
8.1.3.
Integral doble en coordenadas polares
En una seccion mas adelante se estudiara el cambio de variables en general, en esta seccion, se tratara el caso de una integral doble. Considere la
siguiente region R
La teora establece que:

Z bZ
ZZ
f ()
f (r cos , r sen ) |J (r, )| dr d
f (x, y) dx dy =
a
g()
La expresion J (r, ) se llama Jacobiano y se calcula de la manera siguiente:
J (r, ) =
x
r
x
y
r
y
cos
sen
r sen r cos
!
= r cos2 + r sen2 = r
Al utilizar coordenadas polares solo interesa la variacion de y r ya no de

x e y.
Ejemplo 8
Z
2a
2axx2

x2 + y 2 dy dx
0
1. Dibuje la region de integracion.

2. Escriba la integral en coordenadas polares y calc
ulela.
3. Escriba la integral en coordenadas polares generalizadas y calc
ulela.
c
2012
Mynor Chac
on Daz
178
MA1003
8.1
Soluci
on:
1. En este caso
0 x 2a
0 y 2ax x2
por lo tanto la region de integracion es la siguiente:
2. El valor de r va de 0 a la curva, la curva en polares se escribe:
y=
2ax x2
x2 + y 2 = 2ax
r2 = 2ar cos
r = 2a cos
179
c
2012
Mynor Chacon Daz
8.1
MA1003
Entonces:
2a
2axx2
2
x +y
0
2a cos
r2 r dr d
dy dx =
0
2a cos
r3 dr d
=
0
=
0
= 4a
r4
4
2a cos
d
0
cos4 d
= 4a
=a
1 + cos 2
2
2
d

1 + 2 cos 2 + cos2 2 d
"
= a4
=a
= a4
+ {sen sen 0} +
2
1
+
2 2
h
2
Z
0
1 + cos 4
2

1
+ (sen 2 sen 0)
2 4
i
4
3
= a4
4
3. Puesto que la curva se puede escribir en la forma (x a)2 + y 2 = a2 ,

esto da pie a nombrar x a = r cos ; y = r sen , es decir:
c
2012
Mynor Chac
on Daz
180
#
MA1003
8.1
x = a + r cos
y = r sen
Este cambio de variable tiene el efecto geometrico de trasladar el crculo

de manera que quede centrado en el origen (0,0), de manera que va
a variar de 0 a y el radio varia de 0 a a. Note que el Jacobiano queda
igual que en las coordenadas polares |J (r, )| = r. De manera que
2a
2axx2
2
x +y
0

(a + r cos )2 + (r sen )2 r dr d
dy dx =
=
0

r3 + 2ar2 cos + ra2 dr d
3
= a4
4
Las dos imagenes siguientes muesran el calculo realizado con la TI y

con la MuPAD, respectivamente:
181
c
2012
Mynor Chacon Daz
8.1
MA1003
Ejemplo 9
Determine el volumen del solido en el interior de la esfera x2 + y 2 + z 2 = 4
y (x 1)2 + y 2 = 1.
1. Use coordenadas rectangulares y la TI-89..
2. Use coordenadas polares generalizadas y la TI-89.
3. Use coordenadas polares y calcule manualmente.
Soluci
on:
El solido se aprecia en la siguiente imagen:
Debido a la simetra del solido solo se calculara la cuarta parte:
La region R es la siguiente:
c
2012
Mynor Chac
on Daz
182
MA1003
8.1
Por lo tanto se tiene:

V =4
Z Z p
4 x2 y 2 dA
1. La TI no calcula algebraicamente la integral, sin embargo recuerde que

si se digita
antes de precionar ENTER, se obtiene la aproximacion
numerica correspondiente.
V =4
Z Z p
Z
4
x2
y 2 dA
2xx2
=4
4 x2 y 2 9,6440970803
2. Dado que la region R tiene ecuacion (x1)2 +y 2 = 1, entonces usaremos

las coordenadas polares generalizadas
(
x = 1 + r cos
y = r sen
Al usar las coordenadas generalizadas se debe entender que la region

R queda centrada en el origen, por lo que se tiene
0 , 0 r 1
de manera que
183
c
2012
Mynor Chacon Daz
8.1
MA1003
Z Z p
V =4
4 x2 y 2 dA
R
=4
0
4 (1 + r cos )2 (r sen )2 drd
=4
0
r 3 r2 2r cos drd 9,6440970803
La u
ltima integral no se puede calcular facilmente de forma manual
por lo que se obtuvo la aproximacion utilizando la TI:
3. Al usar las coordenadas polares:

(
x = r cos
y = r sen
La ecuacion (x 1)2 + y 2 = 1 se transforma en r = 2 cos y la region

R no sufre ninguna alteracion, por lo tanto
0
,
2
entonces:
c
2012
Mynor Chac
on Daz
184
0 r 2 cos
MA1003
8.1
V =4
Z Z p
4 x2 y 2 dA
2 cos
=4
0
= 2
= 2
32
3
u du d
u du d
u2
i4 sen2
32
=
3
4 sen2
22
=
3
44 cos2
haciendo u = 4 r2
4 r2 r dr d
4
=
3

8 sen3 8 d
3

sen 1 d
"Z
sen3 d
d
0

32 2
9,6440970803
=
3 3 2
185
c
2012
Mynor Chacon Daz
8.1
MA1003
Las coodenadas polares son u

tiles cuando la region de integracion incluye el
2
2
2
crculo x + y = a , note que en tal caso x = r cos , y = r sen por lo tanto
x2 + y 2 = r2 , de manera que r vara entre 0 y a. Si la region de integracion
x2 y 2
incluye la elipse 2 + 2 = 1 resulta natural modificar un poco las polares
a
b
x2 y 2
y escribir x = ar cos , y = br sen , de esta manera se cumple 2 + 2 = r2
a
b
por lo que r de nuevo vara entre 0 y 1. Estas coordenadas se conocen como
coordenadas el
pticas.
En tal caso el Jacobiano J (r, ) queda de la manera siguiente:
J (r, ) =
x
r
x
y
r
y
a cos
b sen
ar sen br cos
!
= abr
La integral se trasforma de la siguiente manera:
Z bZ
ZZ
f ()
f (ar cos , br sen ) abr dr d
f (x, y) dx dy =
R
g()
Ejemplo 10
Utilizar las coordenadas elpticas para determinar el volumen del cuerpo

limitado por z = 5x2 + y 2 y 12 z = x2 + y 2 .
Soluci
on:
Se tiene z = 5x2 + y 2 y z = 12 x2 y 2 , al igualar ambas ecuaciones,

x2 y 2
se obtiene
+
= 1, esta elipse representa la proyeccion (R) en el plano
2
6
XY del volumen.
c
2012
Mynor Chac
on Daz
186
MA1003
8.2
De manera que la integral de volumen es:
Z Z

12 x2 y 2 5x2 + y 2 dx dy
v =
R
Z Z

12 6x2 2y 2 dx dy
=
R
Se utilizaran coordenadas elpticas, x = 2r cos y y = 6r sen , como se

sabe el jacobiano es J = 2 6r = 12r, entonces:
Z Z

12 6x2 2y 2 dx dy
v =
R
12 1 r2
= 4
0
12r dr d
Z
= 48 12
0

r r3 dr d
= 12 3
187
c
2012
Mynor Chacon Daz
8.2
MA1003
8.2.
Integral triple
Recuerde que una funcion continua f : I R R, I = [a, b], permite

definir
Z
f (x) dx
a
si f : R R2 R es continua sobre R = [a, b] [c, d], o R es una region

cerrada y acotada en el plano, se puede entonces definir
ZZ
f (x, y) dxdy
R
Una region plana R es una region acotada si es posible incluirla dentro de un

disco de radio fijo y es cerrada si R contiene todos los puntos de su frontera.
De forma analoga, dada una funcion continua f : B R3 R, siendo
B un solido cerrado y acotado en el espacio, como por ejemplo: una bola
solida, un cuerpo solido hecho con plasticina, un cubo solido. Entonces se
puede definir
ZZZ
f (x, y, z) dxdydz
B
Un solido B es acotado si es posible incluirlo dentro de una esfera de radio

fijo y es cerrada si B contiene todos los puntos de su frontera o contorno, en
tal caso la frontera esta constituida por una o varias superficies.
Suponga que la base de un solido B esta representada por la region R en el

plano XY y el techo es el pedazo de superficie z = g(x, y) 0 encima de R,
entonces:
ZZZ
ZZZ
f (x, y, z)dzdydx =
ZZ Z
dzdydx =
g(x,y)
ZZ
dzdydx =
g(x, y)dA
R
Es decir en tal caso la integral triple de 1 da como resultado el volumen entre

z = g(x, y) y el plano XY, o sea el volumen del solido B; V (B). En general
si la superficie z = g(x, y) 0 y R es la proyeccion de la superficie en XY,
entonces la integral triple de la funcion 1 da como resultado el volumen bajo
z = g(x, y):
ZZZ
ZZ
dzdydx =
g(x, y) dA = volumen bajo g(x, y)
B
c
2012
Mynor Chac
on Daz
188
MA1003
8.2
Los Teorema de Fubinni siguen siendo validos y el calculo algebraico es

analogo al utilizado en integrales dobles, el siguiente ejemplo trata sobre el
cambio de orden de integracion en una integral triple, luego se expondra el
calculo de una integral triple en coordenadas rectangulares y en el tercer
ejemplo se obtendra el volumen de un solido utilizando la intregrales triples.
Ejemplo 11
Cambie el orden de integracion de
Z
4 2x
3
3 x2 3y
4
f (x, y, z) dz dy dx
0
a una integral con el orden dy dx dz

Soluci
on:
Observando el orden original de integracion dz dy dx y haciendo la correspondencia con los lmites de integracion se deduce que
0x6
0y 4
x 3y
0z 3 ,
2
4
De nuevo el orden dz
2x
3
x 3y
z vara del piso al plano z = 3
2
4
dy dx , indica que la proyeccion debe realizarse en

| {z }
proyeccion
el plano XY, esto es, en el proceso de calculo, las dos u

ltimas integrales
se hacen con respecto a Y y a X, seg
un el orden establecido en la integral
original. Las dos primeras desigualdades determinan la region de integracion
en el plano XY, es decir la proyeccion en XY:
189
c
2012
Mynor Chacon Daz
8.2
MA1003
De manera que la integral con el nuevo orede de integracion es:

Z
62z
4 2x
3y
3
4
f (x, y, z) dz dx dz
0
Donde f (x, y, z) vara sobre el solido siguiente:
Ejemplo 12
ZZZ
Calcule
1
3 dz dy dx, siendo B el tetraedro limitado por
B (1 + x + y + z)
los tres planos coordenados y el plano x + y + z = 1
c
2012
Mynor Chac
on Daz
190
MA1003
8.2
Soluci
on:
El solido es similar al del ejemplo anterior:
La proyeccion en el plano XY, es la region triangular siguiente:
El calculo de la integral es el siguiente:

191
c
2012
Mynor Chacon Daz
8.2
MA1003
ZZZ
B
1
dz dy dx =
(1 + x + y + z)3
=
=
=
1x
1xy
1
dz dy dx
(1 + x + y + z)3
0
0
0
1xy
Z 1 Z 1x
1
dy dx
2 (1 + x + y + z)2 0
0
0

Z Z
1
1 1 1x 1
dy dx
2 0 0
4 (1 + x + y)2

1x )
Z (
1 1 1x
1
+
dx
2 0
4
1+x+y 0

Z
1
1 1 1x 1
+
dx
2 0
4
2 1+x
#
"
1
1 1
x2
1
x
+ ln(2)
2 4
2 0 2

1 1 1
+ ln(2)
2 8 2
5
1
ln(2)
2
16
Ejemplo 13
Calcule el volumen del solido en el primer octante limitado por la superficies
z = 1 x2 , y = 2x, x = 3.
Soluci
on:
El solido se puede observar en la siguiente imagen:
La proyeccion en el plano XY, es la region R siguiente:

c
2012
Mynor Chac
on Daz
192
MA1003
8.2
La manera mas facil para realizar la integracion, es elegir la variacion

y
x3
2
El planteo y la integracion es la siguiente:
0 y 1;
1y 2

1 y 2 dx dy
dz dx dy =
0
y
2
y
2
3
1 y 2 x y dy
2
y
dy
2

y y3
2
3 3y +
dy
2
2
1 y2
=
1
=
0

1

y2 y4
3
+
= 3y y
4
8 0

1
1
=
1 +3
8
4
15
8
8.2.1.
Coordenadas cilndricas
Consideremos un punto (x, y, z), si las dos primeras coordenadas son

determinadas usando polares, mientras que z sigue variando de arriba hacia abajo en su eje respectivo, entonces se obtienen lo que se conoce como
Coordenadas cilndricas:
193
c
2012
Mynor Chacon Daz
8.2
MA1003
x = r cos
y = r sen
z=z
En tal caso
J (r, , z) =
x
r
x

x
z
y
r
y
y
z
z
r
z

z
z
cos
sen 0
= r sen r cos 0 = r
0
0
1
de manera que, la integral sobre alg

un solido B se podra trasformar como
sigue:
ZZZ
f2 ()
g2 (r,)
f (x, y, z) dx dy dz =
B
f (r cos , r sen , z)r dz dr d

1
f1 ()
g1 (r,)
Ejemplo 14
Haga uso de las coordenadas cilndricas para encontrar el volumen del cuerpo
limitado por z = 2x2 + y 2 , z = 12 x2 2y 2 .
Soluci
on:
c
2012
Mynor Chac
on Daz
194
MA1003
8.2
El volumen a calcular es el del huevo limitado por los dos paraboliodes,

como se aprecia en la imagen siguiente:
Igualando las ecuaciones de ambas superificies se obtiene x2 + y 2 = 4, esto

significa que la proyeccion en XY es el disco x2 + y 2 = 4, por lo que se
utilizaran coordenadas cilndricas, en tal caso se tiene:
z = 2x2 + y 2 z = r2 1 + cos2
z = 12 x2 2y 2 z = 12 r2 1 + sen2
Entonces el vol
umen esta dado por:
12r 2 (1+sen2 )
v =
r dz dr d
0
r 2 (1+cos2 )
0
2
=
0

12r 3r3 dr d
= 24
Ejemplo 15
Utilice coordenadas cilndricas para calcular el volumen del solido limitado
por el paraboloide 4z = x2 + y 2 y el hiperboloide x2 + y 2 z 2 = 4, z > 0.
Debe dibujar la proyeccion en el plano XY y calcular la integral triple de
forma manual.
(Sugerencia: utilizar el GenGCF para visualizar el solido).
195
c
2012
Mynor Chacon Daz
8.2
MA1003
Soluci
on:
Para determinar la interseccion de las superficies:
x2 + y 2 = 4z
4zz 2 = 4 z 2 4z+4 = 0 (z2)2 = 0 z = 2
x2 + y 2 z 2 = 4
Es decir
lassuperficies se intersecan a la altura de z = 2 en el crculo

2
x2 +y 2 = 2 2 . Por otro lado la interseccion del hiperboloide x2 +y 2 z 2 =
4, z > 0 con el plano z = 0, es el crculo x2 + y 2 = 4, como se aprecia en la
siguiente imagen:
La proyeccion en el plano XY es la siguiente:
La proyeccion en plano YZ es:

c
2012
Mynor Chac
on Daz
196
MA1003
8.2
Usando coordenadas polares x = r cos , x = r sen , |J| = r, se tiene que

r2
4
2
2
2
2
2
x + y z = 4 r z = 4 z = r2 4
x2 + y 2 = 4z 4z = r2 z =
de manera que:
r2
4
V =4
2 2
r dz dr d + 4
0
2 2
r dr d + 4
0
16
=
+
2 4
2 2
r3 dr d 4
Z
0
r dz dr d
r2 4
r3
2
r r 4 dr d
4
2 2
r r2 4 dr d
"
#
Z
(2 2)4 16
4 1
= 2 +
2
u 2 du
2
4
4
2 0
= 2 + 6
r2
4

donde u = r2 4
16
3
8
=
3
197
c
2012
Mynor Chacon Daz
8.2
MA1003
Ejemplo 16
Sea T el solido en el primer octante limitado por los planos:
x = 0, y = 0, z = 0, z = x + y, x + y + z = 2
1. Elabore una grafica con Winplot o el GenGCF que muestre el solido

T, debe verse solamente el primer octante.
2. Elabore una grafica con Winplot o el GenGCF que muestre la proyeccion en el plano XY.
3. Calcule el volumen de T usando coordenadas cilndricas, de forma manual.
4. Calcule el volumen de T usando coordenadas cilndricas, con la TI-89
o con MuPAD, debe incluir las imagenes respectivas.
Soluci
on:
1. Sea T el solido en el primer octante limitado por los planos:
x = 0, y = 0, z = 0, z = x + y, x + y + z = 2
Para determinar la proyeccion en el plano XY de la recta de interseccion:
(
z =x+y
x+y+z =2
Se sustituye la primer ecuacion en la segunda, as z + z = 2 z = 1,

luego al sustituir en cualquiera de las dos se llega a la ecuacion x+y = 1.
La siguiente imagen muestra los planos que constituyen el solido, la
recta x + y = 1 proyeccion de la recta de interseccion de los dos planos
y la recta x+y = 2 que representa la interseccion del plano x+y+z = 2
con el plano XY.
c
2012
Mynor Chac
on Daz
198
MA1003
8.2
2. De forma que la proyeccion sobre el plano XY es:
En coordenadas cilndricas se tiene
x = r cos
y = r sen
z=z
Entonces:
z = x + y z = r [cos + sen ]
x + y + z = 2 z = 2 r [cos + sen ]
1
x + y = 1 r [cos + sen ] = 1 r =
cos + sen
1
x+y =2r =
cos + sen
199
c
2012
Mynor Chacon Daz
8.2
MA1003
De manera que el volumen en coordenadas cilndricas es:
1
cos +sen
r[cos +sen ]
V =
0
Z
rdzdrd+
2
cos +sen
2r[cos +sen ]
rdzdrd
1
cos +sen
3. Seguidamente se calculara cada una de las integrales por separado:
1
cos +sen
r[cos +sen ]
I1 =
rdz dr d
0
0
1
cos +sen
=
0
r2 [cos + sen ] dr d
1
=
3
1
=
3
1
=
3
1
=
3
1
=
3
1
=
3
1
[cos + sen ] d
[cos + sen ]3
1
d
[cos + sen ]2
1
d
1 + sen 2
1
1 sen 2
d
1 + sen 2 1 sen 2
1 sen 2
d
cos2 2
sec2 2 sec 2 tan 2 d

2
1 1
1
=
tan 2 sec 2
3 2
2
0

1 1
1
=
[0] [1 1]
3 2
2
|
{z
}
1
1
3
Analogamente:
c
2012
Mynor Chac
on Daz
200
MA1003
8.2
2
cos +sen
2r[cos +sen ]
I2 =
rdz dr d
1
cos +sen
2
cos +sen
2

cos +sen
3
r
r2 [cos + sen ]
d
3
1
=
0
2
3
cos +sen
Z
0
|
=
2r r2 [cos + sen ] dr d
1
cos +sen
1
d
[cos + sen ]2
{z
}
1
2
3
Por lo tanto
Z
V
1
cos +sen
r[cos +sen ]
2
cos +sen
2r[cos +sen ]
rdz dr d
rdz dr d +
0
1
cos +sen
1 2
+
3 3
= 1
4. Las siguientes imagenes corresponden a las integrales I1 e I2 , hechas
con la TI-89 y MuPAD respectivamente:
201
c
2012
Mynor Chacon Daz
8.2
8.2.2.
MA1003
Coordenadas esf
ericas
Recuerdese que el vector que va de (0, 0, 0) a (x, y, z) esta representado por

el punto P = (x, y, z), la proyeccion de v sobre el plano XY es el vector
u = (x, y). El vector v forma un angulo con el eje vertical Z + y u forma
un angulo con el eje X + . Si se denota r = kvk entonces el punto P queda
determinado por
x = r cos sen
y = r sen sen
z = r cos
Para cubrir cualquier posicion de P en el espacio tridimensional se debe
considerar
0 < 2;
0 < ;
0r
Tales coordenadas se llaman Coordenas esfericas, notese que en este caso se

cumplen las siguientes identidades:
x2 + y 2 + z 2 = r 2 ;
c
2012
Mynor Chac
on Daz
tan =
202
y
;
x
cos =
z
r
MA1003
8.2
De modo que cualquier punto P en la esfera x2 + y 2 + z 2 r2 tiene coordenadas (r, , ). El valor absoluto del Jacobiano respectivo es:
|J(r, , )| = r2 sen
De forma que la integral f (x, y, z) sobre un cuerpo B se trasforma de la
siguiente forma:
ZZZ
ZZZ
f (r cos sen , r sen sen , r cos )r2 sen drdd
f (x, y, z)dxdydz =
B
Siendo B el solido determinado en el espacio tridimensional con ejes r, , .

Ejemplo 17
Utilice coordenadas esfericas para calcular la siguiente integral:
Z
4x2
Z 4x2 y2
x
x2 + y 2 + z 2 dz dy dx
Soluci
on:
Seg
un la integral se tiene
0 x 2;
0y
4 x2 ;
0z
4 x2 y 2
Se deduce de aqu que el cuerpo sobre el cual se esta integrando es la octava

parte de la esfera x2 + y 2 + z 2 = 4, ubicada en el primer octante:
La integral se establece y calcula como sigue:

203
c
2012
Mynor Chacon Daz
8.2
MA1003
0
2
=
0
r cos sen
|
{z
}
I =
r2 r2 sen d dr d
| {z }
J
r4 cos sen2 d dr d
5
como se aprecia en la imagen siguiente obtenida al realizar el calculo con la
TI:
=
Ejemplo 18
Haga uso de las coordenadas esfericas para encontrar el volumen del cuerpo
z2
= x2 + y 2 , x2 + y 2 + (z 1)2 = 1. (El calculo siguiente se
limitado por
3
refiere al cono de helado, limitado superiormente por la esfera e inferiormente
por el cono)
Soluci
on:
El cuerpo limitado se aprecia en la siguiente imagen:
c
2012
Mynor Chac
on Daz
204
MA1003
8.2
Al sustitur la primer ecuacion en la segunda se obtiene 2z(2z 3) = 0

3
3
z = 0 z = , esto significa se intersecan a la altura de z = en el crculo
2 !
2
2
3
x2 + y 2 =
(Proyeccion en el XY)
2
De la figura anterior se tiene tan = , lo cual implica = . Utilizando

6
3
coordenadas esfericas
x = r cos sen
y = r sen sen
z = r cos
J = r2 sen
Se tiene que
x2 +y 2 +(z1)2 = 1 x2 +y 2 +z 2 2z = 0 r2 2r cos = 0 r = 2 cos
Entonces el vol
umen esta dado por:
2 cos
r2 sen dr d d
v = 4
0
32
=
3
Z
0
sen cos3 d d
Z 3

2
32
u3 du
3
2
1
=
(u = cos )
7
12
205
c
2012
Mynor Chacon Daz
8.3
MA1003
El calculo que se refiere al volumen limitado superiormente por la cono e

inferiormente por la esfera el volumen es el siguiente:
2 cos
r2 sen dr d d
v = 4
32
=
3
=
8.3.
Z
0
sen cos3 d d
3
4
Cambio de Variables m
as generales
Una funcion sobreyectiva f : A B se conoce como aplicaci

on sobre B,
esto significa f (A) = B. Si la aplicacion es de A en A, entonces se le llama
Transformaci
on. Ademas recuerde que una funcion definida sobre A es uno
a uno(inyectiva) si a, b A se cumple f (a) = f (b) a = b.
En particular si A R2 y r es una funcion vectorial:
r : A R2 |{z}
| {z }
(u,v)
R2
|{z}
7 r(u,v)=(x,y)
se dice que r es una Transformaci

on o cambio de coordenadas en T
si se cumplen las siguientes condiciones:
X(u, v), Y (u, v) tienen derivadas parciales de primer orden continuas,
siendo r(u, v) = (X(u, v), Y (u, v)) con x = X(u, v), y = Y (u, v).
r es uno a uno, excepto quizas en un conjunto finito de puntos.
!
Y
X
(X, Y )
u
u
J(u, v) = X
=
6= 0, (u, v) A
Y
(u, v)
v
v
Por ejemplo en el caso de las coordenadas polares, la transformacion es
(r, ) = (r cos , r sen ), donde r > 0, 0 < 2 y

(X, Y )
=r
|J(r, )| =
(r, )
De nuevo considerese la transformacion r(u, v) = (X(u, v), Y (u, v)). la cual
transforma una region T en el plano UV en otra region S en el plano XY
c
2012
Mynor Chac
on Daz
206
MA1003
8.3
Si las ecuaciones x = X(u, v), y = Y (u, v) se logran despejar en la forma

u = U (x, y), v = V (x, y) se obtiene la transformacion inversa r1 (x, y) =
(U (x, y), V (x, y)), esto no siempre es facil algebraicamente, sin embargo recuerde que r1 va a existir gracias a la inyectividad de r.
Si consideramos v como constante se obtiene r(u, v) = g(u), esta representa
r
una curva en el plano XY cuya pendiente esta dada por
u
Ahora considere un peque

no rectangulo R T de area A(R) = uv y su
aplicacion R S
207
c
2012
Mynor Chacon Daz
8.3
MA1003
La imagen de R es el paralelogramo curvilneo R , cuya area es aproximadamente

r

r
r
r

u u v v = u v uv
Pero observe que:
Y
u
Y
v
r
r X
= u
u v
X
0 = (0, 0, J(u, v))

0
Entonces

r
r

u v = k(0, 0, J(u, v))k = |J(u, v)|
En conclusion el area de R es aproximadamente |J(u, v)| u v; lo que
significa
a (R ) = a (R) |J(u, v)|
Si J = 1 en T entonces R y R tienen la misma area, es decir la aplicacion
conserva areas. De lo contrario se debe multiplicar el aerea de R por |J|,
siendo entonces |J| un factor de ampliaci
on de
areas.
Si T se particiona en un conjunto de subrectangulos R con correspondientes
aplicaciones R en la region S, entonces
X
X
a R =
a (R ) |J(u, v)|
{z
a(S)
En terminos de integrales esta ecuacion significa:

ZZ
ZZ
dx dy =
|J(u, v)| du dv
S
T
| {z }
a(S)
Al introducir el integrando f (x, y), que varia sobre S, se obtiene la f

ormula
general para cambio de variable en una integral doble
ZZ
ZZ
f (X(u, v), Y (u, v)) |J(u, v)| du dv
f (x, y) dx dy =
S
c
2012
Mynor Chac
on Daz
208
MA1003
8.3
En la practica la integral propuesta es la integral sobre la region del plano XY

y mediante una transformacion debemos trasladarnos al plano UV, donde
se efectuara el calculo de la nueva integral, esto significa que se utilizara la
aplicacion inversa r1 (x, y) = (U (x, y), V (x, y)), es decir tendremos u y v en
funcion de x e y:
u = U (x, y)
v = V (x, y)
(X, Y )
, pues esto
(u, v)
requiere tener x e y en terminos de u y v. Para resolver esto podemos tratar
de despejar x e y, pero esto usualmente no es facil, la otra forma es tener en
cuenta la siguiente formula:
Esto presenta un problema al quere calcular J(u, v) =

1
(U, V )
(X, Y )
=
J(u, v) =
(u, v)
(x, y)
8.3.1.
Trasformaciones con Winplot
Winplot ofrece una aplicacion que permite realizar trasformaciones de R2

en R2 . Como se menciono anteriormente, Winplot utiliza directamente la
transformacion inversa r1 (x, y) = (U (x, y), V (x, y)), consecuentemente el
|
{z
}
(u,v)
sistema de ejes de la izquierda representa el plano XY (dominio1.ma2 en

Winplot) y el sistema de la derecha representa el plano UV (rango en Winplot). Para llegar a esta aplicacion de Winplot se sigue la ruta Ventana
>Mapping >Plano xy, luego se a
naden los ejes presionando Crtl + E en
cada uno de los planos, por u
ltimo se introduce la transformacion inversa siguiendo la ruta Map >Nuevo, de oficio Winplot muestra la siguiente la
transformacion
(
u = U (x, y) = 2x y
v = V (x, y) = x + y

209
c
2012
Mynor Chacon Daz
8.3
MA1003
al darle ok y acercar (zoom) la pantalla de la derecha se aprecia la imagen

siguiente:
La siguiente imagen muestra la transformacion que se produce en un peque

no
segmento:
c
2012
Mynor Chac
on Daz
210
MA1003
8.3
Ejemplo 19
x2 + y 2
dy dx siendo R el triangulo en el plano
2 2
R 1x y
XY con vertices (0, 0); (1, 0); (1, 1). Haga el cambio de variable u = 2xy, v =
x2 y 2 .
Z Z
1. Manualmente obtenga la region en el plano UV, en la que se transforma

la region R del plano XY. Deben aparecer todos los calculos necesarios.
2. Utilize Winplot para verificar que la region hallada en el item anterior

es correcta. Presente la imagen correspondiente generada por Winplot.
3. Calcule el Jacobiano J =
(X, Y )
.
(u, v)
4. Calcule la integral.
Soluci
on:
1. La region R tiene como borde las rectas y = x, y = 0, x = 1, como se

aprecia en la siguiente figura:
211
c
2012
Mynor Chacon Daz
8.3
MA1003
Primero note que u = 2xy, v = x2 y 2 , entonces xy 0 u 0.

Luego
y = x y x = 0 v = x2 y 2 = 0
y=0u=0
x = 1 u = 2y, v = 1 y 2 v = 1
u 2
2
4 u2
4
Por lo tanto la region en el plano UV, esta limitada por las curvas
u 0, v = 0, u = 0, v =
2. La imagen en Winplot es la siguiente:
c
2012
Mynor Chac
on Daz
212
4 u2
4
MA1003
8.3
1
(U, V )
. Pero
(x, y)
!
(U, V )
2y 2x
= det
= 4y 2 4x2
(x, y)
2x 2y

(X, Y )
=
3. El Jacobiano J =
(u, v)
Por lo tanto |J| =
1
4 (x2 + y 2 )
4. La integral
x2 + y 2
dy dx =
1 x2 y 2
Z Z
R
Z
0
4u2
4
4u2
4
=
0
1
=
4
=
Ejemplo 20
Z Z
Sea I =
1
x2 + y 2
dv du
2
2
1 (u/2) 4 (x + y 2 )
1
dv du
4 u2
du
0
1
2
(1 + 2x) sen(x + y) dA, donde la region S es el triangulo
limitado por y = x, x + y = 2, y = 0. Calcule I mediante la transformacion

213
c
2012
Mynor Chacon Daz
8.3
MA1003
u = x + y, v = y x2 . Use Winplot para representar la region S en el plano

XY, la region T en el plano UV. Ademas use Mupad o la TI para obtener
el resultado de la integracion.Debe incluir las imagenes correspondientes de
Winplot y MuPAD.
Soluci
on:
La region S tiene como borde las rectas y = x, x + y = 2, y = 0, como se
aprecia en la siguiente figura:
Primero note que u = x + y, v = y x2 , entonces
y = x u = 2x v = x x2 v =
2u u2
=0
4
y = 0 u = x v = x2 v = u2
x+y =2u=2
Las regiones en el plano XY y en el UV, son:
(X, Y )
El Jacobiano J =
=
(u, v)
(U, V )
(x, y)
1
(U, V )
1 2x
= det
(x, y)
1
1
Por lo tanto |J| =
1
1 + 2x
Entonces
c
2012
Mynor Chac
on Daz
214
. Pero
!
= 1 + 2x
MA1003
8.3
Z Z
2uu2
4
(1 + 2x) sen(x + y) dA =
sen u dv du
0
u2
7 sen 2 5 cos 2 3
2
Ejemplo 21
Mediante el cambio de variable x = u2 v 2 , y = 2uv, u 0, v 0, calcule
Z Z
I=
R
1
p
x2 + y 2
dA
donde la region R es la region en el primer cuadrante limitada por las curvas

1
y 2 = 4 + 4x, y 2 = x + , y 2 = 4 4x, y 2 = 64 16x. Utilice Winplot para
4
verificar que la region transformada es la correcta
Soluci
on:
La imagen siguiente muestra la region R en el I cuadrante
215
c
2012
Mynor Chacon Daz
8.3
MA1003
Observe que no es facil despejar u, v en terminos de x e y, por lo tanto no

se puede utilizar la aplicacion Mapping de Winplot de forma directa.
Ahora cada una de las cuatro curvas frontera de la region R deben escribirse
en terminos de u,v, recordando que x = u2 v 2 , y = 2uv:
Primera curva
y 2 = 4 + 4x 4u2 v 2 = 4u2 4v 2 + 4
u2 v 2 1 = u2 v 2
u2 v 2 u2 + v 2 1 = 0

u2 v 2 1 + v 2 1 = 0

u2 + 1 v 2 1 = 0
v = 1
v = 1,
c
2012
Mynor Chac
on Daz
216
pues
u > 0, v > 0
MA1003
8.3
Segunda curva
y2 = x +
1
1
4u2 v 2 = u2 v 2 +
4
4
4u2 v 2 u2 + v 2
4u2 v 2 u2 +
1
=0
4
4v 2 1
=0
4
1

4v 2 1 = 0
u2 4v 2 1 +
4

1
2
u +
4v 2 1 = 0
4
v=
1
2
1
v= ,
2
pues
u > 0, v > 0
Tercer curva
y 2 = 4 4x 4u2 v 2 + 4u2 4v 2 4 = 0

4u2 v 2 + 1 4 v 2 + 1 = 0

4 u2 1 v 2 + 1
u = 1
u = 1,
pues
u > 0, v > 0
Cuarta curva
217
c
2012
Mynor Chacon Daz
8.3
MA1003
y 2 = 64 16x 4u2 v 2 + 16u2 16v 2 64 = 0

4u2 v 2 + 4 16 v 2 + 4 = 0

4 u2 4 v 2 + 4
u = 2
u = 2,
pues
u > 0, v > 0
entonces las cuatro curvas frontera de la region R se convierten en v = 1, v =

1
, u = 1, u = 2, las cuales forman un rectangulo en el plano UV:
2
Luego
(X, Y )
2u 2v
J=
= det
(u, v)
2v 2v
!
= 4 u2 + v 2
falta transformar el denominador del integrando:

x2 + y 2 = u4 2u2 v 2 + v 4 + 4u2 v 2 = u4 + 2u2 v 2 + v 4 = u2 + v 2
1
1
= 2
, de manera que:
por lo tanto p
u + v2
x2 + y 2
ZZ
R
p
dA =
x2 + y 2
c
2012
Mynor Chac
on Daz
Z
1
Z
1
2
1
4 dv du = 4(2 1) 1
2
218

=2
2
MA1003
8.4
Para verificar que la region transformada es la correcta, se cambia el papel

de las variables, es decir ingresamos en Winplot, la trasformacion
(
u = U (x, y) = x2 y 2
v = V (x, y) = 2xy
la trasformacion debera ir del rectangulo en UV a la region R en XY,en
efecto esto es lo que ocurre, como se aprecia en la siguiente figura:
se concluye que la transformacion realizada manualmente esta correcta.
8.4.
F
ormulas para momentos, masa y centro de masa
Seguidamente se presentan algunas de las formulas mas utilizadas en un

curso de calculo, como lo son las correspondientes la masa y el primer momento para placas delgadas que cubren una region R en el plano XY, siendo
(x, y) la densidad de la placa en cada punto (x,y).
Masa y Momentos en R2
Z Z
Masa=M =
(x, y)dx dy
R
219
c
2012
Mynor Chacon Daz
8.4
MA1003
Z Z
Primer momento con restecto al eje x = Mx =
y (x, y)dx dy
R
Z Z
Primer momento con restecto al eje y = My =
x (x, y)dx dy
R

Centro de masa = (x, y) =
My M x
,
M M
Momentos de Inercia
Z Z
Z Z
y (x, y)dx dy
Ix =
x2 (x, y)dx dy
Iy =
Masa y Momentos en R3 , siendo (x, y, z) la densidad en cada punto

(x,y,z) del solido
Z Z Z
(x, y, z) dx dy dz
Masa=M =
D
Z Z Z
z (x, y, z) dx dy dz
Mxy =
D
Z Z Z
y (x, y, z) dx dy dz
Mxz =
D
Z Z Z
Myz =
x (x, y, z) dx dy dz
D

Centro de masa = (x, y, , z) =
Myz Mxz Mxy

,
M M M
Momentos de Inercia
Z Z Z
Ix =
Z Z Z
(y +z )(x, y, z)dxdy
D
Iy =
(x2 +z 2 )(x, y, z)dx
Z Z Z
Iz =
(x2 + y 2 ) (x, y, z)dx dy
c
2012
Mynor Chac
on Daz
220
MA1003
8.4
Ejemplo 22
Una lamina tiene la forma limitada por la elipse
x2 y 2
+
= 1, 0 y 4, y = 0
4
16
en cada punto la densidad esta dada por (x, y) = |x| y. Calcule el centro de
masa de la lamina.
Soluci
on:
La imagen siguiente representa la lamina (R):
Puesto que (x, y) = (x, y) =, se tiene que la funcion de densidad es

simetricaZcon
Z respecto al eje Y, esto significa que el calculo de la masa,
|x| y dA, se puede hacer sobre la parte derecha de la region y
M =
R
multiplicar por dos, es decir:

ZZ
Z 21 y
|x| y dA = 2
M=
R
xy dx dy = 16
0
Puesto que la lamina es simetrica su centro de masa debe encontrarse sobre

su eje de simetra(eje Y), por lo tanto debe tener la forma (x, y) = (0, y), es
Mx
decir basta determinar y =
, luego
M
ZZ
Z 4 Z 21 y
2
512
Mx =
|x| y 2 dA = 2
xy 2 dx dy =
15
R
0
0
Por lo tanto
221
c
2012
Mynor Chacon Daz
8.4
MA1003
512
15
32
y=
=
centro de masa =
16
15
32
0,
15
La imagen siguiente muestra el centro de masa:
El valor de las integrales se obtuvo con MuPaD:
Ejemplo 23
Determine el volumen y el centroide del solido que yace dentro de la esfera
x2 + y 2 + z 2 = a2 y arriba del cono x2 + y 2 z 2 = 0, z 0. la densidad en
cada punto es constante (x, y, z) = 1.
Soluci
on:
De las ecuaciones x2 + y 2 + z 2 = a2 , x2 + y 2 z 2 = 0, z 0 se obtiene
c
2012
Mynor Chac
on Daz
222
MA1003
8.4
2z 2 = a2 z 2 =
a2
a
a
z = z =
2
2
2
Esto dice que la interseccion de la esfera y el cono ocurre a la altura de

a
a2
z = formando el crculo x2 + y 2 = , de aqu que la proyeccion en el
2
2
a2
plano XY es precisamente el disco x2 + y 2 . Como muestra la siguiente
2
imagen.
Debido a la simetra del solido se calculara el volumen en el I octante y

se multiplicara por cuatro, ademas se utilizaran coordenadas esfericas, el
volumen es el siguiente:
Z
V =4
0

a3 2 2
r sen dr d d =
3
2
La integral se calculo con el emulador de la TI-89:
223
c
2012
Mynor Chacon Daz
8.4
MA1003
A causa de la simetra del solido el centroide tiene que ubicarse en el eje Z,

Mxy
, pero
de manera que solamente se debe calcular z =
M
Z Z Z
Mxy =
r cos r2 sen drdd
z (x, y, z) dxdy dz =
s
olido
Z Z Z
Mxy =
s
olido
r cos |{z}
1 r2 sen dr d d
| {z }
| {z }
=
0
0
2
=
0
|J|
r3 cos sen dr d d
a4
8
Mxy
=
Entonces z =
M
a4
8
a3 [2 2]
3
Por lo tanto el centroide es
3a
.
8 2 2
!
3a
.
0, 0,
8 2 2
Recuerde que al ser = 1, el volumen V es igual a la masa M.

Ejemplo 24
Si la densidad en cada punto es inversamente proporcional a la distancia
del punto al origen de coordenadas, calcule la masa del casquete esferico
limitado por:
2
x + y 2 + z 2 = a2
2
2
2
2
x +y +z =b
x2 + y 2 = z 2
z0
0<a<b
c
2012
Mynor Chac
on Daz
224
MA1003
8.4
Soluci
on:
En la siguiente imagen se bordea con una lnea negra el casquete esferico:
Para calcular la masa se va a utilizar las coordenadas esfericas:
x = cos sen
y = sen sen
z = cos
J(, , ) = 2 sen
Para apreciar mejor la variacion de , , , considere la proyeccion del

casquete esferico en el primer cuadrante del plano YZ:
k
y gracias a la
r
+ +
simetra del casquete esferico, se tiene que la masa viene dada por:
Se sabe que la densidad es (x, y, z) = p
x2
225
y2
z2
c
2012
Mynor Chacon Daz
8.4
MA1003
r2 sen
V =
0
=k
k
dr d d
r
r sen dr d d
0

k 2
2
=
2 2
b a
2
Como se aprecia en la siguiente imagen:
c
2012
Mynor Chac
on Daz
226
CAPITULO
Integral de Lnea e Integral de Superficie
En este captulo se extiende el concepto de integral sobre un intervalo

a una integral sobre una curva, as como la integral doble sobre una region
plana se extiende a una integral sobre una superficie.
9.1.
9.1.1.
Integral de Lnea
Integral de lnea de un campo escalar
Como se ve en el primer curso de calculo, dada una funci

Z b on continua
f (x) dx. Si
definida en el intervalo [a, b], se puede definir la integral
a
ahora f : [a, b] R R3 ( f un Campo escalar) y el intervalo [a, b] se

sustituye por la curva tridimensional C := r(t) = (X(t), Y (t), Z(t)). Entonces se tiene una integral de lnea sobre C.
Para definir tal integral la curva se divide o particiona en n subarcos de
longitud Sk , k = 1, ..., n, en cada subarco se elige un punto (xk , yk , zk ),
siendo f continua desde a = r(t0 ) hasta b = r(t1 ) y X(t), Y (t), Z(t) con
primer derivada continua entonces cuando n o sea si Sk 0 se
tiene:
lm
n
X
!
f (xk , yk , zk ) Sk
Z
=
f (x, y, z) ds
C
k=1
Recuerdese que, si s es la longitud de arco, entonces

227
ds
= kr0 (t)k, por lo
dt
9.1
MA1003
que ds = kr0 (t)k dt, de esta forma

Z
t1
t1
f (X(t), Y (t), Z(t)) kr (t)k dt =
f (x, y, z)ds =
f (r(t)) kr0 (t)k dt
t0
t0
Ejemplo 1
Z
Calcule
(2x + 9xy) ds, sobre la curva x = t, y = t2 , z = t3 , t [0, 1].
C
Soluci
on

Se tiene r(t) = t, t2 , t3 kr0 (t)k = 1 + 4t2 + 9t4 , de manera que
Z
(2x+9xy) ds =
C

1
2t + 9t3
1 + 4t2 + 9t4 dt =
4
Z
1
14
14 14 1
u du =
6
1
2

(Siendo u = 1 + 4t2 + 9t4 , du = 8t + 36t3 dt = 4 2t + 9t3 dt)
Ejemplo 2
Z
Calcule
(x 3y 2 + z) ds, sobre la curva en la imagen siguiente:
C
c
2012
Mynor Chac
on Daz
228
MA1003
9.1
Soluci
on
En este caso
Z
(x 3y + z) ds =
(x 3y 2 + z) ds
(x 3y + z) ds +
C1
C2
El segmento de recta, C1 , desde (0, 0, 0) hasta (1, 1, 1), es una curva que se
parametriza como r(t) = (t, t, t), t [0, 1], entonces r0 (t) = (1, 1, 1)
kr0 (t)k = 3, as:
Z

3 dt = 3
t 3t + t
(x 3y + z) ds =
C1
2t 3t2 dt = 0
El segmento de recta, C2 , desde (1, 1, 1) hasta (1, 1, 0), tiene vector director
(0, 0, 1) la parametrizacion es r(t) = (1, 1, 1 t), t [0, 1], entonces se
tiene r0 (t) = (0, 0, 1) kr0 (t)k = 1, por lo tanto:
(x 3y + z) ds =
C2
Z
(1 3 + 1 t) dt =
1 t dt =
0
3
2
por u
ltimo
Z
(x 3y 2 + z) ds =
3
2
(Este valor es independiente de la parametrizaci

on del camino Ci )
9.1.2.
Aplicaciones fsicas
Las siguiente son las formulas para la masa y el primer momento de resortes, varillas delgadas y alambres a lo largo de una curva suave C, siendo
(x, y, z) la densidad en cada punto (x,y,z) de la curva.
229
c
2012
Mynor Chacon Daz
9.1
MA1003
Z
Masa=M =
(x, y, z)ds
C
Z
Primer momento con restecto al plano yz = Myz =
x (x, y, z) dt
C
Z
Primer momento con restecto al plano xz = Mxz =
y (x, y, z) dt
C
Z
Primer momento con restecto al plano xz = Mxy =
z (x, y, z) dt
C

Centro de masa = (x, y, z) =
Myz Mxz Mxy

,
,
M M M
Momentos de Inercia con respecto a los ejes y otras rectas

Z
Ix =
y +z
Z
(x, y, z) dt
Iy =

x2 + z 2 (x, y, z) dt
Z
Iz =
x +y
Z
(x, y, z) dt
IL =
r2 (x, y, z) dt
Siendo r(x, y, z) = distancia del punto (x, y, z) a larrecta L.

IL
Radio de giro con respecto a una recta L = RL =
M
c
2012
Mynor Chac
on Daz
230
MA1003
9.1
Las siguiente son las formulas para la masa y los momentos para cuerpos en
el espacio, siendo (x, y, z) la densidad en cada punto (x,y,z) del solido
Z Z Z
Masa=M =
(x, y, z) dx dy dz
D
Z Z Z
Mxy =
D
Z Z Z
Mxz =
y (x, y, z) dx dy dz
D
Z Z Z
x (x, y, z) dx dy dz
Myz =
D

Centro de masa = (x, y, , z) =
Myz Mxz Mxy

,
M M M
Momentos de Inercia
Z Z Z
Ix =
Z Z Z
(y + z ) (x, y, z)dx dy
D
Iy =
(x2 + z 2 ) (x, y, z)dx
Z Z Z
Iz =
(x2 + y 2 ) (x, y, z)dx dy
231
c
2012
Mynor Chacon Daz
9.1
MA1003
Ejemplo 3
Calcule la masa de un alambre cuya tarea es la curva de interseccion de la
esfera x2 + y 2 + z 2 = 1 y el plano x + y + z = 0, si la densidad del alambre
es (x, y, z) = x2 .
Soluci
on
En este caso
(
C=
x2 + y 2 + z 2 = 1
x+y+z =0
Despejando z de la segunda ecuacion y reemplazando en la primera se obtiene:
1 = x2 + y 2 + (x y)2
= 2x2 + 2xy + 2y 2
c
2012
Mynor Chac
on Daz
232
MA1003
9.1
1
= x2 + xy + y 2
2
1
1
= x2 + xy + y 2 y 2 + y 2
4
4

2
1
3
=
x + y + y2
2
4
!2

2
1
3
=
x+ y +
y
2
2

la ecuacion
1
x+ y
2
2
+
3
y
2
!2

=
x + 2y =
3
y
2
1
2
2
, sugiere:
cos t
1
2
1
2
cos t 21 y
sen t
x=
y = 2 sen t
3
x=
1
2
cos t
1
6
sen t
0 t 2
y = 2 sen t
3
Puesto que z = (x + y), se tiene:
z = x y

r
2
1
1
sen t
= cos t sen t
3
2
6
1
1
= cos t sen t
2
6
233
c
2012
Mynor Chacon Daz
9.1
MA1003
De manera que la parametrizacion de la curva C es:
x = 12 cos t 16 sen t
y = 26 sen t
0 t 2
z = 1 cos t 1 sen t
2
Las coordenadas de r0 (t) son :
x0 = 12 sen t 16 cos t
y 0 = 26 cos t
z 0 = 1 sen t 1 cos t
2
entonces
(x0 )2 = 12 sen2 t + 2 12 16 sen t cos t + 61 cos2 t
(y 0 )2 = 64 cos2 t
(z 0 )2 = 1 sen2 t 2 1 1 sen t cos t + 1 cos2 t

2
6
2 6
de forma que
(x0 ) + (y 0 ) + (z 0 ) = sen2 t +
4
2
cos2 t + cos2 t = cos2 t + sen2 t = 1
6
6
q
kr (t)k = (x0 )2 + (y 0 )2 + (z 0 )2 = 1
0
la masa del alambre esta dada por la integral de lnea:

Z
Z
ds =
x ds =
2
1
1
cos t sen t
kr0 (t)k dt
2
6
por u
ltimo
Z
x ds =
C

1
1 1
1
2
2
2
cos t 2 sen t cos t + sen t dt =
2
6
3
2 6

c
2012
Mynor Chac
on Daz
234
MA1003
9.1
Ejemplo 4
Z
Calcule

x2 + y 2 z 2 ds, C es la cuarta parte de la circunferencia
(
C=
x2 + y 2 + z 2 = R 2
x=y
Soluci
on
Puesto que x = y se tiene
x2
z2
2x + z = R 2 + 2 = 1
R
R

R
R
Entonces C = f (t) = cos t, cos t, R sen t , definiendo esta funcion
2
2
en la TI, se obtiene kr0 (t)k = |R| = R
235
c
2012
Mynor Chacon Daz
9.1
MA1003
luego el campo escalar g(x, y, z) = x2 + y 2 z 2 = R2 cos(2t), de forma que
x +y z
Z
ds =
R cos(2t) R dt = R
cos(2t) dt = 0
0
9.1.3.
Integral de lnea de un campo vectorial
Sea la familia de curvas y = c x2 + 1, c R
Si se considera los vectores tangentes en cada punto de cada una de las

curvas de la familia:
se obtiene un campo de vectores:

c
2012
Mynor Chac
on Daz
236
MA1003
9.1
En realidad a cada punto del plano,(x, y), se le asigna un vector, v R2 , a

esto se le llama Campo Vectorial. Formalmente este es un Campo Vectorial
F : R2 R2 . Se denotara
F (x, y) = (P (x, y), Q(x, y)) = P i +Q j

siendo P y Q campos escalares de R2 en R.
Si F : R3 R3 , se denotara
F (x, y, z) = (P (x, y), Q(x, y), R(x, y)) = P i +Q j +R k

siendo P , Q y R campos escalares de R2 en R.
La imagen de Campo Vectorial de R3 en R3 es la siguiente:
En general un campo vectorial es una funcion F : Rn Rm , donde

F (x1 , ..., xn ) = (f1 (x1 , ..., xn ), f2 (x1 , ..., xn ), ....fm (x1 , ..., xn ))
siendo fi : Rn R, campos escalares.
Sea r(t) una funcion vectorial, la curva descrita por r es Regular o suave
si r0 (t) es continua y nunca se anula, es decir si sus componentes tienen
237
c
2012
Mynor Chacon Daz
9.1
MA1003
primeras derivadas parciales continuas que no se anulan de forma simultanea.

La curva descrita por r(t) es Regular a trozos o suave a trozos si esta
formada por un n
umero finito de curvas suaves.
Ahora sea C := r(t) una curva suave que comparte dominio con el campo
vectorial F (x, y, z) = (P, Q, R), F : R3 R3 . Suponga ademas que la
curva C se particiona en n tramos de longitud Sk y que F representa
una fuerza que act
ua sobre cada punto de la curva, por ejemplo la fuerza
de gravedad. La imagen siguiente muestra una ampliacion de uno de estos
tramos, junto con la fuerza F aplicada a un punto del tramo, el punto r(tk )
y el vector tangente unitario T (t):
Recuerde que la proyeccion ortogonal del vector Fk sobre el tangente Tk es:

ProyTk (Fk ) =
Fk Tk
= Fk Tk
kTk k2
Puesto que
trabajo = (fuerza constante) (distancia)
Entonces el trabajo realizado en cada segmento de la curva es
Wk = Fk Tk Sk
As el trabajo realizado a lo largo de la curva es C es:
lm
n
X
!
Fk Tk Sk
Z
F T ds
=
C
k=1
Conocida como la integral de lnea de un campo Vectorial.

c
2012
Mynor Chac
on Daz
238
MA1003
9.1
Recuerdese que ds = kr0 (t)k dt, entonces:

Z
F T ds
W =
C
Z
=
C
Z
=
r0 (t)
kr0 (t)k dt
F 0
kr (t)k
F r0 (t) dt
Z
F dr
=
C
Z
(P, Q, R)
=
C
dx dy dz
, ,
dt dt dt

dt
Z
(P, Q, R) (dx, dy, dz)
=
C
Z
P dx + Q dy + R dz
=
C
En el caso en que F : R2 R2 , la integral de lnea se define analogamente

como:
Z
Z
Z
F T ds =
F dr =
P dx + Q dy
C
De forma analoga si C es una curva cerrada suave en el plano, N el vector

normal unitario a C que apunta hacia afuera, entonces el Flujo de F a
trav
es de C es
I
I
F N ds =
P dy Q dx
C
Ejemplo 5

Calcule el trabajo realizado por el campo de fuerzas F (x, y, z) = y x2 i

+ z y 2 j + x z 2 k , a lo largo de la curva r(t) = t, t2 , t3 , desde
(0, 0, 0) hasta (1, 1, 1).
Soluci
on
239
c
2012
Mynor Chacon Daz
9.1
MA1003
La integral a calcular es
Z
F T ds
W =
C
Z
F dr
=
C
Z
=
P dx + Q dy + R dz
C

y x2 dx + z y 2 dy + x z 2 dz
=
C

t2 t2 dt + t3 t4 2t dt + t t6 3t2 dt
=
0
4

2t 2t5 + 3t3 3t8 dt
29
60
Ahora se va a obtener el resultado utilizando Winplot, en la opcion para

graficar en 3D, sigua la ruta Ecua >Ecua Diferencial...
ahora en las casillas x0 , y 0 , z 0 se escriben las componentes P, Q, R y se habilita

la ventanilla campo
c
2012
Mynor Chac
on Daz
240
MA1003
9.1
en la ventana dimensiones de la caja se deja todo igual,
obteniendose una grafica del campo de vectores. Seguidamente haga clic

en caja
y fijemos la posicion 1 < x < 1, 1 < y < 1, 1 < z < 1, con el objetivo
de apreciar mejor la grafica
241
c
2012
Mynor Chacon Daz
9.1
MA1003
ahora se ingresa la curva r(t) y se obtiene la imagen con el campo de vectores

y la curva
por u
ltimo se sigue la ruta Dos >Integral de l
nea...
c
2012
Mynor Chac
on Daz
242
MA1003
9.1
y se obtiene el valor de la integral
Z
C

29
0,48333
y x2 dx + z y 2 dy + x z 2 dz =
60
Ejemplo 6
Calcule el flujo del campo de fuerzas F (x, y) = (x y) i +x j , a traves

de la circunferencia x2 + y 2 = 1.
Soluci
on
La integral a calcular es
I
I
F N ds =
(x y) dy x dx
C
Para efectos de utilizar Winplot se debe ver como sigue:

243
c
2012
Mynor Chacon Daz
9.1
MA1003
I
x dx + (x y) dy
C
Es decir, al momento de ingresar el campo en Winplot se hara como
G(x, y) = x i + (x y) j
Procediendo de forma analoga al ejemplo anterior, se introduce el campo

G(x, y) = (x, x y, 0) y la curva r(t) = (cos t, sen t, 0), de esta manera se
obtiene el campo de fuerzas (visto desde el eje Z) y el valor de la integral de
lnea = 3,14159:
en efecto
I
I
F N ds =
(x y) dy x dx =
C
Ejemplo 7
Calcular el trabajo realizado por el campo de fuerzas F (x, y, z) = y 2 ~i +
z 2 ~j + x2 ~k a lo largo de la esfera x2 + y 2 + z 2 = a2 y el cilindro x2 + y 2 = ax,
siendo a > 0, z 0. El camino es recorrido de modo que observando el plano
XY desde el eje positivo Z el sentido es el de las agujas del reloj.
Soluci
on
Se tiene la curva
c
2012
Mynor Chac
on Daz
244
MA1003
9.1
(
C=
x 2 + y 2 + z 2 = a2
x2 + y 2 = ax

a 2
a2
Recuerde que x2 + y 2 = ax equivale a escribir x
+ y 2 = , esto nos
2
4
sugiere
a
a
x = 2 cos t + 2
0 t 2
y = a sen t
2
Por otro lado sustituyendo la ecuacion x2 + y 2 = ax en x2 + y 2 + z 2 = a2 se

obtiene ax + z 2 = a2 , entonces:
ax + z 2 = a2 z 2 = a(a x)

a
a
z 2 = a a cos t
2
2
z2 = a
z2 =
a
2

a
cos t
2
a2
(1 cos t)
2
a
z = 1 cos t
2
De esta manera la parametrizacion de la curva C es
a
a
x = cos t +
2
2
y = a sen t
2
z
=
1 cos t
0 t 2
por lo que
245
c
2012
Mynor Chacon Daz
9.1
MA1003
dx = a2 sen t dt
dy = a cos t dt
2
sen t
dz = 2a 2 1cos
dt
luego
a

a2
a3
sen2 t sen t dt = sen3 t dt
4
2
8
a
a3

a2
z 2 dy = (1 cos t)
cos t dt =
cos t cos2 t
2
2
4
a
a 2 a
sen t
a3 (cos t + 1)2 sen t

cos t +
dt =
x2 dz =
dt
2
2
2 2 1 cos t
8 2 1 cos t
y 2 dx =
Note que
Z 2
a3
sen3 t dt = 0
y dx =
8 0
C
I
3 Z 2

a
a3
z 2 dy =
cos t cos2 t dt =
4 0
4
C
Z 2
I
2
3
(cos t + 1) sen t
a
dt = 0
x2 dz =
1 cos t
8 2 0
C
I
(use u = 1 cos t)
Entonces
I
I
F dr =
y 2 dx + z 2 dy + x2 dz = 0
a3
+ 0
4
2
Note que la posicion de P = r(0) = (2, 0, 0) y Q = r( ) = (1, 1, ) en la
2
2
curva indica que el recorrido se hace en contra de la orientacion pedida.
c
2012
Mynor Chac
on Daz
246
MA1003
9.1
esto implica que se le debe cambiar el signo al resultado anterior:

a3
4
C
C
Para comprobar el resultado con Winplot, se le debe asignar un valor cona3
= 2 , la
creto a la variable a. Por ejemplo si a = 2, el resultado
4
imagen siguiente muestra el campo de vectores, la curva C y el resultado
2 6,28319
I
F dr =
y 2 dx + z 2 dy + x2 dz =
247
c
2012
Mynor Chacon Daz
9.1
MA1003
El resultado es negativo debido a que Winplot no sabe que la direccion que

se quera era la inversa.
9.1.4.
Campos conservativos.
El termino conservativo proviene de la fsica, por ejemplo en los campos

gravitacionales o electricos el trabajo realizado para mover una masa de un
punto a otro solo depende de estos y no de la trayectoria que une ambos
puntos. En terminos matematicos si el campo F (x, y) = (P, Q) esta definido
en una region abierta D, donde P, Q tienen derivadas parciales continuas
entonces F es conservativo si y solo si
P
Q
=
y
x
Analogamente el campo F (x, y, z) = (P, Q, R) definido en una esfera abierta
B, donde P, Q, R tienen derivadas parciales continuas entonces F es conservativo si y solo si:
rot(F ) = 0
La expresion rot(F ), significa rotacional de F es un operador
rot : R3 R3
definido como sigue

i

rot(F ) = F = x

P

k
P
Q P
R Q P

=
z
y
z z
x x
y

R
Definici
on
Una region abierta D Rn es convexo si todo par de puntos pueden unirse
mediante un camino regular a trozos contenido en D
Definici
on
Una region abierta D Rn es convexo si todo par de puntos pueden unirse
mediante un camino rectilneo contenido en D
c
2012
Mynor Chac
on Daz
248
MA1003
9.1
Primer teorema fundamental del c

alculo para integrales de lnea
Sea F : D Rn Rm continuo, D conexo. Supongase que es conservativo,
o sea la integral de lnea de D no depende del camino, entonces
F = f,
x D,
Z
siendo
f (x) =
F dr
a
al campo escalar f : Rm R se le conoce como Funci

on Potencial de F
Segundo teorema fundamental del c
alculo para integrales de lnea
Sea f : D Rn R diferenciable, tal que f es continuo en D conexo,

C : r(t) una curva regular a trozos en D. Entonces
Z
f dr = f (b) f (a)
a
249
c
2012
Mynor Chacon Daz
9.1
MA1003
Corolario
Si C : r(t) es una curva cerrada entonces
I
f dr = 0
Teorema
Sea F : D Rn Rm continuo, D conexo. Son equivalentes:
Existe una campo escalar f , tal que F = f .
Z
F dr es independiente del camino, sobre cualquier curva regular a
C
trozos C = r(t).
I
F dr = 0, sobre cualquier curva cerrada regular a trozos C = r(t).
C
Ejemplo 8
Z
(3,1)
x + 2y
y
dx
+
dy, sobre un camino que no atraviesa
2
(x + y)2
(1,1) (x + y)
la recta y = x
Calcule
Soluci
on

x + 2y
y
El ejercicio sugiere que el campo F (x, y) =
,
es conser(x + y)2 (x + y)2
vativo, de hecho esto se verifica luego de ver que las derivadas cruzadas son
iguales, tal como se aprecia en la siguiente imagen:
c
2012
Mynor Chac
on Daz
250
MA1003
9.1
P
2y
Q
=
=
2
y
(x + y)
x
por lo tanto f / f = F es decir tal que
f
x
x+2y
(x+y)2
f
y
y
(x+y)2
integrando cada ecuacion tenemos

R x+2y
y
f
(x,
y)
=
dx = ln |x + y| x+y
+ K(y)
(x+y)2
R
y
y
f (x, y) = (x+y)
2 dy = ln |x + y| x+y + H(x)
Se deduce que K(y) = H(x) = C(constante), es decir

f (x, y) = ln |x + y|
y
+ C
x+y
de esta forma se tiene que:
251
c
2012
Mynor Chacon Daz
9.1
MA1003
(3,1)
(1,1)
x + 2y
y
1
dx
+
dy
=
f
(3,
1)
f
(1,
1)
=
ln
2
+
(x + y)2
(x + y)2
4
Ejemplo 9
Z
(2,2,0)
e ln y dx +
Siendo y > 0, calcule

(1,1, 2 )
ex
+ sen z
y

dy + y cos z dz
Soluci
on

ex
Se utilizara MuPaD para calcular el rot(F ), siendo F (x, y, z) = e ln y, + sen z, y cos z .
y
La imagen siguiente muestra el comando que se debe de usar y el resultado
correspondiente:

Puesto que rot (F ) = 0, se tiene un campo conservativo, por lo que existe

un campo escalar f (x, y, z), tal que f = F , esto significa:
f
x
= ex ln y
f
y
f
y
= y cos z
ex
y
+ sen z

c
2012
Mynor Chac
on Daz
252
MA1003
9.1
R x
f
(x,
y,
z)
=
e ln y dx = ex ln y + H(y, z)
R x
f (x, y, z) = ey + sen z dy = ex ln y + y sen z + I(x, z)
f (x, y, z) = R y cos z dz = y sen z + J(x, y)

comparando se deduce que H(y, z) = y sen z; J(x, y) = ex ln y; I(x, z) =
K, Kconstante, de donde
f (x, y, z) = ex ln y + y sen z + K
por lo tanto
Z
(2,2,0)

x

e
+ sen z dy+y cos zdz = f (2, 2, 0)f 1, 1,
e ln ydx+
= e2 ln 21
y
2
x
(1,1, 2 )
Ejemplo 10
Z
Determine
F dr, siendo C el arco de la curva que va de A(1, 0, 1) a
B(0, 1, 1):
siendo
253
c
2012
Mynor Chacon Daz
9.1
MA1003
(
:=
5x2 + 5y 2 + 2z 2 8xy + 2xz + 2yz = 9

x+yz =0
y F (x, y, z) = 2x sen(y) + 5yz, x2 cos(y) 6yez + 5xz, 3y 2 ez + 5xy
Soluci
on
Debido a que la parametrizacion de la curva no parece simple, calculamos

el rotacional del campo con MuPaD:
Con esto se concluye que la integral de lnea es independiente de la trayectoria, determinemos entonces la funcion potencial
f
x
= 2x sen y + 5yz
f
y
= x2 cos y 6yez + 5xz
f
y
= 3y 2 ez + 5xy
f (x, y, z) = x2 sen(y) + 5xyz + H(y, z)
f (x, y, z) = x2 sen(y) 3y 2 ez + 5xyz + I(x, z)
f (x, y, z) = 3y 2 ez + 5xyz + J(x, y)

Por tanto se deduce que H(y, z) = 3y 2 ez ; J(x, y) = x2 sen y; I(x, z) =
K, Kconstante, de donde
f (x, y, z) = x2 sen(y) 3y 2 ez + 5xyz + K
c
2012
Mynor Chac
on Daz
254
MA1003
9.2
Z
F dr =
C
9.2.

(2x sen y + 5yz) dx+ x2 cos y 6yez + 5xz dy+ 3y 2 ez + 5xy dz = 3e
Teorema de Green
Sean P, Q campos escalares, con derivadas parciales de primer orden continuas en un abierto A R2 . Sea C una curva cerrada simple regular a trozos,
sea R la union de C y su interior, ademas suponga que R esta contenida en
A, entonces se cumple
ZZ
R
Q P
x
y
I
dx dy =
P dx + Q dy
C
(La integral de lnea se calcula siguiendo la direccion positiva, esto es la

direccion contraria al movimiento de las agujas de un reloj)
255
c
2012
Mynor Chacon Daz
9.2
MA1003
Ejemplo 11

Sea F (x, y) = 2xy x2 , x + y 2 , C := r(t) la curva cerradaZ que limita la
F dr.
region R, determinada por las curvas y = x2 , x = y 2 calcule
C
1. De forma directa.
2. Utilizando el teorema de Green.
Soluci
on
1. Para el calculo directo, se debe considerar C = C1 C2 y calcular dos
integrales de lnea.
I
Z
2xy x
F dr =
C
dx+ x + y
Z
dy+
C1
C2
La parametrizacion de C1 es:
x=t
C1 :=
y = t2
0t1
as
(
c
2012
Mynor Chac
on Daz

2xy x2 dx+ x + y 2 dy
dx = dt
dy = 2tdt
256
MA1003
9.2
De esta manera:

2(t)(t2 ) t2 dt + t + t4 2tdt
F dr =
0
C1

2t5 + 2t3 + t2 dt
=
0
7
6
La parametrizacion de C2 es:
x=t
C1 :=
0t1
y=t
note que t = 0 genera el punto (0, 0) y t = 1 genera el punto (1, 1), esto
significa que la parametrizacion produce que el recorrido de la curva
sea el inverso al deseado, por lo que se le a
nadira un (-) a la integral
de lnea. Note que se tiene
(
dx = 2tdt
dy = dt
De esta manera:

2(t2 )(t) t4 2tdt + t2 + t2 dt
F dr =
C1
Z
=

4t4 2t5 + 2t2 dt
17
15
finalmente se tiene
I
C

7 17
1
2xy x2 dx + x + y 2 dy =
=
6 15
30
257
c
2012
Mynor Chacon Daz
9.2
MA1003
2. Para usar el teorema de Green, note que P (x, y) = 2xy x2 , Q(x, y) =

x + y 2 , de manera que:
I
2xy x
dx + x + y
ZZ
dy
=
R
Q P
x
y

dx dy
ZZ
(1 2x) dx dy
=
R
(1 2x) dy dx
=
x2
1
30
Ejemplo 12
Calcular
I

arctan x + y 2 dx + ln y x2 dy
siendo C el contorno de la region R de la siguiente imagen:
Soluci
on
Note que C = C1 C2 C3 C4 , de manera que el calculo directo
requiere de 4 integrales no muy faciles, este hecho sugiere utilizar el teorema
de Green, dado que la curva es cerrada simple regular a trozos, por lo que
I
arctan x + y
dx + ln y x
2(x + y) dA
dy =
R
utilizando coordenadas polares, se obtiene:

c
2012
Mynor Chac
on Daz
ZZ
258
MA1003
I
arctan x + y
9.2
dx + ln y x
ZZ
dy
2(x + y) dA
2 (r cos + r sen ) r dr d
=
0
2 r2 (cos + sen ) dr d
=
0
224
3
Ejemplo 13
Z
Use el teorema de Green para calcular
(ex sen y my) dx+(ex cos y m) dy,
siendo C la semicircunferencia superior x2 + y 2 = ax, recorrida en sentido

positivo desde (a, 0) hasta (0, 0).
Soluci
on
a
Observe que la parametrizacion de la curva C es x = [a + cos t] , y =
2
a
sen t, al sustitur en el campo escalar quedan expresiones que no son facil2
mente integrables, este hecho sugiere usar el teorema de Green, sin embargo
259
c
2012
Mynor Chacon Daz
9.2
MA1003
en este caso la curva C no es cerrada, de manera que se va a considerar el segmento desde (0, 0) hasta (a, 0) para cerrar la curva, es decir se aplicara el
teorema de Green a la curva = C
entonces
ZZ
(e sen y my) dx + (e cos y m) dy
ZZ
Q P
x
y
ZZ
=
m dx dy
R
ZZ
= m
dx dy
R
= m a(R)
luego
c
2012
Mynor Chac
on Daz
260
dx dy
ex cos y ex cos y + m dx dy
= m

a 2
2
2
a2 m
8
MA1003
9.2
I
F dr =
Z
Z
a2 m
=
F dr +
F dr
8
C
Z
Z
a2 m
F dr
F dr =
8
Z
F dr +
F dr
Z
de forma que solo hace falta calcular
F dr, para esto se parametriza
:=
x = t dx = dt
0t1
y = 0 dy = 0
entonces
Z
(ex sen y my) dx + (ex cos y m) dy = 0

|{z}
|
{z
}
0
en conclusion:
a2 m
F dr =
Z
C
Z
F dr =
| {z }
a2 m
8
o sea
Z
(ex sen y my) dx + (ex cos y m) dy =
9.2.1.
a2 m
8
Area
de una regi
on plana, usando Teorema de Green
Sea R una region plana limitada por una curva cerrada simple suave a trozos
C. Si se considera el campo F (x, y) = (y/2, x/2) el Teorema de Green
establece que:
I
y
x
dx + dy =
2
2
C
Z Z
dA = a(R)
R
Reescribiendo lo anterior se obtiene que el area de la region R, a(R), puede

obtenerse por una integral de l`nea:
1
a(R) =
2
I
x dy y dx
C
261
c
2012
Mynor Chacon Daz
9.2
MA1003
Use lo anterior para obtener el area de la region indicada en la siguiente

figura, la cual tiene como borde la Hypocycloide de ecuacion x2/3 +y 2/3 = 32/3

y esta representada por la funcion vectorial r(t) = 3 cos3 t, 3 sen3 t :
En este caso la parametrizacion es
3
2
x = 3 cos t dx = 9 cos t sen t dt

C :=
y = 3 sen3 t dy = 9 sen2 t cos t dt
0 t 2
De manera que
1
a(R) =
2
Z
=
I
x dy y dx
C
2

3 cos3 t 9 sen2 t cos t dt 3 sen3 t 9 cos2 t sen t dt

27 cos4 t sen2 t + 27 cos2 t sen4 t dt
=
0

27 cos2 t sen2 t cos2 t + sen2 t dt
=
0
cos2 t sen2 t dt
= 27
0
27
=
4
27
4
=
=
4 cos2 t sen2 t dt
0
2
sen2 (2t) dt
27
4
c
2012
Mynor Chac
on Daz
262
MA1003
9.2.2.
9.2
Teorema de Green para una regi

on con agujero
Considerese la siguiente region R limitada por dos curvas cerradas C1 , C2 ,

denotaremos C = C1 C2
Considerese ademas un campo con derivadas parciales de primer orden continuas en una region abierta que contiene a R. Si partimos la region en dos
subregiones R1 , R2 , como se aprecia en la imagen siguiente:
Intuitivamente es claro que

ZZ
R
Q P
x
y
ZZ
dA =
R1
Q P
x
y
263
ZZ
dA +
R2
Q P
x
y

dA
c
2012
Mynor Chacon Daz
9.2
MA1003
si se considera solamente la parte superior, es posible aplicar el teorema de

Green, siendo 1 la parte superior de C1 y C2 , denotado por C1s y C2s , unido
con los segmentos rectilneos AB y P Q (La flecha indica la orientacion
de la curva, en este caso de izquierda a derecha), es decir
s
s
= C1
C2
AB P Q
as la igualdad de Green establece lo siguiente

ZZ
R1
Q P
x
y
Z
dA =
Z
F dr =
Z
F dr+
s
C2
s
C1
Z
F dr+
AB
Z
F dr+
F dr
P Q
consideraciones similares para la parte inferior
llevan a establecer la igualdad de Green de la siguiente manera

ZZ
R2
Q P
x
y
Z
dA =
c
2012
Mynor Chac
on Daz
Z
F dr =
i
C1
264
Z
F dr+
i
C2
Z
F dr+
AB
Z
F dr+
P Q
F dr
MA1003
9.2
al sumar

ZZ
R1
Q P
x
y
ZZ
dA +
R2
Q P
x
y

dA
se obtiene
Z
F dr+
Z
s
C1
Z
F dr+
s
C2
Z
F dr+
AB
Z
F dr+
Z
F dr+
i
C1
P Q
Z
F dr+
i
C2
Z
F dr+
AB
F dr
P Q
las integrales sobre los segmentos se anulan por tener orientaciones contrarias, de manera que
ZZ
R1
Q P
x
y
ZZ
dA+
R2
Q P
x
y
dA =
s
C1
F dr +
{z
R
C1
ZZ
R
Q P
x
y
F dr
C2
C1+
Z
=
F dr
C
lo que significa que el teorema de Green es valido para la region con agujero,
la notacion C1+ , significa que la curva se recorre en sentido positivo y . C1 ,
significa que la curva se recorre en sentido negativo.
265
F dr +
F dr +
s
C2
} |
{z
R

ZZ
Q P
Q P
=
dA +
x
y
x
y
R1
R2
Z
Z
=
F dr +
F dr
ZZ
dA
i
C1
c
2012
Mynor Chacon Daz
C2
i
C2
F dr
F dr
}
9.2
MA1003
ZZ
R
Q P
x
y
Z
dA =
F dr
C
La igualdad anterior puede escribirse de la siguiente forma:

ZZ
R
Q P
x
y
Z
F dr
dA =
F dr
C2+
C1+
En tal caso si el campo F = (P, Q), sea conservativo en R, se tiene:

ZZ
0=
R
Q P
x
y
Z
0 =
Z
F dr
dA =
C1+
F dr
Z
F dr =
C1+
F dr
C2+
De estos resultados se deduce el siguiente:
c
2012
Mynor Chac
on Daz
F dr
C2+
266
F dr
C2+
C1+
MA1003
9.2
Teorema
Sea F = (P, Q) conservativo en todo R2 excepto quizas en en el punto P .
Entonces si C1 y C2 son dos curvas que rodean a P , se cumple
F dr
F dr =
C2+
C1+
A
un mas si consideramos un crculo, C, alrededor del punto P , entonces
podermos concluir:
Z
F dr =
C1+
Z
F dr =
C2+
F dr
C+
Ejemplo 14

Sea F (x, y) =

y
x
,
, definido en R = R2 {(0, 0)}.
x2 + y 2 x2 + y 2
1. Demuestre que F (x, y) es conservativo en R.

I
y
x
2. Calcule
dx + 2
dy siendo C cualquier curva cerrada
2
2
x + y2
C x +y
simple regular a trozos contenida en R, que contenga a (0, 0) en su
interior y recorrida en el sentido positivo.
267
c
2012
Mynor Chacon Daz
9.2
MA1003
y
x
dx + 2
dy siendo C cualquier curva cerrada
2
+y
x + y2
C
simple regular a trozos contenida en R, que no contenga a (0, 0) en su
interior.
I
3. Calcule
x2
Soluci
on
1. Observe que
P
Q
=
x
y
esto implica que F (x, y) es conservativo en R.

2. La curva C puede ser como la que aparece en la imagen siguiente:
c
2012
Mynor Chac
on Daz
268
MA1003
9.2

Para calcular la integral considere el disco D = (x, y) /x2 + y 2 = 1
y la region R1
de acuerdo al teorema anterior aplicado a la region R1 se tiene

I
C
y
x
dx + 2
dy =
2
2
x +y
x + y2
I
D
x2
y
x
dx + 2
dy
2
+y
x + y2
la parametrizacion del disco es x = cos t, y = sen t, 0 t 2, por

lo tanto
I
C
y
x
dx
+
dy =
x2 + y 2
x2 + y 2
I
D
y
x
dx
+
dy
x2 + y 2
x2 + y 2
sen t ( sen t dt) cos t ( cos t dt)
=
0
sen2 t cos2 t dt
=
0
dt
=
0
= 2
3. La curva C en este caso puede ser como la que aparece en la imagen

siguiente:
269
c
2012
Mynor Chacon Daz
9.3
MA1003
en tal caso el campo F (x, y) es conservativo en R2 de manera que es

aplicable el teorema de Green, de forma que
I
C
y
x
dx + 2
dy =
2
2
x +y
x + y2
I
C

Q P
x
y
|
{z
}
0
= 0
c
2012
Mynor Chac
on Daz
270
MA1003
9.3.
9.3.1.
9.3
Integral de superficie
Consideraciones b
asicas y Producto Fundamental
En una curva espacial los puntos se mueven con un grado de libertad,

si se considera el lugar geometrico de puntos espaciales que se mueven con
dos grados de libertad se obtiene lo que se llama una superficie. Puede ser
representada explcitamente en la forma z = f (x, y) o de forma implcita
F (x, y, z) = 0. Existe una tercer forma en que se puede representar una
superficie, esta es la forma parametrica:
x = X(u, v)
S :=
y = Y (u, v)
z = Z(u, v)
o r(u, v) = (X(u, v), Y (u, v), Z(u, v)) ; r : T R3 , siendo T una region
del plano U V .
Considerese la siguiente imagen
la imagen bajo r del peque

no rectangulo rojo en la region T , de area uv,
es la porcion r(T ) de color rojo en S, la cual se aproxima con el paralelor
r
gramo determinado por los vectores
u y
v, este paralelogramo
u
v
tiene area

r
r

r
r
=
u v

u
u
u
v
v
271
c
2012
Mynor Chacon Daz
9.3
MA1003
Esto sugiere que el area a(S) de la superficie S esta dada por:
ZZ
a(S) =
S

ZZ
r

r

du dv
dS =
u
v
T
Un punto (u, v) en T en el cual
r
r r
r
6= 0 y
,
son continuas, se
u
v
u u
llama punto regular de r. Los puntos en los que no se cumplen alguna de

estas condiciones se llaman puntos singulares de r. Una superficie r(T )
se llama superficie regular si todos sus puntos son regulares..
Al producto
9.3.2.
r
r
se le llama Producto Fundametal.

u
v
Ejemplos de parametrizaciones de algunas superficies

El cono
z 2 = x2 + y 2 ,
0za
x = u cos v
y = u sen v
z=u
0<u<a
r(u, v) = (u cos v, u sen v, u)

i
j
k

r
r

= cos v
sen v 1 = (u cos v, u sen v, u)

u v
u sen v u cos v 0
De manera que la norma del producto fundamental es

r

r

= 2u2 = 2u
u v
c
2012
Mynor Chac
on Daz
272
MA1003
9.3
La esfera
x 2 + y 2 + z 2 = a2
x = a cos u sen v
y = a sen u sen v
z = a cos v
0 < u < 2
0<v<
r(u, v) = (a cos u sen v, a sen u sen v, a cos v)

i
j
k

r r

= a sen u sen v a cos u sen v

0
= a2 cos u sen2 v, sen u sen2 v, sen v cos v

u v
a cos u cos v a sen u cos v a sen v

r
r
= a2 cos2 u sen4 v + sen2 u sen4 v + sen2 v cos2 v

u v
= a2 sen4 v + sen2 v cos2 v

p
= a2 sen2 v (sen2 v + cos2 v)
= a2 sen v
El cilindro
x 2 + y 2 = a2 ,
0zb
x = a cos u
y = a sen u
z=v
273
0 < u < 2
0<v<b
c
2012
Mynor Chacon Daz
9.3
MA1003
r(u, v) = (a cos u, a sen u, v)

i
j
k

r
r

= a sen u a cos u 0 = (a cos u, a sen u, 0)

u v

0
0
1

r

r

u v = a
Ejemplo 15
Determinar dos parametrizaciones distintas del cilindro x2 + (y 3)2 =
9, 0 z 5
Siguiendo lo visto anteriormente se tiene:
x = 3 cos u
y 3 = 3 sen u
0 < u < 2
z=v
0<v<5
O sea:
x = 3 cos u
y = 3 + 3 sen u
z=v
0 < u < 2
0<v<5
De otra manera, usando coordenadas cilindricas se tiene:
x2 + (y 3)2 = 9 x2 + y 2 6y + 9 = 9 r2 6r sen = 0 r = 6 sen

entonces
c
2012
Mynor Chac
on Daz
274
MA1003
9.3
x = r cos = 6 sen cos = 3 sen 2

y = r sen = 6 sen sen = 6 sen2
z=z
Es decir, la nueva parametrizacion es:
x = 3 sen 2u
y = 6 sen2 u
z=v
0 u ,
0v5
Usando GenGCF y la segunda parametrizacion se obtiene la imagen siguiente:
El Toro de revoluci
on
Se obtiene al girar un crculo de radio r > 0 y centro en (R, 0, 0), en el plano

yz, alrededor del eje Z.
275
c
2012
Mynor Chacon Daz
9.3
MA1003
La grafica de esta superficie y su parametrizacion son las siguientes:
x = (R + r cos u) cos v
y = (R + r cos u) sen v
z = r sen u
0 < u < 2
0 < v < 2
r(u, v) = ((R + r cos u) cos v, (R + r cos u) sen v, r sen u)
r r
== ((R + r cos u)r cos u cos v, (R + r cos u)r sen v cos u, , (R + r cos u)r sen u)
u v

r
r

u v = r(R + r cos u)
Superficie z = f (x, y)
Si se utilizan x,y como los parametros u,v, entonces la parametrizacion elemental es
x=u
y=v
z = f (u, v)
c
2012
Mynor Chac
on Daz
276
MA1003
9.3
r(u, v) = (u, v, f (u, v))
r
=
En tal caso tenemos
u
f
1, 0,
u
r
y
=
v
f
0, 1,
v
De manera que el producto fundamental da como resultado:

i j k

r
r
f
f
f
= 1 0 u = , , 1

u v
u v
0 1 f
v
r
r
resulta ser igual al vector normal al plano tangente a

u v
la superficie explcita z = f (x, y).
La norma del producto fundamental es

q

r
r
2
2

u v = 1 + (fu ) + (fv )
Note que el
Superficie z = f (x, y) definida implcitamente por G(x, y, z) = 0
z
Gx
z
Gy
= ,
= ,
x
Gz
y
Gz
como en el caso anterior y utilizando x,y como parametros tenemos:
En tal caso, la derivacion implcita garantiza que
r
r
=
x y
Esto implica que:

z z
Gx Gy
, ,1 =
,
,1
x y
Gz Gz

q
r

r

1 + (fx )2 + (fy )2
x y =
s
2
2
Gx
Gy
=
+
+1
Gz
Gz
s
G2x + G2y + G2z
=
G2z
p 2
Gx + G2y + G2z
=
|Gz |
277
c
2012
Mynor Chacon Daz
9.3
MA1003
Ejemplo 16
Para el caso en que la superficie sea el plano ax + by + cz = d, se tiene que
la ecuacion implcita de la superficie es G(x, y, z) = ax + by + cz d = 0,
as entonces Gx = a, Gy = b, Gz = c, siguiendo las formulas en el anterior
caso se concluye que:
r
r
=
x y

a b
, ,1
c c
r r
resulta ser un vector paralelo al vector normal a plano

x y
(a, b, c), pero no son iguales, a menos que c = 1. Por u
ltimo
Note que el

2
2
2
r

r

= a +b +c
x y
|c|
9.3.3.
C
alculo de
areas de superficies
Resumiendo lo visto hasta ahora tenemos:

ZZ
A(S) =
dS =
S=r(T )
RR r

u
T
r
v
dudv =
RR q
R
1 + gx2 + (gy )2 dA
en el caso S := z = g(x, y)
RR Gx2 +Gy2 +Gz2
dA
|Gz |
R
en el caso S := G(x, y, z) = 0
Ejemplo 17
Dada la superficie S de ecuacion vectorial
r(u, v) = (u cos v, u sen v, u2 )
donde 0 u 4, 0 v 2
1. Haga el dibujo de la superficie.
2. Obtenga el area de S usando la parametrizacion r(u, v).
c
2012
Mynor Chac
on Daz
278
MA1003
9.3
3. Obtenga el area de S, considerandola como una superficie del tipo

z = f (x, y).
Soluci
on
Dado que superficie S tiene ecuacion vectorial:
r(u, v) = (u cos v, u sen v, u2 )
donde 0 u 4, 0 v 2.
entonces el maximo valor de z es z = 16,
1. En la proxima imagen aparecere la superficie limitada por el plano
z = 16.
2. Evidentemente la superficie es una porcion del paraboloide circular,

z = x2 + y 2 . Seg
un la parametrizacion r(u, v) = (u cos v, u sen v, u2 ),
se tiene entonces las variables x = u cos v, y = u sen v, por lo que u, v
juegan el papel de r, en las coordenas polares.
El producto fundamental

i
j
k

r
r

= cos v
sen v 2u = 2u2 cos v, 2u2 sen v, u

u v
u sen v u cos v 0
279
c
2012
Mynor Chacon Daz
9.3
MA1003
Para calcular el area de S se necesita calcular la norma del producto

fundamental

r
r
4
2
4
2
2
4
2
2

u v = 4u cos v + 4u sen v + u = 4u + u = u 4u + 1
Entonces el area de S viene dada por

Z Z
Z 2 Z 4
i
r

r
h
2

65 65 1
u 4u + 1 du dv =
u v du dv =
6
T
0
0
3. Explcitamente la superficie S, es z = g(x, y) = x2 + y 2 , en este caso el
normal exterir es (gx , gy , 1) = (2x, 2y, 1), de manera que

p

r
r
= k(2x, 2y, 1)k = 4x2 + 4y 2 + 1

x y
Z Z
Z Z
A(S) =
p
4x2 + 4y 2 + 1 dA
dS =
R
Dado que el valor mayor de z es 16, la proyeccion (R) en en plano XY

es el disco de radio 4, x2 + y 2 = 16, utilizando polares se tiene
Z Z
A(S) =
Z
p
4x2 + 4y 2 + 1 dA =
Z
0
4r2
i
h
+ 1rdrd =
65 65 1
6
Los calculos se realizaron con la TI-89, note que en los incisos e) y f)

se llega a la misma integral:
c
2012
Mynor Chac
on Daz
280
MA1003
9.3
Ejemplo 18
Calcular el area de la region que en el plano x + y + z = 3 determina el
cilindro x2 + y 2 = 9.
Soluci
on
La imagen siguiente muestra la superficie S y el vector perpendicular al
plano x + y + z = 3, v = (1, 1, 1), como se vio en la seccion anterior en este
caso tomando como parametros x,y el producto fundamental es:

r
r
= k(1, 1, 1)k = 3

x y
y la region T es la proyeccion de S en el plano XY, a saber el disco x2 +y 2 = 9
de esta forma

Z Z
Z Z
r

r

dS =
dxdy
=
3
dxdy
=
3A(T
)
=
9
3

y
S
T x
T
Z Z
A(S) =
281
c
2012
Mynor Chacon Daz
9.3
MA1003
Ejemplo 19
Considere la superficie r(, ) = ( cos , sen , g (, ))
1. Siendo z = g(, ), pruebe que

s
2 2
r r
z
= 2 + 2 z
+

2. Use la opcion Superficie param

etrica en GenGCF o Ecua>Param
etrica
en Winplot para obtener una imagen de la superficie.
x = cos
r(, ) :=
y = sen
z=
[0, 1], [0, ]
3. Obtenga el area de la rampa.
Solucion:
1. Puesto que r(, ) = ( cos , sen , g (, )), entonces

r
=
g
cos , sen ,

,
r
=
g
sen , cos ,
Efectuando el producto cruz se obtiene

r r
= (g sen g cos , g sen + g cos , )

as

s
2 2
r r
z
= 2 + 2 z
+

No olvide que
g
z
g
z
=
y que
=
c
2012
Mynor Chac
on Daz
282
MA1003
9.3
En el caso en que
x = cos
r(, ) :=
y = sen
z=
[0, 1], [0, ]
se tiene z = g(, ) = o sea r(, ) = ( cos , sen , ), la imagen

respectiva generada por GenGCF es la siguiente
Si se ubica la superficie dentro del cilindro x2 + y 2 = 1, se aprecia as:

283
c
2012
Mynor Chacon Daz
9.3
MA1003
2. En este caso entonces se tiene

s
2 2 q
p
r r
z
2 + 2 (0)2 + (1)2 =
= 2 + 2 z
2 + 1
+
=

De manera que

Z Z
Z Z p
Z Z 1p
r r
dA =
A(S) =
2 + 1dA =
2 + 1d d

T
T
0
0
Utilizando la TI

Z Z
h

i
r r
dA =
Se obtiene A(S) =
2
+
ln
1
+
2

2
T
c
2012
Mynor Chac
on Daz
284
MA1003
9.3.4.
9.3
C
alculo de Integrales de superficie
Considerese una funcion f (x, y, z) definida sobre los puntos de la superficie

S, en tal caso se define la integral de superfice
ZZ
ZZ
f (x, y, z) dS =

r
r
du dv
f (r(u, v))
u
v
Por ejemplo, si f (x, y, z) es la carga por unidad de

area en cada punto de S,
esto es la densidad de carga, entonces la integral de superficie representa la
Carga Total en S.
En el caso en que la superficie este dada en forma explcita, z = g(x, y) se
tiene:

r

r

f (r (u, v))
=
u v dudv
T
ZZ
q
f (x, y, g(x, y)) 1 + gx2 + (gy )2 dA
=
ZZ
ZZ
f (x, y, z) dS
S=r(T )
En el caso en que la superficie este dada en forma implcita, G(x, y, z) =

z g(x, y) = 0 se tiene:
ZZ
f (x, y, z) dS
S=r(T )

r

r

=
f (r (u, v))
u v dudv
T
p 2
ZZ
Gx + Gy2 + Gz2
=
f (x, y, g(x, y))
dA
|Gz |
R
ZZ
En ambos casos R es la proyeccion sobre el plano XY

285
c
2012
Mynor Chacon Daz
9.3
MA1003
Ejemplo 20
Siendo S la porcion del primer octante del cilindro y 2 + z 2 = 9, entre los
planos x = 0, x = 4.
Z Z
(x + z) dS, utilizando x e y co-
1. Calcule la integral de superficie

S
mo
Debe aparecer el calculo del producto fundamental

parametros.

r
r

x y y la region donde varan x e y.
2. Calcule la integral anterior, utilizando la parametrizacion r(u, v), que

se obtiene de las coordenadas
polares.
Debe aparecer el calculo del

r

r

producto fundamental
u v y la region donde varan u y v.
Soluci
on
1. La superficie S, limitada por y 2 + z 2 = 9, x = 0, x = 4, en el primer

octante, es la siguiente porcion de cilindro.
En
p este caso la superficie se puede escribir en forma explcita z =
9 y 2 , de esta manera entonces
| {z }
g(x,y)
2y
y
gx = 0, gy = p
=
z
2 9 y2
por lo tanto
c
2012
Mynor Chac
on Daz
286
MA1003
9.3

q
r

r

(gx )2 + (gy )2 + 1
x y =
r
z2 + y2
=
z2
p
z2 + y2
=
z
p
9 y2 + y2
=
z
3
= p
9 y2
la region donde vara x e y, es el rectangulo
as
287
c
2012
Mynor Chacon Daz
9.3
MA1003
Z Z

3
9 y2 p
dy dx
9 y2
0
0
!
Z 3Z 4
x
p
+ 1 dx dy
= 3
9 y2
0
0
!
Z 3
8
p
= 3
+ 4 dy
9 y2
0

0
3
+ 36
= 24 arc sen arc sen
3
3
(x + z) dS =
S
= 24

2
x+
+ 36
= 12 + 36
2. Puesto que y 2 +z 2 = 9, se parametriza y = 3 cos(u), z = 3 sen(u), luego
x = v, es decir
h i
r(u, v) = (v, 3 cos(u), 3 sen(u)) , v [0, 4], u 0,
2
r
= (0, 3 sen(u), 3 cos(u))
u
r
= (1, 0, 0)
v
entonces se tiene que
r
r
= (0, 3 cos(u), 3 sen(u))

u v

r
r

u v = 3
de forma que
Z Z
c
2012
Mynor Chac
on Daz
(v + 3 sen u) 3 dv du = 12 + 36
(x + z) dS =
0
288
MA1003
9.3
Como se aprecia en la siguiente imagen, que muestra la integral hecha

en MuPad
Ejemplo 21
Z Z p
x2 + y 2 dS; donde S es la superficie lateral del cono
Calcule
S
x2 y 2 z 2
+
2 = 0;
a2 a2
b
0zb
Soluci
on
x2 y 2 z 2
En este caso la superficie es es cono G(x, y, z) = 2 + 2 2 = 0, de manera
a
a
b
que:
G2x
G2y
G2z

4 x2 y 2
4z 2
= 2 2+ 2 + 4
a a
a
b

4 z2
4z 2
= 2 2 + 4
a b
b

4z 2 1
1
=
+
b2 a2 b2

4z 2 a2 + b2
=
b2
a2 b2
Entonces tenemos:
G2x + G2y + G2z
b2 2z 2
2z
=
a + b2 =
2
2
2z
ab
b
p
dS =
289
a2 + b2
a
c
2012
Mynor Chacon Daz
9.3
MA1003
Por lo tanto
Z Z p
x2
y2
dS =
a2 + b 2
a
Z Z p
x2 + y 2 dA
Donde T es la proyeccion del cono en el plano XY , a saber el disco x2 + y 2 =

a2 , usando coordenadas polares x = r cos , x = r sen
Por u
ltimo
Z Z p
x2
y2
dS =
a2 + b 2
a
Z Z p
a2 + b 2
a
x2 + y 2 dA
T
2
r2 dr d
a2 + b 2 3
a 2
3a
2a2 a2 + b2
=
3
=
9.3.5.
Superficies orientables
Una superficie S es orientable o tiene dos caras si se puede definir un

vector normal N en todo punto interior(que no pertenezca a la frontera de
S), de forma tal que los vectores normales pueden variar de forma continua
sobre la superficie. El ejemplo clasico de una superficie no orientable es la
Cinta de Mobius.
Suponga que cada hormiga representa un vector normal N , si bien las hormigas avanzan de forma continua sobre la cinta, notese que sin importar el
punto de partida cada hormiga regresara a este pero con las patas y el
cuerpo en posicion invertida. En terminos vectoriales, N luego de caminar
sobre S sin cambiar de cara(exterior o interior) se convierte en N , este
c
2012
Mynor Chac
on Daz
290
MA1003
9.3
hecho no se puede lograr con ning

un campo de vectores se concluye que S
es No-Orientable. Una Superifice entonces debe tener dos caras y una vez
elegida la direccion hacia la que se desea que apunte el vector N , se dira que
S esta orientada.
Si S es orientable y F es un campo definido sobre S se define la integral de
superficie en forma vectorial
ZZ
F N dS
S
Del algebra lineal se conoce que la proyeccion ortogonal del vector F sobre
F N
N es el vector ProyN F =
=F N
kN k2
es decir
ZZ
ZZ
F N dS =
ProyN F dS
Geometricamente representa la integral sobre S de la componente normal

del campo F , a esta u
ltima integral se le conoce como Integral de Flujo.
En el caso en que (x, y, z), es la densidad de un fludo en el punto (x, y, z)
entonces
ZZ
F N dS = masa de fludo que fluye a traves de S por unida de tiempo
S
9.3.6.
C
alculo de Integrales de Flujo
En general dada una superficie S := r(T ), recuerdese que

r
u
N =
r
u
291
r
v
r
v
c
2012
Mynor Chacon Daz
9.3
MA1003
Entonces
ZZ
ZZ
F N dS =
S=r(T )
F (r(u, v))
T
r
u
r
u
|

ZZ
=
F (r(u, v))
T
r
v

r
v
{z }
N
r
r
u v

r
r

u v du dv
|
{z
}
dS

du dv
Dada una superficie S := G(x, y, z) = z g(x, y) = 0, y utilizando x e y

como parametros recuerdese que
r
x

N=
r
r
y
G
=

r
kGk
y
en tal situacion se tienen dos casos:
(g , g , 1)
p x 2y 2
1 + gx + gy
Normal hacia arriba
N=
(g , g , 1)
p x y2
1 + gx + gy2
Normal hacia abajo

puesto que
dS =
q
1 + gx2 + gy2 dx dy N dS = (gx , gy , 1) dx dy
entonces
ZZ
F N dS =
S=r(R)
c
2012
Mynor Chac
on Daz
ZZ
F (x, y, g(x, y)) (gx , gy , 1) dx dy
Normal hacia arriba
ZZ
F (x, y, g(x, y)) (gx , gy , 1) dx dy
Normal hacia abajo

292
MA1003
9.3
Si la superficie z = g(x, y) esta definida de forma implcita por la ecuacion

G(x, y, z) = 0, se tiene en forma analoga las formulas siguientes
ZZ
F N dS =
S=r(R)
ZZ
R

Gx Gy
F (r(x, y))
,
, 1 dx dy
Gz Gz
Normal hacia arriba

ZZ
Gx Gx
F (r(x, y)) , , 1 dx dy
Gz
Gz
Normal hacia abajo
La Region R es la proyeccion de la superficie en el plano XY

Ejemplo 22
Sea S la parte del paraboloide z = 4 x2 y 2 , sobre el plano XY , orientada
por un vector hacia arriba. Un fludo de densidad fluye a traves de S
seg
un el campo de velocidades F (x, y, z) = (x, y, z). Hallar la razon de flujo
de masa a traves de S.
Soluci
on
Entonces
293
c
2012
Mynor Chacon Daz
9.4
MA1003
Z Z
Z Z
F N dS =
F N dS
ZZ
(x, y, z) (gx , gy , 1) dx dy
=
R
ZZ
=
(x, y, , 4 x2 y 2 ) (2x, 2y, 1) dx dy
ZZ
x2 + y 2 + 4 dx dy
=
R
La region R es el disco x2 +y 2 = 4, por lo tanto se utilizaran las coordenadas

polares, para obtener:
ZZ
x + y + 4 dx dy =
R

r2 + 4 r dr d = 24
9.4.
Teoremas de Stokes y de Gauss
9.4.1.
Teorema de Stokes
Teorema de Stokes
Sean S una superficie orientable, con un vector normal unitario N , limitada
por curva cerrada simple regular a trozos, C. Si F (x, y, z) = (P, Q, R) es un
campo vectorial cuyas componentes tienen derivadas parciales continuas en
una region abierta que contiene a S y a C, entonces
I
ZZ
rot(F ) N dS
F dr =
C
equivalente a escribir
I
ZZ
rot(F ) N dS
P dx + Q dy + R dz =
C
c
2012
Mynor Chac
on Daz
294
MA1003
9.4
Visto desde la punta del vector N , la direccion en que se recorre la curva

debe ser la positiva(contraria al movimiento de las agujas del reloj), para una
superficie como la de la figura anterior se pueden dar los siguientes casos:
Ejemplo 23
ZZ
rot(F ) N dS; en una integal de lnea
Transforme la integral de superficie

S
y calcularla. Considerese F (x, y, z) = (y, z, x), S la parte del paraboloide

z = 1 x2 y 2 ; z y N la normal con componente z no negativa.
Soluci
on
El teorema de Stokes, en este caso establece que
ZZ
I
rot(F ) N dS =
y dx + z dy + x dz
C
siendo C la curva x2 + y 2 = 1, recorrida en la direccion que se observa en la

figura siguiente:
295
c
2012
Mynor Chacon Daz
9.4
MA1003
una parametrizacion de C es
x = cos t
y = sen t
z=0
0 t 2
entonces
dx = sen t dt
dy = cos t dt
dz = 0
el calculo de la integral queda de la siguiente manera:
ZZ
I
rot(F ) N dS
y dx + |{z}
z dy + x |{z}
dz
0
sen t [ sen t dt]

0
Z
=
c
2012
Mynor Chac
on Daz
296
sen2 t dt
MA1003
9.4
Ejemplo 24
I
Aplicando el Teorema de Stokes calcule
x dx + (x + y) dy + (x + y + z) dz,
siendo C, la curva
x = 3 sen t
C :=
y = 3 cos t
z = 3 sen t + 3 cos t
0 t 2
Soluci
on
El teorema de Stokes, establece que
I
ZZ
rot(F ) N dS
x dx + (x + y) dy + (x + y + z) dz =
C
Observe que C es la interseccion del cilindro x2 +y 2 = 9 y el plano z = x+y.

Ahora hay que responder a la pregunta C
ual es la superficie S cuyo borde
o frontera es C? En este ejemplo, la respuesta es sencilla, tal superficie es el
pedazo de plano dentro del cilindro. De forma que la superficie a considerar
es S := x + y z = 0, en tal caso el producto fundamental es
r
r
r
r
= (1, 1, 1)
= (1, 1, 1)
x
y
x
y
Antes de elegir uno de ellos, puesto
que r(t) = (3 sen t, 3 cos t, 3 sen t+3 cos t),
se tiene r(0) = (0, 3, 3) y r
= (3, 0, 3). Para que la orientacion de la
2
r
r
curva sea la positiva, se debe entonces elegir el vector
= (1, 1, 1).
x
y
297
c
2012
Mynor Chacon Daz
9.4
MA1003
La imagen siguiente muestra el calculo del rotacional del campo F (x, y, z) =

(x, x + y, x + y + z)
Es decir rot(F ) = (1, 1, 1), por u

ltimo
I
ZZ
x dx + (x + y) dy + (x + y + z) dz
rot(F ) N dS
ZZ
(1, 1, 1) (1, 1, 1)
=
R
ZZ
=
dA
R
A(R)
(3)2
No olvide que R es la proyeccion en el plano XY , esto es el disco x2 +y 2 = 9.

Ejemplo 25

Sea F (x, y, z) = x, z 2 , x2 + y 2 + z y considere la curva
2
2
x +y =2z
C :=
2x + 2y + z = 0
Z
Calcule
F dr
C
c
2012
Mynor Chac
on Daz
298
MA1003
9.4
1. De forma directa
2. Aplicando el teorema de Stokes y tomando C como la frontera del
plano.
3. Aplicando el teorema de Stokes y tomando C como el borde del paraboloide.
Soluci
on
1. Despejando z en la ecuacion del cilindro, sustituyendola en la ecuacion

del plano y completando cuadrados se obtiene

(x 1)
2
2

+
(y 1)
2
2
= 1
De aqu se llega a x = 1 + 2 cos t, y = 1 + 2 sen t, por lo tanto z =

2x 2y = 4 [1 + cos t + sen t]
x = 1 + 2 cos t
y = 1 + 2 sen t
z = 4 [1 + cos t + sen t]
0 t 2
Las coordenadas de r0 (t) son :

0
x = 2 sen t dt
y 0 = 2 cos t dt
0
z = 4 [sen t cos t] dt
luego utilizando MuPaD se obtiene
x dx + z 2 dy + (x2 + y 2 + z) dz =
[1 + 2 cos t] [2 sen t] + 16 [1 + cos t + sen t] [2 cos t] + 8 [sen t cos t] =

24 cos t + 6 sen t + 60 cos t sen t + 64 cos2 t + 32 cos3 t + 32 cos t sen2 t + 64 cos2 t sen t dt
299
c
2012
Mynor Chacon Daz
9.4
MA1003
y luego
24 cos t + 6 sen t + 60 cos t sen t + 64 cos2 t + ...
32 cos3 t + 32 cos t sen2 t + 64 cos2 t sen t
dt = 64
en resumen se tiene
I
F dr
x dx + z 2 dy + (x2 + y 2 + z) dz
64

Note que r(0) = (3, 1, 8) y r
= (1, 3, 8), por lo que la direccion
2
inducida por la parametrizacion es la positiva, vista desde el eje positivo
Z.
c
2012
Mynor Chac
on Daz
300
MA1003
9.4
Esto va a determinar la direccion del vector N (que es la misma del

producto fundamental), ya que vista desde la punta de tal vector, la
direccion del recorrido de la curva, debe coincidir con la la inducida
por la parametrizacion.
2. Ahora se va a aplicar el teorema de Stokes, considerando la curva C,
como la frontera del pedazo de plano que queda denttro del cilindro,
es decir S := 2x + 2y + z = 0, en tal caso el producto fundamental es
r
r
r
r
= (2, 2, 1)
= (2, 2, 1)
x
y
x
y
En este caso desde el vector (2, 2, 1), la direccion de la curva es la
r
r
= (2, 2, 1),
deseada, por lo que se elige
x
y
luego al calcular el rotacional de F se obtiene rot(F ) = (2y2z, 2x, 0),

entonces
301
c
2012
Mynor Chacon Daz
9.4
MA1003
ZZ
ZZ
rot(F ) N dS
(2y 2z, 2x, 0) (2, 2, 1) dA
ZZ
(4y 4z 4x) dA
=
R
ZZ
=
[y x (2x 2y)] dA
4
R
ZZ
=
(x + 3y) dA
R
Procediendo como en la parametrizacion de la curva se obtiene
(x 1)
2
2

+
(y 1)
2
2
= 1
O equivalentemente (x 1)2 + (y 1)2 = 22 , de donde se obtiene que

la proyeccion es

R = (x, y)/(x 1)2 + (y 1)2 22
De manera que para calcular la integral sobre la proyeccion se utilizaran

las coordenadas polares generalizadas
x = 1 + r cos
y = 1 + r sen
|J| = r
0 2
0r2
as que
c
2012
Mynor Chac
on Daz
302
MA1003
9.4
ZZ
ZZ
rot(F ) N dS
(2y 2z, 2x, 0) (2, 2, 1) dA
ZZ
=
(x + 3y) dA
R
Z
=
(4 + r cos + 3r sen ) r dr d
4
0
64
3. Ahora se va a aplicar el teorema de Stokes, considerando la curva C,

como el borde del pedazo de paraboloide que queda por encima del
plano, es decir S := z = 2 x2 y 2 , entonces gx = 2x, gy = 2y, en
| {z }
g(x,y)
tal caso el producto fundamental es
( gx , gy , 1) = ( 2x, 2y, 1)
r
r
x
y
( g , g , 1) = ( 2x, 2y, 1)
x
y
observe que en el punto (0, 0, 2),
(0, 0, 1)
r
r
x
y
(0, 0, 1)
En este caso desde el vector (0, 0, 1), la direccion de la curva es la

r
r
deseada, por lo que se elige
= ( 2x, 2y, 1 ),
x
y
303
c
2012
Mynor Chacon Daz
9.4
MA1003
Evidentemente el valor del rotacional es el mismo y la integral de superficie es la siguiente
ZZ
ZZ
rot(F ) N dS
(2y 2z, 2x, 0) ( 2x, 2y, 1 ) dA
ZZ
=
(4xz) dA
R
ZZ
=

4x 2 x2 y 2 dA
ZZ
=

4xy 2 + 4x3 8x dA
Usando polares generalizadas obtenemos el integrando

8 8r cos + 4(1 + r cos )3 + 4(1 + r cos )(1 + r sen )2 r
Usando MuPaD se obtiene

Z
0
8r2 cos t + 8r2 sen t + 8r3 cos t sen t + 12r3 cos2 t + ...
4r4 cos3 t + 4r3 sen2 t + 4r4 cos t sen2 t dr d = 64

Como se aprecia en la imagen siguiente
c
2012
Mynor Chac
on Daz
304
MA1003
9.4
En resumen se obtuvo
ZZ
I
rot(F ) N dS
F dr = 64 =
S
9.4.2.
Teorema de Gauss o de la divergencia
Teorema de Gauss
Sean una region solida limitada por una superficie orientada por un vector
normal unitario dirigido al exterior de . Si F (x, y, z) = (P, Q, R) es
un campo vectorial cuyas componentes tienen derivadas parciales continuas
en , entonces
ZZZ
ZZ
div (F ) dV
F N dS =
equivalente a escribir
ZZZ
ZZ
F N dS =
S
305
P
Q R
+
+
x
y
z

dV
c
2012
Mynor Chacon Daz
9.4
MA1003
Ejemplo 26
Sea F (x, y, z) = (x, y, 2z) y S la superficie x2 + y 2 + (z 1)2 = 1, utilizar
el teorema de Gauss para determinar el valor de
ZZ
F N dS
S
Soluci
on
Es evidente que

div (F ) = 1 1 + 2 = 2 = (x, y, z)/x2 + y 2 + (z 1)2 1
F N dS =
dV = 2V () = 2
2 dV = 2
ZZZ
ZZZ
ZZ
4 3
1
3

=
8
3
Ejemplo 27
Sea S la superficie limitada por el paraboloide z = x2 + y 2 y el plano z = 4,
N la normal unitaria apuntando hacia afuera de las superficies y F (x, y, z) =
(x, y, 3) un campo definido sobre S.
ZZ
F N dS
1. Calcule
S
ZZ
F N dS utilizando el teorema de la divergencia.
2. Calcule
S
Soluci
on
1. El solido , esta limitado por dos superficies, el plano S1 := z = 4 y el

paraboloide S2 := z = x2 + y 2 , de manera que
ZZ
ZZ
ZZ
F N dS =
S
c
2012
Mynor Chac
on Daz
F N1 dS +
S1
306
F N2 dS
S2
MA1003
9.4
Primero se calculara la integral sobre S1 := z = 4, este caso sabemos

r
r
= ((0, 0, 1) (note que este vector apunta al exterior del

que
x
y
solido), entonces
ZZ
F N1 dS =
S1
ZZ
ZZ
(x, y, 3)(0, 0, 1)dA = 3

dA = 3 22 = 12
Note que la proyeccion R es es disco x2 + y 2 4
Para calcular la integral sobre S2 := z = x2 + y 2 ; 0 z 4, se

pueden utilizar x, y como parametros, sin embargo, aunque resulte un
poco mas laborioso se parametrizara como sigue:
x = u cos v
S2 :=
y = u sen v
z = u2
0 v 2
0 u2 4
0u2
307
c
2012
Mynor Chacon Daz
9.4
MA1003
De manera que
r(u, v) = u cos v, u sen v, u2

i
j
k

r
r

= cos v
sen v 2u = 2u2 cos v, 2u2 sen v, u

u v
u sen v u cos v 0
Notese que al punto de la superficie P = (0, 1, 1), corresponden u =
1; v = , por tanto en P , el producto fundamental es

2
r
r
= (0, 2, 1)
u v
Al graficar este vector sobre el punto P sobre S2 , se obtiene
Es decir el producto fundamental apunta hacia adentro del solido, por

lo que se debe cambiar la orientacion del producto fundamental, es
decir se tomara

r
r
= 2u2 cos v, 2u2 sen v, u

u v
c
2012
Mynor Chac
on Daz
308
MA1003
9.4
ZZ
ZZ
F N2 dS
r
r
du dv
u v
F (r(u, v))
S2

(u cos v, u sen v, 3) 2u2 cos v, 2u2 sen v, u du dv
=
0
0
2

2u3 cos2 v + 2u3 sen2 v 3u du dv
=
0
0
2 Z

2u3 3u du dv
=
0
= 4
En resumen se tiene
ZZ
F N dS = 12 + 4 = 16
S
2. El solido , tiene como proyeccion el disco x2 + y 2 4, por lo tanto se

utilizaran coordenadas cilndricas para el calculo de la integral triple
que se obtiene al aplicar el teorema de la divergencia, a saber:
ZZ
ZZZ
F N dS
div (F ) dV
ZZZ
=
(1 + 1 + 0) dV
ZZZ
=
dV
r dz dr d
0
r2
= 16
309
c
2012
Mynor Chacon Daz
9.4
MA1003
Este resultado verifica el teorema de Gauss

ZZ
ZZZ
F N dS = 16 =
div (F ) dV
Ejemplo 28

Verifique el teorema de Gauss para el campo F (x, y, z) = x2 , y 2 , z 2 y el
solido limitado por las superficies x2 = 1 y; z = 0; z = y.
Soluci
on
El solido es el siguiente
Esta limitado por tres superficies

ZZ
ZZ
F N dS =
ZZ
F N1 dS +
S1
ZZ
F N2 dS : +
S2
Siendo S1 := x2 = 1 y; S2 := z = 0 y S3 := z = y.
Sobre S1 := x2 = 1 y, tenemos
c
2012
Mynor Chac
on Daz
310
F N3 dS
S3
MA1003
9.4
x=u
S2 :=
y = 1 u2 sen v
z=v
1u1
0 v 1 u2
Note que como x = u, entonces 1 u 1, luego sobre la superficie, z

vara de 0(el piso) hasta el borde del plano, es decir |{z}
z = y
|{z}
v
1u2
De manera que
r(u, v) = u cos v, u sen v, u2

i

j
k

r
r

= 1 2u 0 = (2u, 1, 0)

u v
0
0
1
Notese que al punto de la superficie P = (0, 1, 0), corresponden u = 0; v = 0,
por tanto en P , el producto fundamental es
r
r
= (0, 1, 0)
u v
Al graficar este vector sobre el punto P sobre S1 , se obtiene
Es decir el producto fundamental apunta hacia adentro del solido, por lo que
se debe cambiar la orientacion del producto fundamental, es decir se tomara
r
r
= (2u, 1, 0)
u v
Entonces
311
c
2012
Mynor Chacon Daz
9.4
MA1003
ZZ
ZZ
F N1 dS
F (r(u, v))
S1
1u2
=
1
u2 , 1 u2
2
, v2
(2u, 1, 0) dv du
1u2

2u3 + 1 2u2 + u4 dv du
=
1
r
r
du dv
u v
32
35
r
r
r
r
= (0, 1, 1) o
=
x
y
x
y
(0, 1, 1), de estos dos el que apunta hacia a fuera del solido es (0, 1, 1),
entonces
Sobre S2 := y + z = 0, sabemos que
ZZ
ZZ
F N2 dS

x2 , y 2 , z 2 (0, 1, 1) dA
S2
=
1
1x2
z 2 y 2 dy dx
| {z }
0
= 0

Note que la proyeccion es R = (x, y)/y 1 x2
c
2012
Mynor Chac
on Daz
312
MA1003
9.4
Por u
ltimo, sobre S3 := z = 0, sabemos que
r
r
x
y
= (0, 0, 1) o
r
r
= (0, 0, 1), de estos dos el que apunta hacia a fuera del solido
x
y
es (0, 0, 1), entonces
ZZ
ZZ
F N3 dS

x2 , y 2 , z 2 (0, 0, 1) dA
S3
ZZ
=
R
z 2 dA
|{z}
0
= 0
En resumen
ZZ
F N dS =
S
32
32
+ 0+ 0 =
35
35
Para verificar tal resultado el teorema de Gauss establece que
ZZ
ZZZ
F N dS =
div (F ) dV
313
c
2012
Mynor Chacon Daz
9.4
MA1003
ZZ
ZZZ
F N dS
div (F ) dV
ZZZ
=
(2x + 2y + 2z) dV
ZZZ
=
(x + y + z) dV
Z
=
1x2
(x + y + z) dz dy dx
2
1
Z
0
32
35
Este resultado verifica el teorema de Gauss

ZZZ
ZZ
32
div (F ) dV
F N dS =
=
35
c
2012
Mynor Chac
on Daz
314
Indice alfabetico
Campo conservativo, 248

Campo vectorial, 237
Campos escalares
Criterio de las derivadas de segundo orden, 132
extremos absolutos de..., 144
extremos condicionados de..., 148
extremos relativos de..., 129
Formula de Taylor de segundo orden, 130
punto crtico de ..., 129
punto de silla de ..., 129
Conjunto conexo, 248
Conjunto convexo, 248
curva en el espacio, 68
curva plana, 68
Curva Regular o suave, 238
Curva Regular(suave) a trozos, 238
Derivacion Implcita, 117
Derivada direccional, 118
valores estremos de.., 125
Derivadas Parciales, 97
con TI-89 y MuPAD, 99
de orden superior, 100
con TI-89 y MuPAD, 100
funcion vectorial, 63
derivacion e integracion
TI-89, 64
derivacion, integracion, 63
productos, 67
funciones de dos variables, 89
diferencial de.., 101
Diferencial Total de.., 106
Aproximacion lineal de.., 104
continuidad de.., 94
dominio de.., 91
lmite de.., 93
Gauss-Jordan
metodo de, 26
Gradiente, 113
Integral de lnea, 227
aplicaciones fsicas, 229
de un campo conservativo, 248
de un campo escalar, 227
de un campo vectorial, 236
Funcion potencial, 249
teorema fundamental 1, 249
teorema fundamental 2, 249
definicion de trabajo, 239
flujo en una region plana, 239
Integral de superficie, 271
area de una superficie, 272
calculo de , 285
calculo de areas, 278
315
9.4
MA1003
calculo de integral de flujo, 291

en forma vectorial, 291
integral de flujo, 291
parametrizaciones, 272
producto fundamental, 271, 272
punto regular, 272
punto singular, 272
superficie regular, 272
Integral doble, 159
cambio de variables, 206
en coordenadas elpticas, 186
en coordenadas polares, 178
en coordenadas polares generalizadas, 178
en region R, 163
en region rectangular, 159
formula general para cambio de
variables, 208
Teorema de Fubinni, 162
volumen, 160
Integral triple
calculo de un volumen, 192
cambio de orden de integracion,
189
como volumen, 189
en coordenadas cilndricas, 193
en coordenadas esfericas, 202
en coordenadas rectangulares, 188
integracion iterada, 190
Integrales m
ultiples
momentos, masa y centro de masa,
219
matriz
definicion de, 21
determinante, 24
escalonada reducida, 27
producto de, 22
suma de, 22
transpuesta de, 23
valor propio, 31
vector propio, 31
Matriz Hessiana, 131
MuPaD
calculo del rotacional, 252
parametrizacion, 69
Plano Tangente, 122
Regla de la cadena, 108
en sistema de ecuaciones, 116
involucrando funciones no explcitas,
111
utilizando MuPAD, 112
sistema de ecuaciones
MuPAD, 30
TI-89, 28
Superficies orientables, 290
Teorema de Gauss, 305
Teorema de Gauss o de la divergencia, 305
Teorema de Green, 255
area de una region plana, 261
para una region con agujero, 263
Teorema de Stokes, 294
Lagrange
valor propio
Funcion de, 148
MuPAD, 33
metodo para clasificar extremos
TI-89, 33
condicionados, 149
vector propio
Matriz Hessiana de L, 149
MuPAD, 33
Multiplicadores de..., 148
vectores
c
2012
Mynor Chac
on Daz
316
MA1003
9.4
magnitud de , 13
producto escalar, 16
producto por escalar, 12
producto vectorial, 17
suma de, 12
unitarios, 15
Winplot
transformaciones, 209
vectores en R2 , R3 , 10
integral de lnea, 240
representacion de punto, 7
317
c
2012
Mynor Chacon Daz

Calculo Iii - M. Chacon PDF

Cargado por

Copyright:

Formatos disponibles

Calculo Iii - M. Chacon PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Calculo Iii - M. Chacon PDF

Cargado por

Copyright:

Formatos disponibles

Material complementario

Lic. Mynor Chac

2. Operaciones con vectores

2.1. Representacion de un punto en R y R

2.2. Representacion de un vector en R2 y R3 . . . . . . . . . . . . . . .

2.3. Adicion, sustraccion y magnitud(o norma) de vectores en la TI-89

2.4. Vector unitario . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

2.5. Producto Punto o Escalar . . . . . . . . . . . . . . . . . . . . . . .

2.6. Producto Vectorial o Producto Cruz . . . . . . . . . . . . . . . . .

2.7. Uso de la tecla STO> . . . . . . . . . . . . . . . . . . . . . . . . . .

3. Operaciones con matrices

3.1. Representacion, adicion y producto de matrices . . . . . . . . . . .

3.2. Transpuesta de una matriz

3.3. Determinante de una matriz . . . . . . . . . . . . . . . . . . . . . .

3.4. Metodo de Gauss-Jordan

3.5. Uso del comando soluc() en el caso de infinitas soluciones . . . .

3.6. Valores Propios y Vectores Propios de una matriz . . . . . . . . . .

4.1. Representacion de superficies generadas por curvas planas, con

4.2.3. Hiperboloide de dos hojas:

4.2.8. Cilindro con Directriz en forma parametrica . . . . . . . . .

4.2.10. Superficie de Revolucion . . . . . . . . . . . . . . . . . . . .

4.2.11. Un caso particular . . . . . . . . . . . . . . . . . . . . . . .

5.1. Definicion y propiedades . . . . . . . . . . . . . . . . . . . . . . . .

5.2. Derivacion e integracion . . . . . . . . . . . . . . . . . . . . . . . .

5.4.1. Parametrizacion de y = f (x) . . . . . . . . . . . . . . . . .

5.4.2. Parametrizacion de F (x, y) = 0 . . . . . . . . . . . . . . . .

5.4.3. Parametrizacion de curvas en

5.5. Dibujo de curvas con Winplot y GenGCF . . . . . . . . . . . . . .

5.6. El Triedro Intrnseco o Marco de Frenet . . . . . . . . . . . . . . .

5.7. Curvatura, Torsion y Longitud de arco . . . . . . . . . . . . . . . .

6. Funciones de Varias Variables o Campos Escalares

6.1. Graficacion de funciones de dos variables . . . . . . . . . . . . . . .

6.2. Dominio de una funcion de dos variables . . . . . . . . . . . . . . .

6.3. Lmite de una funcion de dos variables . . . . . . . . . . . . . . . .

6.4. Derivadas Parciales . . . . . . . . . . . . . . . . . . . . . . . . . . .

6.5. Derivadas Parciales de orden superior . . . . . . . . . . . . . . . . 100

7. Extremos de Campos Escalares

9. Integral de Lnea e Integral de Superficie

Este material representa un esfuerzo por incorporar programas computacionales

Operaciones con vectores

Para representar el punto A = (2, 4) en el graficador Winplot, siga la

Luego aparece la ventana interior conteniendo el sistema de ejes cartesianos,

Aparece la ventana donde se digitan las coordendas del punto, se elige su

El punto debe verse como sigue.

Analogamente para representar el punto A = (1, 2, 3), se sigue la ruta

Se eligen las caractersticas del punto.

El punto debe verse como sigue.

Sean A = (2, 4) y B = (3, 1) dos puntos, el vector v =AB= B A =

siguiente ventana apareceran representados los puntos y el vector v =AB,

En la ventana segmento, las coordenadas de los puntos (x1 , y1 ) y (x2 , y2 )

Se obtiene as la imagen siguiente:

Para graficar el vector anclado se repite el proceso anterior, dibujando el

La imagen final es la siguiente:

v =AB y el vector anclado en el origen, esto es v = (1, 1, 1). Haga girar el

Seguidamente se va a calcular la magnitud

Luego de presionar Enter, se obtiene:

Una vez marcado el vector al que se le va a calcular la magnitud.

Dado el vector v Rn , un vector con la misma direccion pero de magniv

Producto Punto o Escalar

Si se desea el valor 3u5v+2(uv), simplemente se digita 3u-5v+2*prodVect(u,v),

Luego se digita 3a-5b y se presiona Enter.