Método del punto fijo

Un punto fijo de una función g, es un numero p, tal que g (p) = p, el problema se centra en
encontrar soluciones de una ecuación f(x)=0 y el de encontrar los puntos fijos de una
función h(x) son equivalentes en el siguiente sentido: dado el problema de hallar las
soluciones de una ecuación f(x)=0, podemos definir una función g con un punto fijo p de
muchas formas; por ejemplo, f(x)=x - g(x). En forma inversa, si la función g tiene un punto
fijo en p, entonces la función definida por f(x)=x - g(x) posee un cero en p.

La forma general para el método del punto fijo es x= g(x), a partir de la función se pueden
se pueden conformar diferentes g(x). Por ejemplo:
x−cosx =0

x=cosx x=cos x
−1
x=2 x−cosx
x=√ x + x−cosx
2

En el método del punto fijo resulta importante el teorema del punto dijo, el cual va a ayudar
a saber cuál es el mejor g(x) en una ecuación.
El método del punto fijo dice lo siguiente: si g= [ a , b ] → [ a ,b ] es continua y derivable en
[ a , b ] , con |g ´ ( x )|≤ k <1 ∀ x ∈ [ a , b ] y dado un xo ∈ [ a , b ] , sucesión
Xn=g (Xn−1) , es convergente, tal que lim Xn=a , entonces a es una raíz de la
n→π
ecuación x= g(x)
Ejercicio método del punto fijo:

f ( x )=x 3 +4 x−10 , teniendo en cuenta que x=0

 3 2
x=x + 4 x + x−10
 x=
√ 10
x +4
 10
 x= −4
x4

i g1(x) g2(x) g3(x)

0 3,875 1,3484 0,4444
1 111,12 1,3674 46,8351
2 1459854 1,3649 -3,59

Interpolación
Teniendo en cuenta que los polinomios de Taylor no son adecuados para la interpolación, es
necesario hacer uso de métodos alternos.

Existen tres métodos para interpolación: Lagrange, Newton y Hermite

Lagrange
La fórmula general para este método es la siguiente: L ( x )=L 0 y 0+ L1 y 1 y las
x−x 1 x−x 0
formulas particulares para L0 y L1 son las siguientes: L0= L1=
x 0−x 1 x 1−x 0
n−1 n−1 n−1
x−xj
Formula general L ( x )=∑ ligi=∑ yi ∏ ; i≠ j
i =0 i=0 j=0 xj−xj

Ejemplo:
x y
P0 -4 1
P1 0 2

x−0 −x x+ 4 x +4
L0= = , L 1= =
−4−0 4 0+ 4 4
Por tanto:

L ( x )=L 0+ L 1 y 1= ( −x4 ) ( 1) +( x +44 ) ( 2)

−x x+ 4 −x +2 x +8 x +8 x
¿ + = = = +2
4 2 4 4 4
Ejemplo:
x y
P1 1 6
P0 0 9

x−x 1 x−x 0
L0= L1=
x 0−x 1 x 1−x 0

Por lo tanto,
x−1 x−0
L0= =−x +1 , L1= =x
0−1 1−0

( (−x+1 ) )∗9+ x ( 6 )=−9 x +9+6 x=9−3 x

Newton (Diferencias divididas)

La fórmula para este método es el siguiente: N ( x ) =a0 + a1 ( x−x 0)

Ejercicio: Hallar a0 teniendo en cuenta que a 4 =3 y p0 x = -1

x y
p0 -1 ?
p1 0 0
p2 1 -2
p3 2 38
p4 4 748
 86−a 0
3=
30
a0 =86−90

a0 =−4

Hermite
Para el método de Hermite está la siguiente fórmula:

H ( x ) =Q0,0 +Q1,1 ( x−x 0 ) +Q2,2 (x−x 0)2 +Q3,3 ( x−x 0 )2 ( x−x 1 ) +Q4,4 ( x−x 1 )2 (x−x 1)2 … … … … .

Ejercicio: Hallar a0 , a1 , a2 , a3 , a4 teniendo en cuenta la siguiente tabla

x Tiempo 0 3 5 8 13

y Distancia 0 22 38 62 99
5 6 3 3
y’ Velocidad 75 77 80 74 72
Integración numérica
La necesidad de aproximar numéricamente el valor de una integral surge por dos motivos
fundamentalmente: la dificultad o imposibilidad en el cálculo de una primitiva y la función
a integrar solo se conoce por una tabla de valores.
Existen dos métodos de integración numérica los cuales son Trapecio y Simpson, estos de
desarrollan de manera simple y múltiple. Primeramente se desarrollaran en la metodología
simple.

