Apuntes Sobre Análisis Matemático II
Apuntes Sobre Análisis Matemático II
Apuntes Sobre Análisis Matemático II
Lombardo
Apuntes de Análisis Matemático II.
Nelson G. Lombardo
Primera edición: Abril 2008
Edición:
© Copyright 2008-2009, Nelson G. Lombardo.
Esta documentación esta licenciada bajo los términos de la Licencia GNU para la Documentación
Libre. Ver en http:/www.gnu.org/copyleft/
This documentation is licensed under the terms of the GNU Free Documentation License. See
http:/www.gnu.org/copyleft/
<[email protected]>. Para cualquier consulta sobre este texto en forma electrónica se
pide que el subject o sujeto sea: “amII”. Esta edición electrónica se ha realizado íntegramente con
software libre, mediante el procesador Lyx.
Imagen de tapa:
Función de ejemplo ’Trap’ dentro de las "IsoSurface" del software K3DSurf v0.6 (GPL). Editada
usando GIMP 2.4.0-rc3 ("GNU Image Manipulation Program").
Índice general
iii
CAPíTULO 1
1.1. Objetivo
Este texto intenta resolver ejercicios típicos de la cátedra de Análisis Matemático II para ayudar
en el aprendizaje de esta materia básica para Ingeniería.
Espera ayudar a los alumnos que están en la cátedra o por rendir el final de la materia. Los
ejercicios que se presenta son extraídos directamente de los prácticos que generalmente se presentan.
La resolución de los siguientes ejercicios fue en gran parte realizada por el autor y en algunos casos
especiales se necesito de ayuda del JTP (Jefe de Trabajos Prácticos) puesto que en el momento de
resolverlos como es costumbre el alumno no siempre tiene las “herramientas” de resolución pulidas y
cada tanto algunos casos en esos momentos no pudieron resolverse.
La idea en estas páginas, es ayudar al lector que supongo ha de estar en la misma situación que
estuvo el autor al cursar la materia, facilitar la interpretación de los ejercicios es fundamental puesto
que en los momentos en donde se cursan la materia la visión del alumno es un poco más “reducida”,
puesto que entre los estudiantes en primera instancia no son capaces de contemplar todo el panorama
y como los Profesores y Ayudantes ya poseen un conocimiento más amplio que a veces apabullan, con
lo cual es difícil comprender del todo las técnicas y definiciones.
No se intenta evadir de ninguna forma el aporte del estudio formal en la cátedra, sino ampliar e
intentar llegar a mejores interpretaciones de la teoría y práctica a través de ejercicios resueltos como
“recetas”. ¿Cuantas veces se tiene ese click en la mente al ver un ejercicio resuelto o sólo al ver una
parte?
1.2. Metodología
Los ejercicios están resueltos de la forma más lógica teniendo en cuenta de cómo se suelen resolver
entre varias personas.
Cuando se deba resolver sistemas de ecuaciones se resolverán por el método de Gauss-Jordan, es
el método con que el autor se siente más cómodo, el lector podrá evadir (si quiere) esta parte aunque
la resolución siempre será desarrollada.
En el último capitulo del texto se destinará para los cambios realizados sobre el mismo para poder
llevar un control más preciso y que el lector que utilice este material pueda visualizar más cómodamente
el progreso y mejoras como también los errores que se estén dando.
1.3. Consideraciones
1.3.1. Enunciados de los ejemplos. Se opto por introducir al principio de los enunciados un
código sencillo y con colores para tener una primera aproximación sobre la dificultad del ejercicios y
si se resolvió a través de caminos alternativos.
El código es como sigue:
[1] = sencillo
[2] = intermedio
[3] = difícil
[4] = asistido con profesor
Se le concatenara la letra “a” cuando la resolución sea a través de más de una alternativa:
[3a] = difícil con resolución alternativa
A veces algunos ejercicios se resolverán en parte suponiendo los resultados parciales como generales
y dejando al lector la resolución complementaria para la práctica, la forma de indicar esto es utilizando
la letra “p” como por ejemplo:
[1p] = sencillo con resolución parcial
1
1.3. CONSIDERACIONES 2
Algunos ejercicios son problemas de Física por lo que aunque su resolución en general no variara
pero si se indicará con la letra “F”:
[2F] = intermedio, problema Física
En ciertas ocasiones, algunos de los ejemplos necesitan de cálculo numérico para ser resueltos, en
esos casos el siguiente icono indicara el uso de software de cálculo especializado según el ejercicio
a resolver.
1.3.2. Referencias pasadas. Cuando se necesite realizar referencias sobre resoluciones anteri-
ores se etiquetara la zona a tomar nota con un icono y se hará la referencia a la página donde se
esta requiriendo esta marca.
1.3.3. Función trigonométrica. Por cuestiones técnicas a la función seno se la denominara
“sin” y no como “sen”.
1.3.4. Función trigonométrica inversa. Se adopta la siguiente convención:
Ecuaciones Diferenciales
2.1. Teoría
2.1.1. Ecuaciones Lineales.
(
lineales
ordinarias
diferenciales no lineales
parciales
Integrar
R
P (x)dx
e y=C
Despejamos y
R
y = Ce− P (x)dx
y 0 + P (x) y = Q (x)
R
P (x)dx
e = µ (x)
Realizando la sustitución
µ (x) y 0 + P (x) µ (x) y = Q (x) µ (x) y 0
0
(µ (x) y) = Q (x) µ (x)
Integrar
Z
µ (x) y = Q (x) µ (x) dx
Despejamos y
R
Q (x) µ (x) dx
y=
µ (x)
Por lo que
R Z R
y = e− P (x)dx
Q (x) e P (x)dx
dx
dy
= f(x;y)
dx
h(y) dy = g(x) dx
dy g(x)
y0 = =
dx h(y)
.
Homogéneas. Todos los términos tienen el mismo grado
F(x;y) = y 0
√ 1
m em
dy = √ dm
2 m
√
simplificamos m
1 m
dy = e dm
2
sólo resta integrar
1
Z Z
dy = em dm
2
1
y = em + c
2
Volveremos a sustituir para volver a la variable original para obtener la solución:
1 x2
y= e +c
2
y 2 + x2 = c
2.2. EJEMPLOS DE ECUACIONES DIFERENCIALES 6
y y0 + x = 0
x
p
c− x2 −√ +x=0
c − x2
De la solución particular el valor de c√es siempre≥ 0 y además que c ≥ x2 , por lo que resulta de
este análisis la simplificación del término c − x2 :
−x + x = 0
Por lo que resta verificar que la solución particular dada es correcta. Se necesita hallar el valor de
c que satisfaga las condiciones iniciales
√ √ 2
r
2= c− 2
√ √
2= c−2
2=c−2
c=4
Figura 2.2.1.
Figura 2.2.2.
3y 00 + 3y 0 − 4y = 0
Desarrollo 4. Para poder determinar el valor de r se necesita sustituir el valor de la variable
y, pero en primera instancia se ha de derivar y dos veces
y = erx
y 0 = r erx
y 00 = r2 erx
se pueden realizar las sustituciones correspondientes
3 r2 erx + 3 (rerx ) − 4 (erx ) = 0
factor común erx
2.2. EJEMPLOS DE ECUACIONES DIFERENCIALES 8
erx 3r2 + 3r − 4 = 0
Esta ecuación tendrá solución cuando alguno de los dos factores o ambos sea cero, del análisis
erx 6= 0
por lo que el primer factor nunca puede anularse entonces el segundo factor es el único factor que
podría anularse
3r2 + 3r − 4 = 0
utilizando la resolvente
p
(−3) ± (32 ) − [(4) (3) (−4)]
r1−2 =
(−2) (3)
se obtiene
r1 ≈ −0,7583
r2 ≈ 1,7583
Como existen dos raíces reales para esta ecuación, los valores calculados de r1 y r2 verifican la
ecuación diferencial original.
Ejemplo 5. [1] Mediante separación de variables resolver la siguiente ecuación diferencial teniendo
en cuenta la condición inicial y(0) = 1:
y 0 + 5x4 y 2 = 0
Desarrollo 5. Nuevamente se comenzara agrupando convenientemente las variables
dy
+ 5x4 y 2 = 0
dx
dividiendo por y 2 y multiplicando por dx
dy
+ 5x4 dx = 0
y2
dy
= −5x4 dx
y2
se procederá a integrar
1
Z Z
dy = −5 x4 dx
y2
1
reescribir y2
en forma de exponente
Z Z
−2
y dy = −5 x4 dx
integrar directamente
x5
−y −1 = −5 +c
5
1
= −x5 + c
−
y
Reescribir dejando en función de x como generalmente se trabaja
1
− =y
c − x5
2.2. EJEMPLOS DE ECUACIONES DIFERENCIALES 9
Se analizara este resultado a través de las condiciones iniciales para hallar el valor de c:
1
1=−
c − 05
1
1=−
c
c = −1
y 0 = 4 y + 1 cos (2x)
p
dy p
= 4 y + 1 cos (2x)
dx
1
√ dy = 4 cos (2x) dx
y+1
Para realizar una integración más sencilla se realizara un cambio de variables
m=y+1
como en cualquier cambio de variables también se necesita hallar los diferenciales pertinentes
dm = dy
Se realiza la sustitución, se escribe la raíz en forma de exponente y se integra
Z Z
1
m− 2 dm = 4 cos (2x) dx
1
2m− 2 = −4 sin (2x) 2 + c
1
2m− 2 = −8 sin (2x) + c
Resuelta la integral se vuelve a las variables originales:
1
2 (y + 1)− 2 = −8 sin (2x) + c
reescribir en función de x e y a la constante c, una técnica muy utilizada cuando la separación de
variables no es viable en primera instancia
1
2 (y + 1)− 2 + 8 sin (2x) = c
De la condición inicial se intenta obtener el valor de c:
1
− 12
2 (−1 + 1) + 8 sin 2 π =c
4
1
8 sin π =c
2
Ejemplo 7. [1] Resolver la siguiente ecuación homogénea usando separación de variables susti-
tuyendo u = xy :
x2 y 0 = y 2 + xy + x2
2.2. EJEMPLOS DE ECUACIONES DIFERENCIALES 10
dy
x2 = y 2 + xy + x2
dx
dividimos por x2
dy y2 y
= 2 + +1
dx x x
En esta última expresión notamos que se repite la forma de la sustitución sugerida
y
u=
x
como cualquier sustitución requiere de el cálculo de la derivadas correspondientes
ux = y
derivando factores
du dx dy
x+u =
dx dx dx
u0 x + u = y 0
con este último resultado la sustitución ya es posible
dy y2 y
= 2 + +1
dx x x
xu0 + u = u2 + u + 1
cancelamos los términos u
xu0 = u2 + 1
el próximo paso es la separación de variables
du
x = u2 + 1
dx
1 1
du = dx
u2 + 1 x
ya es posible la integración
1 1
Z Z
2
du = dx
u +1 x
La integración es directa:
y
c = atan − ln (x)
x
Desarrollo 8. Se aplicaran dos caminos para la resolución, se comienza utilizando los despejes
normales:
x2 y 2 = u
r
u
y=
x2
√
u
y=
x
derivar y
√
1
√
2 u
u0 x − u
y0 =
x2
Se realizará la sustitución y se trabajará algebraicamente
√
√1
2 u
u0 x − u √ u
2
2x u = tan (u) − 2
x x
x2 u0 − 2xu u
2
= tan (u) − 2
x x
x2 u0 2xu u
− 2 = tan (u) − 2
x2 x x
u u
u0 − 2 = tan (u) − 2
x x
cancelar lo que resulta en
u0 = tan (u)
agrupar
du
= tan (u)
dx
du
= dx
tan (u)
1
sin(u)
du = dx
cos(u)
cos (u)
du = dx
sin (u)
integrar
cos (u)
Z Z
du = dx
sin (u)
para este cálculo una buena opción es la realización de una nueva sustitución
z = sin (u)
dz = cos (u) du
sustituir e integrar
1
Z Z
dz = dx
z
ln (z) + c = x
se volverá a las variables originales z = sin (u) por lo tanto z = sin x2 y 2
2.2. EJEMPLOS DE ECUACIONES DIFERENCIALES 12
ln sin x2 y 2 +c=x
ln sin x2 y 2 =x−c
x − c = ln sin x2 y 2
con las propiedades de logaritmos se trabajara
2 y2
ex−c = eln(sin(x ))
ex
2 2
= sin x y
ec
1
como ec sigue siendo una constante la llamaremos C
sin x2 y 2
C=
ex
En este ejercicio se puede optar por una resolución diferente que logrará el mismo resultado, empezando
por derivar y 2 :
x2 y 2 = u
u
y2 = 2
x
derivar
u0 x2 − u2x
2yy 0 =
x4
u0 x2 − u2x
y0 =
2yx4
Ahora si se reemplaza
y 0 2x2 y = tan x2 y 2 − 2xy 2
!
u0 x2 − u2x
2yx2 = tan (u) − 2xy 2
2yx4
simplificar el primer miembro
u0 x2 − u2x
= tan (u) − 2xy 2
x2
multiplicar por x2
u0 x2 − 2ux = tan (u) x2 − 2x3 y 2
terminar de sustituir
u0 x2 − 2ux = tan (u) x2 − 2ux
cancelar
u0 x2 = tan (u) x2
dividir por x2
u0 = tan (u)
Notar que se ha llegado a la misma forma que en la página anterior por lo que el desarrollo continua
como se realizo anteriormente.
Ejemplo 9. [1] Resolver usando la sustitución u = y − 3:
q
y0 = 6y − y 2 − 8
Desarrollo 9. Se comienza despejando y de la sustitución
2.2. EJEMPLOS DE ECUACIONES DIFERENCIALES 13
u+3=y
derivamos para poder realizar la sustitución
u0 = y 0
reemplazar
q
y0 = 6y − y 2 − 8
q
u0 = 6 (u + 3) − (u + 3)2 − 8
distributiva
q
u0 = 6u + 18 − (u + 3) (u + 3) − 8
trabajo algebraico
q
u0 = 6u + 10 − (u2 + 3u + 3u + 9)
q
u0 = 6u + 10 − (u2 + 6u + 9)
cancelar 6u
p
u0 = 6u + 10 − u2 − 6u − 9
p
u0 = 1 − u2
realizar la separación de variables
du p
= 1 − u2
dx
1
√ du = dx
1 − u2
integramos directamente
1
Z Z
√ du = dx
1 − u2
asin (u) + c = x
reemplazamos u = y − 3
asin (y − 3) + c = x
por último despejamos c por la dificultad que presenta separar x e y:
c = x − asin (y − 3)
y−x+1
y0 =
y−x+5
Desarrollo 10. Se comienza despejando y como normalmente se realiza
2.2. EJEMPLOS DE ECUACIONES DIFERENCIALES 14
y−x=u
y =u+x
también se deriva correspondientemente
y 0 = u0 + 1
Se comienza con el trabajo de sustitución sobre la ecuación:
y−x+1
y0 =
y−x+5
u+1
u0 + 1 =
u+5
se empieza el trabajo algebraico
u+1
u0 = −1
u+5
común denominador (u + 5)
u+1−u−5
u0 =
u+5
−4
u0 =
u+5
se puede realizar la separación de variables
du −4
=
dx u+5
(u + 5) du = −4dx
Ya es posible la integración:
Z Z
(u + 5) du = −4 dx
Z Z Z
u du + 5 du = −4 dx
u2
+ 5u = −4x + c
2
reemplazar u por y − x
(y − x)2
+ 5 (y − x) = −4x + c
2
puesto que la separación de x e y no es posible se despejará c
(y − x)2
+ 5 (y − x) + 4x = c
2
(y − x)2
+ 5y − 5x + 4x = c
2
(y − x)2
+ 5y − x = c
2
Ejemplo 11. [1] Resolver la siguiente ecuación diferencial, comprobar de que tipo es:
6xy − y 3 dx + 4y + 3x2 − 3xy 2 dy = 0
2.2. EJEMPLOS DE ECUACIONES DIFERENCIALES 15
Desarrollo 11. Sin realizar cálculos es posible ver que la separación de variables no va ha ser
posible, revisemos a ver si la ecuación es del tipo diferencial total:
M (x;y) = 6xy − y 3 → My0 = 6x − 3y 2
N (x;y) = 4y + 3x2 − 3xy 2 → Nx0 = 6x − 3y 2
Como las derivadas parciales cruzadas son iguales entonces la ecuación diferencial dada es exacta, por
lo tanto procedemos a resolver, por simpleza de calculo tomamos M (x;y) :
Z Z
M (x;y) dx = 6xy − y 3 dx
Fy = 3yx2 − y 3 x + g(y)
Derivar en y
Fy0 = 3yx2 − y 3 x + g(y) dy
0 , igualando
Calcular g(y)
Fy0 = N(x;y)
0
g(y) = 4y
integrando con respecto a y Z Z
0
g(y) =4 (y) dy
g(y) = 2y 2 + c
Por lo que la función potencial es:
F(x;y) = 3x2 y − xy 3 + 2y 2 + c
Ejemplo 12. [1] Resolver la siguiente ecuación diferencial, comprobar de que tipo es:
(4y) dx + (x) dy = 0
Desarrollo 12. Se tratara como una ecuación diferencial exacta:
M (x;y) = 4y → My0 = 4
N (x;y) = x → Nx0 = 1
La ecuación no es exacta pero las derivadas parciales cruzadas sólo difieren en una constante, es
posible que la ecuación sea reducible a exacta. Se intentara hallar el factor integrante:
My0 − Nx0
f(x) =
N(x,y)
4−1
f(x) =
x
3
f(x) =
x
Como f(x) sólo esta dependiente de una sola variable y es x entonces se ha encontrado un factor
integrante, procedemos a su cálculo
R
f(x) dx
ñ(x) = e
R 1
dx
ñ(x) = e3 x
2.2. EJEMPLOS DE ECUACIONES DIFERENCIALES 16
la integral es directa
ñ(x) = e3 ln x
por propiedades de logaritmos al factor 3 se lo transforma en un exponente
3
ñ(x) = eln x
nuevamente se aplicamos propiedades de los logaritmos
ñ(x) = x3
Se ha de introducción el factor integrante en la ecuación diferencial original para transformarla en
exacta:
F(x;y) = M(x;y)
0
g(y) =0
integrar
Z Z
0
g(y) dy = (0) dy
g(y) = c
La función potencial es:
F(x;y) = x4 y + c
y 00 + 6y 0 + 9y = 0
Desarrollo 13. Realmente sencilla de resolver, por ser homogénea solo debemos hallar la solución
complementaria:
2.2. EJEMPLOS DE ECUACIONES DIFERENCIALES 17
r2 + 6r + 9 = 0
aplicando resolvente
r = −3
La ecuación anterior tiene raíces reales y coincidentes por lo que la solución complementaria es:
YC = (c1 + c2 x)e−3x
0 = 11:
Ejemplo 14. [1] Resolver con las condiciones iniciales y(0) = 7 e y(0)
y 00 − 4y 0 + 3y = 0
Desarrollo 14. En primera instancia se comenzará calculando la solución complementaria:
r2 − 4r + 3 = 0
aplicando resolvente
(
r1 = 3
r2 = 1
Como las raíces son reales y distintas la solución complementaria es:
YC = c1 e3x + c2 ex
Se trabajara con las condiciones iniciales dadas para hallar el valor de c1 y c2 , primero es necesario
hallar las correspondientes derivadas de la solución complementaria
0
YC0 = c1 e3x + c2 ex
0 = 2:
Ejemplo 15. [2] Resolver con las condiciones iniciales y(0) = 1 y y(0)
y 00 − 3y 0 + 2y = 3e−x
1 0 0
−5 6 3
se toma como pivote la posición (1, 1) y se aplica el método
1 0 0
0 36 18
se divide la segunda fila por 36
1 0 0
0 1 21
las soluciones del sistema son
(
B=0
C = 12
Por lo que la primer solución particular tiene la siguiente forma:
1
YP1 = e−x
2
Se procede con genera la segunda parte de la solución particular, en general al existir funciones
trigonométricas (seno o coseno) se trabajará como sigue:
YP2 = A sin (3x) + D cos (3x)
hallamos la primer derivada
−9A sin (3x) − 9D cos (3x) − 9D cos (3x) + 9D sin (3x) + 2A sin (3x) + 2D cos (3x) = −10 cos (3x)
se procede algebraicamente agrupando convenientemente
(−7A + 9D) sin (3x) + (−9A − 7D) cos (3x) = −10 cos (3x)
puesto que la solución de la ecuación diferencial requiere de la función trigonométrica coseno se
debe hallar los coeficientes de la función seno para anularla, se procede a armar el siguiente sistema
de ecuaciones
(
−7A + 9B =0
−9A − 7B = −10
Resolver el sistema de ecuaciones2, se vuelve a recurrir al método Gauss-Jordan,:
−7 9
0
−9 −7 −10
dividir la primer fila por −7
1 − 79
0
−9 −7 −10
el pivote a usar será el de la posición (1, 1)
2El lector puede utilizar cualquier método para la resolución.
2.2. EJEMPLOS DE ECUACIONES DIFERENCIALES 20
" #
1 − 79 0
0 − 130
7 −10
7
se multiplica la segunda fila por − 130
" #
1 − 79 0
7
0 1 13
la posición (2, 1) será el nuevo pivote
" #
9
1 0 13
7
0 1 13
las soluciones del sistema presentado son:
(
9
A = 13
7
D = 13
La forma que adopta la segunda parte de la solución particular es:
9 7
YP2 = sin (3x) + cos (3x)
13 13
Sólo resta agrupar las soluciones y presentar la solución general:
YG = YC + YP1 + YP2
1 9 7
YG = c1 e2x + c2 ex + e−x + sin (3x) + cos (3x)
2 13 13
Luego de obtener YG restaría hallar los valores de c1 y c2 que estén en concordancia con las
condiciones iniciales:
(
y(0) = 1
0 =2
y(0)
(
YC = c1 e2x + c2 ex
YC0 = 2c1 e2x + c2 ex
reemplazando
(
1 = c1 + c2
2 = 2c1 + c2
Se puede resolver el sistema de ecuaciones, se utiliza el método de Gauss-Jordan:
1 1 1
2 1 2
el pivote a usar es la posición (1, 1)
1 1 1
0 −1 0
dividir la segunda fila por −1
1 1 1
0 1 0
la única posición que se puede utilizar como pivote es (2, 1)
1 0 1
0 1 0
las soluciones son:
2.2. EJEMPLOS DE ECUACIONES DIFERENCIALES 21
(
c1 = 1
c2 = 0
Por lo que la solución general con las condiciones iniciales dadas es:
1 9 7
YG = e2x + e−x + sin (3x) + cos (3x)
2 13 13
x
Ejemplo 16. [1] Verificar que la solución particular YP = 6 ex sea solución de la ecuación y
encontrar la solución general:
ex
y 00 − 2y 0 + y = 12
x
Desarrollo 16. Comenzar calculando la primer derivada de YP :
6ex x − 6ex
YP0 =
x2
también calcular la segunda derivada
6ex x3 − 12ex x2 + 12ex x
YP00 =
x4
Reemplazar en la ecuación diferencial dada:
!
6ex x3 − 12ex x2 + 12ex x 6ex x − 6ex ex ex
−2 +6 = 12
x4 x2 x x
distributiva
ex ex ex ex ex ex ex
6 − 12 2 + 12 3 − 12 + 12 2 + 6 = 12
x x x x x x x
cancelar los términos correspondientes
ex ex
12 = 12
x x
Se puede concluir que la solución particular dada es solución de la ecuación diferencial. Se intentara
hallar la solución general, comenzando por al solución complementaria:
r2 − 2r + 1 = 0
a través de la resolvente
r=1
por existir raíces reales coincidentes la solución complementaria es
YC = (c1 + c2 x) ex
Con este último cálculo se puede proponer la solución general:
ex
YG = (c1 + c2 x) ex + 6
x
Ejemplo 17. [2] Hallar la solución general de la siguiente ecuación diferencial de orden superior
teniendo en cuenta las condiciones iniciales y(0) = 1 y y(0) = −3:
y 00 − 4y 0 + 3y = 10e−2x
Desarrollo 17. La resolución comienza por la búsqueda de la solución complementaria:
2.2. EJEMPLOS DE ECUACIONES DIFERENCIALES 22
r2 − 4r + 3 = 0
resolviendo usando la resolvente
(
r1 = 3
r2 = 1
por lo que la solución complementaria sería
YC = c1 e3x + c2 ex
Se propone una solución particular para la forma 10e−2x :
YP = Bxe−2x + Ce−2x
hallando la primera derivada
y 00 − 4y 0 + 3y = 10e−2x
10 −2x
YG = c1 e3x + c2 ex + e
7
CAPíTULO 3
x4 − y 4
f(x;y) =
x2 − y 2
Desarrollo 18. De la expresión anterior notamos que cuando el denominador se hacer cero es
indeterminada la función por lo que tenemos que descubrir para que valores sucede ello:
x2 − y 2 6= 0
Es una diferencia de cuadrados:
(x − y) (x + y) 6= 0
Por lo tanto:
x 6= y
x 6= −y
Ya podemos dibujar el dominio:
Figura 3.1.1.
Nota 1. Para realizarlo utilizando como en la Figura el software Maxima y usando como graficador
el GnuPlot el código sería:
F(x):=x;
H(x):=-x;
plot2d([F,H],[x,-2,2],[y,-2,2],[style,[linespoints]]);
24
3.1. EJEMPLOS DE FUNCIONES DE DOS VARIABLES 25
x4 − y 4
0=
x2 − y 2
x2 − y 2 x2 + y 2
0=
x2 − y 2
0 = x2 + y 2
y 2 = −x2
z = y2
Figura 3.1.2.
Nota 2. Para usar con el maxima:
plot2d(x^2,[x,-2,2],[y,0,2],[ylabel,"z=y^2"],[xlabel,"y"],[legend,"Traza x=0"]);
Ahora cuando y = 0:
x4
z=
x2
z = x2
3.1. EJEMPLOS DE FUNCIONES DE DOS VARIABLES 26
Figura 3.1.3.
Ahora procederemos a calcular algunas curvas de nivel para mayor comprensión de la función,
empecemos con z = 1:
x4 − y 4
1=
x2 − y 2
1 = x2 + y 2
1 − x2 = y 2
p
± 1 − x2 = y
Con z = 2 :
x4 − y 4
2=
x2 − y 2
2 = x2 + y 2
p
± 2 − x2 = y
Con z = 3:
x4 − y 4
3=
x2 − y 2
p
± 3 − x2 = y
3.1. EJEMPLOS DE FUNCIONES DE DOS VARIABLES 27
Figura 3.1.4.
Nota 3. Para usar con el Maxima y GnuPlot:
plot2d([sqrt(1-x^2),sqrt(2-x^2),sqrt(3-x^2),-sqrt(1-x^2),-sqrt(2-x^2),-sqrt(3-x^2)],[x,-2,2],[y,-2,2]);
Figura 3.1.5.
Nota 4. Para reproducirlo usando Maxima con OpenMath usar:
plot3d((x^4-y^4)/(x^2-y^2),[x,-2,2],[y,-1,1],[plot_format, openmath]);
x2 − y 2
f(x;y) =
x2 + y 2
3.1. EJEMPLOS DE FUNCIONES DE DOS VARIABLES 28
Desarrollo 19. En primera instancia se revisa lo que sucede con la sustitución directa:
lı́m x2 y 2 0
=
(x;y)→(0;0) x3 + y 3 0
Puesto que la sustitución directa presenta una forma indeterminada como solución se comenzará
acercándose a través de infinitas rectas de la forma y = mx:
lı́m x2 m 2 x2
(x;mx)→(0;0) x3 + m3 x3
Factor común x2 :
lı́m x4 m2
(x;mx)→(0;0) x3 (1 + m3 )
simplificamos x puesto que x → 0 por lo que nunca vale 0 así que este paso es válido
lı́m m2
(x;mx)→(0;0) 1 + m3
notemos que mx → 0, por lo que cualquiera de los dos factores pueden acercarse a 0 tanto como
se quiera por lo que:
lı́m m2
=0
(x;mx)→(0;0) 1 + m3
Ahora nos acercaremos a 0 a través de x2 :
lı́m x2 x4
(x;x2 )→(0;0) x3 + x6
lı́m x6
(x;x2 )→(0;0) x3 + x6
lı́m x6
(x;x2 )→(0;0) x3 (1 + x3 )
simplificamos x
lı́m x2
(x;x2 )→(0;0) 1 + x3
como x → 0 entonces:
lı́m x2
=0
(x;x)→(0;0) 1 + x3
Concluiremos que:
lı́m x2 y 2
=0
(x;y)→(0;0) x3 + y 3
El lector podrá tomar la tarea de resolver con las sustituciones correspondientes que quedan como
práctica.
~ = cos π π
Ejemplo 20. [1] Hallar la derivada direccional en la dirección del vector U 4 î + sin 4 ĵ:
f(x,y) = 2x2 + 5y 2
3.1. EJEMPLOS DE FUNCIONES DE DOS VARIABLES 29
kU ~ k= 1
Ya es posible armar el gradiente de la función
~ (x,y) = 4x î + 10y ĵ
∇f
→
−
Por último se puede calcular la derivada direccional en la dirección del vector U
√ √ !
2 2
DU~ f(x;y) = 4x î + 10y ĵ · î + ĵ
2 2
√ √
DU~ f(x;y) = 2 2 x + 5 2 y =
√
DU~ f(x;y) = 2 (2x + 5y)
Ejemplo 21. [2F] La temperatura en cualquier punto (x; y) de una placa rectangular está dada
por
T(x;y) = x2 + y 2
Calcular la tasa de variación de la temperatura en el punto P = (3; 4) en la dirección que forma
un ángulo θ de valor
π
θ= = 60ř
3
con el eje de las abscisas.
Determinar el ángulo de la dirección en que la tasa de variación de la temperatura en el punto
(−3; 1) sea un máximo.
Desarrollo 21. Se comienza calculando las derivadas parciales directamente sin ser necesario
cálculos previos:
∂T(x;y) ∂T(x,y)
= 2x ; = 2y
∂x ∂y
Esto permite hallar la forma del gradiente:
→
−
∇T(x;y) = 2x î + 2y ĵ
Ahora se valorizará el gradiente en el punto (3; 4)
~ (3;4) = 6 î + 8 ĵ
∇T
~ en función del ángulo que se nos da como dato
Se calculara el vector U
~ = cos π π
U î + sin ĵ
3 3
√
~ 1 3
U = î + ĵ
2 2
3.1. EJEMPLOS DE FUNCIONES DE DOS VARIABLES 30
~ kU
DU~ T(x;y) = ∇· ~ k · cos α
primero, notemos que el valor del módulo de U ~ es
kU~ k= 1
segundo, el factor que varía es el cos α pero por se una función trigonométrica básica esta varía
entre −1 ≤ cos α ≤ 1, para que el factor permite un valor máximo:
cos α = 1
entonces se puede proceder al cálculo del gradiente valorizado en el punto (−3; 1):
~ (−3;1) = −6 î + 2 ĵ
∇T
Sólo resta hallar el ángulo que se forma con el eje de las x, por ser un vector podemos utilizar las
componentes para hallarlo:
2 1
atan − = atan − = −18, 4ř
6 3
Ejemplo 22. [1] Calcular el gradiente en el punto (2; 5) y también la derivada direccional en la
dirección del vector V = (3; −3) de la siguiente función:
f(x;y) = x2 − 5xy
Desarrollo 22. Se calcularan las derivadas parciales directamente por la forma polinómica que
presenta:
∂f(x;y) ∂f(x;y)
= 2x − 5y ; = −5x
∂x ∂y
~ (gradiente) de la función:
Ahora se presentará el ∇
~ (2;5) = −21 î − 10 ĵ
∇f
3.1. EJEMPLOS DE FUNCIONES DE DOS VARIABLES 31
Ahora se calculará el valor del vector que nace en el punto (2; 5) y se dirige a (3; −3):
~ ) = (3 − 2) î + (−3 − 5) ĵ = î − 8 ĵ
V(RQ
Calcularemos el módulo de este vector para transformar el vector U en unitario o en mejores
términos en un versor
p √
M od = 12 + 82 = 65
~ = √1 î − √8 ĵ =⇒k U
U ~ k= 1
65 65
Por lo que ya es posible el cálculo de la derivada direccional requerida
~ · ∇f
DU~ f(2;5) = U ~ (2;5)
1 8 21 80 101
DU~ f(2;5) = √ î − √ ĵ · −21 î − 10 ĵ = − √ − √ = − √
65 65 65 65 65
≈ 12, 5275
Ejemplo 23. [1] Determinar los extremos relativos de la siguiente función:
f(x;y) = 2x4 + y 2 − x2 − 2y
Desarrollo 23. Iniciamos calculando las derivadas parciales, como la función es tipo polinómica
no es necesario ningún trabajo algebraico previo:
∂f(x;y) ∂f(x;y)
= 8x3 − 2x ; = 2y − 2
∂x ∂y
∂ 2 f(x;y) 2 ∂ 2 f(x;y)
= 24x − 2 ; =2
∂x2 ∂y 2
∂f(x;y) ∂f(x;y)
=0; =0
∂x∂y ∂y∂x
Como las primeras derivadas parciales son exclusivas cada una de la variable derivada entonces se
puede calcular los ceros por separado, con respecto a x:
8x3 − 2x = 0 =⇒ x 8x2 − 2 = 0
solución
x = 0
1
x= 2
x = − 1
2
luego la derivada con respecto a y
2
2y − 2 = 0 =⇒ y =
2
solución
n
y=1
De las combinaciones de las soluciones anteriores obtendremos los puntos críticos
1 1 1
(0; 1) ;
; ; − ;1
2 2 2
Ahora se calculan el Helessiano para cada punto crítico.
Calculamos | H(0;1) |:
3.1. EJEMPLOS DE FUNCIONES DE DOS VARIABLES 32
−2 0
| H(0;1) | = = −4 < 0 ∴
0 2
4 0
| H( 1 ;1) |= =8>0
2 0 2
1
y como 4 > 0 =⇒ ∃ mR en 2
;1 .
Calculamos | H(− 1 ;1) |:
2
4 0
| H(− 1 ;1) |= =8>0
2 0 2
y como 4 > 0 =⇒ ∃ mR en − 12 ; 1 .
Figura 3.1.6.
Nota 5. Usando Maxima:
plot3d((2*(x^4))+(y^2)-(x^2)-(2*y), [x,-1,1], [y,0.5,1.5], [gnuplot_pm3d, true]);
3.1. EJEMPLOS DE FUNCIONES DE DOS VARIABLES 33
Figura 3.1.7.
Nota 6. Usando Maxima:
plot3d((2*(x^4))+(y^2)-(x^2)-(2*y), [x,-1,1], [y,0.5,1.5], [gnuplot_pm3d, true]);
Figura 3.1.8.
Nota 7. Usando Maxima:
plot3d((2*(x^4))+(y^2)-(x^2)-(2*y), [x,-1,1], [y,0.5,1.5]);
f(x;y) = x3 y 2 (1 − x − y)
Desarrollo 24. Se comenzará con las derivadas parciales de la función, convenientemente en
primer lugar aplicaremos distributiva y luego calcularemos:
f(x;y) = x3 y 2 − x4 y 2 − x3 y 3
∂f(x;y) ∂f(x;y)
= 3x2 y 2 − 4x3 y 2 − 3x2 y 3 ; = 2yx3 − 2yx4 − 3y 2 x3
∂x ∂y
3.1. EJEMPLOS DE FUNCIONES DE DOS VARIABLES 34
∂ 2 f(x;y) 2 2 2 3 ∂ 2 f(x;y)
= 6xy − 12x y − 6xy ; = 2x3 − 2x4 − 6yx3
∂x2 ∂y 2
∂f(x;y) ∂f(x;y)
= 6x2 y − 8x3 y − 9x2 y 2 ; = 6x2 y − 8x3 y − 9x2 y 2
∂x∂y ∂y∂x
Armamos el sistema de ecuaciones para hallar los puntos críticos:
∂f(x;y) = 0
( (
∂x 3x2 y 2 − 4x3 y 2 − 3x2 y 3 = 0 x2 y 2 (3 − 4x − 3y) = 0
→ →
∂f(x;y) = 0 2yx3 − 2yx4 − 3y 2 x3 = 0 yx3 (2 − 2x − 3y) = 0
∂y
Notamos dos grupos de factores en cada ecuación, con lo cuál deducimos por ahora que:
(
x=0
y=0
ahora el segundo sistema
( (
3 − 4x − 3y = 0 4x + 3y = 3
→
2 − 2x − 3y = 0 2x + 3y = 2
Se resuelve el sistema de ecuaciones:
4 3 3
2 3 2
dividimos la primer fila por 4
1 43 3
4
2 3 2
el pivote será la posición (1; 1)
" #
3 3
1 4 4
3 1
0 2 2
2
multiplicamos por 3 la segunda fila
" #
1 43 3
4
1
0 1 3
el nuevo pivote será (2; 2)
" #
1
1 0 2
1
0 1 3
1 1
∴x= ; y=
2 3
Los puntos críticos que se hallan las siguientes combinaciones:
1 1 1 1
(0; 0) ; 0; ; ;0 ; ;
3 2 2 3
Ahora se calculará el Helessiano para cada punto crítico.
Se calcula | H(0;0) |:
0 0
| H(0;0) |= =0
0 0
no aporta información.
Se calcula | H(0; 1 ) |:
3
0 0
| H(0; 1 ) |= =0
3 0 0
no aporta información.
3.1. EJEMPLOS DE FUNCIONES DE DOS VARIABLES 35
Se calcula | H( 1 ;0) |:
2
0 0
| H( 1 ;0) |= =0
2 0 0
no aporta información.
Se calcula | H( 1 ; 1 ) |:
2 3
Figura 3.1.9.
Nota 8. Usando Maxima:
plot3d((x^3)*(y^2)*(1-x-y), [x,-1,1], [y,-1,1], [gnuplot_pm3d,true]);
3.1. EJEMPLOS DE FUNCIONES DE DOS VARIABLES 36
Figura 3.1.10.
Nota 9. Usando Maxima:
plot3d((x^3)*(y^2)*(1-x-y), [x,-1,1], [y,-1,1], [gnuplot_pm3d,true]);
Figura 3.1.11.
Nota 10. Usando Maxima:
plot3d((x^3)*(y^2)*(1-x-y), [x,-5,5], [y,-5,5], [gnuplot_pm3d,true]);
Así del estudio del gráfico se deduce que la función no posee puntos extremos.
Ejemplo 25. [1] Calcular si las derivadas parciales cumplen con
∂2u ∂2u
+ 2 =0
∂x2 ∂y
de la siguiente función
∂u
= −e−x cos y − −e−y sin x = e−y sin x − e−x cos y
∂x
∂u
= −e−x sin y − −e−y cos x = e−y cos x − e−x sin y
∂y
∂2u
= −e−x cos x − −e−y cos y = e−y cos x + e−x cos y
∂x2
∂2u −x −y
= −e−y cos x − e−x cos y
= −e cos x − e cos y
∂y 2
Por lo que:
∂2u ∂2u
+ 2 == 0
∂x2 ∂y
Ejemplo 26. [1] Aplicando la regla de la cadena donde f(x;y) que es dependiente de la tercera
variable t:
f(x;y) = ln x2 + y
con
∂f(x;y) 1 ∂f(x;y) 1
= 2 2x ; = 2
∂x x +y ∂y x +y
dx dy
= − sin t ; = 2e2t
dt dt
df(x;y) 2
2t
= − cos t sin t + e
dt cos2 t + e2t
Ejemplo 27. [1] Hallar la derivada implícita de la expresión:
1
2y 2 + (xy) 3 = 3x2 + 17
Desarrollo 27. Se trabaja algebraicamente
1
2y 2 + (xy) 3 − 3x2 − 17 = 0
Se hallan las derivadas parciales
∂f 1 2 ∂f 1 2
= (xy)− 3 y − 6x ; = 4y + (xy)− 3 x
∂x 3 ∂y 3
Aplicamos la regla de la cadena
3.1. EJEMPLOS DE FUNCIONES DE DOS VARIABLES 38
2
dy ∂f
∂x
1
(xy)− 3 y − 6x
= − ∂f =−3 2
dx ∂y 4y + 1 (xy)− 3 x3
Ejemplo 28. [1] Hallar tres números naturales cuyo producto sea máximo y su suma sea igual a
90.
Desarrollo 28. Presentemos las restricciones:
xyz = M ax = C
y la idea es llegar a
x + y + z = 90
así que podemos plantear los multiplicadores de Lagrange
x=y=z
Por lo que
x + y + z = 90
si tomamos
x=y=z=η
entonces
η + η + η = 90
3η = 90
η = 30
Concluimos que si
η = 30 =⇒ x = y = z = 30
CAPíTULO 4
Integrales
x3 = 2 − 2y
1
y = 1 − x2
2
Figura 4.1.1.
Nota 11. Usando Maxima:
plot2d([1-(0.5*x^2)], [x,0,2], [y,0,1], [gnuplot_preamble, "set grid"])$
como se integra en primer lugar con respecto a y la variable x trabaja como un valor constante
√
Z2" #1− 1 x2
y2 2
I= x dx
2 0
0
√
Z2" 2 #
1 1 2
I= x 1 − x2 − x (0) dx
2 2
0
√
Z2 2
1 1
I= x 1 − x2 dx
2 2
0
√
Z2
1 1
I= x 1 − x2 + x4 dx
2 4
0
√
Z2
1 1
I= x − x3 + x5 dx
2 4
0
distribuimos la integral
√ √ √
Z2 Z2 Z2
1
3 1 5
I= x dx − x dx + x dx
2 4
0 0 0
volvemos a integrar y sólo aparece la variable x, esto es lógico puesto que el diferencial es con
respecto a esta variable y si esto no sucede algún error se ha cometido. Como la función que queda
representada es polinómica la integración es sencilla
√ 2 √ 4 √ 6
2 4 6
1 2 (0) 2 (0) 1 2
(0)
I= − − − + −
2 2 2 4 2 4 6 6
1 1 16 1 8 4
I= 1−1+ = =
2 4 6 2 3 3
4
I=
3
La gráfica de f(x;y) = xy es:
4.1. EJEMPLOS DE INTEGRALES DOBLES 41
Figura 4.1.2.
Nota 12. Usando Maxima y Gnuplot:
plot3d(x*y, [x,-10,10], [y,-10,10], [gnuplot_pm3d,true]);
Figura 4.1.3.
Nota 13. Usando Maxima y Gnuplot:
plot3d(x*y, [x,-10,10], [y,-10,10], [gnuplot_pm3d,true]);
y2 + 1 = x
√
y =± x−1
por lo que
4.1. EJEMPLOS DE INTEGRALES DOBLES 42
Figura 4.1.4.
Nota 14. Utilizando Maxima:
plot2d([sqrt(x-1),-sqrt(x-1), [parametric, 4, t, [t, -3, 3]]], [x,0,5]);
Se comenzará el cálculo dejando constante al dy, por lo que los límites de integración para dx
serán: (
x1 = 4
x0 = y 2 + 1
ahora bien los límites para dy
( √ √
y1 = 4 − 1 = 3
√ √
y0 = − 4 − 1 = − 3
por lo que se puede presentar la forma de la integral
√
Z3 Z4
IA = dx dy
√
− 3 y 2 +1
√
Z3
IA = 3 − y 2 dy
√
− 3
" # √3
y3
IA = 3y −
3 √
− 3
√ 3 √ 3
√ 3 √ − 3
IA = 3 3 − − −3 3 −
3 3
√ 3 √ 3
√ 3 √ 3
IA = 3 3 − +3 3−
3 3
√ 2 √
IA = 6 3 − 3 3
3
√
IA = 4 3
Ahora en cambio dx será constante entonces los límites de dy
4.1. EJEMPLOS DE INTEGRALES DOBLES 43
( √
y1 = x − 1
√
y0 = − x − 1
y los límites de dx serán
(
x1 = 4
x0 = 1
Presentamos y calculamos la integral:
√
Z4 Zx−1
IB = dy dx
√
1 − x−1
Z4 √
x−1
IB = [y]−√ x−1
dx
1
Z4 √
IB = 2 x − 1 dx
1
se realiza la sustitución
µ=x−1
dµ = dx
reemplazamos
Zµ1
1
IB = 2 µ 2 dµ
µ0
" 3 #µ1
µ2
IB = 2 2
3 µ0
se retorna a las variables originales y extremos también
3
4
(x − 1) 2
IB = 2
2
3 1
√
33
IB = 2 2
3
√
3 3
IB = 2 2
3
√
IB = 4 3
Notar que IA = IB , por lo que el cálculo de la integral se realizo correctamente.
CAPíTULO 5
Transformada de Laplace
5.1. Requisitos
5.2. Transformadas indirectas
5.2.1. Trigonométricas.
L {cos (ω0 t + θ)} = L {cos (ω0 t) cos (θ) − sin (ω0 t) sin (θ)}
$ ω0
L {cos (ω0 t + θ)} = cos (θ) 2 − sin (θ) 2
$2 + ω0 $ + ω02
44
CAPíTULO 6
Tablas
6.1. Derivadas
6.1.1. Básicas.
1. Dx (un ) = nun−1 Dx u
2. Dx (u + v) = Dx u + Dx v
3. Dx (uv) = uDx v + vDx u
4. Dx uv = vDx u+uD xv
v2
u u
5. Dx (e ) = e Dx u
6. Dx (au ) = au ln a Dx u
7. Dx (ln u) = u1 Dx u
6.1.2. Trigonométricas.
1. Dx (sin u) = cos u Dx u
2. Dx (cos u) = − sin u Dx u
3. Dx (tan u) = sec2 u Dx u
4. Dx (cot u) = − csc2 u Dx u
5. Dx (sec u) = sec u tan u Dx u
6. Dx (csc u) = − csc u cot u Dx u
6.1.3. Trigonométricas inversas.
1
1. Dx (asinu) = √1−u2
Dx u
1
2. Dx (acosu) = − √1−u2
Dx u
1
3. Dx (atanu) = 1+u 2 Dx u
1
4. Dx (acotu) = − 1+u2 Dx u
5. Dx (asecu) = u√u12 −1 Dx u
6. Dx (acscu) = − u√u12 −1 Dx u
6.1.4. Hiperbólicas.
1. Dx (sinh u) = cosh u Dx u
2. Dx (cosh u) = sinh u Dx u
3. Dx (tanh u) = sinh2 u Dx u
6.2. Integrales
6.2.1. Básicas.
R
1. R du = u + C
2. a
h du = au +i C
R R R
3. f(u) + g(u) du = f(u) du + g(u) du
n+1
4. un du = un+1 + C con (n 6= −1)
R
5. u1 du = ln | u | +C
R
Historial
27/03/10
• Se realiza una revisión ortográfica usando el diccionario aspell.
• Se escribe una mínima parte sobre transformadas indirectas de Laplace.
• Se comienza a experimentar con una presentación mucho más atractiva.
18/01/09
• Se revisan usando ispell 4598 palabras del texto, en general errores de acentuación.
13/08/08
• El texto se licencia bajos los términos de Licencia1 GNU para la Documentación Libre.
Ver en http:/www.gnu.org/copyleft/
24/07/08
• Se adjuntan más gráficos para proveer mayor aporte visual a las interpretaciones.
22/07/08
• Se comienza capítulo sobre integrales.
• Se han adherido casi todas las derivadas más usada en el capítulo de tablas.
• Una mala elección del identificador vectorial (→−
u ) es resuelto por el identificador sopor-
tado para tal fin (~u).
• Información sobre la imagen de tapa se adjunta.
• Se adopta las notas numeradas en cada imagen realizada por sistemas CAD’s para poder
llegar a realizar referencias rápidas cuando se quiera aportar información técnica sobre
el tema.
• Más desarrollo sobre el sistema indicativo de los enunciados de los ejemplos.
• Comienza ha escribirse la bibliografía consultada y el trabajo de recolección de datos.
18/07/08
• Mejor redacción de los textos con mayor coherencia y cohesión.
• Se destacan con estilo negrita (BoldStyle) los resultados parciales y generales.
• Mayor cantidad de aclaraciones sobre los pasos algebraicos.
• Revisión de la ortografía y gramática, corrección de muchos errores.
• Diferente presentación de los sistemas de ecuaciones para mayor compresión.
• Se ha cambiado la presentación de las matrices por un modelo más común.
• Ingreso de un capítulo sobre tablas de derivación e integración.
• Solucionado un problema con el tipo de documento “amsbook” de lyx que producía
errores al traducir los layout de ejemplo y nota al español y los dejaba en ingles como
así una redeferenciación del “Proof” como Desarrollo, se escribe el siguiente código en el
“LaTex Preamble” (sin utilizar babel):
\newtheorem{exa}{Ejemplo}
\renewenvironment{example}{\begin{exa}}{\end{exa}}
\newtheorem{nota}{Nota}
\renewenvironment{note}{\begin{nota}}{\end{nota}}
\newtheorem{des}{Desarrollo}
\renewenvironment{proof}{\begin{des}}{\end{des}}
• Se eliminó el token “*” como operador de la multiplicación y se introdujo una escritura
normalizada alejando los factores del mismo tipo numérico con un pequeño espacio de
separación (Thinspace) y si fuesen signado se utilizan paréntesis.
• Se introduce un código en los enunciados de los ejemplos para reconocer la dificultad y
caminos alternativos de resolución (KeyLevelFeature).
47
Bibliografía
[1] EC7 Leithold - El cálculo séptima edición. De Louis Leithold. ISBN 970-613-182-5.
[2] Maxima - Common Lisp version of MACSYMA symbolic mathematics package released under the GNU Public
License - http://maxima.sourceforge.net
[3] LYX - http://www.lyx.org/
48