Capítulo 1

Interpolación.

Interpolar signi…ca estimar el valor desconocido de una función en un punto.

Interpolación lineal: Dados los puntos (x0 ; y0 ) y (x1 ; y1 ) ; el polinomio interpolador
es
y1 y0
P (x) = y0 + (x x0 )
x1 x0
Matemático Francés, Joseph Lagrange, descubrío una forma ligeramente distinto
y1 y0
P1 (x) = y0 + (x x0 ) (x x0 )
x1 x0 x1 x0
x x1 x x0
= y0 + y1
x0 x1 x1 x0

Denotando
x x1 x x0
L1;0 (x) = ; L1;1 (x) =
x0 x1 x1 x0
Escribimos,

P1 (x) = y0 L1;0 (x) + y1 L1;1 (x)

X1
= yk L1;k (x)
k=0

Los L1;k se llaman polinomios coe…cientes de Lagrange. y P1 (x) polinomio interpo-

lador de Lagrange.

Teorema 1 Si x0 ; x1 ; ; xn son n + 1 números distintos y si f es una función cuyos

valores están dados en xi ; entonces existe un único polinomio de grado n; Pn (x) ; tal
que f (x) = Pn (xk ) ; k = 0; 1; ; n:

Este polinomio es dada por

Pn (x) = f (x0 ) Ln;0 (x) + f (x1 ) Ln;1 (x) + + f (xk ) Ln;n (x)
Xn
= f (xk ) Ln;k (x) ;
k=0

1
2 CAPÍTULO 1. INTERPOLACIÓN.

donde
(x x0 ) (x x1 ) (x xk 1 ) (x xk+1 ) (x xn )
Ln;k (x) =
(xk x0 ) (xk x1 ) (xk xk 1 ) (xk xk+1 ) (xk xn )
Yn
(x xj )
= :
j=0
xk xj
j6=k

Nótese que
0 j 6= k
Ln;k (xj ) =
1 j=k

Observación 1 En adelante denotaremos Ln;k (x) = Lk (x) :

Ejemplo 1 Sea f (x) = x1 ; x 2 [2; 4] : Use tres nodos x0 = 2; x1 = 2;5 y x2 = 4; para

construir el polinomio interpolador cuadrático P2 (x) :

Teorema 2 (Fórmula del error) suponga que f 2 C n+1 [a; b] y que x0 ; x1 ; ; xn 2

[a; b] son n + 1 nodos de interpolación. Si x 2 [a; b] entonces

f (x) = Pn (x) + En (x)

donde Pn (x) es el polinomio inrterpolador de Lagrange y el error En (x) es

(x x0 ) (x x1 ) (x xn ) f (n+1) ( (x))
En (x) =
(n + 1)!
para algún (x) en ha; bi :

Cabe señalar que los polinomios de Lagrange se emplean para derivar diferenciación
numérica y los métodos de integración.

Teorema 3 (Cotas de error para la interpolación de Lagrange) Suponga que f

de…nida en [a; b], que contiene los nodos equiespaciados xk = x0 + hk: Además, suónga
que f y sus derivadas hasta de orden 4 son continuas y por tanto acotados en [x0 ; xn ] ;
es decir
f (n+1) (x) < Mn+1 ; x0 x xn ; n = 1; 2; 3
Entonces
h2 M2
jE1 (x)j ; x 2 [x0 ; x1 ]
8
h3 M3
jE2 (x)j p ; x 2 [x0 ; x2 ]
9 3
h4 M4
jE3 (x)j ; x 2 [x0 ; x3 ]
24
Observación 2 Al usar polinomio interpolador de Lagrange uno de los inconvenoientes
es que Pn 1 y Pn no se relacionan. Cada polinomio se construye individualmente y por
tanto el trabajo es arduo si el grado es alto.
Capítulo 2

Integración Numérica
Rb
En este capítulo pretendemos determinar a f (x) dx:
El método básico con que se aproxima esta integral recibe el nombre de cuadratura
numérica y emplea una suma del tipo
X
n
ai f (xi )
i=0
Rb
para aproximar a f (x) dx:
El método de la cuadratura que presentamos se basa en los polinomios interpolantes
de Lagrange.
Primero seleccionamos los nodos fx0 ; x1 ; ; xn g de [a; b] y luego el polinomio inter-
polante
X n
Pn (x) = f (xi) Li (x)
i=0
y su término de error de truncamiento en [a; b]. Entonces,
Z b Z b Z b
f (x) dx = Pn (x) dx + En (x) dx
a a a
Z bX
n Z bY
n
f (n+1) ( (x))
= f (xi) Li (x) dx + (x xi ) dx
a i=0 a i=0 (n + 1)!

