Calculo Teoria Trimestre 1 PDF
Calculo Teoria Trimestre 1 PDF
Calculo Teoria Trimestre 1 PDF
Edoardo Provenzi
Pr
ologo
Estos apuntes quieren ser una gua al curso de calculo del primer trimestre
y un soporte a las clases magistrales.
En escribir las notas he intentado poner el enfasis sobre las motiva-
ciones y los problemas pr aticos que llevan a la exigencia de construir cada
modelo matem atico que examinaremos. De hecho, historicamente, es el re-
conocimiento de los lmites de los instrumentos matematicos a nuestra dis-
posicion que ha permitido el desarrollo de la ciencia hasta llegar al nivel de
evolucion de nuestra epoca.
Espero que esta forma de presentar los conceptos matematicos resulte
mas comprensible y que los estudiantes puedan apreciar ya a partir de este
primer curso de c alculo la fundamental importancia de las matematicas en
la ciencia, tanto teorica como aplicada.
Para reforzar esta idea de conexion entre diferentes disciplinas cientificas,
he introducido en los apuntes aplicaciones a la fsica, a la mecanica y a la
teora de los se
nales.
Adem as, de acuerdo con la importancia cada vez mayor del ordenador
en la ciencia, examinaremos un algoritmo numerico extremadamente u til en
las aplicaciones: el metodo de Newton para aproximar los ceros de funciones
de una variable real, que extenderemos a funciones de mas variables en el
segundo trimestre.
Quiero agradecer, en orden alfabetico, a X`enia Alba, Pablo Arias, Felipe
Calderero, Juan Calvo y Gloria Haro, cuyas sugerencias y correcciones han
contribuido a mejorar sensiblemente estos apuntes.
El autor.
1
Indice general
3. Derivadas 64
3.1. Definici
on de derivada de una funcion y su interpretacion ge-
ometrica y mec
anica . . . . . . . . . . . . . . . . . . . . . . . 64
3.2. F
ormulas para el calculo de derivadas . . . . . . . . . . . . . 68
2
3.3. Derivabilidad y continuidad . . . . . . . . . . . . . . . . . . . 70
3.4.
Algebra de las derivadas . . . . . . . . . . . . . . . . . . . . . 71
3.5. Derivada de la funcion compuesta (regla de la cadena) y de
la funci
on inversa . . . . . . . . . . . . . . . . . . . . . . . . . 72
3.6. Diferenciabilidad y derivabilidad en R: el concepto de lineal-
izaci
on local de una funcion . . . . . . . . . . . . . . . . . . . 75
3.6.1. El diferencial de una funcion . . . . . . . . . . . . . . 78
3.7. La f ormula de Taylor con resto de Peano . . . . . . . . . . . . 80
3.8. El teorema de lH opital y sus aplicaciones al calculo de lmites
en forma indeterminada . . . . . . . . . . . . . . . . . . . . . 86
3.9. Jerarqua de infinitos e infinitesimos . . . . . . . . . . . . . . 88
3.10. Teoremas fundamentales sobre derivadas . . . . . . . . . . . . 91
3.11. El teorema del valor medio de Lagrange y sus aplicaciones al
estudio de la monotona de una funcion . . . . . . . . . . . . 93
3.11.1. La formula de Taylor con resto de Lagrange . . . . . . 102
3.12. La optimizaci on . . . . . . . . . . . . . . . . . . . . . . . . . . 103
3.13. El problema del fontanero . . . . . . . . . . . . . . . . . . . . 103
5. Integrales 116
5.1. Definici
on de la integral de Riemann . . . . . . . . . . . . . . 116
5.2. Propiedades de la integral . . . . . . . . . . . . . . . . . . . . 120
5.3. Valor medio y valor eficaz de una funcion . . . . . . . . . . . 121
5.4. El teorema fundamental del calculo diferencial . . . . . . . . 122
5.5. Integrales y dinamica de las partculas . . . . . . . . . . . . . 131
5.6. Calculo de la longitud del grafico de una funcion . . . . . . . 132
5.7. Integrales generalizadas: integracion de funciones no contin-
uas o en intervalos ilimitados . . . . . . . . . . . . . . . . . . 133
5.8. La funcion integral . . . . . . . . . . . . . . . . . . . . . . . . 136
5.8.1. Una aplicacion de la funcion integral: la ecualizacion
de histogramas . . . . . . . . . . . . . . . . . . . . . . 138
5.9. Relaciones de ortogonalidad de las funciones sin y cos . . . . 140
3
5.10. Aproximaci on numerica de integrales . . . . . . . . . . . . . . 141
5.10.1. Integraci
on por aproximacion rectangular . . . . . . . 142
5.10.2. Integraci
on por aproximacion trapezoidal . . . . . . . 143
5.10.3. Integraci
on por aproximacion parabolica (formula de
Cavalieri-Simpson) . . . . . . . . . . . . . . . . . . . . 145
5.10.4. Comandos Matlab para aproximar numericamente in-
tegrales . . . . . . . . . . . . . . . . . . . . . . . . . . 148
4
Captulo 1
1.1. Introducci
on
En este captulo profundizaremos en el conocimiento de los objetos basicos
del calculo: los conjuntos y las relaciones entre ellos, que se representan a
traves de funciones.
X = {x X : x satisface la propiedad P}
5
se lee como pertenece a y : se lee como tal que.
Veamos unas definiciones u tiles para el resto del curso:
Numeros enteros: Z = {. . . , 2, 1, 0, 1, 2, . . .} ;
p
Numeros racionales: Q = x : p, q Z, M CD(p, q) = 1 : x = .
q
6
La u on significa que para cada x Q existen ()1 dos n
ltima definici umeros
enteros p, q coprimos, o sea cuyo maximo com un denominador M CD(p, q) es
1, tales que x se puede escribir como el cociente entre p y q. La eleccion de p
y q coprimos sirve para evitar repetir la misma fraccion escrita con n umeros
que son m ultiplos de p e q, por ejemplo 32 y 46 definen el mismo n umero
racional, con la diferencia de que 2 y 3 son coprimos mientras que 4 y 6 no.
Una fraccion entre numeros coprimos se denomina fracci on irreducible. Ob-
viamente dos n umeros pares no pueden ser coprimos, porque, como mnimo,
el MCD entre ellos es 2. N otese tambien que si elegimos q = 1, obtenemos
que todos los n umeros enteros se pueden considerar n umeros racionales con
denominador igual a 1. Resultan evidentes las siguientes inclusiones estric-
tas:
N Z Q.
N tiene dos subconjuntos notables, que son el conjunto de los n
umeros pares
e impares:
P = {p N : n : p = 2n} , I = {i N : n : i = 2n + 1} .
7
ver que cualquier numero natural n restado por un numero natural m mayor
que n devuelve un n umero entero negativo que, por definicion, no pertenece
a N. Por ejemplo, n = 4, m = 7, n m = 3 6 N. De aqu sigue que
N no es un conjunto apto para resolver las ecuaciones del tipo x = a b,
con a, b N, porque el resultado x de la resta entre a y b podra no ser un
n
umero natural.
Hemos visto que para demostrar que un conjunto no tiene una propiedad
es suficiente ense
nar un contraejemplo. Viceversa, es decir,
8
dos lados, o sea , eso significa que las dos propiedades son equivalentes,
porque se implican mutuamente. se lee: si y solo si.
Veamos ahora la prueba de una proposicion en la cual se usa tanto el
concepto de elemento generico de un conjunto, como el de doble implicacion
y que, adem as, nos sirvir
a despues para demostrar un teorema importante
sobre los n
umeros racionales.
on: p P p2 P y i I i2 I.
Proposici
Prueba. Para demostrar una doble implicacion hay que demostrar, gen-
eralmente por separado, que vale la implicacion hacia la derecha y tambien
la implicacion hacia la izquierda. Comencemos con las implicaciones hacia
la derecha.
(p P = p2 P): escribimos el elemento generico del conjunto P y de-
mostramos que esta propiedad vale para ese elemento. La demostracion se
extender a automaticamente a todos los otros numeros pares. Sea entonces
p = 2n, con n N, elevamos p al cuadrado: p2 = (2n)2 = 2(2n2 ), entonces
tambien p2 se puede escribir como el producto de 2 por un n umero natural,
2n2 , eso prueba que p2 es par.
(i I = i2 I): an alogamente a lo que hicimos antes, escribimos el ele-
mento generico del conjunto I y demostramos que esta propiedad vale para
ese elemento, la demostracion se extendera automaticamente a todos los
otros numeros impares. Sea entonces i = 2n + 1, con n N, elevamos i al
cuadrado: i2 = (2n + 1)2 = 4n2 + 4n + 1 = 2(2n2 + 2n) + 1, entonces tambien
i2 se puede escribir como el producto de 2 por un numero natural, 2n2 + 2n,
mas 1, eso prueba que i2 es impar.
(p P = p2 P): sea p2 un n umero par obtenido elevando al cuadrado
un n umero natural p. Dado que N es la union disjunta del conjunto de los
numeros pares e impares, solo quedan dos posibilidades p es par o impar. Si
p fuera impar, por lo que acabamos de demostrar, p2 sera impar! Entonces
s
olo queda la posibilidad que tambien p sea par.
(i I = i2 P): se usa el mismo razonamiento del punto anterior. 2
9
Para hacer eso, necesitamos introducir una tecnica de demostraccion muy
potente que se llama
T
ecnica de demostracci
on por reducci
on al absurdo.
10
Este resultado es muy sorprendente y contraintuitivo porque, siendo Q
cerrado con respecto a la division, dados dos n umeros racionales, siempre
podemos construir otro n umero racional que este entre ellos, por ejemp-
lo su media aritmetica: p+q2 . Intuitivamente, podr amos pensar que pode-
mos seguir tomando promedios de n umeros racionales llegando a identificar
cualquier numero como un n umero racional, pero el teorema que acabamos
de demostrar ense na que nuestra intuicion no es correcta: no podemos repre-
sentar gr
aficamente los numeros racionales como una recta continua porque
en Qexisten lagunas, como por ejemplo 2. Para dibujar la laguna dada
por 2 se suele construir (vease la figura siguiente) un cuadrado de lado
1 con el vertice de suroeste en el punto 0 de una recta horizontal:
por el
teorema de Pit agoras, la diagonal del cuadrado tiene longitud 1 + 1 = 2,
si rotamos ladiagonal sobre la recta encontramos un punto que tiene una
distancia de 2 de 0 (porque la rotacion mantiene constante la longitud de
la diagonal del cuadrado). Como acabamos de ver, ese punto no puedeser
representado como un n umero racional. Existen infinitos puntos como 2,
por ejemplo . Esos puntos se llaman irracionales y se puede demostrar
que son muchos m as que los numeros racionales, pero eso va mas alla del
fin de este curso.
Proposicion: los n
umeros racionales coinciden con los numeros decimales
finitos o peri
odicos (con periodo diferente de 9) y los n umeros irracionales
coinciden con los n umeros decimales infinitos no periodicos.
11
R = Q {Irracionales}.
Se puede demostrar que cada punto sobre la recta se puede iden-
tificar con un n umero real, entonces R no tiene las lagunas de los
racionales. Por esta razon, a partir de ahora, llamaremos R el conjuntos
de los n
umeros reales o la recta real o el eje real. A pesar de tener lagu-
nas, Q es denso en R, una afirmacion que se puede formalizar en el siguiente
teorema.
El teorema de densidad nos dice que cada n umero real puede ser aprox-
imado con precisi on arbitraria por un numero racional, ya que, si queremos
aproximar por ejemplo r1 , podemos elegir r2 muy cercano a r1 , y por el
teorema anterior, existir
a un numero racional q a um mas cercano a r1 que
r2 ! Este hecho se usa cada da en los ordenadores.
Otra propiedad de R que resultara u til en el curso es que R es un con-
junto totalmente ordenado, es decir, para cada x, y R se puede definir una
on binaria < entre ellos con este significado: x < y x y < 0.
relaci
1.2.1. Intervalos en R
Dados a, b R, definimos estos subconjuntos notables de R:
12
(, b] = {x R : x b}, intervalo semiabierto e ilimitado inferior-
mente;
En particular, identificamos
R = (, +);
R+ = (0, +);
R+
0 = [0, +);
R = (, 0);
R
0 = (, 0].
13
C = {z : x, y R : z = x + iy} .
C se llama conjunto de los n
umeros complejos, x se llama parte real e
y se llama parte imaginaria de z. La teora de los n umeros complejos se
desarrolla en el curso de Algebra, al cual hacemos referencia. Aqu decimos
simplemente que C no es un conjunto ordenado y que R C porque cada
z C con parte imaginaria nula se puede identificar con el n
umero real que
da su parte real. La cadena de inclusiones entre conjuntos numericos es:
NZQRC.
14
1.3. Funciones
El concepto de funcion es central en matematica y, sorprendentemente,
s
olo en una fecha relativamente reciente (1829) ha sido formalizado.
La idea basica que est
a detras del concepto de funcion es la conexi
on
entre input y output de un sistema:
r Area r2 .
Resulta muy util abstraer la definicion y no atarse a un ejemplo concreto y
tener una definici
on general, que debemos a Johann Dirichlet, matematico
aleman (1805-1859) y Nikolai Lobachevsky, matematico ruso (1792-1856).
15
El ejemplo mas trivial de funcion esta dado por la funci
on identidad
sobre un conjunto cualquiera X:
idX : X X
x 7 idX (x) = x,
o sea, la funci
on que deja cualquier punto del conjunto invariado. Veremos
la utilidad del concepto de funcion identidad mas adelante cuando hablemos
de la relaci
on entre composicion e inversion de funciones.
Cuando X R se dice que f es una funcion de variable real, cuando
Y R se dice que f toma valores reales, si tanto X como Y son subconjuntos
de R, se dice que f es una funci on real de variable real:
f : X R Y R
x 7 y = f (x),
llamaremos f (x) a la expresi on analitica o correspondencia funci onal
de f . Como dominio y codominio de una funcion real de variable real son
siempre subconjuntos de R, se suelen denotar estas funciones simplemente
a traves de la correspondencia funcional, escribiendo y = f (x), expresando
a parte el dominio.
Aunque tambien nosotros muchas veces usaremos esta comoda conven-
ci
on, invitamos a los estudiantes a practicar la notacion completa, porque
resultarautil a la hora de analizar la composicion de funciones y la teora
de las funciones de m as variables.
Remarcamos adem as que una funci on esta completamente e unvo-
camente determinada a trav es de su dominio, codominio y corre-
spondencia funci onal. As pues, dos funciones que tienen la misma cor-
respondencia funci onal, pero dos dominios o codominios diferentes, tienen
que considerarse dos funciones diferentes.
Las funciones reales de variable real mas sencillas son las funciones
constantes.
k : R {k} R
x 7 k(x) = k.
16
todos los n
umeros reales en una constante k (que es la constante que define la
funci
on). En particular, tenemos la funci on id enticamente nula 0(x) 0
y la funcion unidad 1(x) 1.
El hecho de que la imagen de una funcion constante no contenga la vari-
able x no tiene que confundir: el dominio y el codominio de una funcion no
tienen porque coincidir, los datos de input y los de output de una funcion
pueden ser muy diferentes entre ellos. Como ejemplo concreto podemos pen-
sar a la funci
on que calcula el area de un rectangulo: los datos de input son
dos, la longitud de los dos lados del rectangulo, pero el dato de output es
s
olo uno, el
area del rect
angulo.
Otros ejemplos notables de funciones reales de variable real son la funci
on
valor absoluto y signo, que nos permiten mostrar tambien que las funciones
se pueden definir por trozos:
Funcion signo: devuelve el signo de un n
umero real multiplicado por
1. Se puede definir en dos formas, una definida sobre R quitando el
cero y la otra definida sobre todo R:
| | : R R+
0 (
x si x 0
x 7 | |(x) = |x| = .
x si x < 0
|x| = x sign(x) , x R;
x |x|
sign(x) = = , x R \ {0}.
|x| x
17
Aunque las funciones reales de variable real sean funciones de enorme
importancia en la ciencia, remarcamos que estas son solo una clase particular
de funciones y que es importante recordar la definicion general de funci on.
A cada funci
on se asocia un grafico, cuya definicion necesita recordar lo
que es el producto cartesiano entre dos conjuntos X, Y :
X n = X . . }X
| .{z potencia cartesiana n-esima de X,
n veces
X n est
a dado por n-uplas ordenadas de elementos de un solo conjunto. En
particular, si X = R, tenemos:
18
on f : X Y , x 7 f (x), su gr
Def.: dada la funci afico es el subcon-
junto de X Y definido por
Si f es una funci
on real de variable real, su grafico es una curva en el
plano real, como muestra la figura de izquierda de la siguiente imagen.
19
afico de la funcion identidad sobre R, o sea id(x) = x x R es
nula). El gr
la recta bisectriz al primer y tercer cuadrante:
El la p
agina siguiente introducimos unas definiciones sobre propiedades
notables de las funciones reales de variable real.
20
Propiedades de las funciones reales de variable real
on f : D R R real de variable real se dice:
Una funci
21
existen tres puntos con abscisas diferentes x1 , x2 , x3 , que tienen la misma
ordenada y = f (x1 ) = f (x2 ) = f (x3 ).
Def.: dadas f, g : D R R y c R,
f g : D R R, (f g)(x) = f (x)g(x), x D;
f n : D R R, f n (x) = (f (x))n , x D, n N;
c f : D R R, (c f )(x) = cf (x), x D;
f f f (x)
g : D R R, g (x) = g(x) , x D tal que g(x) 6= 0.
(f + g)(x) = sin(x) + x2 ;
(f g)(x) = x2 sin(x);
22
(3 f )(x) = 3 sin(x), (7 g)(x) = 7x2 ;
f sin(x)
g (x) = x2
, con x 6= 0.
1.5. Composici
on de funciones
Pensemos otra vez en una funcion como en una caja negra que conecta
un dato de input con un dato de output (y uno solo) de un cierto sistema.
Supongamos que en nuestro sistema existen dos cajas negras representadas
por dos funciones puestas en cascada (vease la proxima figura), en este caso
la relaci
on entre input y output esta dada por la funci
on compuesta entre
las dos. Obviamente el razonamiento tiene sentido con la condicion que la
segunda funci on pueda interpretar los datos de output de la primera, es
decir, si las imagenes de la primera funcion pertenecen al dominio de la
segunda.
f2 : W Z
w 7 z = f2 (w)
tales que Y W , podemos definir la composicion entre f1 y f2 , en este
on f2 f1 definida as
orden, como la funci
f2 f1 : X Z
x 7 (f2 f1 )(x) = f2 (f1 (x)).
Puede resultar u
til visualizar los pasos de la composicion como sigue:
f1 f2
X Y W Z
x 7 y = f1 (x) 7 z = f2 (y),
23
entonces f2 f1 puede ser pensada, en conjunto, como la funcion que conecta
x directamente a z, aunque intrnsecamente lo haga a traves de dos pasos: el
primero transformando x en y con la funcion f1 y el segundo transformando
y en z con f2 .
Observemos que:
Ejemplo:
f1 = log : R+ R
x 7 y = log(x),
f2 = ()3 : R R
y 7 y 3 ,
on f2 f1 es posible: f2 f1 (x) = f2 (log(x)) = (log(x))3 . La
la composici
composicion f1 f2 no es posible porque el logaritmo esta bien definido solo
+
sobre R , mientras que el codominio de f2 es todo R.
Terminamos esta secci on enunciando un resultado que relaciona mono-
tonia y composici on:
La composici
on de dos funciones monotonas crecientes es creciente;
La composici
on de dos funciones monotonas decrecientes es creciente;
24
1.6. Inversi
on de funciones
Volviendo a la analoga entre una funcion y una caja negra que asocia un
dato de input a un dato de output, es interesante preguntarse si es posible,
una vez que los datos hayan sido transformadso, volver atras y determinar
a que dato de input corresponde cada dato de output.
Es f
acil darse cuenta que esto es posible solo si la funcion es inyectiva.
Consideremos, por ejemplo, la situacion que muestra la figura siguiente, en
la cual la funcion f asocia a dos datos de input distintos, x3 y x4 , el mismo
dato de output y3 . Si queremos asociar un valor unvoco al dato de input
que f ha transformado en y3 , nos enfrentamos a una ambig uedad: tenemos
que elegir x3 entre x4 ? Si queremos asociar y3 a un dato de input a traves
de una correspondencia funcional no podemos asociar y3 a ambos porque
perderamos la caracterstica fundamental de una funcion: la unicidad de la
asociacion entre input y output! Por otro lado, tanto x3 como x4 tienen el
mismo derecho de ser elegidos. No hay forma de obviar este problema, de
ah que la falta de inyectividad haga que no sea posible invertir una funcion.
Def.: sea
f : X f (X) Y
x 7 y = f (x),
una funci
on inyectiva. Definimos la funcion inversa de f como sigue:
f 1 : f (X) X
y 7 f 1 (y) = x,
donde x es la u
nica pre-imagen de y, o sea el u
nico elemento de X tal que
f (x) = y.
25
Ejemplos:
3 : R R
x 7 y = 3 x,
x R.
Este es un hecho general: la composici
on entre una funci
on f y su
inversa f 1 es la funci
on identidad sobre el dominio de f , mientras
que la composici on entre f 1 y f es la funci
on identidad sobre el
codominio de f :
26
1.7. Los ceros de las funciones reales de variable
real y sus caracterizaci
on abstracta
En esta seccion definiremos los ceros de las funciones reales de variable
real y examinaremos una caracterizacion que tendra la doble utilidad de
practicar la tecnica de demostraccion por absurdo y de adelantar unos con-
ceptos b
asicos para desarrollare la teora de los lmites.
Una funcion generica puede tener ningun cero (por ejemplo la funcion
exponencial), un solo cero (por ejemplo el logaritmo, que tiene un cero en x =
1), un numero finito de ceros (por ejemplo un polinomio) o hasta un n umero
infinitos de ceros (por ejemplo las funciones trigometricas seno, coseno y
tangente).
27
1.8. Los conjuntos funcionales
Existen conjuntos cuyos elementos son funciones. Estos conjuntos se lla-
man, por razones obvias, conjuntos funcionales. Podemos dar unos ejemplos:
F = {f : D R R}
FI F = {f : D R R, f inyectiva}
FE F = {f : D R R, f exhaustiva}
etc.
El estudiante podra comenzar a ver la importancia de los espacios fun-
cionales durante el curso de Ecuaciones Diferenciales del tercer trimestre.
28
1.9. Las funciones elementales de la matem
atica y
sus gr
aficos
Durante las clases de practicas se presentaran las funciones elementales
del c
alculo. En esta secci
on mostramos sus graficos.
29
30
31
32
33
34
35
Captulo 2
Introduccion al calculo
diferencial: el concepto de
lmite
36
basta para describir movimientos mas complejos como, por ejemplo, la
cada de un cuerpo bajo la accion de la gravedad, en el cual la veloci-
dad cambia en cada instante. En estos casos tenemos que considerar
la velocidad instant anea, o sea la velocidad media en un instante de
tiempo infinitamente breve. Pero, que significa eso? Matematicamente
esa frase no tiene sentido, porque es evidente que si hacemos tender
t a 0, tambien s tendra a 0, entonces llegaremos a una expresion
que, por ahora, no tiene ning un significado: 00 . La primera persona
que atac o este problema fue uno de los genios mas grandes de la his-
toria, el fsico-matematico ingles Isaac Newton (1642-1727), nacido
en el a no en el cual se murio Galileo Galilei! Con su metodo de las
fluxiones, Newton dio la primera definicion de velocidad instantanea
como derivada de la funcion s(t) que asocia a cada instante de tiempo
t el punto en el espacio donde se encuentra la partcula que se esta de-
splazando. Como veremos, la derivada de s(t) es el lmite, para t
que tiende a 0, de la razon entre s y t. Este ejemplo ense na como
la exigencia de la formalizacion del concepto de lmite y de derivada
esta ntimamente relacionada con el nacimiento de la ciencia moderna.
37
como la u nica recta que pasa por P y que es perpendicular al radio
de C, o bien como aquella recta que pasa por P y que no intersecta C
en ning un otro punto. Pero si intentamos extender estas definiciones
a curvas m as generales podemos ver enseguida que pierden su signifi-
cado: el radio es una caracterstica de las circunferencias, no existe
el radio de parabolas o hiperbolas, etc. La segunda definicion parece
m as facil de generalizar, y en cierta medida lo es porque sigue valida
para las c onicas (circunferencias, elipses, parabolas, hiperbolas), pero
si consideramos, por ejemplo, la curva dada por el grafico de la funcion
y = x3 , podemos ver que (aparte de en el origen) la recta tangente al
grafico en un punto intersecta siempre otro punto del grafico msmo,
como se ve en la figura aqu debajo.
Nos damos cuenta de que, a pesar de que nos parezca obvio e intuitivo
38
definir la recta tangente a una curva en una cierta forma, no podemos
dar una definicion rigurosa y general.
Geometricamente, la recta tangente en P se puede aproximar con una
recta secante en dos puntos al grafico de f , uno fijo, P , y uno movil, Q,
que puede acercarse o alejarse de P . A medida que Q se acerca a P , sin
coincidir con P , la (
unica) recta que pasa por P y Q se acerca a la recta
tangente al grafico de f en P . Es importante observar que este proceso
no termina con Q que viene a coincidir con P , porque en ese caso no
tendramos la posibilidad de definir una sola recta, dado que por un
punto pasan infinitas rectas! En cambio, para poder definir la recta
tangente, Q tiene que acercarse a P indefinidamente. Encontramos de
nuevo un problema que necesita la definicion rigurosa de un proceso
infinito.
Como veremos en breve, la definicion de recta tangente a una curva
est
a relacionada con la derivada, y entonces con la velocidad, pero ob-
servamos que su utilidad va mas alla: considerando esta figura es facil
39
computar cu antas veces el perfil de ese objeto contiene una unidad
de medida de longitud canonica (un milmetro, centmetro, metro,
kil
ometro, etc.), pero si un objeto tiene un perfil dado por una curva
no rectilinea este metodo no se puede utilizar. Sin embargo, si con-
sideramos trozos peque nos de esa curva, estos seran mas semejantes
a una recta y podremos aproximar sus longitudes con la tecnica de
arriba. Obviamente en matematicas el concepto de trozos peque nos
de una curva no tiene sentido y, de nuevo, aparece en un problema
real la necesidad de definir rigurosamente que quiere decir acercar in-
definidamente dos puntos;
Como comentario final a estas motivaciones podemos decir que los prob-
lemas de la geometra curvilnea, sin el calculo diferencial, no solo no podran
ser resueltos, sino tampoco definidos!
40
estudiantes de las universidades politecnicas de Turn y Pars, se dio cuen-
ta de que no saba contestar con precision a las dudas de sus alumnos con
respecto a conceptos como acercarse indefinidamente a un punto o de
cantidad infinitesima. Cuando Lagrange puso en evidencia el problema de
la formalizaci on del concepto de lmite, su fama internacional y el respeto
de la comunidad cientifica hicieron que un n umero cada vez mayor de ci-
entificos prestigiosos se dedicaran a un problema que hasta aquel momento
haba sido considerado de interes secundario: la formalizacion del concepto
de lmite.
Una primera respuesta a la exigencia de rigor formulada por Lagrange
fue dada por el matem atico frances Jean-Baptiste dAlembert (1717-1783),
pero fue Augustin-Louis Cauchy (1789-1857), tambien frances, el que re-
solvio definitivamente el problema de la definicion formal de lmite, y lo hizo
con instrumentos algebraicos de sorprendente sencillez y que ya estaban a
disposicion de los matem aticos desde hace siglos: las inecuaciones!
La definici on de Cauchy ha pasado a la historia como la definicion -
y para ayudar a su comprension es u til introducir el concepto de entorno de
un numero real y entorno del infinito:
Def.: sea x0 R,
porque |f (x)| < corresponde a la union de las soluciones del los dos sis-
temas (
f (x) g(x) <
f (x) g(x)
41
(
(f (x) g(x)) <
f (x) < g(x).
on del primer sistema es g(x) f (x) < g(x) + , la solucion del
La soluci
segundo sistema es g(x) < f (x) < g(x), la union de las dos soluciones es
g(x) < f (x) < g(x) + .
En particular,
42
on (Cauchy): sea f : D R R, con D
Def. - de lmite de funci
subconjunto abierto de R, y sea x0 D o bien sea x0 un punto de la
frontera de D. Se dice que f converge a ` cuando x tiende a x0 o que el
lmite de f para x que tiende a x0 es `, y se escribe:
lm f (x) = ` ,
xx0
43
Cauchy ha popularizado su teora sobre los lmites en su Cours dAnalyse
(1821) para la Ecole Polytechnique de Paris, la eleccion de las letras griegas
y se motiva con el hecho de que el sola abreviar la palabra error con
y la palabra distancia con . Como debera resultar claro de su definicion,
y representan, respectivamente, el error que se hace considerando el valor
de f (x) en vez de ` y la distancia entre x y x0 .
N otese que Cauchy logro dar una definicion rigurosa de distancia in-
finitesima (y entonces de lmite) a traves de una implicacion entre desigual-
dades que tiene que valer para cada eleccion de un numero real positivo. El
hecho de que el eje real no tenga los vacos que caracterizan los conjuntos
discretos N, Z y Q, es fundamental para estar seguros de que el proceso de
regresi on de y hacia cero (sin llegar nunca a cero) este bien definido.
on: si f es una funcion constante k(x) k x R, entonces
Observaci
f (x) toma siempre el mismo valor, por lo tanto es evidente que el lmite
para x x0 de una funci on constante es igual al valor constante k para
cada x0 R. Entonces la operaci on de lmite aplicada a una funci on
constante es trivial.
El lmite de una funci
on en un punto de la frontera de su dominio puede
no ser finito, como vamos a ver en las siguientes definiciones.
1. el lmite de f es +, y se escribe:
lm f (x) = + ,
xx0
44
si para cada entorno UM (+), existe un entorno UM (x0 ) D tal
que, para cada x UM (x0 ), x 6= x0 , entonces f (x) UM (+);
2. el lmite de f es , y se escribe:
lm f (x) = ,
xx0
45
1.
lm f (x) = ` ,
x+
2.
lm f (x) = + ,
x+
si para cada entorno UK (+), existe un entorno UMK (+) tal que,
para cada x UMK (+), se tiene que f (x) UK (+);
3.
lm f (x) = ,
x+
si para cada entorno Um (), existe un entorno UMm (+) tal que,
para cada x UMm (+), se tiene que f (x) Um ().
46
2.3. Lmites laterales
Las definiciones dadas arriba se pueden modificar ligeramente para definir
los lmites laterales. En sustancia, lo u
nico que cambia en estas definiciones
es que se permite a x acercarse a x0 , o a f (x) acercarse a `, solo desde la
izquierda o desde la derecha, respectivamente.
lm f (x) = `
xx0
47
2. f converge a `+ R, o que `+ es el lmite derecho de f cuando x
tiende a x0 , y se escribe
lm f (x) = `+
xx0
Ejemplos.
48
Dado que | x+1 1
x 1| = | x |, () se puede reescribir como
1
> 0 M > 0 tal que x > M < (?)
x
donde hemos quitado el valor absoluto porque estamos considerando
x > M > 0, entonces tambien x1 > 0. Se ve muy facilmente que
si tomamos M 1 , entonces la implicacion contenida en la (?) es
verdadera, porque si invertimos la desigualdad x > M 1 obtenemos
1
x < , lo cu al demuestra que lmx+ x+1x = 1 seg un la definicion de
Cauchy. Dejamos como ejercicio mostrar que tambien lmx x+1 x =
1. Entonces, la recta y = 1 es una as ntota horizontal para la funcion
f (x) = x+1
x cuando x , como se ve en la figura aqu debajo. Un
til ejercicio es probar que lmx0 x+1
u x = , as
que la recta x = 0
resulta una asntota vertical para f (x) = x+1
x .
En http://www.math24.net/definition-of-limit-of-function.html se
pueden encontrar m as ejemplos de lmites calculados a traves de la defini-
cion . Se puede observar que la definicion de Cauchy, a pesar de haber
dado rigor a un concepto extremadamente dficil como el de una cantidad
indefinidamente cercana a otra (y con herramientas muy sencillas como de-
sigualdades e implicaciones logicas), no es un instrumento muy comodo
para el c alculo efectivo de los lmites. La dificultad de la tecnica de la defini-
cion de Cauchy crece al aumentar la complejidad de la expresion analtica
de f , hasta llegar a ser pr acticamente in util si f es una funcion compuesta
muy complicada. En las secciones siguientes desarrollaremos unas tecnicas
m as eficientes para calcular lmites. A traves de estas tecnicas, el calculo de
lmites como los dos presentados arriba sera inmediato.
49
Ejemplos de funciones que no tienen lmite son las funciones trigonometri-
cas cuando x :
o sea, si el valor absoluto de h(x) esta limitado por una funcion que
tiende a 0, entonces tambien h tiende a 0 (porque |h(x)| g(x) equiv-
ale a g(x) h(x) g(x), y, si g(x) 0, tambien g(x) 0,
entonces, por el teorema citado arriba, tambien h(x) 0). Veremos
aplicaciones de este resultado en las clases de practicas.
50
Cambio de variable en los lmites: Si f g esta bien definida en
un entorno de x0 y g(x) t0 , entonces, poniendo t = g(x), vale que
xx0
lm f (g(x)) = lm f (t).
xx0 tt0
1
Ejemplo: si queremos calcular lmx0 e x , podemos poner t = x1 , que
1
tiende a cuando x 0 , por tanto: lmx0 e x = lmt et =
0;
1. f (x) g(x) `1 `2 ;
2. cf (x) c`1 , c R ;
3. f (x) g(x) `1 `2 ;
f (x) `1
4. , siempre y cuando g(x), `2 6= 0 para que la fraccion
g(x) `2
tenga sentido;
5. f (x)g(x) ``12 , siempre y cuando f (x), `1 > 0 para que la op-
eraci
on de exponencial tenga sentido.
1.
` = ,
o sea, la suma de los limites de una funcion convergente a ` R
y una funcion divergente es una funcion con lmite divergente;
2.
+ + = + ;
51
3.
= ;
4.
`= , ` 6= 0;
donde ` es el lmite de una funcion, y es diferente de cero. El signo
de depende del valor de ` y del en el lado de izquierda;
5.
`
=,
0
donde 0 es el lmite de una funcion que tiende a cero, no el
n
umero cero!
6.
`
=0.
+ ;
0 ;
0
0;
;
1 ;
00 ;
(+)0 .
52
2.5. La extensi
on de R
Ya hemos dicho que y + no son n umeros reales, sino el compor-
tamiento o tendencia de funciones divergentes. Sin embargo, teniendo en
on la representacion de R como una recta, podemos pensar
consideraci
como aquel punto que esta siempre a la izquierda de cada numero real y
+ como aquel punto que esta siempre a la derecha de cada numero real.
El conjunto dado por la union de R y de los dos puntos y + se
llama extensi
on de R (o recta extendida/ampliada) y se suele escribir como
sigue:
R = R {} {+} .
La interpretacion grafica canonica de R se obtiene considerando la figura
debajo: se ve que podemos poner en correspondencia biyectiva cada punto
del eje real con los puntos de la semicircunferencia simplemente conectando
el centro C con la recta horizontal. Solo los puntos A y B no tienen corre-
spondencia con puntos de la recta, sin embargo, resulta natural pensar que
el punto A esta en correspondencia con y que B esta en correspondencia
con +.
lm f (x) = f (x0 ) ;
xx0
53
Si escribimos explcitamente la definicion de lmite para una funcion
continua, es decir:
nos damos cuenta de que, para una funcion continua, variaciones peque nas
de x alrededor de x0 corresponden a variaciones peque nas de f (x)
alrededor de f (x0 ).
Si escribimos x = x0 + h, para un cierto h R oportuno, tenemos que
x x0 si y solo si h 0, entonces, gracias al teorema de cambio de variable
en los lmites, podemos decir tambien que f es continua en x0 si
lm f (x0 + h) = f (x0 ) .
h0
54
por la derecha e izquierda, lmxx f (x) = `1 , lmxx+ f (x) = `2 , pero
0 0
`1 6= `2 , se dice que f tiene una discontinuidad de salto; el salto es |`2 `1 |
y puede ser finito o infinito. Las funciones que tienen discontinuidades de
salto como el signo resultan muy utiles para modelar un cambio abrupto de
una cantidad.
Hay casos en los que una funcion f no esta definida en un punto x0 ,
pero existe su lmite para x x0 . En ese caso es posible prolongar por
continuidad la funci on f como sigue:
(
f (x) si x D
f(x) =
lmxx0 f (x) si x = x0 ,
55
Entonces el c
alculo de lmites de funciones continuas en los pun-
tos de sus dominio es trivial. Los problemas surgen a la hora de calcular
el lmite en los puntos de la frontera del dominio o cuando aparecen formas
de indeterminaci on, como veremos.
2. f toma valores con signo opuestos en los extremos de [a, b], o sea,
f (a) f (b) < 0,
2
Bernard Bolzano, matem atico ceco (1781-1848).
3
Karl Weierstrass, importantsimo matem atico aleman (1815-1897).
56
Teorema de los valores intermedios (Darboux4 ): con las mismas no-
taciones del teorema de Weierstrass, entonces, [m, M ] x [a, b] tal
que f (x) = .
Los dos u
ltimos teoremas se pueden resumir diciendo que una funci on f
continua sobre un intervalo compacto toma todos los valores entre
el mnimo y el m aximo, con notacion matematica: f ([a, b]) = [m, M ].
La figura debajo muestra que la hipotesis de continuidad es fundamental
para la validez del teorema: en la imagen de la derecha todos los valores
contenidos entre y1 e y2 no corresponden a ning un valor de x en [a, b] porque
la funci
on dibujada tiene una discontinuidad de salto en x = c.
Terminamos esta secci on con una observacion que nos ayudara entender
mejor la importancia de la continuidad de los n umeros reales. Considere-
mos la funci on f : Q Q, x 7 x2 , o sea, la funcion que aplica un numero
racional cualquiera en su cuadrado (bien definida porque sabemos que Q es
cerrado con respecto a esa operacion). Se puede demostrar muy facilmente
que f no posee la propiedad de los valores intermedios con un contraejemplo:
f (1) = 1 y f (2) = 4, el numero 2 esta entre 1 y 4, sin embargo no existe
ningun numero racional que elevado al cuadrado devuelva 2, porque solo 2
elevado al cuadrado genera 2 y sabemos que 2 6 Q! Este ejemplo muestra
claramente como las propiedades de los numeros reales y el calculo diferen-
cial est
an ntimamente relacionados y cuan util es el conjunto de los n
umeros
reales como ambiente de trabajo para desarrollar el calculo diferencial.
4
Jean Gaston Darboux (1842-1917), matem
atico frances.
57
2.7. Infinitos e infinit
esimos de orden inferior, su-
perior y del mismo orden, la relaci on de ser
asint
otico a
En esta secci
on vamos a examinar la rapidez relativa de convergencia
a cero y de divergencia al infinito. Ya hemos observado que la operacion
de lmite define el comportamiento o tendencia de una funcion cuando su
variable se acerca a un determinado valor. Esta tendencia puede ser mas
o menos marcada seg un la funcion y seg
un el punto donde se eval ua el
lmite. Saber comparar la rapidez de convergencia o divergencia de funciones
es extremadamente importante en muchas aplicaciones del calculo, como
veremos durante el curso y como el estudiante podra apreciar durante el
desarrollo de sus estudios.
Comenzamos con una definicion preliminar.
f es infinit
esima (o un infinitesimo) en x0 si lmxx0 f (x) = 0;
58
En el caso (3) se dice que f es un infinito de orden superior a g
en x0 , o bien que f diverge mas rapidamente que g (de nuevo, intu-
itivamente, el lmite puede ser infinito solo si f llega a tomar valores
grandes antes que g).
`
habiendo usado la regla de aritmetizacion parcial del dada por: =0
` R.
Por tanto, podemos decir que f (x) = 2000x2 es un infinito de orden
inferior a g(x) = 5x3 para x + o bien que g(x) = 5x3 es un infinito de
orden superior a f (x) = 2000x2 para x +.
Sean ahora f, g : D R R funciones infinitesimas en x0 D. Ahora
la cuestion es saber si podemos medir la rapidez relativa de convergencia
a 0 entre f y g. Como antes, vamos a evaluar el lmite: lmxx0 fg(x)
(x)
. Si este
lmite no existe, se dice que f y g son infinitesimos no comparables en
x0 . Si el lmite existe, distinguimos los casos siguientes:
0 (1);
f (x)
lm = k R (2);
xx0 g(x)
(3).
59
Ejemplo: consideremos otra vez f (x) = 2000x2 y g(x) = 5x3 pero exam-
inemos sus comportamiento cuando x 0 y ambas funciones convergen a
0.
2000x2 400
lm 3
= lm = +,
x0 5x x+ x
habiendo usado la regla de aritmetizacion parcial del dada por: 0` = +
` R+ .
Por tanto, podemos decir que f (x) = 2000x2 es un infinitesimo de orden
inferior a g(x) = 5x3 para x 0 o bien que g(x) = 5x3 es un infinitesimo
de orden superior a f (x) = 2000x2 para x +.
Entonces, para cada pareja de coeficientes a, b 6= 0, ax3 es tanto un
infinito de orden superior a bx2 cuando x + como un infinitesimo de
orden superior a bx3 cuando x 0, como confirma la figura de aqu debajo.
f (x)
f (x) g(x) lm =1.
xx0 xx0 g(x)
60
La relacion de ser asintotico a es extremadamente u til a la hora de
calcular lmites, por eso es importante destacar sus propiedades:
1. es una relaci
on de equivalencia, o sea, son ciertas las siguientes
propiedades:
f (x) f(x)
f (x)g(x) f(x)
g (x) y
xx0 g(x) xx0 g(x)
61
Para mostrar un contraejemplo sobre la funcion exponencial vamos a
2 2
comprobar si ef (x) = e3x x es asintotica a ef (x) = e3x para x +
con la definici
on:
2
e3x x
lm = lm ex = 0,
x+ e3x2 x+
entonces ef (x) 6 ef (x) , en general.
lm f (x) = lm an xn
xx0 x
62
lmx f (x) lmx f (x) f (x)
entonces lmx an xn = 1, pero de lmx an xn = lmx an xn se sigue que
an xn + an1 xn1 + . . . + a1 x + a0 an xn ,
x
an xn + an1 xn1 + . . . + a1 x + a0 an xn an nm
= x .
bm xm + bm1 xm1 + . . . + b1 x + b0 x bm x m bm
30x4 +2000x 2
+80 .
Ejemplo: calcular lmx 7x8 2x
Tenemos que
30x4 + 2000x2 + 80 30 4
x 0,
7x8 2x x 7 x
30x4 +2000x 2
+80 = 0.
entonces lmx 7x8 2x
5 +x3 +8x
Ejemplo: calcular lmx x 3x5 19x2
. Tenemos que
x5 + x3 + 8x 1 1 1
x0 = .
3x5 19x2 x 3 3 x 3
la suma de infinit
esimos es asintotica al infinit
esimo de orden
inferior.
63
Captulo 3
Derivadas
3.1. Definici
on de derivada de una funci on y su
interpretaci
on geometrica y mec
anica
Como ya dijimos, el concepto de derivada esta relacionado con el proble-
ma de encontrar la recta tangente al grafico de una funcion o, de forma mas
general, a una curva en un cierto punto. Para entender mejor este problema,
consideremos la figura aqu debajo.
64
de A a B con la raz
on incremental definida por
f (x0 + h) f (x0 )
lm ,
h0 h
f 0 (x0 ) , (notaci
on de Lagrange), usado sobre todo por los matematicos;
f(x0 ) , (notaci
on de Newton), usado sobre todo por los fsicos,
cuando f depende del tiempo;
df
(x0 ) , (notaci
on de Leibnitz), que resulta muy util a la hora
dx
de hacer operaciones formales con las derivadas, como veremos mas
adelante;
65
on de Arbogast1 ), que resulta u
Df (x0 ) , (notaci til cuando se ex-
tende la derivaci
on a funciones de mas variables o a operadores difer-
enciales en espacios funcionales.
Resulta evidente que f 0 (x0 ) expresa la rapidez de variaci on de la
funcion f en un entorno infinitesimo de x0 .
Por construcci on, f 0 (x0 ) es la pendiente de la recta tangente al
gr
afico de f en (x0 , f (x0 )). Para determinar la ecuacion de esta recta es
suficiente utilizar la f
ormula que nos proporciona la recta que pasa por un
punto con una pendiente conocida. En nuestro caso se obtiene:
Recta tangente al gr
afico de f en (x0 , f (x0 )) :
f 0 : (a, b) R
x 7 f 0 (x),
df
otros smbolos para la funcion derivada son f, dx .
0
Si f tambien es derivable, podemos calcular la derivada de la funcion
derivada de f , que llamaremos derivada segunda y escribiremos f 00 o f o
d2 f
dx2
o D2 f :
f 0 (x0 + h) f 0 (x0 )
f 00 (x) = lm ,
h0 h
otese que ahora en la razon incremental aparece f 0 y no f como antes.
n
Veremos m as adelante el significado geometrico de la derivada segunda en
relacion a la curvatura del gr afico de una funcion.
Podemos iterar este razonamiento, obteniendo la derivada de orden
n
n 1 en x0 , que denotamos con los smbolos f (n) (x0 ), ddxnf (x0 ) o Dn f .
La derivada (que se suele llamar tambien derivada primera) y la deriva-
da segunda tienen significados fsicos importantes: dado un cuerpo en
movimiento a lo largo de un eje espacial, podemos definir la funcion s que
asocia a cada instante temporal t el lugar s(t) donde se encuentra el cuerpo
en el eje:
s : [0, +) R
t 7 s(t)
1
Antoine Arbogast (1759-1803), matem
atico frances.
66
s: funci
on posici
on espacial.
v(t) = s(t)
: velocidad instant anea del cuerpo en movimiento en el
instante t (rapidez de variacion de la funcion posicion);
a(t) = v(t)
= s(t) : aceleraci
on instant anea del cuerpo en movimien-
to en el instante t (rapidez de varicion de la funcion velocidad);
donde m es la masa del cuerpo, que Newton identifico como aquel atrib-
uto de un cuerpo que se opone a variaciones de velocidades provocadas
por una fuerza externa (inercia).
67
3.2. F
ormulas para el c
alculo de derivadas
Utilizando el lmite que define las derivadas se pueden calcular las derivadas
de las funciones elementales, que son los ladrillos fundamentales para calcu-
lar todas las otras derivadas. Veamos un ejemplo, calculemos la derivada de
on y = x2 con la f
la funci ormula de la definicion:
(x + h)2 x2 x2 + 2xh + h2 x2 2xh + h2
= = = 2x + h 2x
h h h h0
68
N
otese ademas que ex es la u nica funci on que coincide con su
derivada y que (sin x)00 = sin x y (cos x)00 = cos x, estas dos observa-
ciones ser
an importantes para el curso de Ecuaciones Diferenciales.
Vamos ahora a examinar la derivada de la funcion valor absoluto:
(
0 (x)0 = 1 si x > 0
(|x|) = 0
(x) = 1 si x < 0
0 es un punto donde |x| cambia definicion, entonces hay que examinarlo con
cuidado con el lmite que define la derivada en x = 0:
|0 + h| |0| |h|
lm = lm = lm sign(h),
h0 h h0 h h0
este lmite no existe, dado que los lmites por la izquierda y por la derecha
valen -1 y +1, respectivamente.
Entonces la funci on x 7 |x| no es derivable es 0 y su derivada para
valores x 6= 0 coincide con el signo de x:
(
0 1 si x > 0
|x| = sign(x) =
1 si x < 0.
Mas all
a de los puntos angulosos existen otras dos singularidades de las
funciones con respecto a las derivadas: los puntos de inflexion con tangente
vertical y las c
uspides.
69
Def.: una funci
on f tiene una c uspide en x0 si es continua en x0 , no
derivable en x0 pero existen los lmites de las razones incrementales por la
izquierda y derecha de f en x0 y son infinitos con signos opuestos.
p
3
Un ejemplo de funci
on con una c
uspide en x = 0 es f (x) = |x|:
70
on.: f : (a, b) R R derivable en x0 (a, b) = f continua
Proposici
en x0 .
f (x0 + h) f (x0 )
f (x0 + h) f (x0 ) = h.
h
Ahora operamos los lmites para h 0 a ambos lados:
f (x0 + h) f (x0 )
lm [f (x0 + h) f (x0 )] = lm h. (3.1)
h0 h0 h
Gracias a la hip
otesis de derivabilidad de f y a las propiedades de los lmites,
el lado de la derecha se puede reescribir como:
f (x0 + h) f (x0 )
lm lm h = f 0 (x0 ) lm h = 0.
h0 h h0 h0
3.4.
Algebra de las derivadas
Utilizando la f
ormula definitoria, se puede demostrar que la derivada de
la suma o resta, del producto y de la division entre funciones se calculan
como dice la proposici
on siguiente.
(f g)0 = f 0 g 0
(k f )0 = k f 0 k R
71
(f g)0 = f 0 g + f g 0 Regla de Leibnitz2
(f /g)0 = (f 0 g f g 0 )/g 2
en particular, de la u
ltima propiedad y del hecho de que la derivada de una
funci
on constante es nula, sigue que
0
1 g0
= 2
g g
dz dz dy
= ,
dx dy dx
72
Para usar esta regla hay que aprender a interpretar una funcion complicada
como composici on sucesiva de funciones elementales. Para ello es u
til pensar
c
omo se calcula una funcion compuesta a traves de una calculadora.
sin ()3
x 7 y = sin x 7 z = (sin x)3 .
1
(f 1 )0 (y) = F
ormula de derivaci
on de la funci
on inversa.
f 0 (x)
Es a traves de esta f
ormula que se pueden calcular las derivadas de
las funciones trigonometricas inversas. Newton descrubrio esta formula ge-
ometricamente: como se recuerda en la figura siguiente, el grafico de f y de
su inversa son simetricos con respecto a la bisectriz y = x, sigue que los
angulos y son complementarios, o sea + = 2 .
73
Sigue que f 0 (x) = tan() = tan( 2 ) = 1 1
tan() = f 1 (y) , habiendo
usado la identidad de la tangente: tan( 2 ) = 1
tan() y el hecho de que
(f 1 )0 (y) = tan().
Con la notaci on de Leibnitz podemos escribir la regla de derivacion de
la funci on inversa con una expresion casi tautologica:
dx 1
= dy .
dy
dx
En la secci
on siguiente vamos a analizar el origen de la notacion de
Leibnitz examinando el concepto de diferencial.
74
3.6. Diferenciabilidad y derivabilidad en R: el con-
cepto de linealizaci
on local de una funcion
El c
alculo diferencial fue inventado contemporaneamente por Newton y
Leibnitz durante la segunda mitad del 1600. Mientras Newton estaba mas
interesado en las aplicaciones del concepto de derivada a la fsica, obtenien-
do resultados excepcionales como la ley de gravitacion universal o la de-
mostraci on matem atica de que los planetas del sistema solar siguen orbitas
elpticas, Leibnitz estaba mas interesados en la componente algebraica del
c
alculo diferencial a traves del estudio del diferencial, que ahora vamos a
introducir volviendo a examinar la figura:
Una operaci on muy com un tanto en matematica pura como en las apli-
caciones es la linealizaci on, o sea la sustitucion de una funcion que depende
no linealmente de sus variables con una funcion lineal. Como en cada aprox-
imaci on, tambien en el caso de la linealizacion hay que proporcionar una
estimaci on cuantitativa o cualitativa del error que se cumple aproximando.
Si una funci
on f es derivable en un punto x0 , la funcion lineal que pode-
mos asociar de forma m as natural a f es la recta tangente al grafico de f en
x0 : y(x) = f (x0 ) + f 0 (x0 )(x x0 ).
En la figura, la variaci on de la ordenada del punto sobre el grafico de
f cuando x vara de x0 a x0 + h esta dibujada en verde, mientras que la
variacion de la ordenada del punto sobre la recta tangente cuando x vara
de x0 a x0 + h est a dibujada en azul.
Es razonable considerar el error aceptable si la diferencia entre las or-
denadas del punto sobre el grafico de f y sobre la recta tangente (o sea la
75
diferencia de longitud del segmento verde y azul, respectivamente) tiende a
0 mas rapidamente que el desplazamiento h. La diferencia entre la ordenada
de f en x0 + h y la ordenada del punto sobre la recta tangente y(x) en
x = x0 + h se suele escribir como o(h) o-pequena y vale:
La condici
on de aceptabilidad del error, o sea que o(h) tienda a 0 mas
rapidamente que h 0 tiene un nombre particular:
Dado que la aproximaci on de f se hace con una funcion lineal y dado que
esta aproximaci on produce un error despreciable solo en puntos cercanos a
x0 , se suele decir que la diferenciaci on es un proceso de linealizaci on
local de una funci on.
Por lo tanto, las funciones no diferenciables son aquellas que no se pueden
linealizar localmente de forma eficaz porque no tinen recta tangente o
porque el error que se cumple aproximando la funcion con su recta tangente
es demasiado grande en el sentido especificado arriba.
Para funciones reales de variable real, y solo para estas, resulta que ser
derivable y ser diferenciable son condiciones equivalentes, como vamos a
demostrar.
76
Teorema: La funci on f : (a, b) R R es diferenciable en x0 (a, b) si
y solo si es derivable en x0 .
f (x0 + h) f (x0 )
lm = f 0 (x0 ),
h0 h
que es la definici
on de derivabilidad.
f (x0 + h) f (x0 )
lm = lm f 0 (x0 ),
h0 h h0
77
si pasamos el lmite que aparece en el lado de la derecha al lado de
la izquierda y utilizamos las propiedades algebraicas de los lmites
podemos escribir:
f (x0 + h) f (x0 ) 0
lm f (x0 ) = 0,
h0 h
78
Def.: Sea f : (a, b) R R diferenciable en x0 (a, b). El valor del
diferencial (de orden 1) de f en x0 es
df (x0 ) = f 0 (x0 )h .
df
df (x0 ) = f 0 (x0 )dx (x0 ) = f 0 (x0 ) ,
dx
79
3.7. La f
ormula de Taylor con resto de Peano
En la secci
on antecedente hemos comentado que, si f es derivable en x0 ,
podemos aproximar f localmente con una funcion lineal, la recta tangente a
su gr
afico en x0 , con menos de un error que resulta un infinitesimo de orden
superior con respecto a x x0 :
f (x) ' f (x0 ) + f 0 (x0 )(x x0 ) + f 00 (x0 )(x x0 )2 + . . . + f (n) (x0 )(x x0 )n NO!
(3.3)
donde ' es el smbolo de egualdad aproximada. Esta formula no puede ser
correcta por una raz on muy sencilla. Al fin de entender esta razon, conviene
suponer por el momento que n sea igual a 2. En este caso, la formula (incor-
recta) (3.3) se escribe as: f (x) ' f (x0 ) + f 0 (x0 )(x x0 ) + f 00 (x0 )(x x0 )2
y entonces, si calculamos las derivadas con respecto a x de las funciones
que aparecen en ambos lados obtenemos: f 0 (x) ' f 0 (x0 ) + 2f 00 (x0 )(x x0 )
y f 00 (x) ' 2f 00 (x0 ) (observamos que f (x0 ) y f 0 (x0 ) son numeros reales y se
comportan como constantes, o sea se anulan, si las derivamos con respecto
a la variable x). Si calculamos la derivada en x = x0 , la u ltima relacion
00 00
que acabamos de escribir se vuelve f (x0 ) ' 2f (x0 ) lo cual es evidente-
mente incorrecto (a menos que f 00 (x0 ) sea nula, caso en el cual tenemos una
aproximaci on lineal!). Si queremos modificar la formula (3.3) para eliminar
80
la incoherencia que acabamos de subrallar, tenemos que introducir un co-
eficiente que anule el factor 2 que aparece delante de la derivada segunda,
es decir: f (x) ' f (x0 ) + f 0 (x0 )(x x0 ) + 21 f 00 (x0 )(x x0 )2 . Cuales son los
coeficientes de los terminos de orden superior al segundo? Se podra pensar
1/3, 1/4, . . . , 1/n, en realidad no es as, y para darse cuenta de eso es sufi-
ciente repetir el razonamiento de arriba para n = 3 obteniendo la f ormula
incorrecta f 000 (x0 ) ' 3 2 f 000 (x0 ), entonces esta vez el coeficiente que ten-
1
emos que introducir en la formula (3.3) para anular el factor 3 2 es 32 .
Iterando el razonamiento hasta el orden n, llegamos a obtener la f ormula
incorrecta f (n) (x0 ) ' n (n 1) (n 2) 2 f (n) (x0 ) y entonces habra que
corregir la f ormula (3.3) como sigue:
81
on que significa que la egualdad vale x U (x0 ) y
donde = es una notaci
xx0
donde o((x x0 )n ) es un infinitesimo superior con respecto a (x x0 )n :
o((x x0 )n )
0.
(x x0 )n xx0
es decir, explicitamente,
Observaciones:
1. La f
ormula de Taylor, independientemente del resto, es una de las
f
ormulas mas importantes de la matematica: permite aproximar una
funcion complicada con la funci on no lineal m as sencilla, un
polinomio, por menos de infinit esimos de orden superior. Ver-
emos que esta propiedad resultara extremadamente u til tambien a la
hora de resolver lmites en forma indeterminada;
82
Como muestra la figura superior, estos detalles finos permiten aprox-
imar cada vez mejor la funcion tambien lejos del punto central del
desarrollo de Taylor. Cerca del centro del desarrollo, 0 en este caso,
todos los polinomios dan una aproximacion excelente de la funcion
y = sin(x) (derivable infinitas veces con continuidad), a medida que
nos alejamos de 0, los polinomios de orden mas grande aproximan la
funci
on con un error cada vez mas pequeno.
4. La f
ormula de Taylor de orden 1 con resto de Peano coincide con la
(1)
on de diferenciabilidad: f (x) = Tx0 + o((x x0 )), con
definici
xx0
Tx(1)
0
= f (x0 )+f 0 (x0 )(xx0 ) Polinomio de Taylor de orden 1.
83
f (x) Desarrollo de Taylor-Mac Laurin de orden n
2 n xk
ex 1 + x + x2! + . . . + xn! + o(xn ) = nk=0 + o(xn )
P
k!
3 5 x2n+1 x2k+1
x x3! + x5! + . . . + (1)n (2n+1)! + o(x2n+1 ) = nk=0 (1)k + o(x2n+1 )
P
sin x
(2k + 1)!
2 4 x2n x2k
1 x2! + x4! + . . . + (1)n (2n)! + o(x2n ) = nk=0 (1)k + o(x2n )
P
cos x
(2k)!
3 5 x2n+1 x2k+1
x + x3! + x5! + . . . + (2n+1)! + o(x2n+1 ) = nk=0 + o(x2n+1 )
P
sinh x
(2k + 1)!
2 4 x2n x2k
1 + x2! + x4! + . . . (2n)! + o(x2n ) = nk=0 + o(x2n )
P
cosh x
(2k)!
2 3 n xn
x x2 + x3 . . . + (1)n1 xn + o(xn ) = nn=1 (1)n1 + o(xn )
P
log(1 + x)
n
(1 + x) 1 + x + (1) x2 + . . . + n xn + o(xn ) = nk=0 k xk + o(xn )
P
2
Observamos que:
1. En el desarrollo de Taylor-Mac Laurin de funciones pares aparecen solo
potencias pares, en el desarrollo de Taylor-Mac Laurin de funciones
impares aparecen solo potencias impares;
2. El desarrollo de Taylor-Mac Laurin de sin x y sinh x son identicos,
exceptuado el signo alterno
(
1 si n es par
(1)n =
1 si n es impar.
que est
a presente en el desarrollo de sin x;
3. El desarrollo de Taylor-Mac Laurin de cos x y cosh x son identicos,
exceptuado el signo alterno que esta presente en el desarrollo de cos x;
4. En el desarrollo de log(1 + x) aparece
(
1 si n es impar
(1)n1 = (1)n+1 =
1 si n es par.
84
En particular:
( 1)
= = 1, = = 1, = .
n 0 1 n1 2 2
(x)2
e(x) 1 + (x) , e(x) 1 + (x) +
xx0 xx0 2
(x)3
sin (x) (x) , sin (x) (x)
xx0 xx0 3!
(x)2
cos (x) 1 , cos (x) 1
xx0 xx0 2!
(x)3
sinh (x) (x) , sinh (x) (x) +
xx0 xx0 3!
(x)2
cosh (x) 1 , cosh (x) 1 +
xx0 xx0 2
(x)2
log(1 + (x)) (x) , log(1 + (x)) (x)
xx0 xx0 2
2
(1 + (x)) 1 + (x) , (1 + (x)) 1 + (x) + (x)2 .
xx0 xx0 2
Ejemplo: f (x) = e5 x2 , es evidente que (x) = 5 x 2 0, entonces
x2
el derarrollo de
Taylor de f (x) alrededor de x = 2 hasta el orden 2 es
f (x) 1 + 5 x 2 + 252 |x 2|.
x2
85
3.8. El teorema de lH opital y sus aplicaciones al
c
alculo de lmites en forma indeterminada
El desarrollo de Taylor y el teorema de lHopital, que vamos a examinar
en esta secci
on, son los dos metodos mas potentes para resolver los lmites
que se presentan en forma indeterminada.
Teorema de lH
opital :
Hip
otesis:
2. c (a, b);
f 0 (x)
lm = L R .
xc g 0 (x)
86
El teorema dice que, bajo sus hipotesis, el lmite del cociente entre las
derivadas coincide con el lmite del cociente de las funciones de partida.
Sin embargo, si el lmite del cociente de las derivadas no existe, no
se puede decir nada sobre el lmite del cociente de las funciones de
partida!
lm x log x = 0 .
x0+
porque cos x no tiene lmite para x , como ya hemos visto. Por lo tanto
no todas las hip otesis del teorema de lHopital son validas y tenemos que
usar otra tecnica para calcular el lmite:
sin x
x sin x x(1 x ) x
lm = lm sin x
= lm = 1.
x+ x + sin x x+ x(1 + x+ x
x )
Como se ve, hubiera sido incorrecto inferir que el lmite original no existe
solo porque el lmite del cociente de las derivadas no existe!
87
Observaci on: el teorema de lHopital se puede aplicar tambien de forma
iterada (si siguen valiendo sus hipotesis, obviamente), o sea se puede intentar
calcular el lmite indeterminado con la derivada segunda, tercera, etc.
A veces estas iteraciones pueden generar un bucle infinito, como en el
caso siguiente:
x2 + 1 f (x)
lm = ,
x+ x g(x)
otesis del teorema de lHopital, pero f 0 (x) =
valen las hip x
x2 +1
y g 0 (x) = 1,
f 0 (x)
entonces = x sigue teniendo la misma forma de indeterminacion.
g 0 (x) x2 +1
00
Si derivamos otra vez obtenemos fg00 (x)
(x) 2
= xx+1 y as hemos vuelto al lmite
inicial. Se genera un bucle infinito que no permite resolver el lmite. Hay
que usar otra tecnica (el desarrollo de Taylor, como veremos en las clases de
pr
acticas y seminarios).
Una de las aplicaciones mas u tiles del teorema de lHopital consiste en
la determinacion de la jerarqua de infinitos e infinitesimos, vamos a ver de
que se trata.
88
Teorema:
Jerarqua de infinitos:
loga x (loga x) x x x bx xx , x +
Jerarqua de infinit
esimos:
1 1 1 1 1 1 1
x x , x +
loga x (loga x) x x x b x
donde es el smbolo que expresa que el infinito o infinitesimo a su derecha
es de orden superior al infinito o infinitesimo que aparece a su izquierda y
donde:
a > 1;
> 1;
(0, 1);
> 1;
b > 1;
> 0.
Observaciones importantes:
1
x
1. Observando que bx = 1b , b > 1 as que 1b = c (0, 1), podemos
1
sustituir bx en el orden de infinitesimos crecientes con cx , 0 < c < 1;
1
2. x
= x1 , (0, 1), x1 = x1 , 1, entonces podemos reescribir
la jerarqua de esos infinitesimos considerando el lmite para x 0
como sigue:
x x , 0 < < , x 0,
es decir: potencias m as grandes de x tienden a cero m as rap-
idamente que potencias m as peque nas de x para x 0 (volver
afico de la seccion 2.7 que compara x2 y x3 alrededor de
a mirar el gr
0!).
89
Veamos unos ejemplos de aplicacion de la jerarquas de infinitos e infinitesimos:
Calcular
85 log x 10x f (x)
lm = lm .
x+ 5 x + 19x 14 x+ g(x)
Calcular
lm (log 2)x x40 .
x+
Dado que la base de log es e ' 2,7 y que 2 < e, log 2 (0, 1), entonces
(log 2)x 0 para x , mientras que x40 para x , por
tanto el lmite se presenta con la forma de indeterminacion [0 ].
Podemos reescribir el lmite as:
(log 2)x (log 2)x
lm (log 2)x x40 = lm 1 = lm =0
1 40
x+ x+ x+
x40 x
x40 x40
lm (log 2)x x40 = lm 1 = lm x = 0
x+ x+ x+ 1
(log 2)x log 2
En las clases de pr
acticas y seminarios se vera que combinando opor-
tunamente desarrollo de Taylor, criterio asint otico y jerarqua de
infinitos e infinitesimos, se pueden resolver todos los lmites que
se presentan en forma indeterminada.
90
3.10. Teoremas fundamentales sobre derivadas
Vamos ahora a examinar unos teoremas fundamentales sobre las derivadas
de funciones reales de variable real. Con estos teoremas completaremos el
analisis de todos los resultados basicos que sirven para estudiar el compor-
tamiento gr afico de las funciones y los aplicaremos a problemas de opti-
mizacion.
Comencemos definiendo los extremos de una funcion.
Los puntos de m aximos y mnimos se llaman, con una palabra, puntos ex-
tremos de f . Una funcion puede tener mas extremos locales, como se ve en
las figuras siguientes.
91
Prueba. Haremos la prueba suponiendo que x0 sea un punto de mnimo
local, si es un m
aximo la prueba es analoga.
Si x0 es un mnimo local, existe un entorno abierto U (x0 ) tal que, para cada
x U (x0 ) podemos escribir f (x0 ) f (x), entonces f (x) f (x0 ) 0.
Observamos tambien que podemos representar x como x = x0 + h, donde h
es un n umero oportuno 0. Seg un el signo de h tenemos:
h > 0, x x0 = x0 + h x0 = h > 0, entonces
f (x) f (x0 ) f (x0 + h) f (x0 ) f (x0 + h) f (x0 )
= 0 = lm 0
x x0 h h0+ h
por el teorema de permanencia del signo;
92
3.11. El teorema del valor medio de Lagrange y
sus aplicaciones al estudio de la monotona
de una funcion
Supongamos haber recorrido en coche los 160 km que separan dos puntos
A y B en 2 horas. La velocidad media con la cual hemos viajado es 80 km/h.
Si hemos viajado de forma continua, sin parar o tener otras irregularidades
durante el trayecto, entonces en algun instante t0 (0, 2horas) la velocidad
con la cual hemos viajado ha coincidico con la velocidad media de 80 km/h,
o sea:
f (b) f (a)
= f 0 (c).
ba
93
en el punto c (puede haber mas de un punto!) la recta tangente al grafico
de f es paralela a la recta secante al grafico de f que conecta A (a, f (a))
y B (b, f (b)).
f (b) f (a)
y = f (a) + (x a)
ba
e introducir la funci
on auxiliar
f (b) f (a)
g(x) = f (x) f (a) + (x a) .
ba
f (b) f (a)
g 0 (x) = f 0 (x) ,
ba
entonces demostrar que f satisface la tesis del teorema de Lagrange equivale
a demostrar que existe por lo menos un punto c (a, b) tal que g 0 (c) = 0,
porque, en ese caso, 0 = f 0 (c) f (b)f
ba
(a)
, o sea f 0 (c) = f (b)f
ba
(a)
.
Siendo g continua en un conjunto compacto [a, b], por el teorema de Weier-
strass g tiene mnimo y maximo absolutos en [a, b]: x1 , x2 [a, b] tales
que
m = g(x1 ) = mn {g(x)}
x[a,b]
94
Vamos a ver un corolario (consecuencia inmediata) del teorema de La-
grange.
f (x2 ) f (x1 )
= f 0 (c).
x2 x1
Dado que hemos hecho da hipotesis de que f 0 (c) 0 y que x2 > x1 , sigue
que f (x2 ) f (x1 ) 0, o sea f (x2 ) f (x1 ), o sea f es creciente.
Las otras dos implicaciones inversas se demuestran con esta misma tecnica.
2
Los resultados de esta seccion nos proporcionan una tecnica para buscar
extremos de funciones reales de variable real.
95
Busqueda de extremos:
Sea f : D R R, para buscar sus extremos hay que:
f 0 (x) = 0.
96
El estudio del signo de f 0 alrededor de x0 podra resultar complicado,
por eso resulta u til tener un metodo alternativo para averiguar si un punto
estacionario es un punto extremo. En un punto estacionario x0 de f , la
f
ormula de Taylor de f con centro en x0 no tiene el termino lineal porque
f 0 (x0 ) = 0. Entonces, si f es derivable dos veces en un punto estacionario,
el polinomio aproximante de Taylor alrededor de x0 se escribe as:
f 00 (x0 )
Tx(2)
0
(x) = f (x0 ) + (x x0 )2 , x0 punto estacionario .
xx0 2
Como ense
na la figura siguiente
97
f 00 (x0 ) > 0: concavidad hacia arriba;
98
Tenemos todas las informaciones para poder desarrollar el esquema para
el estudio grafico de una funcion. Aunque los ordenadores pueden propor-
cionar f
acilmente gr aficos de funciones, es importante saber estudiar cualita-
tivamente y cuantitativamente el grafico de una funcion, dado que a menudo
los gr
aficos contienen informaciones extremadamente importantes, como el
estudiante podr a apreciar tanto en las practicas y seminarios de este curso,
como en otras asignaturas de su carrera.
99
Estudio gr
afico de funciones
1. Determinar el dominio D de f ;
f (x)
lm =m, lm [f (x) mx] = q .
x x x
100
tiene puntos extremos y si son maximos o mnimos. Si el estudio del
signo de f 0 es muy complicado, puede resultar u
til el test de monotona
de la derivada segunda;
9. Se aconseja dibujar el grafico paso por paso, introduciendo cada vez las
informaciones que cada punto proporciona. Si se ve que hay contradic-
ciones entre estas informaciones hay que verificar que las cuentas sean
correctas. Por eso insistimos sobre el punto 5. del analisis asintotico en
los extremos, que a menudo ayuda a mantener una coherencia global
en el estudio gr
afico de la funcion.
101
3.11.1. La f
ormula de Taylor con resto de Lagrange
Otra consecuencia notable del teorema del valor medio de Lagrange es
la posibilidad de escribir el desarrollo de Taylor con un resto diferente con
respecto a lo de Peano. Para entender porque es suficiente re-escribir la tesis
del teorema de Lagrange as: f (b) = f (a) + f 0 (c)(b a), pero si ponemos
a = x0 y b = x = x0 + h, obtenemos la siguente aproximacion lineal de f en
x0 + h:
f (x) = f (x0 ) + f 0 ()(x x0 ),
con (x0 , x0 + h). Esta formula establece una linealizacion de f (bajo
las hip
otesis del teorema de Lagrange obviamente) alternativa al desarrollo
de Taylor de orden 1 con resto de Peano y se llama formula de Taylor con
centro en x0 y resto de Lagrange.
Tambien esta formula se puede extender a una aproximacion polinomial
mas precisa, puesto que la funcion f satisfaga las condiciones de derivabili-
dad contenidas en el siguiente teorema.
102
3.12. La optimizaci
on
Se llama optimizaci on un campo muy amplio de las matematicas ded-
icado a la busqueda de extremos de funciones u objetos mas complicados
llamados funcionales. La optimizacion tiene un radio de aplicacion practi-
camente universal.
En el curso de c alculo del segundo trimestre examinaremos las fun-
ciones reales de varias variables y eso nos dara la oportunidad de mostrar
aplicaciones de optimizaci on muy interesantes. Sin embargo, ya con fun-
ciones reales de una sola variable real podemos discutir un problema de
optimizacion que viene de la vida practica y que dara una idea de la im-
portancia de las tecnicas de optimizacion en matematica aplicada. Veremos
otros ejemplos en las clases de practicas y seminarios.
103
gracias al teorema de Pit
agoras, se puede escribir como
p
L(a, b) = a2 + b2 ,
que tiene sentido cuando a > 2 (vease la figura arriba). Derivando L con
respecto a a tenemos:
2
2a 1 (a2) 3
L0 (a) = r .
2
a
2 a2 + a2
entonces a = 3 2 + 2 es un punto de mmimo para L(a). Si L(a ) es mayor
de 4, el palo puede pasar por la ca
nera, si no, el palo se queda atrapado.
Afortunadamente para el fontanero L(a ) ' 4, 16 y el palo podra pasar por
la ca
nera.
104
Captulo 4
El m
etodo de Newton 1-D
4.1. Construcci
on del m
etodo de Newton 1-D
La idea principal del metodo de Newton consiste en observar que los ceros
de las funciones lineales son extremadamente faciles de calcular: f (x) =
ax + b = 0 si y solo si x = b/a. Si una funcion es derivable, sabemos
que, en cada punto x0 de su dominio, puede ser linealizada a traves del
proceso de diferenciaci on: en un entorno suficientemente peque no de x0 la
funci
on f coincide con su recta tangente (o polinomio de Taylor de orden 1)
a menos de infinitesimos de orden superior (con un resto que tiende a cero
cuando nos aproximamos a x0 ). Dado que la aproximacion de f por su recta
tangente es valida solo localmente, sera esperar demasiado que el cero de f
coincidiera con el cero de su recta tangente; sin embargo, veremos que si se
itera la linealizaci
on de f en los puntos determinados por los ceros de sus
rectas tangentes, se aproximara cada vez mejor el valor del cero de f . Esto
explica por que el metodo de Newton se suele llamar tambien metodo de las
tangentes o metodo de las linealizaciones iteradas.
Comencemos presentando la idea geometrica subyacente al metodo con
105
la figura 4.1.
f (xk1 )
xk = xk1 ,
f 0 (xk1 )
106
4.2. Aplicacion al calculo de los valores decimales
de las races cuadradas
Presentamos una aplicacion notable del metodo de Newton 1-D: el calcu-
lo de races cuadradas.
p Vamos a ver c omo se puede determinar por ejemplo el valor decimal de
p 3/4, que es 0,8660 . . . con el esquema de Newton. Tenemosp que interpretar
3/4 como el cero de una funci pon, para ello definimos x = 3/4 y definimos
unap funcion de x que tenga 3/4 como cero. No podemos definir f (x) =
x 3/4 porque su derivada vale 1, entonces la recta tangente tiene siemprep
la misma pendiente y no hay manera de converger iterativamente a 3/4.
Para evitar este problema podemos elevar al cuadrado la definicion de x
obteniendo x2 = 3/4, o sea 4x2 3p= 0. Si definimos f (x) = 4x2 3,
los ceros de f corresponden a x = 3/4. Vamos a aplicar el metodo de
Newton: necesitamos la derivada primera de f , que es f 0 (x) = 8x, y un
punto inicial. Dado que 3/4 esta entre 0 y 1 y dado que la raiz cuadrada
de un n umero que est a entre 0 y 1 sigue en este intervalo, elegimos por
0
ejemplo x = 0,5 para no estar demasiado lejos de la solucion y converger
mas rapidamente.
Iteraci
on 1:
f (0,5) 2
x1 = 0,5 = 0,5 = 1;
f 0 (0,5) 4
Iteraci
on 2:
f (1) 1 7
x2 = 1 0
= 1 = ' 0,875;
f (1) 8 8
Iteraci
on 3:
7 f(7)
x3 = 0 87 ' 0,8660.
8 f (8)
107
on {xk }k0 se suele utilizar el error relativo entre dos
los puntos de la sucesi
iteraciones sucesivas:
|xk xk1 |
Ek .
|xk |
Se suelen parar las iteraciones del algoritmo de Newton despues de un
umero de iteraciones k tales que xk y xk1 sean cercanas. Formalmente
n
se puede expresar como sigue:
Ek T,
siendo T > 0 un valor de tolerancia que el usuario del algoritmo elige segun
sus exigencias de precisi
on. Otra opcion sera parar el algoritmo de Newton
umero de iteraciones k tales que
despues de un n
|f (xk )| T,
108
Figura 4.2: Oscilaciones en el metodo de Newton.
Est
a claro que el algoritmo de Newton converge si y solo si se cumple que
porque en ese caso el valor absoluto del error disminuye en cada iteracion,
hasta llegar a determinar el cero x con la precision deseada.
Queremos utilizar la condicion |ek+1 | < |ek | para determinar una restric-
on sobre el valor inicial x0 que nos garantize la convergencia del algoritmo
ci
de Newton.
Empecemos definiendo el intervalo de b usqueda del cero x de f como
sigue: I = [x x0 , x + x0 ] y supongamos que f sea dos veces derivable con
derivadas continuas en I. Bajo esta u ltima condicion podemos escribir, para
cada x I, el desarrollo de Taylor de f centrado en xk I hasta el orden
1 con resto de Lagrange para cada paso k del algoritmo:
1
f (x) = f (xk ) + f 0 (xk )(x xk ) + f 00 ()(x xk )2 , con | xk | < |x xk |.
2
Dado que este desarrollo vale para cada x I, podemos seleccionar en
particular x = x , en este caso el desarrollo de Taylor se escribe:
1
f (x ) = f (xk )+f 0 (xk )(x xk )+ f 00 ()(x xk )2 , con | xk | < |x xk |.
2
on que f (x ) = 0 y dividimos ambos miembros por
Si tenemos en consideraci
f (x ) (suponiendo que f 0 (xk ) 6= 0), obtenemos
0 k
f (xk ) k 1 f 00 ()
0= + x x + (x xk )2 ,
f 0 (xk ) 2 f 0 (xk )
es decir
f (xk ) 1 f 00 ()
x xk + = (x xk )2 .
f 0 (xk ) 2 f 0 (xk )
109
k
Ahora usamos el hecho de que xk+1 = xk ff0(x )
(xk )
(dado por el paso basico
de las iteraciones del metodo de Newton) para reescribir la u ltima ecuacion
como:
1 f 00 ()
x xk+1 = 0 k (x xk )2 .
2 f (x )
La ventaja de esta ultima f ormula es que pone en evidencia la relacion entre
el error ek+1 y el error ek :
1 f 00 () 2
ek+1 = e ,
2 f 0 (xk ) k
tomando el valor absoluto en ambos lados:
00
f ()
|ek+1 | = 0 k |ek |2 ,
(4.1)
2f (x )
lo cual expresa el hecho de que el error absoluto e en el paso k + 1 es
proporcional al cuadrado del error
00 absoluto
e en el paso k, con coeficiente
f ()
de proporcionalidad dado por 2f 0 (xk ) . El valor maximo que puede tomar
este coeficiente es1 :
1 f 00 (t)
M max 0 ,
tI 2 f (t)
donde, por sencillez, hemos usado la variable generica t I en lugar de las
dos variables , xk I. Entonces la ecuacion (4.1) implica que:
|ek+1 | M |ek |2 = M |ek ||ek |, (4.2)
por lo que se cumplir a la condicion necesaria y suficiente de convergencia,
|ek+1 | < |ek |, si y solo si M |ek | < 1. Para que M sea finito necesitamos que
f 0 (t) 6= 0 t I (en particular es necesario que f 0 (x ) 6= 0). Si se verifica
esta condici on, podemos garantizar que la magnitud del error absoluto e
decrece en cada iteraci on (es decir, que |ek+1 | < |ek |) siempre que:
1
M< ,
|ek |
pero como |ek+1 | < |ek |, es suficiente pedir que
1 1
M< = , (4.3)
|e0 | |x x0 |
(aqu es importante que el smbolo de desigualdad < sea estricto). Compro-
bemos esta afirmacion para las primeras dos iteraciones:
1
Recordar que, para maximizar una fracci
on, hay que maximizar el numerador y min-
imizar el denominador.
110
Iteracion 1: k = 0. Sustituyendo k = 0 en la ecuacion (4.2) tenemos
que |e1 | M |e0 |2 , entonces, si M cumple con la propiedad (4.3),
2
|e1 | < |e|e00|| = |e0 |;
Iteraci
on 2: k = 1. Sustituyendo k = 1 en la ecuacion (4.2) tenemos
que |e2 | M |e1 |2 entonces, si M cumple con la propiedad (4.3) y
2 2
teniendo en consideracion que |e1 | < |e0 |: |e2 | < |e|e10|| < |e|e11|| = |e1 |;
2
|x x0 | < 00 .
maxtI ff 0 (t)
(t)
suficientemente pr
oximo a la solucion).
111
Dejamos como ejercicio mental razonar sobre el valor del maximo cuando
f 0 y f 00 sean mon otonas y negativas.
El algoritmo de Newton, cuando converge, lo hace cuadr aticamente, es
decir, la distancia entre el cero de f y la solucion proporcionada en el paso
k decrece cuadr aticamente a medida que aumentamos k. Esa es una de las
propiedades m as importantes del metodo de Newton: si converge, lo hace
muy r apidamente. La contrapartida es que tienen que cumplirse las condi-
ciones de f especificadas anteriormente y el punto inicial x0 no tiene que
estar demasiado lejos del cero x . Esto hace que sea mas difcil poder aplicar
el metodo de Newton a problemas donde no se conoce el intervalo de valores
probables donde se puede encontrar x . Por esa razon, un estudio cualitativo
cuidadoso antes de aplicar el metodo de Newton puede resultar extremada-
mente u til, como veremos en el ejemplo siguiente.
4.5. Aplicacion al c
alculo aproximado de los pun-
tos de intersecci
on entre gr
aficos de funciones
El metodo de Newton resulta particularmente u til a la hora de aproximar
la posici
on de los puntos de interseccion entre graficos de dos funciones. Sean
f1 y f2 dos funciones derivables, en los valores x por los cuales los graficos
de f1 y f2 se intersectan tendremos obviamente f1 (x) = f2 (x), por lo tanto
encontrar las intersecciones entre f1 y f2 es equivalente a encontrar los ceros
de la funcion F (x) = f1 (x) f2 (x), que se puede hacer con el metodo de
Newton.
Veamos una aplicaci on: supongamos de querer estudiar graficamente la
funcion f (x) = e x+2 13 x3 . Sabemos que los intervalos de crecimiento y
decrecimiento de f se determinan a traves del signo de la derivada primera
de f , o sea f 0 (x) = ex+2 x2 . La expresion analitica de f 0 (x) no es trivial
y la inecuaci on f 0 (x) 0, o sea ex+2 x2 , que determina los valores de
x por los cuales f es creciente o decreciente, respectivamente, no tiene una
solucion analtica explicita. En estos casos es aconsejable actuar as:
112
x [2, 0] donde ex+2 = x2 ;
113
F 0 tiene un mnimo en x = log 22 en el cual toma el valor de 4, 61 > 0,
entonces: 00
F 00 (0) e2 1
F (x)
max 0
= = .
x[2,0] F (x) F 0 (log 2 2) 4, 16
Entonces la condici
on de convergencia se escribe como:
2 8, 32
|x0 x | < e2 1
= ' 1, 3,
e2 1
4,16
eso nos dice que si escogemos un valor inicial x0 que esta alejado del
on x por menos de 1,3, seguramente el algoritmo
punto de intersecci
114
de Newton convergera a x . Dado que el intervalo de busqueda es
I = [2, 0], es suficiente escoger x0 = 1 para estar seguros que la
condicion valga. Inicializando el algoritmo de Newton con x0 = 1 y
tomando como valor de tolerancia 0, 01, encontramos como salida final
x = 1, 37. Dado que ex+2 > x2 solo cuando x > 1, 37 sigue que
f (x) es creciente para cada x (1, 37, +), decreciente para cada
x (, 1, 37) y que x = 1, 37 es un punto de mnimo absoluto
para f .
115
Captulo 5
Integrales
5.1. Definici
on de la integral de Riemann
Consideramos una funcion continua f : [a, b] R y operamos una par-
tici
on de [a, b] en n sub-intervalos:
116
tambien los extremos), y evaluamos la funcion f en j , como en la figura
siguiente. .
Si operamos una particion mas fina de [a, b], por ejemplo esta
la posici
on de los j variara, y con ella el valor de f calculada en esos
puntos y la distancia x entre los extremos de los sub-intervalos. El siguiente
teorema, seg un Riemann, dice que, si f es continua, existe el lmite de n
117
tendiendo a infinito, independientemente de la eleccion de los j , y coincide
con el
area por debajo del grafico de f .
areas elementales de rect angulos con base infinitesima y altura dada por el
valor de la funci
on f en todos los puntos x que recorren el intervalo [a, b].
La genialidad de la notacion de Leibnitz esta en el hecho de que el smbolo
que denota la integral sintetiza toda la construcci on con la cual se llega a la
integral misma!
Obviamente Z a
f (x)dx = 0 a R, f
a
Z a
f impar = f (x)dx = 0 .
a
Nomenclatura:
f (x): funci
on integranda, x: variable de integraci
on.
118
Observaci
on importante: si f cambia de signo (o sea pasa de ser positiva
Rb
a negativa, o viceversa), a f (x)dx representa una suma de areas con
signo, por eso la integral puede tambien ser negativa o nula, como muestra
la figura siguiente (autoexplicativa).
119
5.2. Propiedades de la integral
Las propiedades de las integrales siguen de las analogas propiedades de
los lmites:
Linealidad: dados , R,
Z b Z b Z b
[f (x) + g(x)]dx = f (x)dx + g(x)dx ;
a a a
Si a c b,
Z b Z c Z b
f (x)dx = f (x)dx + f (x)dx ;
a a c
120
5.3. Valor medio y valor eficaz de una funci
on
La media aritmetica de un conjunto finito de valores viene dada por
la suma de esos valores dividido por el n umero de valores. Si queremos
calcular la media de una magnitud que vara de forma continua tenemos
que generalizar la definicion antecedente con el uso de las integrales como
en la definici
on siguiente.
Z b 1/2
1 2
fE = valor eficaz de f en [a, b] = [f (x)] dx .
ba a
1
Rb
entonces ba a f (x)dx est
a contenido entre el mnimo y el maximo valor
que f toma en [a, b]. Siendo f continua, en virtud del teorema de los valores
Rb
1
intermedios tiene que existir c [a, b] tal que ba a f (x)dx = f (c). 2
121
5.4. El teorema fundamental del c
alculo diferen-
cial
Como veremos m as adelante, la definicion de integral que hemos dado
se presta muy bien para encontrar formulas de aproximacion numerica de
las integrales. Sin embargo, el calculo analtico (o sea a traves de funciones
elementales) de las integrales se hace de manera completamente diferente y
necesita el concepto de primitiva o anti-derivada.
122
Observaci on importante: aunque f sea continua, no siempre es posible es-
cribir su primitiva a traves de una composicion de funciones elementales, en
partcular, se demuestra rigurosamente (como corolario de un teorema de-
bido al matem atico frances Joseph Liouville (1809-1882)) que la primitiva
2
de las funciones del tipo f (x) = eax con a 6= 0, en partcular de la funci
on
x 2
Gaussiana f (x) = e , no se puede escribir a trav es de la composi-
ci
on de funciones elementales. Veremos despues que una primitiva de la
Gaussiana se puede escribir de forma implcita.
123
existe j (xj1 , xj ) tal que
124
con f
acilidad problemas que, sin esas, necesitaran calculos mucho mas com-
plicados.
Observamos ahora que, si derivamos ambos lados de la ecuacion (5.1),
Rb
obtenemos2 a f 0 (x)dx = G0 (b) G0 (a) pero G0 (b) G0 (a) = f (b) f (a),
entonces
Z b
f 0 (x)dx = f (b) f (a) = [f (x)]ba . (5.3)
a
125
estudiantes a tener siempre presente que la forma correcta de interpretarlo
es a traves de la ecuaci
on (5.3).
R Rb
Observaci on: f (x)dx y a f (x)dx tienen dos significados completa-
mente distintos:
R
f (x)dx es un conjunto de funciones (las primitivas de f );
Rb
a f (x)dx es un n
umero real!
Leyendo al reves la tabla de derivacion de funciones elementales, podemos
obtener gratuitamente la tabla de las integrales inmediatas:
R
f (x) f (x)dx
k kx + c
x+1
x +1 + c, 6= 1
1
x log |x| + c
sin x cos x + c
cos x sin x + c
1
cos2 (x)
= 1 + tan2 x tan x + c
ex ex + c
ax
ax log a , a > 0, a 6= 1
1
1+x2
arctan x + c
1 arcsin x + c
1x2
sinh x cosh x + c
cosh x sinh x + c
126
Observese que en el segundo caso hay que poner -1 al numerador porque la
derivada de x es -1. La primitiva log |x| tiene en realidad que interpretarse
como la pareja de primitivas log(x) para x < 0 y log(x) para x > 0.
Para encontrar primitivas de funciones mas complicadas resultan ex-
tremadamente u tiles estos dos teoremas.
127
f (x)g 0 (x) g 0 (x)dx
R
f (x) g(x) =
xn+1
xn log(mx + q), log(mx + q) n+1
1
xn sin(mx) xn m cos(mx)
1
xn cos(mx) xn m sin(mx)
1 mx
xn emx , n > 0 xn me
1
enx sin(mx) enx m cos(mx)
1
enx cos(mx) enx m sin(mx)
R R
Ejemplo: log x dx. Interpretamos esta integral como 1log x dx, tomamos
como factor finito f (x) R= log x, f 0 (x)R = x1 y como factor diferencial g 0 (x)dx =
1 dx, entonces g(x) = g 0 (x)dx = 1 dx = x + c. Seleccionamos c = 0 por
sencillez3 y usamos la f ormula de integracion por partes:
Z Z Z
1
log x dx = x log x x dx = x log x 1 dx = x log x x + c,
x
entonces Z
log x dx = x log x x + c .
Se ver
an muchos otros ejemplos en las clases de practicas y seminarios.
128
por tanto, G(x) es una primitiva de f , si y solo si, G((t)) es una primitiva de
f ((t))0 (t). Usando la notacion de integral indefinida este hecho se traduce
escribiendo Z Z
f (x)dx = f ((t))0 (t)dt ,
129
La f
ormula de integraci on por sustitucion resulta partcularmente u til y
f
acil de implementar cuando x aparece en una combinacion lineal: x + ,
, R. En este caso, poniendo t = x + y despejando con respecto a
x obtenemos x = t b, entonces dx = 1 dt. Se ve que el diferencial dx y
dt en este caso difieren s olo por una constante multiplicativa que tiene el
significado de cambio de escala.
Vamos a ver porqueR 1 x hablamos de cambio de escala considerando por
ejemplo la integral 0 100 dx y resolviendola sin sacar afuera 1/100 sino us-
x
ando la f ormula de sustitucion de variable, t = 100 , x = 100t, dx = 100dt,
eso significa que, si medimos la distancia infinitesima dx en metros, para que
valga la igualdad dx = 100dt, tenemos que medir la distancia infinitesima
dt en centmetros. La consecuencia es que cuando x vara de 0 a 1 metro, t
variara de 0 a 100cm, que es la misma longitud, pero medida con una escala
diferente!
En general, un cambio de variable lineal t = x + genera un cambio de
escala homogeneo en todo el intervalo de integraci on, que se realiza a traves
del factor de conversi on 1/.
Un cambio de variable no lineal t = 1 (x) genera un cambio de escala
que vara, en general, punto por punto en el intervalo de integraci on y que
se realiza a traves del factor de conversi on 1/0 (t).
En las clases de practicas y seminarios se veran las tecnicas de sustitucion
m as utilizadas en la resolucion de las integrales. Aqu analizamos solo una de
esta tecnicas: consideramos la integral f (x)f 0 (x)dx y operamos
R
el cambio
0 entonces f (x)f 0 (x)dx =
R
de
R variable t = f (x), dt = d(f (x)) = f (x)dx,
tdt, que es una integral inmediata y vale t /2 + c = [f (x)]2 /2+c, c R.
2
[f (x)]2
Z
f (x)f 0 (x)dx = +c ;
2
[f (x)]+1
Z
f (x)f 0 (x)dx = +c , 6= 1;
+1
Z 0
f (x)
dx = log |f (x)| + c .
f (x)
130
5.5. Integrales y din
amica de las partculas
En esta seccion pondremos en evidencia la importancia de las integrales
en la dinamica de las partculas. Consideramos la funcion posicion de una
partcula que se mueve en un eje, supongamos el eje horizontal: s : [0, T ]
R, t 7 s(t) = x, posici on de la partcula en el instante t, con s(0) = a
y s(T ) = b. Como sabemos, la longitud del recorrido de la partcula es
Rb
a dx, pero dado que x = Rb
s(t), podemos usar la formula de integracion
RT RT
por sustitucion y escribir a dx = 0 ds(t) = 0 s(t)dt,
donde s(t)
= v(t)
es la velocidad instantanea de desplazamiento de la partcula, quedando
demostrado que
Z T
Espacio recorrido en [0, T ] = v(t)dt
0
131
5.6. Calculo de la longitud del gr
afico de una fun-
ci
on
En esta seccion mostraremos como se puede calcular la longitud del
gr
afico de una funci
on a traves de las integrales.
Sea y = f (x) una funci
on continua con derivada continua en [a, b]. Vamos
a considerar la contribucion a la longitud del grafico dada por la parte del
gr
afico entre dos puntos x y x + dx (dx infinitesimo positivo), tal y como se
observa en la figura siguiente:
= x3/2 para
Ejemplo: Calculamos la longitud del grafico de la funcion y q
a
x [0, a], a R. Dado que y 0 = 23 x obtenemos L = 0 1 + 94 xdx,
R
132
Observaci
on importante: la forma de razonar que acabamos de ver es
comun en muchas aplicaciones. A menudo es demasiado difcil resolver un
problema no lineal globalmente, entonces lo que se hace es considerar la
linealizaci
on local del problema en un entorno infinit esimo de un
punto, este problema suele ser m as facil y resoluble, la soluci on
que se obtiene se extiende del infinitamente peque no al finito con
el proceso de integraci on.
133
Sea ahora f : [a, b) una funcion continua que tiene una asintonta ver-
tical en x = b, o sea lmxb f (x) = + como en la figura debajo
134
Sea ahora f : [a, ) una funcion continua sobre un intervalo ilimitado,
como en la figura siguiente
Ejemplos notables:
(
Z b + si 1
1
dx = (ba)1
a (b x) 1 si < 1
(
+
+ si 1
Z
1
dx =
1 x 1
1 si > 1
Z +
Z + (x)2
2 1
ex dx = , e 22 dx = 1 ,
2
Podemos observar que el hecho de que la u ltima integral valga 1 explica el
p
orque de los coeficientes en la funcion Gaussiana: el area debajo de una
distribuci
on de probabilidad debe dar 1 (probabilidad del 100 %!).
135
5.8. La funci
on integral
Dada una funci on f : [a, b] R continua, sabemos que su integral defini-
Rb
da4 a f (t)dt es un numero real. Si fijamos el extremo inferior de integracion
a y dejamos libre el valor del extremo superior, llamandolo x, obtenemos una
funci
on llamada funci on integral:
F : [a, b] R Rx
x 7 F (x) = a f (t)dt,
su funci
on integral, entonces F es derivable y
136
Este teorema tiene consecuencias interesantes:
funci
on que se puede prolongar con continuidad a la funcion valor
absoluto, que es continua en todo R pero no derivable en x = 0. Pode-
mos resumir estas consideraciones diciendo que integrar una funcion
aumenta su grado de regularidad y puede arreglar sus discon-
tinuidades, interpretando esta frase en el sentido que le hemos dado
arriba;
137
5.8.1. Una aplicaci
on de la funci
on integral: la ecualizaci
on
de histogramas
La funci
on integral aparece en una de las transformaciones mas comunes
del procesado de se nales digitales: la ecualizacion de histogramas. Aunque
las variables que aparecen en las se nales digitales sean discretas (o sea tienen
un paso entero), la teora que sigue resultara mas facil si suponemos que la
variable r que describe un nivel de intensidad de la se nal f es continua y vara
en el intervalo normalizado [0, 1] (cada intervalo finito se puede normalizar
al intervalo [0, 1], as que esta hipotesis no es restrictiva). Indicamos con |f |
el numero de muestras de f (o su longitud, area o volumen, si consideramos
f como una se nal continua 1-D, 2-D o 3-D, respectivamente).
El histograma h de la se nal es una funcion h : [0, 1] R+ 0 que asocia
a cada nivel r [0, 1] el n umero de veces que el se nal toma ese mismo
valor (por ejemplo, si se mide la intensidad 0,5 de f 20 veces, entonces
h(0,5) = 20).
Muchas veces un se nal tiene un histograma no balanceado, en el sentido
que unos pocos valores aparecen muchas veces (donde el histograma presenta
un pico) y otros pocas o nunca (donde el histograma vale 0 o casi 0). En
cambio, si una se nal tiene un histograma homogeneo, significa que todos los
niveles aparecen con la misma frecuencia. Es facil enteder que esta u ltima
se
nal tiene un contenido informativo superior a la primera.
La ecualizaci on de histograma es la transformacion : [0, 1] [0, 1],
r 7 (r) = s que modifica la distribucion de los niveles de una se nal f
para hacer que su histograma sea lo mas homogeneo posble. Para evitar la
inversion del orden entre los niveles de f se suele asumir que 0 (r) 0 r,
o sea, es una funci on monotona no decreciente.
Para determinar analticamente resulta u til observar que el n umero
total de pxeles tienen que conservarse durante la transformacion, es decir,
el n umero de veces que aparecen los niveles en el rango de intensidades
[r, r + dr] en la se nal original tiene que coincidir con el n umero de veces
que aparecen los niveles en el rango de intensidades [s, s + ds] en la se nal
transformada.
Si llamamos hi el histograma de la se nal original (i=input) y ho el
histograma de la se nal transformada (o=output) lo que acabamos de de-
cir se traduce escribiendo: hi (r)dr = ho (s)ds. Si queremos que ho sea ho-
mogeneo, o sea ho (s) |f |, tenemos que escribir la relacion precedente como
hi (r)dr = |f |ds.
Integrando el lado de izquierda de 0 a r y el lado de derecha de (0) a
138
(r) = s obtenemos:
Z r Z (r)=s
hi (t)dt = |f | ds = |f |(s (0)).
0 (0)
lo cual indica que es 1/|f | por la funcion integral del histograma original de
f , llamado histograma acumulado normalizado de f y representado por
H. Entonces, realiza ecualizacion de histograma mapeando cada nivel r >
0 en el histograma acumulado normalizado de f en r (suma, normalizada,
del n umero de veces que se han detectado los niveles de 0 hasta r):
139
5.9. Relaciones de ortogonalidad de las funciones
sin y cos
En esta secci
on demostraremos unas formulas que el estudiante volvera a
encontrar en el cursos de teora de se
nales en relacion a la transformada
de Fourier.
140
5.10. Aproximaci
on num
erica de integrales
Ya dijimos que no siempre se puede calcular explcitamente la primitiva
de una funci
on a traves de la composicion de funciones elementales, como en
x2
el caso de la Gaussiana f (x) = e 2 . De hecho son muchas mas las funciones
continuas de las que no sabemos calcular explcitamente una primitiva que
las que s se puede!
Por esta raz on hace falta desarrollar tecnicas alternativas para calcular
integrales. Algunas de estas tecnicas se basan en el analisis de funciones
de variable compleja (an alisis complejo) pero son demasiado complicadas
para ser examinadas en un curso de calculo del primer a no. En esta seccion
presentamos tecnicas numericas que se basan en la definicion misma de la
integral y en su relaci on con el calculo del area. Pero en este caso, la definicion
de integral como lmite de sumas de Cauchy-Riemann no es comoda para el
c
alculo de primitivas, pero resulta u til tanto para desarrollar teoremas sobre
integrales como para su aproximacion numerica.
La idea b asica detr
as de las aproximaciones numericas de integrales con-
siste en la sustituci on de la funcion f por una funcion f mas simple para
la cual el c alculo de la integral resulte simple, rapido y eficaz (o sea que el
error de aproximaci on resulte peque no).
En esta secci on veremos tres metodos: integracion por rectangulos, por
trapecios y por par abolas. Como estos tres metodos estan basados en las
sumas de Cauchy-Riemann, conviene recordar la notacion antes de entrar
en detalle: dado el intervalo de integracion [a, b] realizamos una partici on
de [a, b] en n sub-intervalos de extremos a = x0 , x1 , x2 , . . . , xn1 , xn = b.
La distancia entre estos puntos es constante y vale x = ba n , entonces
xj = a + jx, j = 0, . . . , n. Si en cada uno de los sub-intervalos [xj , xj+1 ],
j = 0, . . . , n 1 elegimos un punto arbitrario j y evaluamos la funcion f en
j , podemos definir las sumas de Cauchy-Riemann como sigue:
n
X n
X
Sn = f (j )(xj xj1 ) = f (j )x.
j=1 j=1
141
5.10.1. Integraci
on por aproximaci
on rectangular
En el metodo de integracion numerica por rectangulos se sustituye f
con una funci on f constante por trozos. Los trozos son los intervalos
[xj , xj+1 ], j = 0, . . . , n 1, y la constante puede ser el valor de f calculado
en el extremo izquierdo de [xj , xj+1 ] o en su punto medio:
o bien
xj + xj+1
f2 (x) = f , x [xj , xj+1 ], j = 0, . . . , n 1.
2
La figura siguiente simboliza la aproximacion por rectangulos cuando se
escoge el extremo izquierdo o el punto medio, respectivamente.
n1
baX xj + xj+1
Sn (f2 ) = f .
n 2
j=0
142
Z b
f (x)dx Sn (f2 ) C(x)2 .
m
ax
a
ba
Dado que x = ny que, generalmente n b a, es preferible una
f
ormula que esta acotada por potencias mas altas de x, por eso la aproxi-
macion rectangular de punto medio es mas precisa (si se considera la misma
partici
on) que la aproximacion rectangular de extremo izquierdo.
5.10.2. Integraci
on por aproximaci
on trapezoidal
En el metodo de integracion numerica por trapecios se sustituye f con
una funci on f lineal por trozos. El segmento rectilneo en cada intervalo
[xj , xj+1 ], j = 0, . . . , n1, puede ser el segmento de recta tangente al grafico
de f en el punto medio de [xj , xj+1 ] (figura debajo a la izquierda), el segmen-
to que conecta (xj , f (xj )) con (xj+1 , f (xj+1 )) (figura debajo a la derecha),
o bien el segmento que pasa por (xj , f (xj )) con una pendiente asignada, por
ejemplo el valor de la tangente en el punto medio. En todos estos casos ya
no se aproxima el area por debajo del grafico de f con rectangulos, sino con
trapecios.
143
presiones5 de la funci on aproximada f:
0 xj + xj+1
f1 (x) = f (xj ) + f (x xj ), x [xj , xj+1 ], j = 0, . . . , n 1,
2
o bien,
f (xj+1 ) f (xj )
f2 (x) = f (xj ) + (x xj ), x [xj , xj+1 ], j = 0, . . . , n 1.
xj+1 xj
144
resultado esperado por la geometra del trapecio aproximada en este caso.
Las sumas de Cauchy-Riemann para la funcion aproximada f2 son:
n1
b a X f (xj ) + f (xj+1 )
Sn (f2 ) =
n 2
j=0
ba
= [f (a) + 2f (x1 ) + 2f (x2 ) + . . . + 2f (xn1 ) + f (b)]
2n
b a f (a) + f (b)
= + f (x1 ) + f (x2 ) + . . . + f (xn1 )
n 2
n1
b a f (a) + f (b) X
= + f (xj ) .
n 2
j=1
5.10.3. Integraci
on por aproximaci
on parab
olica (f
ormula de
Cavalieri-Simpson)
En el metodo de integracion numerica de Cavalieri-Simpson se sustituye
f con una funcion f parab olica por trozos, como en la figura siguiente.
Como se ha visto antes, pasar de una funcion constante a una funcion
lineal ha permitido una mejora sustancial de la aproximacion. Considerar
una funcion cuadratica aumenta a un mas la precision de la aproximacion.
Una parabola necesita 3 puntos para ser unvocamente determinada, por
tanto hay que refinar la particion anterior cortando los n intervalos [xj , xj+1 ],
145
x +x
por ejemplo, en el punto medio j = j 2 j+1 , y obteniendo 2n intervalos de
longitud = x ba
2 = 2n cada uno, entonces xj+1 j = y xj j = .
La funci
on parabolica aproximada f tiene ecuacion:
146
area debajo de f en el intervalo [xj , xj+1 ] es
El
Z xj+1 Z xj+1
f (x)dx = [(x j )2 + (x j ) + ]dx
xj xj
xj+1
(x j )3 (x j )2
= + + x
3 2 xj
(xj+1 j )3 (xj+1 j )2
= + + xj+1
3 2
(xj j )3 (xj j )2
xj
3 2
3 2 ()3 ()2
= + + xj+1 xj
3 2 3 2
3
= 2 + (xj+1 xj )
3
2 3
= + 2.
3
R xj+1
Para explicitar xj f(x)dx necesitamos el valor de , dado que ya conoce-
mos = f (j ) y = x ba
2 = 2n .
Podemos determinar a partir de las condiciones que aparecen en el
sistema anterior, (1)-(2) y (3)-(2) devuelven
f (xj ) f (j ) = 2
f (xj+1 ) f (j ) = 2 + ,
sumando estas dos ultimas ecuaciones y despejando con respecto a obten-
emos:
f (xj ) + f (xj+1 2f (j ))
= .
2 2
R xj+1
on de en xj f(x)dx llegamos a escribir:
Introduciendo esta expresi
Z xj+1
b a f (xj ) + f (xj+1 ) + 4f (j )
f(x)dx =
xj n 6
y sumando la contribuci on de todos los intervalos de la particion para j =
ormula de Cavalieri-Simpson6
0, . . . , n obtenemos la f
6
La f
ormula de Cavalieri-Simpson fue usada por primera vez por el fsico y astr
onomo
alem
an Johannes Kepler (1571-1630) y llamada Keplersche Fassregel, o sea regla de
Kepler. El matem atico italiano Bonaventura Cavalieri (1598-1647) fue el primero en for-
malizarla en un contexto matem atico y el matem
atico brit
anico Thomas Simpson (1710-
1761) la populariz
o en la comunidad angl ofona.
147
n1 n
" #
b a X X
Sn (f) = f (a) + f (b) + 2 f (xk ) + 4 f (i ) .
6n
k=1 i=0
148
Captulo 6
En este u
ltimo captulo volveremos a examinar las propiedades de los
n
umeros enteros y las relacionaremos con unos de los conceptos que hemos
visto hasta ahora, en particular lo de lmite de una funcion y de su desar-
rollo de Taylor. Comenzamos con una seccion dedicada a una tecnica de
demostraccion peculiar de los n
umeros naturales.
para que todas las piezas del domino se caigan es necesario y suficiente
que se cumplan estas dos condiciones:
149
1. Que se caiga la primera pieza;
2. Que cualquier pieza este posicionada de tal forma que, si se cae, provo-
ca la caida de la pieza siguiente.
El principio de inducci
on extiende esta idea al caso de un n
umero infinito
de piezas y a propiedades matematicas mas generales que la caida de una
pieza de domin o.
T
ecnica de demostracci
on por inducci
on
Supongamos que queremos demostrar una propiedad P para todos los
umeros enteros n n0 , con n0 N;
n
1. Comenzamos demostrando que P vale cuando n = n0 ;
2. Supongamos que la propiedad valga para un valor generico de n > n0
(hip on) y demostramos que eso implca que P vale
otesis de inducci
para n + 1.
Como por el domino, eso garantiza que la validez de P, comenzando por n0 ,
va propag
andose de un entero al siguiente hasta el infinto.
Es imprescindible presentar un ejemplo (se veran mas en los seminarios).
150
Para entender la importancia de la tecnica de induccion, citamos el hecho
de que Peano logr
o demostrar formalmente la validez de la formula de Taylor
con una prueba muy simple y elegante basada en la induccion.
Ejemplos:
n 0: (an ) = (n2 ) = (0, 1, 4, 9, 16, 25, . . .);
n 0: (an ) = ((1)n ) = (1, 1, 1, 1, 1, . . .);
1
n 1: (an ) = (2 n ) = (2, 2, 3 2, 4 2, . . .).
Def.: Se dice que una sucesion (an ) tiene una propiedad P definitiva-
mente si existe un n0 0 tal que an satisface P para cada n n0 .
Ejemplos:
151
on (n 10 n) es definitivamente positiva porque lo es para
La sucesi
cada n 100 (dejamos como ejercicio averiguar que es as);
on ( n12 ) es definitivamente menor de
La sucesi 1
10 porque lo es para cada
n 4.
Estamos listos para dar la definicion de lmite de una sucesion.
lm an = `
n+
su gr
afico (puntos discretos!), a partir de un cierto n0 , queda contenido en
una banda horizontal de ancho 2 centrada en `. Esto sucedera no importa
que tan peque na sea . Cuanto mas peque na, mas habra que esperar hasta
que la sucesi
on entre en la banda, pero tarde o temprano entrara.
lm an = + (resp. lm an = )
n+ n+
152
Teorema de Heine: Sean f : (a, b) R R y x0 (a, b). f converge
a ` R para x x0 si y solo si para cada sucesion (an ) tal que an
(a, b) definitivamente, an 6= x0 definitivamente y (an ) x0 , la sucesion
n+
numerica (f (an )) converge a ` para n +.
loga n (loga n) n n n bn n! nn , n +
153
a > 1;
> 1;
(0, 1);
> 1;
b > 1;
> 0.
xn
En particular, observamos que, fijado cada x R, n! 0 para n +.
Utilizaremos esta observaci
on en breve.
X = N;
: D N F
n 7 (n) fn .
Ejemplos:
n 0: fn : R \ {2} R,
2xn 2x2
(fn (x) = 2+n
) = (f0 (x) = 2 , f1 (x) = 2x
, f (x) = , . . .);
x+2 x+2 3
x+2 2 4
x+2
n 1: fn : (0, 1) (1, +) R,
3x1 3x1 3x1 3x1
(fn (x) = log(x n ) ) = (f1 (x) = log(x) , f2 (x) = , f (x)
log(x2 ) 3
= log(x3 )
, . . .).
154
Def.: Se dice que la sucesi
on de funciones (fn ) converge puntualmente
en x si la sucesi
on numerica (fn (x)) es convergente. Si (fn ) converge pun-
tualmente en cada punto de un intervalo E R, resulta bien definida la
on f : E R dada por
funci
f (x) = lm fn (x)
n+
y llamada funci
on lmite puntual de la sucesion de funciones (fn ).
155
6.3.1. El teorema de aproximaci
on de Stone-Weierstrass
El teorema de aproximacion de Stone-Weierstrass2 es uno de los mas
importantes de toda la matematica. Fue obtenido por primera vez en el
1885 por Weierstrass y sucesivamente generalizado por Stone en el 1937.
El teorema dice que cada funcion continua en un intervalo compacto puede
ser aproximada con un nivel de precision cualquiera a traves de una funcion
polinomial. Dado que los polinomios son las funciones no lineales mas simples
de la matem atica y dado que los ordenadores pueden computarlos de forma
muy rapida, resulta evidente la imporancia tanto teorica como practica de
este resultado. Rigurosamente:
2
Marshall Stone, 1903-1989, matem
atico americano.
156
6.4. Series n
umericas
La teora de las series numericas permite expresar de forma matematica-
mente rigurosa el concepto de suma finita de un n umero infinito de terminos,
un concepto que en la antig uedad pareca no tener sentido, como se deduce
de la famosa paradoja de Zen on (-489 -431) de Aquiles y la tortuga.
El filosofo Zenon, imagino Aquiles y una tortuga compitiendo en una carrera.
Como Aquiles es m as rapido, los jueces de la carrera deciden que Aquiles
comienze en la posici on x0 y la tortuga en la posicion x1 > x0 . Obviamente
en un tiempo t0 > 0 Aquiles alcanzara la posicion de la tortuga y despues la
superarar a. Sin embargo Zenon consideraba esa una paradoja logica, veamos
por que. En el instante de tiempo t1 = t20 , obviamente Aquiles estara detras
de la tortuga, y tambien en el instante t2 = t21 = t40 , etc. siguiendo has-
ta el infinito. La paradoja, seg un Zenon, es que Aquiles ha alcanzado la
tortuga en un instante de tiempo finito obtenido por la suma de una can-
tidad infinita de instantes de tiempo mas peque nos pero siempre positivos.
Zenon, considerando sus elucubraciones filosoficas mas logicas que la reali-
dad de los hechos, crey o haber demostrado logicamente la imposibilidad del
movimiento!
Vamos a ver como se puede dar rigor al concepto de suma de infinitos
terminos. Consideramos una sucesion de n umeros reales (an ) y construimos
una sucesi on auxiliar (sn ) definida como sigue:
s0 = a0
s1 = a0 + a1
s2 = a0 + a1 + a2
..
.
157
n
X
sn = a0 + a1 + a2 + . . . + an = ak ,
k=0
Def.: Se dice serie numerica con termino general (an ) la sucesion nu-
merica de las sumas parciales sn . Si (sn ) es convergente, se dice que la serie
es convergente, si (sn ) diverge o es irregular (no tiene lmite), se dice que
la serie diverge o es irregular, respectivamente. Si (sn ) coverge a S para
n + llamamos S la suma de la serie y escribimos:
X n
X
an = lm ak = S .
n+
n=0 k=0
P
Si la sucesi
on (an ) comienza en n0 > 0 entonces escribiremos n=n0 an .
158
P 2 2
por: n=1 n , a pesare de que
n n+ 0, tenemos que
X 2
= 2 + 1 + 2/3 + 2/4 + 2/5 + 2/6 + 2/7 + 2/8 + . . .
n
n=1
> 2 + 1 + 1/2 + 1/2 + 1/4 + 1/4 + 1/4 + 1/4 + . . .
= 2 + 1 + 1 + 1 + ... = +.
Como
P 1consecuencia inmediata tenemos que tambien la serie harmonica
n=1 n es divergente.
Terminamos con un ejemplo practico para entender que una sumaPinfinita
puede, efectivamente, dar un resultado finito es el siguiente: la serie 1
n=1 2n
modeliza coger la mitad de un pastel, despues la mitad de lo que queda (un
cuarto), despues la mitad de lo que queda, etc. Claramente, en cada paso
siempre quedar a algo y en el lmite se coge todo el pastel (entonces la suma
infinita da 1).
s0 (x) = f0 (x), x D
Def.: Se dice serie de funciones con termino general (fn (x)) la sucesion
de funciones dada por las funciones sumas parciales sn (x). Si (sn (x)) es una
on de funciones que converge puntualmente en x0 D, diremos que
sucesi
la serie de funciones converge puntualmente en x0 D, si (sn (x)) es una
159
on de funciones divergente o irregular en x0 D, diremos que la serie
sucesi
de funciones diverge o es irregular en x0 D, respectivamente.
160
Todas las funciones elementales del calculo tienen esta propiedad.
Consideramos ahora un intervalo real abierto I tal que x0 I. En gen-
eral, cuando el error de la aproximacion de Taylor de una funcion f infini-
tamente derivable tiende a cero para cada x I, se dice que f se puede
desarrollar en serie de Taylor en I con centro en x0 y vale la formula:
X f (n) (x0 )
f (x) = (x x0 )n x, x0 I,
n!
n=0
en particular, si x0 = 0, tenemos
X f (n) (0)
f (x) = xn x I.
n!
n=0
Se puede observar que las series de Taylor son series de funciones partic-
ulares en las cuales las funciones que aparecen son potencias de la variable
x multiplicadas por coeficientes proporcionales a las derivadas de la funcion
que se est
a desarrollando. En general, una expresion del tipo
X
an (x x0 )n
n=0
161
Funci
on Serie de potencias correspondiente Radio de convergencia
P xn
ex n=0 R = +
n! n
P n+1 x
log(1 + x) n=1 (1) R = +
n
P n x2n+1
sin x n=0 (1) R = +
(2n + 1)!
2n
n x
P
cos x n=0 (1) R = +
(2n)!
P x2n+1
sinh x n=0 R = +
(2n + 1)!
P x2n
cosh x n=0 R = +
P (2n)! n
(1 + x)
n=0 n x R=1
1 P n
n=0 x R=1
1x
( 1) ( n + 1)
con R y n = es el coeficiente binomial.
n!
Observese, en particular, que si ponemos x = 1 en la serie exponencial
obtenemos
X 1
e=
n!
n=0
162
6.6. La formula de la exponencial compleja de Eu-
ler
Presentamos en esta u ltima seccion una formula de grande importancia
debida a Leonhard Euler (1707-1783) matematico y fsico suizo. Euler es
seguramente el matem atico mas importante del siglo XVIII y muchos lo
consideran el matem atico mas importante de la historia junto con Gauss.
Es el matem atico que m as resultados ha publicado y a el debemos la gran
mayora de la notacion matematica que utilizamos hoy en da.
Euler estudi
o las series con terminos complejos y encontro unas formulas
particularmente significativas para la serie exponencial en el campo comple-
jo. Para entender su
analisis, tenemos que dar sentido al concepto de serie de
terminos complejos. Dado que una serie no es nada mas que la sucesion de
las sumas parciales de una sucesion, el problema de definir lo que es una serie
convergente en el campo complejo se reduce al problema de la definicion de
una sucesion convergente de n umeros complejos.
163
|sn `| 0 para n .
Antes que nada observamos que esta serie de funciones converge puntual-
mente en Ptodo C ndado que, fijando cada z, la serie de terminos reales no
|z |
negativos n=0 n! es convergente en virtud de la f
ormula de De Moivre -
Stirling.
Analizamos, en particular, el caso de un numero complejo imaginario
puro, o sea z = ix, x R:
X (ix)n x2 x3 x4 x5 x6 x7
eix = = 1 + ix i + +i i + ...
n! 2 3! 4! 5! 6! 7!
n=0
X x2n X x2n+1
= (1)n +i (1)n = cos x + i sin x.
(2n)! (2n + 1)!
n=0 n=0
164
A traves de las f
ormulas de Euler podemos escribir la forma polar de los
n
umeros complejos en la forma exponencial:
z = (cos + i sin ) = ei .
ei + 1 = 0 .
Se considera la f
ormula m
as bella porque mezcla los cinco n umeros y las
cuatro operaciones m
as importantes de la matematica: 0,1,e,, i e igualdad,
suma, producto y exponenciacion, respectivamente.
THE END.
165