Trapecio simple
Esta es la fórmula general para el método de trapecio simple:
f ( x 0+2 f ( x 1 ) +2 f ( x 2 ) + f ( x 3 ) )
h
¿
2
Nota: Este método se realiza para n par e impar
Dicha fórmula se puede demostrar cuando n=2 de la siguiente manera:
x2 x1 x2

∫ f ( x ) dx=∫ f ( x ) dx +∫ f ( x ) dx
x0 x0 x1

x1 x2

Parte 1. ∫ ( L0 y 0 + L1 y 1 ) dx Parte 2. ∫ ( L1 y 1 + L2 y 2 )dx

x0 x1
La parte 2, se realiza igual que la parte 1, por efectos prácticos se mostrará el desarrollo de
la parte 1.
x−x 0
dx=¿
x1− x0
x1 x1 x1
x−x
∫ ( L0 y 0 + L1 y 1 ) dx= y 0 ∫ x −x1 dx+ y1∫ ¿
x 0 x 0 1 0 x 0

[ ] [ ]
x1 x1 x1 x1
− y0 y 2

h
∫ xdx −∫ x1 dx + h1 ∫ x2 −∫ x 0 xdx
x 0 x 0 x0 x 0

[ ] [ ]
2 2 2 2
− y0 x 1 x 0 2 y1 x 1 x 0 2
¿ − −x 1−x 0 x 1 ¿ − −x 0 x 1+ x 0
h 2 2 h 2 2

[ ]
2 2
y0 x 1 x y1 2
¿
h 2
−x 0 x 1 + 0
2
¿ [ x −2 x 0 x 1 + x 20 ]
2h 1

y0 2 y1 y h
¿
2h
[ x 1−2 x 0 x 1 + x 20 ] ¿
2h
[ ( x 1−x 0 )2 ]= 1
2
y0
¿ [( x −x )2 ]
2h 1 0
y0 2 y0 h
¿ [ h ]=
2h 2

x
f (¿¿ 1)
x
f (¿¿ 0)
f ( x1 )+ ¿
h
f ( x 0)+¿ → ¿
2
h
∴ ¿
2
x 0 + 2 f ( x1 ) + f ( x 2 )
f¿
h
¿
2
Ejercicio
Resolver la siguiente integral para n = 1, n = 3, n = 5
2

∫ x2 dx
−1

Para n = 1
b−a 2+1
h= = =3
n 1

-1 2

x
x
f (¿¿ 1)
f (¿¿ o)+¿
3 15
¿= ( 1+ 4 )=
2 2
h
¿
2

Para n = 3
b−a 2+1
h= = =1
n 3

-1 0 1 2
 −1
¿
¿
0
¿
¿
1
¿
¿
¿
1
f ( x 0 +2 f ( x 1 ) +2 f ( x 2 ) +f ( x 3 ) ) = ¿
2
h
¿
2
1
¿ [ 1+0+ 2+ 4 ] =3,5
2

Para n = 5
b−a 2+1
h= = =0,6
n 5

-1 -0,4 0,2 0,8 1,4 2

−1
¿
¿
−0,4
¿
¿
0,2
¿
¿
0,8
¿
¿
1,4
¿
¿
2
¿
¿
¿
0,6
¿
2
 ¿ 0,3 ( 1+ 0,32+0,08+1,28+3,92+ 4 )=3,18

Simpson (1/3)
Para este método solo se puede trabar si el número de particiones es n = par. La fórmula
general es la siguiente:
h
3[ 0
I= f ( x ) + 4 ∑ f ( x i)+2 ∑ f ( xi ) + f (x n ) ]

Y se demuestra se la siguiente forma:

x1

∫ f ( x ) dx
x0

x1 x1

Entonces: ∫ f ( x ) dx ∫ f 2 ( x ) dx
x0 x0

(x−xm )( x−b) ( x−a)(x−b) (x−a)( x−xm)

Donde: f 2 ( x )=f (a ) + f ( xm ) + f (b)
(a−xm )(a−b) ( xm−a)(xm−b) ( b−a ) ( b−xm )
b−a
h= =xm −a=b−xm
2
(x−xm )( x−b) ( x−a)(x−b) (x−a)( x− xm)
f 2 ( x )=f (a ) + f ( xm ) + f ( b)
(−h)(−2 h) (h)(−h) ( 2h )( h )

f (a ) f ( xm ) f ( b)
f 2 ( x )= 2
( x−xm )( x−b )− 2 ( x−a )( x−b ) + 2 ( x−a)( x −xm)
2h h 2h

b b b b
f (a) f ( xm ) f (b)
∫ f 2 ( x ) dx= 2 h2 ∫ ( x−xm )( x−b ) dx− h2 ∫ ( x−a ) ( x−b ) dx + 2 h2 ∫ ( x−a )( x−xm ) dx
a a a a

