Métodos Numéricos

Blanca Rivera Muñoz – Ingeniería en Logística y Transporte

FORMULARIO

Series de Taylor:
%
𝑓 (") (𝑥$ )
𝑓(𝑥) = ' (𝑥 − 𝑥$ )"
𝑘!
"&$

𝑓(𝑥) = 𝑎$ + 𝑎' (𝑥 − 𝑥$ ) + 𝑎( (𝑥 − 𝑥$ )( + 𝑎) (𝑥 − 𝑥$ )) + ⋯

𝑓 * (𝑥$ ) 𝑓 ** (𝑥$ ) 𝑓 *** (𝑥$ )

𝑓(𝑥) = 𝑓(𝑥$ ) + (𝑥 − 𝑥$ ) + (𝑥 − 𝑥$ )( + (𝑥 − 𝑥$ )) + ⋯
1! 2! 3!

Orden:
𝑇$ (𝑥) = 1

𝑇' (𝑥) = 1 + 𝑥
𝑥(
𝑇( (𝑥) = 1 + 𝑥 +
2!
𝑥(
𝑇( (𝑥) = 1 + 𝑥 +
2!
𝑥( 𝑥)
𝑇) (𝑥) = 1 + 𝑥 + +
2! 3!
𝑥( 𝑥) 𝑥,
𝑇, (𝑥) = 1 + 𝑥 + + +
2! 3! 4!

Polinomios de Taylor: La enésima suma parcial de la serie de Taylor para una función 𝑓 en 𝑎 se
conoce como el enésimo polinomio de Taylor. Por ejemplo, las sumas parciales 0, 1, 2 y 3 de la
serie de Taylor están dadas por:
𝑃$ (𝑥) = 𝑓(𝑎)

𝑃' (𝑥) = 𝑓(𝑎) + 𝑓 * (𝑎)(𝑥 − 𝑎)

𝑓 ** (𝑎)
𝑃( (𝑥) = 𝑓(𝑎) + 𝑓 * (𝑎)(𝑥 − 𝑎) + (𝑥 − 𝑎)( ,
2!
𝑓 ** (𝑎) 𝑓 *** (𝑎)
𝑃) (𝑥) = 𝑓(𝑎) + 𝑓 * (𝑎)(𝑥 − 𝑎) + (𝑥 − 𝑎)( + (𝑥 − 𝑎))
2! 3!
Error Absoluto: |𝐸𝑟𝑟𝑜𝑟 𝑅𝑒𝑎𝑙 − 𝐸𝑟𝑟𝑜𝑟 𝐴𝑝𝑟𝑜𝑥𝑖𝑚𝑎𝑑𝑜|
Error Relativo Porcentual (ER%): ER x 100=
!"
ER=
!##$# &'()

Series de Maclaurin:
%
𝑓 (")(0) "
𝑓(𝑥) = ' 𝑥
𝑘!
"&$

𝑓 * (0) 𝑓 ** (0) ( 𝑓 *** (0) )

𝑓(𝑥) = 𝑓(0) + 𝑥+ 𝑥 + 𝑥 +⋯
1! 2! 3!
 Método de la Bisección:

Supongamos que 𝑓 es una funcion continua y definida dentro del intervalo [𝑎, 𝑏] con 𝑓(𝑎) y 𝑓(𝑏)
de signos opuestos.
El teorema del valor intermedio implica que existe un numero 𝑝 en (𝑎, 𝑏) tal que 𝑓(𝑝) = 0.

𝑏' − 𝑎' 𝑎' + 𝑏'

𝑝' = 𝑎' + =
2 2

• Si (𝑝' ) = 0, entonces 𝑝 = 𝑝' y tenemos la solución.

• Si 𝑓(𝑝' ) ≠ 0 y tiene el mismo signo que 𝑓(𝑎' ), entonces 𝑎( = 𝑝' y 𝑏( = 𝑏'
• Si 𝑓(𝑝' ) ≠ 0 y tiene signos opuestos a 𝑓(𝑎' ), entonces será 𝑎( = 𝑎' y 𝑏( = 𝑝'

Volvemos a aplicar el mismo procedimiento al intervalo [𝑎( , 𝑏( ]

El intervalo [𝑎, 𝑏] será donde ocurra un cambio de signo al tomar valores y evaluar la función en
ellos mismos.
El intervalo serán los valores de la columna 𝑥, [1,2]
𝑥 𝑓(𝑥)
-2 -2
-1 -7
0 -10
1 -5
2 14

*+(
Error Absoluto:
,

Ejemplo de Tabla de iteraciones:

𝒊 𝒂𝒊 𝒃𝒊 𝒑𝒊 𝒇(𝒂𝒊) 𝒇(𝒑𝒊) 𝑬𝑨
 Método de la Regla Falsa:

𝑎' (𝑏' − 𝑎' )𝑓(𝑎' ) 𝑎' 𝑓(𝑏' ) − 𝑏' 𝑓(𝑎' )

𝑝' = =
𝑓(𝑏' ) − 𝑓(𝑎' ) 𝑓(𝑏' ) − 𝑓(𝑎' )

• Si (𝑝' ) = 0, entonces 𝑝 = 𝑝' y tenemos la solución.

• Si 𝑓(𝑝' ) ≠ 0 y tiene el mismo signo que 𝑓(𝑎' ), entonces 𝑎( = 𝑝' y 𝑏( = 𝑏'
• Si 𝑓(𝑝' ) ≠ 0 y tiene signos opuestos a 𝑓(𝑎' ), entonces será 𝑎( = 𝑎' y 𝑏( = 𝑝'

Volvemos a aplicar el mismo procedimiento al intervalo [𝑎( , 𝑏( ]

El intervalo [𝑎, 𝑏] será donde ocurra un cambio de signo al tomar valores y evaluar la función en
ellos mismos.
El intervalo serán los valores de la columna 𝑥, [1,2]
𝑥 𝑓(𝑥)
-2 -2
-1 -7
0 -10
1 -5
2 14

*+(
Error Absoluto:
,

Fórmula para el Numero de Iteraciones:

ln(𝑏 − 𝑎) − ln(𝜀 )
𝑛≥
ln(2)

Tabla para las iteraciones:

𝒊 𝒂𝒊 𝒃𝒊 𝒑𝒊 𝒇(𝒂𝒊) 𝒇(𝒑𝒊) 𝑬𝑨
 Newton-Raphson:

𝑓 (𝑥𝑖 )
𝑥𝑖 + 1 = 𝑥𝑖 −
𝑓 - (𝑥𝑖)

Se escoge un numero inicial para aplicar la formula y así realizar las iteraciones.

𝑥$ = 𝑖
Tabla de iteraciones:
𝑖 𝑥𝑖 + 1 𝐸𝐴

JACOBI – Método de Sistema de Ecuaciones:

𝐴𝑥 = 𝑏

Se despeja 𝑥' de la primera, 𝑥( de la segunda y así sucesivamente.

Elementos de la diagonal principal de 𝐴 sean distintos de 0.
La matriz debe ser diagonalmente dominante.

Norma ℓ' : ‖𝑥‖' = ∑-.&$|𝑥' |, Se realiza la suma del valor absoluto de los números del vector.

Norma ℓ( : ‖𝑥‖( Y∑-.&'|𝑥|( , Se eleva al cuadrado cada número del vector, se suman entre ellos y
después raíz cuadrada.

Norma ℓ% : ‖𝑥‖ = max |𝑥. | , 1 ≤ 𝑖 ≤ 𝑛. El numero más grande del valor absoluto de cada número
del vector.

!
Norma ℓ/ : ‖𝑥‖/ = Y∑-.&'|𝑥|/ , Se eleva a la potencia elegida el valor absoluto de cada numero.

El vector inicial debe ser de ceros: 𝑋 ($) = [0,0,0,0]

Gauss-Seidel:
No cambia la formula general, solo que las incognitas calculadas se utilizan para calcular las
siguientes. Las nuevas incognitas se sustituyen en el arreglo siguiente, para conocer el valor de la
siguiente incognita. Y asi sucesivamente.
 METODO DE PUNTO FIJO MULTIVARIABLE:

Se comienza con un vector inical de 𝑋 ($) , se calcula nuevos valores de 𝑋 (') y se repite el proceso,
esperando que despues de cada iteracion 𝑋 (") se aproxime a la raiz buscada.

El vector inical tiene que ser donde ambas ecuaciones cruzen. 𝑋 ($) = [1,1]

𝑥' = 𝑔' (𝑥' , 𝑥( … . . 𝑥- )

𝑥( = 𝑔( (𝑥' , 𝑥( … . . 𝑥- )

Se busca despejar la 𝑥 y 𝑦 de las ecuaciones dadas.

 METODO NEWTON-RAPHSON MULTIVARIABLE:

𝑓' (𝑥, 𝑦) = 0
c
𝑓( (𝑥, 𝑦) = 0

⎧𝜕𝑓' ℎ + 𝜕𝑓' 𝑗 = −𝑓 (𝑥 (") , 𝑦 (") )

⎪ 𝜕𝑥 𝜕𝑦 '

⎨𝜕𝑓( ℎ + 𝜕𝑓( 𝑗 = −𝑓 (𝑥 (") , 𝑦 (") )

⎪ 𝜕𝑥 𝜕𝑦 (
⎩

Donde ℎ = 𝑥 ("0') − 𝑥 (") y 𝑗 = 𝑦 ("0') − 𝑦 (")

𝜕𝑓' 𝜕𝑓'
𝜕𝑥 𝜕𝑦
|𝐴| = k ≠0
𝜕𝑓( 𝜕𝑓( k
𝜕𝑥 𝜕𝑦

El sistema de ecuaciones lineales tiene solucion unica, siempre que el determinante de la matriz
de coeficientes o matriz jacobiana J sea diferente de cero.

Para resolver el sistema de ecuaciones lineales se puede utilizar el Metodo de la Inversa o La

Regla de Cramer, con la solucion se calculan los valores de (𝑥 ("0') , 𝑦 ("0') ).

𝑥 ( + 𝑦 ( − 10𝑥 = 8
c (
𝑥𝑦 + 𝑥 − 8𝑦 + 5 = 0

𝑓' = 𝑥 ( + 𝑦 ( − 10𝑥 = 8
𝑓( = 𝑥𝑦 ( + 𝑥 − 8𝑦 + 5 = 0

1. Sacar la derivada de 𝑓' con respecto a 𝑥 y con respecto a 𝑦

2. Sacar la derivada de 𝑓( con respecto a 𝑥 y con respecto a 𝑦

𝑋 ($) = [1,1]

−8 2 ℎ 0
n no o = n n
2 −6 𝑗 1

|𝐴| = n−8 2 n = 44
2 −6
Se calcula el determinante de la matriz inical.
 |𝐴' | = n0 2 n = −2
1 −6
Determinante del vector ℎ y 𝑗 junto la segunda columna de la matriz.

|𝐴( | = n−8 0n = −8
2 1
Determinante de la primer columna junto a la columna ℎ y 𝑗 de la matriz.

−2 −1
ℎ= =
44 22
−8 −2
𝑗= =
44 11

𝑋 (') =

𝑋 (') = Evaluar la funcion con los resultados de la derivadas hechas anteriormente y formar una
matriz. Ademas, se multiplica por los valores de ℎ y 𝑗. Y por ultimo, se multiplica (−1).