donde ( (x)) se encuentra en [a; b] para cada x:

Z b Xn Z bYn
1
f (x) dx = ai f (xi) + (x xi ) f (n+1) ( (x)) dx
a i=0
(n + 1)! a i=0

donde Z b
ai = Li (x) dx
a
Despreciando el error, tenemos
Z b X
n
f (x) dx ai f (xi)
a i=0

y se denomina fórmulas de Newton Cotés.

3
4 CAPÍTULO 2. INTEGRACIÓN NUMÉRICA

Teorema 4 (Fórmulas de cuadratura cerradas de Nuewton Cotes) Supongamos

que xk = x0 + kh, son nodos equiespaciados y sea fk = f (xk ) para k = 0; 1; 2; : : : ; m:
Tenemos los siguientes:
1. La regla del Trapecio
Z x1
h h3 (2)
f (x) dx = (f0 + f1 ) f ( )
x0 2 12
2. La regla de Simpson
Z x2
h h5 (4)
f (x) dx = (f0 + 4f1 + f2 ) f ( )
x0 3 90

3. La regla 3/8 de Simpson

Z x3
3h 3h5 (4)
f (x) dx = (f0 + 3f1 + 3f2 + f3 ) f ( )
x0 8 80
4. La regla de Boole
Z x4
2h 8h7 (6)
f (x) dx = (7f0 + 32f1 + 12f2 + 32f3 + 7f4 ) f ( )
x0 45 945
Cabe señalar que la regla del trapecio es exacta para P1 y la de Simpson es exacta
para P2 y P3 :
Ejemplo 2 Aproxime las siguientes integrales aplicando: a) la regla del trapecio b) La
regla de RSimpson. R1
1 2
a) 0;5 x4 dx c) 0 e x dx
R1
b) 0
x2 e x dx

Ejemplo 3 Usando las fórmulas del error, estime el error en cada caso del ejemplo
anterior.

2.1. La regla compuesta del trapecio y de simpson

2.1.1. Trapecio compuesto
Si en [a; b] se hace una partición como la siguiente
a = x0 < x 1 < x 2 < < xn 1 < xn = b
entonces se puede aplicar la regla del trapecio a cada uno de los subintervalos. Nótese
que los nodos no se encuentran necesariamente uniformemente espaciados
Z b X n Z xi
f (x) dx = f (x) dx
a i=1 xi 1

1X
n
(xi xi 1 ) [fi 1 + fi ]
2 i=1
2.1. LA REGLA COMPUESTA DEL TRAPECIO Y DE SIMPSON 5

Teorema 5 Sea f 2 C 2 [a; b] , h = b na y si xk = a+kh; para cada k = 0; 1; ; n; existe

2 (a; b) tal que la regla compuesta del trapecio para n subintervalos puede escribirse
con su término de error como
Z b
hX
n
b a 2 00
f (x) dx = [fk 1 + fk ] hf ( )
a 2 i=1 12
h
= [f0 + 2f1 + 2f2 + + 2fn 1 + fn ] + Et (f; h)
2
h X
n 1
= [f0 + fn ] + h fk + Et (f; h)
2 i=1

2.1.2. Simpson Compuesto

b a
Si en [a; b] se hace una partición regular en la que h = con n par, se tiene que
n
xi = a + ih; i = 0; 2; ; n: Entonces
Z b Z x2 Z x4 Z xn
f (x) dx = f (x) dx + f (x) dx + : : : + f (x) dx
a x0 x2 xn 2

n=2 Z
X x2i
= f (x) dx
i=1 x2i 2

n=2
hX
[f2i 2 + 4f2i 1 + f2i ]
3 i=1