b 2 3 2 3 2 3
( x−b ) ( x−b ) ( a−b ) ( a−b ) (−2h ) (−2 h ) 3 4
I 1 =∫ ( x−xm )( x−b ) dx=( x−xm ) − =a−xm + =−(−h ) =2 h −
a 2 6 2 6 2 6 3
b 2
( x −b ) ( x−b)3 (a−b)3 (−2 h)3 −4 3
I 2 =∫ ( x−a ) ( x−b ) dx =( x−a ) − = = = h
a 2 6 b 6 3
 b 2 3 2 3 2 33
( x−xm ) ( x−xm ) ( b−xm ) ( b−xm ) ( a−xm ) h h (−
I 3 =∫ ( x−a ) ( x−xm ) dx=( x−a ) − =( b−a ) − + =( 2 h ) − +
a 2 6 2 6 6 2 6
b
f ( a) f ( xm) −4 3 f ( b ) 2 3
∫ f 2 ( x ) dx= 2 h2 ( 23 h3)− h
2
3 (2h )
3
h
3
4
( )
3
h h
h + 2 h =f ( a ) + + f ( xm ) + f ( b ) = [ f ( a ) + 4 f ( xm ) + f (b) ]
3 3
a

∫ f ( x ) dx ≅(b−a)
a
[ f ( a )+ 4 f ( xm ) + f (b)
6 ]
b−a
h= , xn ∈ [ x i−1 , x 1 ]
n
b x1 x2 xn

∫ f ( x ) dx=∫ f ( x ) dx+∫ f ( x ) dx +… …+ ∫ f ( x ) dx
a x0 x1 xn−1

[ ] [ ]
b
f ( x 0 )+ 4 ( xm ) + f ( x 1 ) f + 4 f ( xmn) + f ( xn)
∫ f ( x ) dx ≅ ( x 1−x 0 ) 6
+… …+(x n−x n−1) (xn−1 )
6
a

[ ]
n−1 n−2

b f (x 0) + 4 ∑ f ( xi)+ 2 ∑ f ( xj ) + f ( xn)
∫ f ( x ) dx ≅(b−a) i=impar
3n
I= par

Ejercicio
Resolver la siguiente integral para n = 4
4,5

∫ esin ⁡(x) dx
−π

b−a 4,5−(−π )
h= = =1,91
n 4

-3,14 -1,23 0,68 2,59 4,5

1,91
I=
3
[ f (−3,14 ) + 4 [ f (−1,23 ) +f (2,59) ]+2 f ( 0,68 )+ f (4,5)]=8,556
Ejercicio
Resolver la siguiente integral para n = 6
7

∫ ( x 2 +1 ) dx
2

b−a 7−2 5
h= = = =0,83
n 6 6

2 2,83 3,67 4,5 5,33 6,17 7

f (¿6,17 )
f (¿ 4,5)+¿
f (¿5,33 )
f (¿3,67 )+ ¿
f (¿7 )
¿+¿
¿
f (¿2,83 )+ ¿+2¿
¿
f (¿ 2)+ 4 ¿
¿
0,83
I= ¿
3

Trapecio múltiple

Para el trapecio múltiple la fórmula es la siguiente: AT =( hx2 )( hy2 ) [ malla ]

Ejercicio: Resolver la siguiente integral múltiple por el método de trapecio n x =2 ,n y =4
1 3

∫∫ 2 xydydx
0 −2
 bx −a x 1−0
hx= = =0,5
nx 2

0 0,5 1

b y −a y 3+2
hy= = =1,25
ny 4

-2 -0,75 0,5 1,75 3

Malla

Malla= [ 1 f ( 0,−2 ) +2 f ( 0.5, 0.2 )+ 1 f (1,−2 ) +2 f ( 0,−0.75 )+ 4 f ( 0.5,−0.75 )+2 f ( 1,−0.75 ) +2 f ( 0, 0.5 ) +4 f ( 0.5,

AT = ( 0,52 )( 1,252 )[ 16 ]=2,5

Simpson (1/3) múltiple

Para Simpson múltiple la fórmula es la siguiente: AS= ( hx3 )( hy3 ) [malla ]

Ejercicio: Resolver la siguiente integral múltiple por el método de trapecio n x =2 ,n y =4
1 2
2 2

∫∫ e x + y dxdy
−1 0

bx −a x 2−0
hx= = =1
nx 2
 0 1 2

b y −a y 1+1
hy= = =0,5
ny 4

-1 -0,5 0 0,5 1

Malla

Malla= [ 1e 0 +1 + 4 e 1 +1 + 1e 2 +1 +4 e 0 +0,5 +16 e 1 +0,5 + 4 e 2 +0,5 +2 e 0 +0 +8 e 1 +0 +2 e2 +0 + 4 e0 +(−0,5) +16 e 1 +(−0,5) +

2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2

AT = ( 13 )( 0,53 )[ 1176,983 ]=65,387