b a
Teorema 6 Sea f 2 C 4 [a; b] , n = 2m; h = y si xk = a + kh; para cada
2m
k = 0; 1; ; 2m: Existe 2 (a; b) tal que la regla compuesta de Simpson para 2m
subintervalos puede escribirse con su término de error como
Z b
hX
n
b a 4 (4)
f (x) dx = [f2k 2 + 4f2k 1 + f2k ] hf ( )
a 3 k=1 180
h
= [f0 + 4f1 + 2f2 + 4f3 + + 2f2m 2 + 4f2m 1 + f2m ] + Es (f; h)
3
2h X 4h X
m 1 m
h
[f0 + f2m ] + f2k + f2k 1 + Es (f; h)
3 3 k=1 3 k=1

Ejemplo 4 Aplique la regla compuesta del trapecio con los valores indicados de n para
aproximar
R 4 las siguientes integrales:R
2
a) 0 x e x dx; n = 4
2
c) 2 x21+4 dx; n = 6
R2 R5
b) x ln xdx n=4 d) p 1 dx n=8
1 3 x2 +4

Ejemplo 5 Aplique la regla compuesta de Simpson con los valores indicados de n para
aproximar las ntegrales del ejemplo 1.
6 CAPÍTULO 2. INTEGRACIÓN NUMÉRICA

Ejemplo 6 Determine, en cada uno de los siguientes casos, el número n y el tamaño de

los subintervalos h de manera que la regla del trapecio con n subintervalos nos permita
obtener la integral dada con una precisión de 5x10 9 :
R7 R2 R =6 R3 1
a) 2 x1 dx b) 0 xe x dx c) =6
cos xdx d) 2 5 x
dx

Ejemplo 7 Realice el ejemplo anterior con la regla de simpson.

Capítulo 3

Funciones ortogonales y series de

fourier

3.1. Preliminares
Tenemos las siguientes de…niciones:
El producto interno de dos funciones f1 y f2 en un intervalo [a; b] es el número
Z b
(f1 ; f2 ) = f1 (x) f2 (x) dx:
a

Dos funciones f1 y f2 son ortogonales en un intervalo [a; b] si (f1 ; f2 ) = 0:

Un conjunto de funciones de valor real f 0 (x) ; 1 (x) ; 2 (x) ; g es ortogonal en
un intervalo [a; b] si
Z b
( m; n) = m (x) n (x) dx = 0; m 6= n:
a

La norma,o longitud generalizada, de una función n es

1=2
k n (x)k = ( n; n)
Z b 1=2
2
= n (x) dx
a

Si f n (x)g es un conjunto ortogonal de funciones en un intervalo [a; b] y k n (x)k = 1

para n = 1; 2; se dice que f n (x)g es un conjunto ortonormal en el intervalo.dado.

Ejemplo 8 El conjunto f1; cos x; cos 2x; g es ortogonal en [ ; ] :

n o
p x ; cos
Ejemplo 9 El conjunto p12 ; cos p 2x es ortonormal en [ ; ] :

Ejemplo 10 Demuestre que el conjunto dado de funciones es ortogonal en el intervalo

indicado. Calcule la norma de cada función del conjunto.
a) fsenx; sen3x; sen5x; g es ortogonal en [0; =2] :
b) fcos x; cos 3x; cos 5x; g es ortogonal en [0; =2] :

7
8 CAPÍTULO 3. FUNCIONES ORTOGONALES Y SERIES DE FOURIER

c) fsennxg
n ;on = 1; 2; 3 [0; ] :
d) sen np x ; n = 1; 2; 3 [0; p] :
n o
e) 1; cos np x ; n = 1; 2; 3 [0; p] :
n o
f) 1; cos np x; sen mp x s; n = 1; 2; 3 ; n = 1; 2; 3 ; [ p; p] :
Capítulo 4

Serie de Fourier

La serie de Fourier de una función f de…nida en el intervalo ( p; p) es de…nida por

a0 X
1
n x n x
f (x) = + an cos + bn sen (1)
2 n=1
p p

donde
Z p
1
a0 = f (x) dx
p p
Z p
1 n
an = f (x) cos xdx
p p p
Z p
1 n
bn = f (x) sen xdx
p p p
Ejemplo 11 Desarrolle f (x) en una serie de Fourier.
0 ; <x<0
f (x) =
x ; 0 x<
Solución.
Tenemos:
a0 = 2 ;
1 ( 1)n
an = cosn2n +1
= n2
bn = n1
Por tanto,
X
1
1 ( 1)n 1
f (x) = + cos nx + sennx
4 n=1
n2 n

Teorema 7 (Condicones de convergencia) Sean f y f 0 continuas en tramos en el

intervalo ( p; p) ; esto es excepto en un número …nito de puntos en el intervalo y con
discontinuidad solo …nitas en esos puntos. Entonces, la serie de fourier de f en el in-
tervalo converge hacia f (x) en un punto de continuidad. En un punto de discontinuidad
la serie de fourier converge hacia el promedio, es decir,
f (x+ ) + f (x ) 1h i
= l m f (x + h) + l m f (x h)
2 2 h!0 h!0

9
10 CAPÍTULO 4. SERIE DE FOURIER

Así por ejemplo, en el ejemplo anterior la función es discontinua en x = 0 y por tanto

la serie converge a
f (0+ ) + f (0 ) +0
= = :
2 2 2

Observación 3 . n o
1. El conjunto de funciones 1; cos p x; cos 2p x; ; sen p x; sen 2p x; tiene peri-
odo común 2p; por consiguiente la ecuación (1) es periódico.
2. La serie de fourier da la extensión periódica de f fuera del intervalo.

f (x + 2p) = f (x)

3. Cuando f es continua por tramos y existen f 0 y f+0 en x = p y x = p; respecti-

vamente, la serie (1) converge al valor promedio

f (p ) + f ( p+ )
:
2
en esos extremos y hacia valores extendidos periódicamente en x = 3p; 5p; 7p; etc

Asi por ejemplo, la serie de fourier del ejemplo (1) converge hacia la extensión per-
iódica de f (x) dado en todo el eje x: En x = 0; 2p; 4p; 6p; la serie converge
a
f (0+ ) + f (0 ) f (2 + ) + f (2 )
= = = :
2 2 2
y en x = ; 3p; 5 ; la serie converge a
+
f( )+f( )
=0
2

Ejemplo 12 Determine la serie de fourier de cada f en el intervalo dado.

a) c)
0 ; <x<0 0 ; =2 < x < 0
f (x) = f (x) =
ex 1 ; 0 x < cos x ; 0 x < =2

b) 8
< 0 ; 2<x<0
f (x) = x ; 0 x<1
:
1 ; 1 x<2

4.0.1. Serie de fourier de cosenos y senos

Tenemos las de…niciones siguientes:
1. La serie de fourier de una función par en el intervalo ( p; p) es la serie de cosenos

a0 X
1
n
f (x) = + an cos x
2 n=1
p
 11

donde
Z
2 p
a0 = f (x) dx
p 0
Z
2 p n
an = f (x) cos xdx
p 0 p
2. La serie de fourier de una función impar en el intervalo ( p; p) es la serie de senos
X
1
n
f (x) = bn sen x
n=1
p
Z p
2 n
bn = f (x) sen xdx
p 0 p
Ejemplo 13 Desarrolle en forma de una serie de fourier las funciones:
a) f (x) = x2 ; 1<x<1 b) f (x) = x; <x<

Desarrollo en mitad de intervalo

Si f está de…nida en intervalo 0 < x < L; como representamos por una serie
trigonométrica? Esto podemos hacer de diferentes formas sustituyendo la de…nición
para f (x) en L < x < 0: Por brevedad consideremos tres casos mas importantanrtes:
i) Re‡ejando la grá…ca de f con respecto al eje Y en ( L; 0) : De esta manera la
función f es par en ( L; L) :
ii) Re‡ejando la grá…ca de f con respecto al origen en ( L; 0) : De esta manera la
función f es impar( L; L) :
iii) De…na f en ( L; 0) por y = f (x + L)
Ejemplo 14 Desarrolle f (x) = x2 ; 0 < x < L: a) en una serie cosenos b) en una serie
desenos c) en una serie de fourier.
Solución.
a) Tenemos
2 2 4L2 ( 1)n
a0 = L; an =
3 n2 2
f (x) =
b)
Ejemplo 15 Desarrolle en una serie de: cosenos y senos, en medio intervalo de las
siguietes funciones.
a) c)
x ; 0 < x < =2 x ; 0<x<1
f (x) = f (x) =
x ; =2 x < 1 ; 1 x<2
Solución.
Ejemplo 16 Desarrolle la función dada como serie de fourier
a) f (x) = x + 1; 0 < x < 1
b) f (x) = 2 x; 0 < x < 2
c) f (x) = x (2 x) ; 0 < x < 2