Metodos Numericos - 2018

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 165

Cálculo de Varias Variables

MÉTODOS
NUMÉRICOS
INGENIERIA DE SISTEMAS
Métodos Numéricos
sfs

© Corporación Universitaria
Remington
Tercera edición
2018
MÉTODOS NUMÉRICOS
Elkin Ceballos Gómez.
FACULTAD DE CIENCIAS BÁSICAS E INGENIERÍA

Editorial Uniremington
Medellín, Colombia
Derechos Reservados ©2011

Primera edición: 2001


Segunda edición: 2012
Tercera edición: 2018

Responsables
Jorge Mauricio Sepúlveda Castaño
Decano de la Facultad de Ciencias Básicas e Ingeniería
[email protected]

Francisco Javier Álvarez Gómez


Coordinador CUR-Virtual
[email protected]

Edición y Montaje
Vicerrectoría de Educación a Distancia y Virtual
Equipo de diseño gráfico

www.uniremington.edu.co
[email protected]

Derechos reservados: El módulo de estudio del curso de MÉTODOS


NUMÉRICOS es propiedad de la Corporación Universitaria
Remington; las imágenes fueron tomadas de diferentes fuentes que
se relacionan en los derechos de autor y las citas en la bibliografía. El
contenido del módulo está protegido por las leyes de derechos de autor
que rigen al país. Este material tiene fines educativos y no puede
usarse con propósitos económicos o comerciales. El autor(es) certificó
(de manera verbal o escrita) No haber incurrido en fraude científico,
plagio o vicios de autoría; en caso contrario eximió de toda
responsabilidad a la Corporación Universitaria Remington y se declaró
como el único responsable.
Métodos Numéricos
sfs

3
TABLA DE CONTENIDO
Pág.

1 UNIDAD 1 SOLUCIÓN NUMÉRICA DE ECUACIONES NO LINEALES 7


1.1.1 OBJETIVO GENERAL 7
1.1.2 OBJETIVOS ESPECÍFICOS 7

1.2 TEMA 1 GENERALIDADES SOBRE LOS MÉTODOS NUMÉRICOS 7


1.2.1 DEFINICIÓN MÉTODOS NUMÉRICOS 7
1.2.2 EJERCICIOS DE APRENDIZAJE 9
1.2.3 TALLER DE ENTRENAMIENTO 11

1.3 TEMA 2 CONTINUIDAD DE FUNCIONES 11


1.3.1 CONTINUIDAD Y DISCONTINUIDAD 12
1.3.2 EJERCICIOS DE APRENDIZAJE 15
1.3.3 TALLER DE ENTRENAMIENTO 27
1.3.4 EJERCICIO DE APRENDIZAJE 27
1.3.5 TALLER DE ENTRENAMIENTO 33

1.4 TEMA 3 MÉTODOS NUMÉRICOS PARA SOLUCIONAR ECUACIONES NO LINEALES 34


1.4.1 TALLER DE APRENDIZAJE 35
1.4.2 TALLER DE APRENDIZAJE 39
1.4.3 MÉTODO DE LA POSICIÓN FALSA O MÉTODO DE LA REGLA FALSA (REGULA FALSI) 44
1.4.4 MÉTODO DE NEWTON 53
1.4.5 EJERCICIOS DE APRENDIZAJE 55
1.4.6 MÉTODO DE LA SECANTE 63
1.4.7 EJERCICIOS DE APRENDIZAJE 71

2 UNIDAD 2 REGRESIÓN E INTERPOLACIÓN 76


2.1.1 OBJETIVO GENERAL 76
2.1.2 OBJETIVOS ESPECÍFICOS 76

2.2 TEMA 1 SOLUCIÓN DE SISTEMAS DE ECUACIONES LINEALES UTILIZANDO DETERMINANTES 76


2.2.2 EJERCICIO DE ENTRENAMIENTO 82
2.2.3 EJERCICIO DE APRENDIZAJE 82
2.2.4 EJERCICIO DE APRENDIZAJE 84
2.2.5 TALLER DE ENTRENAMIENTO 87
2.2.6 EJERCICIO DE APRENDIZAJE 88

2.3 TALLER DE ENTRENAMIENTO 91

2.4 TEMA 2 REGRESIÓN POLINÓMICA 92


2.4.2 EJERCICIOS DE APRENDIZAJE 99
2.4.3 TALLER DE ENTRENAMIENTO 103

2.5 TEMA 3 INTERPOLACIÓN POLINOMIAL 104


2.5.1 TALLER DE APRENDIZAJE 106
2.5.2 EJERCICIOS DE APRENDIZAJE 114
Métodos Numéricos
sfs

4
2.5.3 EJERCICIOS DE APRENDIZAJE 121
2.5.4 TALLER DE ENTRENAMIENTO 123
2.5.5 TALLER DE APRENDIZAJE 125
2.5.6 ACTIVIDAD: REALIZAR LA PRUEBA DEL EJERCICIO 129
2.5.7 ACTIVIDAD: REALIZAR LA PRUEBA DEL EJERCICIO. 130
2.5.8 ACTIVIDAD: REALIZAR LA PRUEBA DEL EJERCICIO 132
2.5.9 TALLER DE ENTRENAMIENTO 132

3 UNIDAD 3 CÁLCULO INFINITESIMAL Y MÉTODOS NUMÉRICOS 133


3.1.1 OBJETIVO GENERAL 133
3.1.2 OBJETIVOS ESPECÍFICOS 133

3.2 TEMA 1 DIFERENCIACIÓN NUMÉRICA 133


3.2.1 EJERCICIO DE APRENDIZAJE 134
3.2.2 EJERCICIOS DE APRENDIZAJE 140
3.2.3 APLICACIONES EN GEOMETRÍA 141
3.2.4 TALLER DE ENTRENAMIENTO 143

3.3 TEMA 2 INTEGRACIÓN NUMÉRICA 144


3.3.1 TALLER DE APRENDIZAJE 147
3.3.2 TALLER DE ENTRENAMIENTO 153

3.4 TEMA 3 SOLUCIÓN NUMÉRICA DE ECUACIONES DIFERENCIALES 153


3.4.1 MÉTODO DE EULER (O MÉTODO DE LA RECTA TANGENTE) 154
3.4.2 EJERCICIO DE APRENDIZAJE 154
3.4.3 MÉTODO MEJORADO (O RUG KUTA DE SEGUNDO ORDEN) (O MÉTODO PREDICTOR – CORRECTOR DE
EULER) 156
3.4.4 EJERCICIOS DE APRENDIZAJE 156
3.4.5 TALLER DE ENTRENAMIENTO 159

4 PISTAS DE APRENDIZAJE 162

5 GLOSARIO 164

6 BIBLIOGRAFÍA 165
Métodos Numéricos
sfs

PROPÓSITO GENERAL

MÉTODOS NUMÉRICOS

El propósito de este módulo es abordar algunos temas de matemáticas y cálculo


utilizando el análisis numérico, ya que, los métodos numéricos proporcionan una serie
de herramientas que le brindan al estudiante y futuro profesional, elementos para
solucionar, situaciones, donde
Métodos Numéricos
sfs

MÉTODOS NUMÉRICOS

OBJETIVO GENERAL
Aplicar los métodos de manera eficiente en la solución de problemas que
involucran modelos matemáticos, procurando que la solución obtenida
mediante la aplicación de los diferentes algoritmos sea óptima, precisa y exacta.

OBJETIVOS ESPECÍFICOS
 Utilizar algunos métodos para la solución numérica de ecuaciones no
lineales.

 Conocer y aplicar métodos para la interpolación numérica de conjuntos de


datos y para la regresión polinómica.

 Conocer y aplicar algunos métodos numéricos en el cálculo infinitesimal


(derivadas, integrales y solución del problema de valor inicia

UNIDAD1 UNIDAD 2 UNIDAD 3

Solución Numérica De
Cálculo infinitesimal y
Ecuaciones No E Interpolación
Métodos numéricos
Lineales
Métodos Numéricos
sfs

7
1 UNIDAD 1 SOLUCIÓN NUMÉRICA DE ECUACIONES NO
LINEALES
 Video: Método de Newton para ecuaciones
 ”Ver Video”: https://www.youtube.com/watch?v=5boNNInwkhk

1.1.1 OBJETIVO GENERAL


Utilizar algunos métodos para la solución numérica de ecuaciones no lineales.

1.1.2 OBJETIVOS ESPECÍFICOS


 Reconocer los conceptos necesarios para iniciar un curso de métodos

 Determinar los puntos de discontinuidad de algunas funciones.

 Reconocer algunos métodos para la solución numérica de ecuaciones no lineales.

1.2 TEMA 1 GENERALIDADES SOBRE LOS MÉTODOS


NUMÉRICOS
1.2.1 DEFINICIÓN MÉTODOS NUMÉRICOS
Los métodos numéricos son técnicas que permiten solucionar problemas
utilizando operaciones repetitivas de suma, resta, multiplicación y división.
Todo método numérico presenta dos características:

 Muchas operaciones aritméticas.

 Presenta un margen de error. El error deseado es inversamente proporcional al número


de cálculos.
Métodos Numéricos
sfs

8
1.2.1.1 MANEJO DE DECIMALES
Cuando se utilizan métodos numéricos, los cálculos realizados, por lo general contienen decimales
infinitos, estos decimales se pueden tratar de la siguiente manera:

 TRUNCAMIENTO:

Consiste en cortar un número a las cifras significativas indicadas, sin


hacer aproximaciones.

 REDONDEO:

Consiste en aproximar un número a las cifras significativas indicadas.

 REGLAS PARA EL REDONDEO

 Si el digito siguiente a la última cifra significativa es 5 o mayor, la última cifra


significativa se aumenta en 1.

 Si el dígito siguiente a la última cifra significativa es menor que 5, la última cifra


significativa se deja igual.

 Enlaces para cifras significativas.

 Video: Cifras significativas


 ”Ver Video”: https://www.youtube.com/watch?v=nZLlbs44yaw

 Video: CIFRAS SIGNIFICATIVAS


 ”Ver Video”: https://www.youtube.com/watch?v=YDcud6J7jTE
Métodos Numéricos
sfs

9
1.2.2 EJERCICIOS DE APRENDIZAJE
1. Redondee el número 3,469823... con cuatro cifras decimales.
Solución

La quinta cifra es 2, por lo tanto, la cuarta cifra, que es 8, se deja igual, entonces:

3,469823...  3,4698
2. Redondee el número 0,324686... con cuatro cifras decimales.

Solución

La quinta cifra decimal es ocho, por lo tanto, se debe sumar uno a la cuarta cifra, que es 6, entonces:

0,324686...  0,3247
3. Redondee el número 2,356342 con una cifra decimal.

Solución

La segunda cifra decimal es cinco, se debe sumar uno a la primera cifra decimal, que es 3, entonces:

2,356342...  2,4
4. Trunque el número 63,45678...con dos cifras decimales.
Solución

Como el número se trunca, simplemente se escribe con las dos primeras cifras decimales:

63,45678...  63,45

1.2.2.1 CONVERGENCIA
Se dice que un método converge cuando se llega al valor buscado, esto se
puede identificar por qué se puede observar que cada cálculo efectuado lleva a
un valor cada vez más parecido al anterior.
Métodos Numéricos
sfs

10
1.2.2.2 DIVERGENCIA
UN MÉTODO DIVERGE CUANDO NO SE LLEGA AL VALOR BUSCADO, ESTO SE IDENTIFICA,
PORQUE CADA CÁLCULO EFECTUADO LLEVA A UN VALOR CADA VEZ DIFERENTE.

1.2.2.3 TEORÍA DE ERRORES.


ERROR ES LA CANTIDAD DE FALLA O DE IMPRECISIÓN QUE SE ESTÁ
DISPUESTOS A ACEPTAR CUANDO SE EFECTÚA UN CÁLCULO O CUANDO SE
HACE ALGUNA MEDICIÓN.
EL ERROR SE OBTIENE COMO LA DIFERENCIA ENTRE EL VALOR VERDADERO Y
EL VALOR CALCULADO.

E  Vverdadero Vcalculado
 Error absoluto: Es el valor absoluto de la diferencia entre el valor verdadero y el
valor calculado o medido. Posee las mismas unidades que la magnitud que se mide.
E= Vverdadero Vcalculado

Otra manera de calcular este error es restando los dos últimos cálculos.
𝑬𝒂 = 𝑬𝒗 − 𝑬𝑪
𝑬𝒂 = 𝒙𝒏 − 𝒙𝒏−𝟏
 Error relativo: Es el cociente entre el error absoluto y el valor absoluto del valor
verdadero, es un número a dimensional y en algunos casos se expresa en tanto por
ciento.

|𝑽𝒗𝒆𝒓𝒅𝒂𝒅𝒆𝒓𝒐 − 𝑽𝒄𝒂𝒍𝒄𝒖𝒍𝒂𝒅𝒐 |
𝑬𝒓 =
|𝑽𝒗𝒆𝒓𝒅𝒂𝒅𝒆𝒓𝒐 |

También se puede expresar de la siguiente manera:

|𝒙𝒏 − 𝒙𝒏−𝟏 |
𝑬𝒓 =
𝒙𝒏
Métodos Numéricos
sfs

11
1.2.3 TALLER DE ENTRENAMIENTO
1) De acuerdo a lo trabajado en el desarrollo del tema anterior, realice los siguientes
ejercicios:

 Redondee el número 0,03456743.... con tres cifras decimales.

 Redondee el número 12,5678342.... con cinco cifras decimales.

 Redondee el número 304,45666666


... con dos cifras decimales.

 Redondee el número 3,445566776.7. con cuatro cifras decimales.

 Trunque el número 4,6780123... con tres cifras decimales.

 Trunque el número 103,6789675... con una cifra decimal.

 Trunque el número - 0,245764... con tres cifras decimales.

 Trunque el número 20,4567349... con cinco cifras decimales.

1.3 TEMA 2 CONTINUIDAD DE FUNCIONES


Los métodos numéricos exigen que la función que se va a utilizar sea continua
en algún intervalo.

 Video: Continuidad de una funcion


 ”Ver Video”: https://www.youtube.com/watch?v=wuAdn84VSoc

A continuación, se recordará la forma de determinar cuándo una función es continua o


discontinua.
Métodos Numéricos
sfs

12
1.3.1 CONTINUIDAD Y DISCONTINUIDAD
 Video: Tipos de discontinuidad AINTE Mat 1 Bach
 ”Ver Video”: https://www.youtube.com/watch?v=RU4v-7u4blI

 Video: UDEM CD 2.6 discontinuidades


 ”Ver Video”: https://www.youtube.com/watch?v=NmSH6LJbVs8

Se dice que una función es continua en todo su dominio, cuando se puede recorrer toda la gráfica
sin tener que levantar la mano, cuando no hay huecos o espacios entre sus gráficas, si algo de
esto se llega a presentar se dice que la función es discontinua.

La función y = f (x) mostrada en la figura 1 es continua, porque se puede recorrer toda su gráfica
sin necesidad de levantar la mano

La función y = g(x) mostrada en la figura 2 es discontinua (no es continua), porque al recorrer


su gráfica hay que levantar la mano para continuar, o porque hay un espacio entre sus gráficas.
Esta función es discontinua en el punto x=1.

Nota: Para indicar puntos de discontinuidad se debe nombrar la x del punto.

Figura 1. Función Continua (Autor Elkin Ceballos Gómez abril 10 de 2011)


Métodos Numéricos
sfs

13

Figura 2. Función discontinua. (Autor Elkin Ceballos Gómez abril 10 de 2011)

1.3.1.1 CONTINUIDAD DE ALGUNAS FUNCIONES


A continuación, se recordará la forma de determinar cuándo una función es continua, para ello
es necesario conocer su dominio.

 Video: DOMINIOS DE FUNCIONES


 ”Ver Video”: https://www.youtube.com/watch?v=qOCMPXoxJyg

 Video: EJERCICIOS RESUELTOS DE DOMINIO DE UNA FUNCION|Graciela Alvarez


 ”Ver Video”: https://www.youtube.com/watch?v=c4TeoNGlcM0
Métodos Numéricos
sfs

14
Función Polinómica.

El dominio de las funciones polinómicas son todos los números reales, esto se escribe así:

Dom. x lR
Las funciones polinómicas son siempre continuas

 Función racional.

Es una función de la forma:

𝒈(𝒙)
𝒚 = 𝒇(𝒙) = 𝒉(𝒙), con 𝒉(𝒙) ≠ 𝟎

 Se identifica porque la función tiene 𝒙 en el denominador.


 La función racional es la función cociente de dos funciones polinómicas.

Su dominio está formado por todos los números reales menos las asíntotas verticales de la función.

 Video: Problema sobre el cálculo del dominio y asíntotas de una función que es
un cociente de polinomios
 ”Ver Video”: https://www.youtube.com/watch?v=gSDOXzMx_h8

Una asíntota vertical es un valor de 𝒙 donde el denominador se hace cero.

Para determinar las asíntotas verticales de una función racional (si tiene) se procede de la siguiente
manera:

1) guale el denominador a cero.

2) Solucione la ecuación resultante. Si la ecuación no tiene solución, quiere decir que la


función racional no tiene asíntotas verticales.

3) Los valores de x obtenidos son las asíntotas verticales de la función.


Métodos Numéricos
sfs

15

DISCONTINUIDAD DE LAS FUNCIONES RACIONALES

1.3.2 EJERCICIOS DE APRENDIZAJE


Determine el dominio y la continuidad de cada una de las siguientes funciones:
𝟓𝒙
1. 𝒚 = 𝒘(𝒙) =
𝟐𝒙−𝟕

SOLUCIÓN:

 DOMINIO:

- Iguale el denominador a cero:


𝟕
𝟐𝒙 − 𝟕 = 𝟎 → 𝒙 = , Es la asíntota vertical; se debe excluir del dominio.
𝟐

Por lo tanto, el dominio está dado por:


𝟕
𝒙 ∈ ℝ 𝒚 𝒙 ≠ , esto es:
𝟐

𝟕 𝟕 𝟕
𝒙 ≠ → 𝑫𝒇 = (−∞, ) ∪ ( , +∞)
𝟐 𝟐 𝟐

 CONTINUIDAD:

𝟓𝒙 𝟕
La función: 𝒚 = 𝒘(𝒙) = es discontinua en 𝒙 =
𝟐𝒙−𝟕 𝟐

𝟒𝒙 − 𝟕
𝟐. 𝒚 = 𝒘(𝒙) =
𝟐𝒙𝟐 − 𝒙 − 𝟔

Solución:

 Dominio:
Métodos Numéricos
sfs

16
Se iguala el denominador a cero:

𝟐𝒙𝟐 − 𝒙 − 𝟔 = 𝟎
Factorizando el polinomio se tiene:

𝟐 𝟒𝒙𝟐 −𝟏(𝟐)𝒙−𝟔 (𝟐𝒙−𝟒)(𝟐𝒙+𝟑)


(𝟐𝒙𝟐 − 𝒙 − 𝟔 = 𝟎) → =𝟎→ =𝟎→
𝟐 𝟐 𝟐

𝟐(𝒙 − 𝟐)(𝟐𝒙 + 𝟑)
= 𝟎 → (𝒙 − 𝟐)(𝟐𝒙 + 𝟑) = 𝟎 →
𝟐
(𝒙 − 𝟐) = 𝟎 → 𝒙 = 𝟐
𝟑
𝟐𝒙 + 𝟑 = 𝟎 → 𝒙 = −
𝟐
Por lo tanto, el Dominio de la función está dado por:

𝟑
𝑫𝑭 = ℝ − {− , 𝟐}
𝟐
También se puede escribir con el siguiente intervalo:

𝟑 𝟑
𝑫𝑭 = (−∞, − ) ∪ (− , 𝟐) ∪ (𝟐, +∞)
𝟐 𝟐

 Continuidad:

𝟒𝒙−𝟕
𝑒𝑠 𝒅𝒊𝒔𝒄𝒐𝒏𝒕𝒊𝒏𝒖𝒂 𝑒𝑛 𝑙𝑜𝑠 𝑝𝑢𝑛𝑡𝑜𝑠: 𝒙 = − 𝟐 𝒚 𝒙 = 𝟐
𝟑
La función: y  w (x) 
𝟐𝒙𝟐 −𝒙−𝟔

 Determinar el dominio y la continuidad de la siguiente función:

𝟑
𝒚 = 𝒇(𝒙) =
𝒙𝟐 − 𝟏
Métodos Numéricos
sfs

17
 Dominio: Se toma el denominador y se hace igual a cero:

𝒙𝟐 − 𝟏 = 𝟎 → (𝒙 + 𝟏)(𝒙 − 𝟏) = 𝟎 →

𝒙 + 𝟏 = 𝟎 → 𝒙 = −𝟏

𝒙−𝟏=𝟎→𝒙=𝟏
El dominio de la función está dado por:

𝑫𝑭 = ℝ − {−𝟏, 𝟏}
También se puede escribir con el siguiente intervalo:

𝑫𝑭 = (−∞, −𝟏) ∪ (−𝟏, 𝟏) ∪ (𝟏, +∞)


 Continuidad:

𝟑
La función: 𝒚 = 𝒇(𝒙) = :
𝒙𝟐 −𝟏

𝐸𝑠 𝒅𝒊𝒔𝒄𝒐𝒏𝒕𝒊𝒏𝒖𝒂 𝑒𝑛 𝑙𝑜𝑠 𝑝𝑢𝑛𝑡𝑜𝑠: 𝒙 = −𝟏 𝒚 𝒙 = 𝟏

Función irracional. Es una función de la forma:

𝒏
𝒚 = 𝒇(𝒙) = √𝒈(𝒙)

Se identifica porque tiene variable dentro de un radical.

 DOMINIO:

Para determinar su dominio se presentan dos casos:

Caso #1:

Si la raíz es impar el dominio corresponde a todos los números reales: Dom.𝐷𝑓 =


{𝑥/ x lR}
Métodos Numéricos
sfs

18

Caso #2:

Si la raíz es par, se debe garantizar que los valores que se le asignen a “x”
hagan positivo todo dentro de la raíz.

Procedimiento para hallar el dominio de funciones irracionales de raíz o índice par:

1. Plantee la inecuación:𝒈(𝒙) ≥ 𝟎, donde 𝒈(𝒙)𝑒𝑠 𝑙𝑎 𝑝𝑎𝑟𝑡𝑒 𝑠𝑢𝑏𝑟𝑎𝑑𝑖𝑐𝑎𝑙.


2. Se soluciona la inecuación.
3. La solución de la inecuación es el dominio de la función irracional.
4. Si la inecuación planteada no tiene solución o si toda dentro de la raíz par es positiva,
quiere decir que el dominio de la función es todos los reales.
𝐷𝑓 = {𝑥/𝑥 ∈ ℝ}= (−∞, +∞)
Enlaces para función irracional.

 Video: Dominio de Una Funcion Irracional Video 1


 ”Ver Video”: https://www.youtube.com/watch?v=BCjQrBEBGdU

 Video: Hallar dominio


 ”Ver Video”: https://www.youtube.com/watch?v=GwkbhPJiHDk

CONTINUIDAD DE FUNCIONES IRRACIONALES:

Una función irracional es continua en todo su dominio.

Para indicar la continuidad de una función irracional, se nombra


los mismos intervalos del dominio.

Ejercicios de Aprendizaje
Métodos Numéricos
sfs

19
1. Determine el dominio y la continuidad de:

𝒚 = 𝒇(𝒙) = √𝒙𝟐 + 𝟐𝒙 − 𝟏𝟓

SOLUCIÓN:

Se tiene una raíz par.

DOMINIO: Para calcularlo se procede de la siguiente manera:

PASOS:

1. Todo lo que está dentro de la raíz debe ser mayor o igual que cero, por lo tanto, se debe
plantear la inecuación:

𝒙𝟐 + 𝟐𝒙 − 𝟏𝟓 ≥ 𝟎
1. Encuentre los números ceros de la expresión resultante. Esto es, iguale a cero y
resuelva la ecuación resultante, los valores obtenidos son los ceros de la expresión.
Un cero es un número donde posiblemente hay cambio de signo
en la expresión.
Para el ejemplo:

x2  2x 15  0  x  5x  3 x  5  0  x  3  0 

x  5  x  3
Estos son los ceros o números críticos.

1. Ubique los ceros en la recta numérica. Véase figura 3.

2. Evalúe el signo de la expresión obtenida en el paso uno. Para ello se toma un número que
se encuentre a la izquierda del primer cero, se toma un número que se encuentre entre
ambos ceros y se toma un número que se encuentre a la derecha del segundo cero. Estos
números se reemplazan en la expresión obtenida en el paso uno y el signo del resultado
se escribe en la recta numérica. Véase la figura 3.
Métodos Numéricos
sfs

20
3. La respuesta o solución de la inecuación, se da tomando los intervalos que cumplan con
el sentido de la desigualdad. Para ello observe el sentido de la desigualdad de la expresión
obtenida en el paso dos.

Si dice: > 0 Se toman los + + + + +, sin incluir los números críticos.

Si dice: ≥ 0 Se toman los + + + + +, incluyendo los números críticos.

Figura 3. Dominio de una función irracional.


(Autor Elkin Ceballos Gómez)

La solución de la inecuación es: x   ,5 3, 


Este es el mismo dominio de la función: 𝑫𝒇 = {𝒙/ x   ,5 3, +
CONTINUIDAD:

La función es continua en:  ,53, 



Ejemplo2: Determine el dominio y la continuidad de la siguiente función:

𝒚 = 𝒇(𝒙) = √𝟑𝒙 − 𝟓
Métodos Numéricos
sfs

21
SOLUCIÓN

Se debe solucionar la siguiente inecuación:

𝟓
3x  5  0  3x  5  0  3x  5  x 
𝟑

Ubicando este valor en la recta numérica, se tiene que:

Figura 4. Dominio de una función


racional. (Autor Elkin Ceballos
Gómez)
5
𝐷𝑓 = [ , +∞)
3

Continuidad: La función es continua en [5 / 3, )

Función Exponencial.
Es una función de la forma:

𝒚 = 𝒇(𝒙) = 𝒃𝒈(𝒙)
Métodos Numéricos
sfs

22
Con: 𝒃 >𝟎𝒚𝒃≠𝟏
 DOMINIO:

El dominio depende de g(x).

 CONTINUIDAD:

Una función exponencial es continua en todo su dominio.

Ejercicios de Aprendizaje
Determine el dominio y la continuidad de las siguientes funciones:

1. 𝒚 = 𝒇(𝒙) = 𝟓𝒙
Es una función exponencial que tiene en el exponente un polinomio, su dominio
corresponde a todos los números reales.

𝒙
𝑫𝒇 = { ∈ ℝ}
𝒙
La función es siempre continua.

2. y  f (x)  e2x1
SOLUCIÓN:

 Dominio: Como en el exponente hay una función lineal, entonces su dominio está dado por todos

los números reales:

𝑫𝒇 = {𝒙/𝒙 ∈ ℝ}
 Continuidad: La función es continua en todo su dominio.

𝟓𝒙
3. y 𝒇(𝒙) = 
𝒆𝒙−𝟒
Métodos Numéricos
sfs

23
SOLUCIÓN

Como en el exponente hay una función racional, el dominio será todos los reales menos las asíntotas
verticales.

x40x4

 Dominio:
𝑫𝒇 = ℝ − {𝟒} →= (−∞, 𝟒) ∪ (𝟒, +∞)

 Continuidad: La función es discontinua en x 4

Enlaces para funciones exponenciales

 Video: FUNCION EXPONENCIAL (MATEMATICA)


 ”Ver Video”: https://www.youtube.com/watch?v=cXnw6kzqASI

 Video: FUNCION EXPONENCIAL GRAFICA 1


 ”Ver Video”: https://www.youtube.com/watch?v=pogo0SmmSbE

FUNCIÓN LOGARÍTMICA: Es una función de la forma:

𝒚 = 𝒇(𝒙) = 𝐥𝐨𝐠 𝒃 𝒈(𝒙)


Con b > o y b≠ 𝟏
 Dominio: Para hallar el dominio se de
be resolver la inecuación:

𝒈(𝒙) > 𝟎

 Continuidad:
Las funciones logarítmicas son continuas en todo su dominio.
Métodos Numéricos
sfs

24

Taller de Entrenamiento

1. Determine el dominio y la continuidad de las siguientes funciones:

2. y  f (x)  log2 (2x  3)

SOLUCIÓN

DOMINIO:
𝟑
2x  3  0  2x  3  0  2x  3  𝒙 =
𝟐

Ubicando en la recta numérica se tiene:

Figura 5. Dominio de una función


logarítmica. (Autor Elkin Ceballos
Gómez)
𝟑
 El dominio de la función está dado por: 𝑫𝒇 = (𝟐 , +∞)

𝟑
 Continuidad: La función es continua en el mismo intervalo del dominio, (𝟐 , +∞)

𝟐. 𝒚 = 𝒇(𝒙) = 𝐥𝐨𝐠(𝒙𝟐 − 𝟒)

SOLUCIÓN:
Métodos Numéricos
sfs

25
DOMINIO: Para calcular el dominio se hace: x2 40
Factorizando se tiene:

x2  4  0  x  2x  2 0  x  2  0  x  2  0

x  2  x  2

Ubicando estos puntos en la recta numérica, se tiene:

Figura 6. Dominio de una función logarítmica. (Autor Elkin Ceballos Gómez)

Dom. x  ,22, 


Continuidad: La función es continua en  ,2 2, 

 Funciones trigonométricas.
Se analizará el dominio y la continuidad de las funciones trigonométricas para el
caso en que el ángulo es x.

 Función seno. (http://www.youtube.com/watch?v=XmAzV3ubl_U)

Es la función de la forma : 𝒚 = 𝒇(𝒙) = 𝐬𝐢𝐧 𝒙


Y su gráfica está determinada por:
Métodos Numéricos
sfs

26

El dominio de esta función es todos los reales, por lo tanto, esta


función es siempre continua.

 Función coseno. (http://www.youtube.com/watch?v=neiQSmqWI8w):


Es la función que se representa como:

𝒚 = 𝒇(𝒙) = 𝒄𝒐𝒔𝒙
Y su gráfica está determinada por:

El dominio de esta función es todos los reales, por lo tanto, esta función
es siempre continua.
Métodos Numéricos
sfs

27
1.3.3 TALLER DE ENTRENAMIENTO
Para cada una de las siguientes funciones determine: Dominio, Continuidad o discontinuidad
(intervalos donde es discontinua) y su respectiva gráfica:

Nota: Recuerde que estas funciones se pueden expresar en función de sen 𝑥 𝑦 cos 𝑥.

 𝑦 = sec 𝑥
 𝑦 = csc 𝑥
 𝑦 = 𝑡𝑎𝑛𝑥
 𝑦 = cot 𝑥

1.3.3.1 SEMILLA
Es una estimación inicial de la raíz. Es el valor inicial con el cual se empieza a operar el
método numérico que se va a aplicar. La estimación de la semilla se puede hacer de dos
formas diferentes:
 Haciendo un bosquejo de la gráfica y el (los) punto(s) donde la gráfica corte el eje x es
(son) la(s) raíz (o raíces) de la ecuación, como el ojo no nos puede dar buena exactitud,
este valor sólo sirve como semilla para iniciar los cálculos.
 Otra forma puede ser evaluando la función en diferentes valores de x, entre los puntos
donde se observe cambio de signo se encuentra una raíz, por lo tanto, se escoge como
semilla un valor de x que se encuentre entre dichos puntos.

1.3.4 EJERCICIO DE APRENDIZAJE

f (x)  3x3  9x2  3x  2 es siempre continua por ser un polinomio

1. Dada la siguiente ecuación:

3x3  9x2  3x  2  0
Su gráfica se observa en la figura 8.
Métodos Numéricos
sfs

28
La ecuación tiene raíces en los siguientes intervalos:
1, 0, 0, 1 2, 3,
Esto es cierto, ya que en estos intervalos la gráfica de la función corta
el eje x.

Esto se puede confirmar, dando valores a x de la siguiente manera:

x  1  f (1)  313  912  31  2  3  9  3  2  13

x  0  f (0)  303  902  30  2  0  0  0  2  2

x  1  f (1)  313  912  31  2  3  9  3  2  1

x  2  f (2)  323  922  32  2  24 36 6  2  4

x  3  f (3)  333  932  33  2  81 81 9  2  11

Ubicando estos puntos en la siguiente tabla de valores:

x -1 0 1 2 3

y -13 2 -1 -4 11

Se puede observar que de:


-1 y 0 la y cambia de signo (de signo – a signo +)
0 a 1 la y cambia de signo + a signo menos.

2 a 3 la y cambia de signo – a signo +.

En resumen:
La ecuación 3x3  9x2  3x  2  0 tiene raíces en:

1, 0 Porque:
Métodos Numéricos
sfs

29
f (1) f (0)  0 (o tienen signo diferente) y

f (x) es continua en este intervalo.


La semilla para la iniciar el cálculo de la raíz en este intervalo es:

−𝟏 + 𝟎
𝒙𝟏 = = −𝟎, 𝟓
𝟐
0,1 Porque f (0) f (1)  0 (o tienen signo diferente) yf (x) es continua en este intervalo.

La semilla para la iniciar el cálculo de la raíz en este intervalo es:

𝟏+𝟎
𝒙𝟏 = = 𝟎, 𝟓
𝟐
2, 3 Porque
f (2) f (3)  0 (o tienen signo diferente) y f (x) es continua en este intervalo

La semilla para la iniciar el cálculo de la raíz en este intervalo es:

𝟐+𝟑
𝒙𝟐 = = 𝟐, 𝟓
𝟐
Métodos Numéricos
sfs

30

Figura 8. Gráfica de una función.


(Autor Elkin Ceballos Gómez abril 10 de 2011)

Generalizando se tiene que:

Si una función f (x) es continua en un intervalo a, b y además, f a y f b tienen:

diferente signo (esto se simboliza como f af b 0 ) entonces, la ecuación f (x)  0 tiene
una raíz en a, b ; la semilla para determinar esta raíz se calcula como:

𝒂+𝒃
𝒙=
𝟐
2. Determine en que intervalos se encuentran las raíces de la ecuación:

4x4  9x  7  0.
Determine la semilla para cada raíz.

SOLUCIÓN

Por ser
Métodos Numéricos
sfs

31
f (x)  4x4  9x  7 una función polinómica, es continua para cualquier valor de x:

Dando diferentes valores x, obteniendo la respectiva y;

x  3  y  4 3  9(3)  7  y  344


4

x  2  y  4 2  9(2)  7  y  75
4

x 1  y  41  9(1)  7  y  6
4

x  0  y  40  9(0)  7  y  7
4

x  1  y  41  9(1)  7  y  12


4

x  2  y  42  9(2)  7  y  39
4

x  3  y  43  9(3)  7  y  290


4

Se llena la siguiente tabla.

x -3 -2 -1 0 1 2 3

y 344 75 6 -7 -12 39 290

1. Se observa cambio de signo entre -1 y cero, se pasa de + a -.

2. Se observa cambio de signo entre 1 y 2, pasa de – a +.

Se concluye que:

La ecuación 4x4 - 9x- 7 = 0 tiene dos raíces reales.

 Raíz número 1 está en el intervalo [1, 0, la semilla es:

−𝟏+𝟎
𝒙𝟏 = = −𝟎, 𝟓
𝟐
Métodos Numéricos
sfs

32

Raíz número 2 está en el intervalo [1, 2], la semilla es:

𝟏+𝟐
𝒙𝟐 = = 𝟏, 𝟓
𝟐
3. Determine en que intervalos se encuentran las raíces de la ecuación 5 ln x 3x 4 0
Determine la semilla para cada raíz.

SOLUCIÓN:

El dominio de: f (x)  5ln x  3x  4 es el dominio de 𝐥𝐧 𝒙

Un logaritmo no acepta ni números negativos, ni el cero, entonces el dominio de la


función está dado por:
𝑫𝒇 = {𝒙/𝒙 ∈ (𝟎, +∞)}
f (x)  5ln x  3x  4,

Dando diferentes valores x, obteniendo la respectiva y;

x  0.01 y  5ln0,01 3(0,01)  4  19,0558...

x  0.5  y  5ln0,5 3(0,5)  4  0,9657...

x  1  y  5ln1 3(1)  4  1

x  2  y  5ln2 3(2)  4  1,4657... x

 3  y  5ln3 3(3)  4  0,4930...

x  4  y  5ln4 3(4)  4  1,0685...


Métodos Numéricos
sfs

33
Se llena la siguiente tabla.

x -0,01 0,5 1 2 3 4

y -19,0558… -0,9657… 1 1,4657… 0,4930… -1,0685

Se observa cambio de signo entre 0,01 y 0,5, se pasa de - a +.

Se observa cambio de signo entre 3 y 4, pasa de + a -


Se concluye que:
La ecuación 5 ln x - 3x +4 = 0 tiene dos raíces reales.

Raíz número 1 está en el intervalo [ 0.01, 0.5, la semilla es:

−𝟎. , 𝟏 + 𝟎, 𝟓
𝒙𝟏 = = 𝟎, 𝟐𝟒𝟓
𝟐
Raíz número 2 está en el intervalo [3, 4, la semilla es:

𝟑+𝟒
𝒙𝟐 = = 𝟑. 𝟓
𝟐

1.3.5 TALLER DE ENTRENAMIENTO


Determine semillas para las raíces de cada una de las siguientes ecuaciones:

1. 5x3 10x2  8x  3  0

2. 2 x  3x 2  5  0

3. x2  3 x  6  0

4. 3x 2  5senx 1  0 Solo dos semillas.


5. 3x  tan𝛼1  0 Solo dos semillas.
Métodos Numéricos
sfs

34
1.4 TEMA 3 MÉTODOS NUMÉRICOS PARA SOLUCIONAR
ECUACIONES NO LINEALES
En los cinco métodos que se explica a continuación, se parte del mismo problema:

Problema: Solucionar la ecuación:

f (x)  0

Se sabe que hay una raíz en el intervalo a, b porque:

f (x) es continua en el intervalo a, b

Porque también se sabe que fa f b 0

Método Bisección

Este un método simple, confiable que siempre converge, pero requiere de


muchos pasos para lograr la exactitud deseada.

 Video: Calculo Numérico - Bisección


 ”Ver Video”: https://www.youtube.com/watch?v=5e1GUYfnc9I

Para aplicar este método tenga en cuenta el siguiente proceso.

PASO UNO: Lleve la ecuación a la forma f (x) 0

PASO DOS: Para la primera raíz halle su correspondiente semilla:

𝒂 𝟏 + 𝒃𝟏
𝒙𝟏 =
𝟐
Cerciórese que se cumpla que:

 f (x) Sea continua en a1 , b1 


 f a1  f b1  0
Métodos Numéricos
sfs

35
PASO TRES: Calcule f x1 . Pueden resultar tres posibilidades:

POSIBILIDAD CARACTERÍSTICA

f x1   0 En este caso la raíz es x1 se termina el proceso.

f (x1) f a1   0 Esto quiere decir que la raíz se encuentra en el intervalo a1 , x1 .
o Para este caso haga a2 =a1  b2  x1

f (x1 )  f b1   0
Esto quiere decir que la raíz se encuentra en el intervalo x1 ,
b1 .

o Para este caso haga a2  x1  b2  b1

PASO CUATRO: Calcule:

𝒂𝟐 + 𝒃𝟐
𝒙𝟐 =
𝟐
Regrese al paso 3.

El proceso termina cuando el error sea menos o igual al error deseado.

Para ilustrar el método de bisección, observe los siguientes dos ejemplos.

1.4.1 TALLER DE APRENDIZAJE

Encuentre una raíz real de la ecuación: f (x)  5x3 16x2  6x  4  0 con una precisión de
0.01,

Es decir, entre dos cálculos consecutivos se debe repetir dos decimales


SOLUCIÓN

Si x  3, y  5 3 16 3  6 3  4  13


3 2

Dando valores a x entre – 2 y 3 se llena la siguiente tabla:


Métodos Numéricos
sfs

36
x -3 -2 -1 0 1 2 3
y -13 8 1 -4 23 112 293

Al observar la tabla se puede ver que la ecuación tiene tres aíces, así: Raíz número 1

Está en el intervalo [ 3,  2 esto se puede afirmar, ya que, f (x)  5x3 16x2  6x  4 es

continua en el intervalo  3,  2 y porque f  3 f  2 tiene signos contrarios o lo que es


lo

mismo: f  3 f  2 0

Raíz número 2

Está en el intervalo 1, 0esto se puede afirmar, ya que f (x)  5x3 16x2  6x  4
continua en el intervalo 1, 0 y porque:

f 1 f 0 tiene signos contrarios o lo que es lo mismo f 1 f 0 0

Raíz número 3

Está en el intervalo 0, 1 esto se puede afirmar, ya que,

f (x)  5x3 16x2  6x  4 es continua, en el intervalo 0, 1 y porque f 0 f 1 tiene


signos contrarios o lo que es lo mismo

f 0 f 1 0
 Procedimiento para hallar la raíz número 3 con una precisión de 0,01 o 10-2.

Paso uno. La ecuación ya está en la forma f (x)  0

Paso dos. Cálculo de la semilla:

𝟎+𝟏
𝒙𝟏 = = 𝟎, 𝟓
𝟐
Paso tres. Se calcula f (x1 )
Métodos Numéricos
sfs

37
Como f 0,5 f 10

f 0,5  50,5 160,5  60,5  4  3,626


3 2

La raíz se encuentra en: 0, 0.5

Paso cuatro. Cálculo de x2

𝟎 + 𝟎, 𝟓
𝒙𝟏 = = 𝟎, 𝟐𝟓
𝟐
Comparando x2  0,25 con x1  0,5 , no se repite ninguna cifra decimal, por lo tanto, se debe
continuar el proceso en el paso tres:

Cálculo de f x2 

f 0.25  50.25 160.25  60.25  4  -1.421875.


3 2

f 0,25 f 0,5 0 , la raíz está en:

[0,25, 0,5]
Cálculo de x3
𝟎, 𝟐𝟓 + 𝟎, 𝟓
𝒙𝟑 = = 𝟎, 𝟑𝟕𝟓
𝟐
Observando 𝑥3 ∧ 𝑥2 no se repite decimal

f 0.375  50.375 160.375  60.375  4  0.763671857...


3 2

Como

f 0,25 f 0.375 0 , la raíz está en [0.25, 0.375], se calcula 𝒙𝟒

Cálculo de

𝟎, 𝟐𝟓 + 𝟎, 𝟑𝟕𝟓
𝒙𝟒 = = 𝟎, 𝟑𝟏𝟐𝟓
𝟐
Métodos Numéricos
sfs

38
Observando x4  x3 se repite un decimal (deben ser dos). Se continúa el proceso.

Se calcula: f x4 

f 0.3125  50.3125 160.3125  60.3125  4  -0.4099120193...


3 2

Como: f 0,3125 f 0.375 0 , la raíz está en [0.3125, 0.375]

Cálculo de x5
𝟎, 𝟑𝟏𝟐𝟓 + 𝟎, 𝟑𝟕𝟓
𝒙𝟓 = = 𝟎, 𝟑𝟒𝟑𝟕𝟓
𝟐

Observando x5  x4 se repite un decimal (deben ser dos). Se continúa el proceso.


Se calcula:

f 0.34375  50.34375 160.34375  60.34375  4  0.1562194284...


3 2

Como: f 0,3125 f 0.34375 0 , la raíz está en [ 0


.3125, 0.34375]
Cálculo de x6

𝟎, 𝟑𝟏𝟐𝟓 + 𝟎, 𝟑𝟒𝟑𝟕𝟓
𝒙𝟔 = = 𝟎, 𝟑𝟐𝟖𝟏𝟐𝟓
𝟐
Observando x6  x5 se repite un decimal (deben ser dos). Se continúa el proceso.

Se calcula f x6 

f 0.328125  50.328125 160.328125  60.328125  4  -0.1319549131..


3 2

Como: f 0,328125 f 0.34375 0 , la raíz está en: [0


.328125, 0.34375]

Cálculo de x7
𝟎, 𝟑𝟐𝟖𝟏𝟐𝟓 + 𝟎, 𝟑𝟒𝟑𝟕𝟓
𝒙𝟕 = = 𝟎, 𝟑𝟑𝟓𝟗𝟑𝟕𝟓
𝟐
Métodos Numéricos
sfs

39
Observando x7  x6 se repite un decimal (deben ser dos). Se continúa el proceso.

Se calcula f x7 

f 0.3359375
  50.3359375 160.3359375  60.3359375  4  0.0108485218..
3 2

f 0,3359375 f 0.34375  0 , la raíz está en:


[0
.328125, 0.3359375] 

Cálculo de x8
𝟎, 𝟑𝟐𝟖𝟏𝟐𝟓 + 𝟎, 𝟑𝟑𝟓𝟗𝟑𝟕𝟓
𝒙𝟖 = = 𝟎, 𝟑𝟑𝟐𝟎𝟑𝟏𝟐𝟓
𝟐

Observando x8  x7 se repiten dos decimales, por lo tanto el proceso ya termina.

La raíz es 𝒙𝟖  0.33203125

1.4.2 TALLER DE APRENDIZAJE

Queda como ejercicio hallar las otras dos raíces utilizando el mismo error.

2. Con una precisión de 10-3 encuentre una raíz de real de la ecuación ln x  cos x  0
Solución:

El dominio de la función cosx son todos los reales. Dom. x  lR

El dominio de la función ln x es Dom: x  0,  

El dominio de la función f x  ln x  cosx se obtiene mediante la intersección del


dominio de
Métodos Numéricos
sfs

40
cada función que la conforman, es decir, Dom. x  0,  
Para determinar en qué intervalo hay una raíz, se debe dar valore a x desde cero
(sin incluirlo) hasta que se obtenga un cambio de signo.

La función f x  ln x  cosx es continua en 0,  


Antes de dar valores a x, no olvide que la calculadora debe estar en modo radianes.

x  0.01 y  ln 0.01 cos0.01 -5.605120816...

x  0.5  y  ln 0.5  cos0.5  -1. 570729472..

x  1  y  ln1 cos1  -0. 5403020358..

x  1.5  y  ln1.5  cos1.5  0.3347279604...

Se puede observar cambio de signo para x  1  x  1.5


Una raíz está en el intervalo [𝟏, 𝟏, 𝟓]. Esto es cierto ya que 𝒇(𝒙)es continua
en [𝟏, 𝟏, 𝟓] y además se cumple que: 𝒇(𝟏). 𝒇(−𝟏, 𝟓) < 𝟎 o lo que es lo mismo, tienen
signos contrarios.

Cálculo de la semilla

𝟏 + 𝟏, 𝟓
𝒙𝟏 = = 𝟏, 𝟐𝟓
𝟐
Cálculo de f x1 f x1   ln1,25  cos1,25  -0. 09217818108..

Como f 1.25 f 1,5 0 , la raíz está en el intervalo [𝟏, 𝟐𝟓, 𝟏, 𝟓] 

Cálculo de x2
𝟏, 𝟐𝟓 + 𝟏, 𝟓
𝒙𝟐 = = 𝟏, 𝟑𝟕𝟓
𝟐
Métodos Numéricos
sfs

41
Comparando x2  1,375con x1  1,25, no se repite ninguna cifra decimal, por lo tanto, se debe
continuar el proceso.

Cálculo de f x2 f 1,375  ln (1.375) - cos1.375  0.1239060321

Como: f 1,25 f 1,375 0 , la raíz está en el intervalo [𝟏, 𝟐𝟓, 𝟏, 𝟑𝟕𝟓]

Cálculo de x3
𝟏, 𝟐𝟓 + 𝟏, 𝟑𝟕𝟓
𝒙𝟑 = = 𝟏, 𝟑𝟏𝟐𝟓
𝟐

Observando x3  x2 se repite decimal (deben ser tres) se continua el proceso

Se calcula f x3 f 1.3125   ln1.3125  cos1.3125   0. 01649998458..

Como f 1,25 f 1.3125 0 , la raíz está en el intervalo: [𝟏, 𝟐𝟓, 𝟏, 𝟑𝟏𝟐𝟓]

Cálculo de 𝒙𝟒
𝟏, 𝟐𝟓 + 𝟏, 𝟑𝟏𝟐𝟓
𝒙𝟒 = = 𝟏, 𝟐𝟖𝟏𝟐𝟓
𝟐

Observando x4  x3 no se repite ningún decimal (deben ser tres). Se continúa el


proceso.

Se calcula f x4 f 1.28125   ln1.28125   cos 1.28125  -0. 03768103221..

f 1.28125 f 1.3125  0, la raíz está en el intervalo [1


. ,28125, 1.3125]
Cálculo de x5
Métodos Numéricos
sfs

42
𝟏, 𝟐𝟖𝟏𝟐𝟓 + 𝟏, 𝟑𝟏𝟐𝟓
𝒙𝟓 = = 𝟏, 𝟐𝟗𝟔𝟖𝟕𝟓
𝟐

Observando x5  x4 se repite un decimal (deben ser tres). Se continúa el proceso.

Se calcula f x5 f 1.296875   ln1.296875   cos1.296875   -0. 0105511247..

Como f 1.296875 f 1.3125 0, la raíz está en el intervalo[𝟏, 𝟐𝟗𝟔𝟖𝟕𝟓, 𝟏, 𝟑𝟏𝟐𝟓 ]

Cálculo de x6

𝟏, 𝟐𝟗𝟔𝟖𝟕𝟓 + 𝟏, 𝟑𝟏𝟐𝟓
𝒙𝟔 = = 𝟏, 𝟑𝟎𝟒𝟔𝟖𝟕𝟓
𝟐

Observando x6  x5 no se repite decimal (deben ser tres). Se continúa el proceso.

Se calcula f x6 f 1.3046875  ln1.3046875   cos1.3046875   0.


002984312150..

Como f 1.296875   f 1.3046875   0 , la raíz está en el intervalo:


[𝟏, 𝟐𝟗𝟔𝟖𝟕𝟓, 𝟏, 𝟑𝟎𝟒𝟔𝟖𝟕𝟓]

Cálculo de x7
𝟏, 𝟐𝟗𝟔𝟖𝟕𝟓 + 𝟏, 𝟑𝟎𝟒𝟔𝟖𝟕𝟓
𝒙𝟕 = = 𝟏, 𝟑𝟎𝟎𝟕𝟖𝟏𝟐𝟓
𝟐

Observando x7  x6 se repiten dos decimales (deben ser tres). Se continúa el


proceso.

Se calcula f x7  f 1.30078125) ln 1.30078125) cos 1.30078125) =


Métodos Numéricos
sfs

43
-0. 003780291734..

Como 𝑓(1,30078125). 𝑓 1,3046875 < 0, por lo tanto la raíz está en el intervalo:

[𝟏. 𝟑𝟎𝟎𝟕𝟖𝟏𝟐𝟓, 𝟏. 𝟑𝟎𝟒𝟔𝟖𝟕𝟓 ]

(-) (+)
Cálculo de x8

𝟏. 𝟑𝟎𝟎𝟕𝟖𝟏𝟐𝟓 + 𝟏. 𝟑𝟎𝟒𝟔𝟖𝟕𝟓
𝒙𝟖 = = 𝟏, 𝟑𝟎𝟐𝟕𝟑𝟒𝟑𝟕𝟓
𝟐

Observando x8  x7 se repiten dos decimales (deben ser tres), por lo tanto el


proceso continua.

f 1.302734357  ln 1.302734357  cos 1.302734357  -0. 0003978616062..

Como f 1.302734375) f 1.30468755  0 , la raíz está en el intervalo:

[𝟏. 𝟑𝟎𝟐𝟕𝟑𝟒𝟑𝟕𝟓, 𝟏. 𝟑𝟎𝟒𝟔𝟖𝟕𝟓 ]

(-) (+)

Cálculo de x9

𝟏, 𝟑𝟎𝟐𝟕𝟑𝟒𝟑𝟕𝟓 + 𝟏. 𝟑𝟎𝟒𝟔𝟖𝟕𝟓
𝒙𝟗 = = 𝟏. 𝟑𝟎𝟑𝟕𝟏𝟎𝟗𝟑 𝟕
𝟐

Observando x8  x9 se repiten dos decimales (deben ser tres), por lo tanto el


proceso continua.

Se calcula f x9  f 1.303710973  ln 1.303710973  cos 1.303710973  0.001293437606.


Métodos Numéricos
sfs

44
Como f 1.302734375 f 1.303710937 0 , la raíz está en el intervalo:

[𝟏. 𝟑𝟎𝟐𝟕𝟑𝟒𝟑𝟕𝟓, 𝟏. 𝟑𝟎𝟑𝟕𝟏𝟎𝟗𝟑 ]

(-) (+)

Cálculo de x10

𝟏, 𝟑𝟎𝟐𝟕𝟑𝟒𝟑𝟕𝟓 + 𝟏. 𝟑𝟎𝟑𝟕𝟏𝟎𝟗𝟑 𝟕
𝒙𝟏𝟎 = = 𝟏, 𝟑𝟎𝟑𝟐𝟐𝟐𝟔𝟓𝟓
𝟐
Observando x9  x10

se repiten tres decimales, por lo tanto, el proceso termino. La raíz es:

𝒙𝟏𝟎 == 𝟏, 𝟑𝟎𝟑𝟐𝟐𝟐𝟔𝟓𝟓

1.4.3 MÉTODO DE LA POSICIÓN FALSA O MÉTODO DE LA REGLA FALSA


(REGULA FALSI)

 Video: Metodo de falsa posicion xvid xvid


 ”Ver Video”: https://www.youtube.com/watch?v=pSyOwt7PkvQ

 Video: Metodo de la Falsa Posición


 ”Ver Video”: https://www.youtube.com/watch?v=7qIPoyDgMM8

Si la función 𝑓(𝑥) es continua en algún intervalo que contenga la raíz buscada, este método
siempre converge.

Sea f (x) una función continua en el intervalo a, bcon f af b 0, sea
r una raíz de f x  0 en el intervalo a, b.
Véase la figura 9
Métodos Numéricos
sfs

45

Figura 9. Deducción de la Regla falsa (Autor Elkin Ceballos Gómez)

El método consiste en encontrar la ecuación de la recta que pasa por los puntos:

(𝒂, 𝒇(𝒂)) ∧ (𝒃, 𝒇(𝒃))

El punto x1 que es el punto donde la recta corta el eje x es la primera aproximación


para la raíz r en el intervalo a, b.

DEDUCCIÓN DE LA FÓRMULA DE LA REGULA FALSI


PASO UNO: Haga a1  a  b1  b
PASO DOS: Encuentre la ecuación de la recta que une los puntos:
Métodos Numéricos
sfs

46
((𝒂𝟏, 𝒇 (𝒂𝟏 ) ∧ (𝒃𝟏,𝒇( 𝒃𝟏 ))
𝑓(𝑏)1 −𝑓(𝑎1 )
La pendiente de la recta está dada por: 𝑚 = 𝑏1 −𝑎1

La ecuación es:

𝒇 (𝒃𝟏 )−𝒇 (𝒂𝟏 )


y  f a   (𝒙 − 𝒂)  

𝒃𝟏 −𝒂𝟏

Cuando:

𝒙 = 𝒙𝟏 , 𝒚 = 𝟎
La ecuación queda:

𝑓 (𝑏1 ) − 𝑓 (𝑎1 )
0 − 𝑓(𝑎1 ) (𝑥1 − 𝑎1 )
𝑏1 − 𝑎1
Despejando x1
𝑏1 − 𝑎1
𝑥1 = 𝑎1 − 𝑓(𝑎1 ) [ ]
𝑓 (𝑏1 ) − 𝑓 (𝑎1 )
𝑎1 𝑓 (𝑏1 ) − b1 𝑓 (𝑎1 )
𝑥1 =
𝑓 (𝑏1 ) − 𝑓 (𝑎1 )
PASO TRES: Halle

f x1  . Puede suceder que

 Si f x1   0 , quiere decir que el proceso terminó y la raíz de la ecuación es x1

 Si f a1 f x1  0 , quiere decir que la raíz se encuentra en a1 , x1  en este caso,
haga:

 a2  a1  b2  x1
Métodos Numéricos
sfs

47
 De lo contrario tome a2  x1  b2  b1

PASO CUATRO: Se aplica nuevamente el proceso anterior al intervalo a2 , b2  , es decir, haga.

𝑎2 𝑓 (𝑏2 ) − b2 𝑓 (𝑎2 )
𝑥=
𝑓 (𝑏2 ) − 𝑓 (𝑎2 )
PASO CINCO: Después de varias iteraciones se tiene que:

𝑎n 𝑓 (𝑏n ) − bn 𝑓 (𝑎n )
𝑥=
𝑓 (𝑏n ) − 𝑓 (𝑎n )
El proceso se termina cuando:

|𝑥𝑛 − 𝑥𝑛−1 | Error deseado


EJERCICIO DE APRENDIZAJE

1. Utilizando el método “Regula Falsi” con un error de 10-2 obtenga una de las
raíces reales de la ecuación:

3x4  2x  9  0
PROCEDIMIENTO

La función polinómica. f x  3x4  2x  9 es continua para cualquier valor de x, por ser una
función polinómica.

Dando valores a x, reemplazando en la función para determinar la y, se obtiene los valores de la


siguiente tabla:

x -2 -1 0 1 2
y 35 -8 -9 -4 43

Se observa cambio de signo entre 𝑥 = −2 𝑦 𝑥 = −1 ; y entre 𝑥 = 1 𝑦 𝑥 = 2

NOTA:
Métodos Numéricos
sfs

48
La ecuación no tiene más raíces reales, ya que, es posible comprobar, dando
más valores a x, que a la izquierda de x = - 2 la función asignará a y solo valores
positivos (no hay más cambio de signos hacia la izquierda) y a la derecha de x =
2 tampoco habrá más cambio de signos. Esto se puede asegurar, porque el signo
que predomina en ambos costados es el de 3x 4 que es positivo.

Las raíces reales de la ecuación 3x4  2x  9  0 se encuentran en:

Raíz número 1: Esta en [−2 − 2.] ya que f  2 f 1 0  f x es continua en  2, 1
+ −

Raíz número 2: Esta en [1; 2.] ya que f 1 f 2 0  f x es continua en 1, 2
− +

CÁCULO DE LA RAÍZ NÚMERO 2

𝑎1 𝑓 (𝑏1 ) − b1 𝑓 (𝑎1 ) 1𝑓 (2) − 2𝑓 (1) 1(3(2)4 + 2(2) − 9) − 2(3(1)4 + 2(1) − 9)


𝑥1 = ⇒ 𝑥1 = =
𝑓 (𝑏1 ) − 𝑓 (𝑎1 ) 𝑓 (2) − 𝑓 (1) (3(2)4 + 2(2) − 9) − (3(1)4 + 2(1) − 9)

Es conveniente escribirlo de la siguiente manera para poder utilizar la calculadora correctamente.

    
x  1 32  22  9  2 31  21  9 / 32  22  9  31  21  9
4 4 4 4

𝒙𝟏  1.08510638
1

Cálculo de f (x1 )

f x   31.085106328...  21.085106328...  9  -2.670580130


4

La raíz se encuentra en:

[1.085106382 … , 2]

− +

Cálculo de x2
Métodos Numéricos
sfs

49
   
x  1 .085106328... 32  22  9  2 31.085106328...  21.085106328...  9
4 4

  
2
/ 32  22  9  31.085106328...  21.085106328...  9
4 4

x2  1.138604649...

Comparando x2  x1 No se repiten decimales, el proceso continúa.

Cálculo de f x2 

f x2  31.138604694...  21.138604694...  9  -1.68067200


4

La raíz está en:

[1.138604649 … , 2]

− +

Cálculo de x3

   
x  1 .138604694... 32  22  9  2 31.138604694...  21.295463593...  9
4 4

3
 
/ 32  22  9  31.138604694...  21.138604694...  9
4 4

x3  1.17100619...

Comparando 𝑥3 ∧ 𝑥2 Se repite un decimal, el proceso continúa.


Cálculo de 𝑓(𝑥3 ):
Métodos Numéricos
sfs

50
1.1 1.1
𝑓(𝑥3 ) = 3 ( ) + 2( ) − 9 = −1.01696214
1.00619 … 1.00619 …

1.1
La raíz está en: [ … 2]
1.006199…

− +
Cálculo de x4

   
𝑥4 =  1 .17100619... 32  22   9  2 31.17100619...  21.17100619...  9
4 4

  
4
/ 32  22  9  31.17100619...  21.17100619...  9
4 4

𝒙𝟒 = 𝟏. 𝟏𝟗𝟎𝟏𝟓𝟗𝟏𝟕𝟎 …
Comparando 𝑥4 ∧ 𝑥3 Se repite un decimal, el proceso continúa.
Cálculo de 𝑓(𝑥4 ):
𝑓(𝑥4 ) = 3(1.190159170)4 + 2(1.190159170) − 9 = −0.600444630

La raíz está en: [1.190159170 … 2]

− +
Cálculo de x5

   
x  1 .190159107... 32  22  9  2 31.190159107...  21.190159107...  9
4 4

  
5
/ 32  22  9  31.190159107...  21.190159107...  9
4 4
Métodos Numéricos
sfs

51

𝑥5  1.20131194...
Comparando 𝑥5 ∧ 𝑥4 Se repite un decimal, el proceso continúa.
Cálculo de 𝑓(𝑥5 ):
𝑓(𝑥5 ) = 3(1.20131194 … ) + 2(1.20131194 … ) − 9 = 1.20131194 …

La raíz está en: [1.201311914 … … 2]

− +
Cálculo de x6

   
x  1 .201311941... 32  22  9  2 31.201311941...  21.201311941...  9
4 4

 
6
/ 32  22  9  31.201311941...  21.201311941...  9
4 4

x6  1.20774808…
Comparando 𝑥6 ∧ 𝑥5 Se repite dos decimales, el proceso terminó.
Con una precisión de 10-2, la raíz de la ecuación 3x4  2x  9  0 es x6  1.20774808..
Queda como ejercicio encontrar la raíz número 1 utilizando el mismo método y la misma
precisión.

1. Teorema del valor medio:

El teorema del valor medio establece que:

Si f (x) es una función continua en un intervalo a, b, existe un número x en a, b tal
que:
Métodos Numéricos
sfs

52
𝑓(𝑏) − 𝑓(𝑎)
𝑓 (𝑥 ) =
𝑏−𝑎
Utilizando el método de la regla falsa, con una precisión de 10-5, encuentre el número

x garantizado por el teorema valor medio para la función:


5𝑥 − 3
𝑓 (𝑥 ) = 2 𝑒𝑛 [3, 6]
𝑥 −4
Procedimiento: Calculando la primera derivada de la función, se tiene que:

( 2
´( )
5𝑥 𝑥 − 4) − 2𝑥(5𝑥 − 3)
𝑓 𝑥 =
(𝑥 2 − 4)2
Realizando las operaciones indicadas y reduciendo términos semejantes, se tiene:

2
´( )
−5𝑥 + 6𝑥 − 20
𝑓 𝑥 =
(𝑥 2 − 4)2

Se sabe que:

𝒂 = 𝟑 𝒚 𝒃 = 𝟔, reemplazando en la ecuación original se tiene que:


12 27
𝑓(𝑎) = y 𝑓(𝑏) = , por lo tanto: f (b)  f (a)  249/160
5 32

𝑓(𝑏) − 𝑓(𝑎) 83
=−
𝑏−𝑎 160
Métodos Numéricos
sfs

53
se debe plantear la ecuación:

−5𝑥 2 + 6𝑥 − 20 83
= −
(𝑥 2 − 4)2 160

Simplificando esta ecuación, queda:

83𝑥 4 − 1464𝑥 2 + 960𝑥 − 1872 = 0

Ejercicio de Entrenamiento:

Encontrar la raíz de la ecuación anterior en el intervalo [3, 6]


Recuerde: para hallar esta raíz, tenga como modelo los ejercicios desarrollados

La solución para el ejercicio está dada por: 𝑥 = 4.01957

1.4.4 MÉTODO DE NEWTON

 Video: Cálculo de raíces por el Método de newton raphson


 ”Ver Video”: https://www.youtube.com/watch?v=PrJsNAR-rhA

 Video: El metodo de Newton


 ”Ver Video”: https://www.youtube.com/watch?v=73McLogMkoA

 Video: Método numérico de Newton


 ”Ver Video”: https://www.youtube.com/watch?v=WCvFJG7ouXQ
Métodos Numéricos
sfs

54
 Video: MÉTODO DE NEWTON
 ”Ver Video”: https://www.youtube.com/watch?v=_qNllgwlXwc

Se tiene que 𝒙𝟏 es una aproximación inicial para la raíz 𝒓 como se ve en la figura 10.
Observe que la línea tangente a la curva en
(𝒙𝟏 , f (𝒙𝟏 )) corta el eje x en el punto (𝒙𝟐 , 𝟎) ,0) y que 𝒙𝟐 es una mejor aproximación a 𝒓
que 𝒙𝟏 .
La ecuación de esta línea tangente es: y  f (𝒙𝟏 )  f '(𝒙𝟏 )(x  𝒙𝟏 ) Se puede encontrar
𝒙𝟐 utilizando la ecuación de la línea tangente.
Como 𝒙𝟐 , 𝟎 está en la línea tangente, sus coordenadas deben satisfacer la
ecuación anterior, por lo tanto:

𝒇 (𝒙𝟏 )
𝟎 − f (𝒙𝟏 )= f '(𝒙𝟏 )( 𝒙𝟐  𝒙𝟏 )− =𝒙𝟐 − 𝒙𝟏 ⟹
𝒇 ′(𝒙𝟏 )

𝒇 (𝒙𝟏 )
𝒙𝟐 = 𝒙𝟏 −
𝒇 ′(𝒙𝟏 )

Una mejor aproximación para r es el valor 𝒙𝟐 .


Repitiendo lo hecho anteriormente se obtiene que:

𝒇 (𝒙𝟐 )
𝒙𝟑 = 𝒙𝟐 −
𝒇 ′(𝒙𝟐 )
Generalizando el proceso se puede decir que:

𝒇 (𝒙𝒏 )
𝒙𝒏+𝟏 = 𝒙𝒏 − , donde 𝒏 = 𝟏, 𝟐, 𝟑, …
𝒇 ′(𝒙𝒏 )

Nota: El proceso se termina cuando se alcanza una exactitud


deseada.
Métodos Numéricos
sfs

55

Figura 10. Determinación del método de Newton (Autor Elkin Ceballos Gómez)

1.4.5 EJERCICIOS DE APRENDIZAJE


1. Utilizando el método de Newton, con una precisión de 10-7 encuentre una raíz
real de la ecuación:

La función 𝑓(𝑥 ) = 𝑠𝑒𝑛𝑥 − 4𝑥 + 20 = 0 es continua en todo su dominio, ya que


𝑠𝑒𝑛𝑥 tiene como dominio todos los números reales y la función  4x tiene como dominio
todos los números reales, es decir la función es continua.

Para determinar intervalos para las raíces, se da valores a x. De esta manera se llena la siguiente tabla:
x 0 1 2 3 4 5
y f 0  20 f 1  16.8 f 2  12.9 f (3)  8.1 f 4  3.2 f 5 -0.9
Métodos Numéricos
sfs

56
4, 5
Se observa el cambio de signo en el intervalo: [ ]
+ −
4, 5
La ecuación 𝑓 (𝑥 ) = 𝑠𝑒𝑛𝑥 − 4𝑥 + 20 = 0 tiene una raíz en el intervalo[ ]
+ −
ya que:

4, 5
𝑓(4). 𝑓(5) < 0 𝑦 𝑓 (𝑥 )𝑒𝑠 𝑐𝑜𝑛𝑡𝑖𝑛𝑢𝑎 𝑒𝑛 [ ]
+ −

 Se calcula la derivada de la Función dada:

Si 𝑓 (𝑥 ) = 𝑠𝑒𝑛𝑥 − 4𝑥 + 20, su derivada está dada por:

𝑓 1 (𝑥 ) = 𝑐𝑜𝑠𝑥 − 4
 Se realiza el cálculo de la semilla:

𝟒+𝟓 𝟗
𝒙𝟏 = = = 𝟒. 𝟓
𝟐 𝟐

 Cálculo de 𝒙𝟐 :
𝒇 (𝒙𝟏 ) 𝒔𝒆𝒏(𝟒. 𝟓) − 𝟒(𝟒. 𝟓) + 𝟐𝟎
𝒙𝟐 = 𝒙𝟏 − = 𝟒. 𝟓 − = 𝟒. 𝟕𝟒𝟐𝟖𝟐𝟏𝟎𝟓𝟔 …
𝒇𝟏 (𝒙𝟏 ) 𝒄𝒐𝒔(𝟒. 𝟓) − 𝟒

Comparando 𝑥2 ∧ 𝑥1 Se repite un decimal, el proceso continúa.


 Cálculo de 𝒙𝟑 :
Métodos Numéricos
sfs

57
𝟒. 𝟕𝟒𝟐𝟖𝟐𝟏𝟎𝟓𝟔 … − (𝒔𝒆𝒏(𝟒. 𝟕𝟒𝟐𝟖𝟐𝟏𝟎𝟓𝟔 … ) − 𝟒(𝟒. 𝟕𝟒𝟐𝟖𝟐𝟏𝟎𝟓𝟔 … ) + 𝟐𝟎
𝑥3 =
(𝒄𝒐𝒔(𝟒. 𝟕𝟒𝟐𝟖𝟐𝟏𝟎𝟓𝟔 … ) − 𝟒)

𝑥3 = 4.75017167 …

Comparando 𝑥3 ∧ 𝑥2 Se repite un decimal, el proceso continúa.


 Cálculo de 𝒙𝟒 :
𝟒. 𝟕𝟓𝟎𝟏𝟕𝟏𝟔𝟕 … − (𝒔𝒆𝒏(𝟒. 𝟒. 𝟕𝟓𝟎𝟏𝟕𝟏𝟔𝟕 … ) − 𝟒(𝟒. 𝟕𝟓𝟎𝟏𝟕𝟏𝟔𝟕 … ) + 𝟐𝟎
𝑥4 =
(𝒄𝒐𝒔(𝟒. 𝟕𝟓𝟎𝟏𝟕𝟏𝟔𝟕 … ) − 𝟒)

𝑥4 = 4.7501784845

Comparando 𝑥4 ∧ 𝑥3 Se repiten cinco decimales, el proceso continúa.


 Cálculo de 𝒙𝟓 :
𝟒.𝟕𝟓𝟎𝟏𝟕𝟖𝟒𝟖𝟒𝟓.(𝒔𝒆𝒏(𝟒.𝟕𝟓𝟎𝟏𝟕𝟖𝟒𝟖𝟒𝟓.)−𝟒(𝟒.𝟕𝟓𝟎𝟏𝟕𝟖𝟒𝟖𝟒𝟓.)+𝟐𝟎
 𝑥5 = (𝒄𝒐𝒔(𝟒.𝟕𝟓𝟎𝟏𝟕𝟖𝟒𝟖𝟒𝟓.)−𝟒)

𝑥5 = 4.75017848.5…

Comparando 𝑥5 ∧ 𝑥4 Se repiten todos los decimales, el proceso terminó.


La raíz es: 𝑥5 = 4.75017848.5

Utilizando el método de Newton, con una precisión de 10-8 Encuentre una raíz real de la ecuación:

3x2  e x  6  0
Por tener una combinación de función polinómica y función exponencial (con polinomio en el
exponente la función f (x)  3x2  ex  6 es continua en todo su dominio.

Es fácil comprobar que la función tiene una raíz en el intervalo 2, 3, ya que:
Métodos Numéricos
sfs

58
f 2  32  e 2  6  -1.38905698
2

f 3  33  e3  6  0.9144630678


2

 Determinación de la derivada:

La función es: f (x)  3x2  ex  6

Su derivada es: f '(x)  6x  ex

 Determinación de la semilla:

2+3 5
𝑥1 = = = 2.5
2 2
Nota: Para escribir la expresión e 2.5 en la calculadora, debe digitar las
teclas:
SHIFT ln 2.5

 Cálculo de 𝑥2
𝒇 (𝒙𝟏 ) 𝟑(𝟐. 𝟓)𝟐 − 𝒆𝟐.𝟓 − 𝟔
= 𝒙𝟏 − = 𝟐. 𝟓 −
𝒇𝟏 (𝒙𝟏 ) 𝟔(𝟐. 𝟓) − 𝒆𝟐.𝟓

𝒙𝟐 = 𝟐. 𝟐𝟗𝟖𝟓𝟕𝟖𝟓𝟗 …

Nota: Con la calculadora debe escribir:

 
x  2.5  32.5  SHIFTln2.5  6 /62.5  SHIFTln2.5
2

2
x2  2.29857859...
Métodos Numéricos
sfs

59
 Cálculo de 𝑥3 :
𝑥3 = = 
2.29857859... 32.29857859...  SHIFTln2.29857859...  6
2

/ 62.29857859...  SHIFTln2.29857859...

𝑥3 = 2.3271903...

 Cálculo de 𝑥4 :
𝑥4 = = 
2.3271903... 32.3271903...)  SHIFTln2.3271903…) 6
2

/ 62.3271903…   SHIFTln2.3271903…) 

𝑥4 = 2.327637268...

 Cálculo de 𝑥5 :
𝑥5 = 2.327637268... 32.327637268...) 2
 SHIFTln2.327637268... ) 6 
/ 62.327637268...   SHIFTln2.327637268... ) 

𝑥5 = 2.3276373282..

 Cálculo de 𝑥6 :
2
 𝑥6 = 2.327637282... 32.327637282...)  SHIFTln2.327637282...) 6

 / 62.327637282...  SHIFTln2.327637282... ) 

 𝑥6 = 2.327637382..
Métodos Numéricos
sfs

60
Observando los dos últimos cálculos, se repiten todos los decimales.

La raíz es: 𝑥6  2.32763738.2..

3. Problemas de cajas:

Se desea construir una caja sin tapa y de forma rectangular con una lámina de cartón d e 50 cm *
35 cm. Para ello se cortarán cuadrados idénticos en las cuatro esquinas y se doblarán los lados hacia
arriba

Si la caja será hecha para contener un volumen 𝟑𝟓𝟎𝟎𝒄𝒎𝟑 con una precisión de 8 decimales
determine.
a. Las dimensiones de la caja.
b. La cantidad de material utilizado.

Nota: Cada raíz se debe aproximar con tres decimales.

Las raíces son: x  2.63064..., x  11.8864...,  x  27.9829...

PROCEDIMIENTO:

Sea x el lado del cuadrado a quitar.

x: Altura.
50  2x: Largo de la caja

 2x: Ancho de la caja.

Se plantea la ecuación:

3500 x35 2x50  2x

Que corresponde al volumen de la caja.

La ecuación a resolver queda: 4x3 170x2 1750x  3500 0


Métodos Numéricos
sfs

61
Solucionando la ecuación presentada se tiene que las raíces se encuentran en 2, 3, 11,
12 27, 28
Queda como ejercicio de entrenamiento encontrar las raíces con el método de
Newton.

Las raíces son: x  2.63064..., x  11.8864...,  x  27.9829...

El valor de x  27.9829 queda descartado, ya que daría longitudes negativas para el largo y
para el ancho.

Con x  2.631
x: Altura.  2.631

50  2x: Largo de la caja 50  2.631 * 2  44.738

35  2x: Ancho de la caja.  35  2.631 * 2  29.738

La cantidad de material es:

44.738* 29.738 2* 22.631* 29.378 2* 44.738* 2.631 1722.311 cm2

Con x  11.8864
x: Altura.  11.886

50  2x: Largo de la caja  50  2 *11.886  26.228 cm.


35  2x: Ancho de la caja.  35  2 *11.886  11.228 cm.

3. Problemas de punto de equilibrio.

La función de ingreso en cierta empresa es: rq 150q  0.2q2 $ y la función de costo
total es: cq q3  80q  600 $ donde q es el número de docenas de cierto artículo
Métodos Numéricos
sfs

62
producidas y vendidas. Determine el punto de equilibrio con una precisión de 4
decimales.

El punto de equilibrio se presenta cuando:

rq  cq
Esto es:

150q  0.2q2 $  q3  80q  600 $

La ecuación a resolver es:

q3  0.2q2  230q  600  0


Actividad: Queda como ejercicio encontrar las raíces de esta ecuación.

Respuesta; Las raíces están en:

17, 16. q  16.4259...Descartada.


(Explica el por qué se descarta)
2, 3. q  2.70068...

13, 14. q  13.5252...


NOTA: Si su elección para la aproximación inicial 𝑥1 , hace que la primera derivada sea igual a cero,
escoja un número diferente que sea cercano a la raíz deseada.
Métodos Numéricos
sfs

63
1.4.6 MÉTODO DE LA SECANTE

 Video: Método de la Secante


 ”Ver Video”: https://www.youtube.com/watch?v=X6zixq7_SZQ

 Video: Metodo de la secante


 ”Ver Video”: https://www.youtube.com/watch?v=aRveHW4cLg4

Se sabe por definición que:

1 𝐿𝑖𝑚 𝑓 (𝑥 ) − 𝑓((𝑥𝑛+1 )
𝑓 (𝑥𝑛+1 ) =
𝑥 → (𝑥𝑛+1 ) 𝑥 − 𝑥𝑛+1
Haciendo x  xn y evaluando el límite, tenemos que:
𝒇(𝒙𝒏 ) − 𝒇((𝒙𝒏+𝟏 )
𝒇𝟏 (𝒙𝒏+𝟏 ) ≈
𝒙𝒏 − 𝒙

Reemplazamos esta expresión en la fórmula de Newton, queda:

𝒇(𝒙𝒏+𝟏 )
𝒙𝒏+𝟐 = 𝒙𝒏+𝟏 − 𝟏
𝒇 (𝒙𝒏+𝟏 )
𝟏
Reemplazando 𝒇 (𝒙𝒏+𝟏 ) Se tiene que:
𝒇(𝒙𝒏+𝟏 )
𝒙𝒏+𝟐 = 𝒙𝒏+𝟏 −
𝒇(𝒙𝒏 ) − 𝒇((𝒙𝒏+𝟏 )
𝒙𝒏 − 𝒙𝒏+𝟏
Métodos Numéricos
sfs

64

Realizando la operación indicada:

𝒇(𝒙𝒏+𝟏 )(𝒙𝒏 − 𝒙𝒏+𝟏 )


𝒙𝒏+𝟐 = 𝒙𝒏+𝟏 −
𝒇(𝒙𝒏 ) − 𝒇((𝒙𝒏+𝟏 )
Se obtiene:

(𝒙𝒏+𝟏 ). 𝒇(𝒙𝒏 ) − 𝒙𝒏 . 𝒇(𝒙𝒏+𝟏 )


𝒙𝒏+𝟐 =
𝒇(𝒙𝒏 ) − 𝒇((𝒙𝒏+𝟏 )

Nota: Para utilizar esta fórmula se deben tener dos


aproximaciones iniciales 𝒙𝟎  𝒙𝟏
Ejercicios de Aprendizaje

Utilizando el método de la secante, con una precisión de 10−6 encuentre una raíz real para la
ecuación:

32x  35x2  3  0
Procedimiento:

La función 𝒇(𝒙) =32x  35x2  3  0 es continua en todo su dominio.


Dando valores a x, y reemplazando en la función, se obtiene la siguiente tabla:
x -3 -2 -1 0 1 2 3
y - - - 4 -23 - 417
311.9… 136.9… 31.8… 56

Se puede observar cambio de signo entre:

x = -1 y x = 0,
Métodos Numéricos
sfs

65
x=0yx=1,y

x = 2 y x = 3.

Las raíces de la ecuación están en:

]
−1, 0
 Raíz número 1: Está en [
− +

0, 1
 Raíz número 2: Está en [ ]
+ −

]
2, 3
 Raíz número 3: Está en [
− +

CÁLCULO DE LA RAÍZ NÚMERO 3.

𝒙𝒐  2  𝒙𝟏  3
 Cálculo de 𝒙𝟐

(𝒙𝟏 ). 𝒇(𝒙𝟎 ) − 𝒙𝟎 . 𝒇(𝒙𝟏 )


𝒙𝟐 =
𝒇(𝒙𝟎 ) − 𝒇(𝒙𝟏 )

Reemplazando por los respectivos valores se tiene:

𝟑. 𝒇(𝟐) − 𝟐. 𝒇(𝟑)
𝒙𝟐 =
𝒇(𝟐) − 𝒇(𝟑)
𝒙𝟐  2.118393243...
Actividad: Verifique este resultado calculando 𝒇(𝟐)⋀𝒇(𝟑)

 Cálculo de 𝒙𝟑
Métodos Numéricos
sfs

66
(𝒙𝟐 ). 𝒇(𝒙𝟏 ) − 𝒙𝟏 . 𝒇(𝒙𝟐 )
𝒙𝟑 =
𝒇(𝒙𝟏 ) − 𝒇(𝒙𝟐 )
Reemplazando los valores correspondientes, se obtiene que:

𝒙𝟑 = 𝟐. 𝟐𝟏𝟏𝟎𝟗𝟒𝟑𝟕𝟏. ..
Actividad: Verifique este resultado realizando las operaciones con los valores indicados.

Comparando 𝒙𝟑  𝒙𝟐 se repiten tres decimales, el proceso continua.

 Cálculo de 𝒙𝟒

(𝒙𝟑 ). 𝒇(𝒙𝟐 ) − 𝒙𝟐 . 𝒇(𝒙𝟑 )


𝒙𝟒 =
𝒇(𝒙𝟐 ) − 𝒇(𝒙𝟑 )
𝒙𝟒 = 𝟐. 𝟓𝟖𝟕𝟐𝟑𝟔𝟔. ..
Actividad: Verifique este resultado realizando las operaciones con los valores indicados.

Comparando 𝒙𝟒  𝒙𝟑 no se repiten decimales, el proceso continua.


 Cálculo de 𝒙𝟓

(𝒙𝟒 ). 𝒇(𝒙𝟑 ) − 𝒙𝟑 . 𝒇(𝒙𝟒 )


𝒙𝟓 =
𝒇(𝒙𝟑 ) − 𝒇(𝒙𝟒 )

𝒙𝟓 = 𝟐. 𝟑𝟓𝟓𝟓𝟑𝟖𝟐. 𝟗..
Actividad: Verifique este resultado realizando las operaciones con los valores indicados.
Métodos Numéricos
sfs

67
Comparando 𝒙𝟓  𝒙𝟒 no se repiten decimales, el proceso continua.
 Cálculo de 𝒙𝟔

(𝒙𝟓 ). 𝒇(𝒙𝟒 ) − 𝒙𝟒 . 𝒇(𝒙𝟓 )


𝒙𝟔 =
𝒇(𝒙𝟒 ) − 𝒇(𝒙𝟓 )

𝒙𝟔 = 𝟐. 𝟑𝟗𝟖𝟑𝟒𝟎𝟖. 𝟓..

Actividad: Verifique este resultado realizando las operaciones con los valores indicados.

Comparando 𝒙𝟔  𝒙𝟓 se repite un decimal, el proceso continua.


 Cálculo de 𝒙𝟕

(𝒙𝟔 ). 𝒇(𝒙𝟓 ) − 𝒙𝟓 . 𝒇(𝒙𝟔 )


𝒙𝟕 =
𝒇(𝒙𝟓 ) − 𝒇(𝒙𝟔 )
𝒙𝟕 = 𝟐. 𝟒𝟏𝟒𝟕𝟕𝟗𝟗. 𝟐..
Actividad: Verifique este resultado realizando las operaciones con los valores indicados.

Comparando 𝒙𝟕  𝒙𝟔 No se repite ningún decimal, el proceso continua.


 Cálculo de 𝒙𝟖

(𝒙𝟕 ). 𝒇(𝒙𝟔 ) − 𝒙𝟔 . 𝒇(𝒙𝟕 )


𝒙𝟖 =
𝒇(𝒙𝟔 ) − 𝒇(𝒙𝟕 )

𝒙𝟖 =2.4132248.9..
Métodos Numéricos
sfs

68
Actividad: Verifique este resultado realizando las operaciones con los valores indicados.

Comparando 𝒙𝟖  𝒙𝟕 Se repiten dos decimales, el proceso continua.


 Cálculo de 𝒙𝟗

(𝒙𝟖 ). 𝒇(𝒙𝟕 ) − 𝒙𝟕 . 𝒇(𝒙𝟖 )


𝒙𝟗 =
𝒇(𝒙𝟕 ) − 𝒇(𝒙𝟖 )

𝒙𝟗 = 𝟐. 𝟒𝟏𝟑𝟐𝟔𝟐𝟔. 𝟐..
Actividad: Verifique este resultado realizando las operaciones con los valores indicados.

Comparando 𝒙𝟗  𝒙𝟖 Se repiten cuatro decimales, el proceso continua.


 Cálculo de 𝒙𝟏𝟎

(𝒙𝟗 ). 𝒇(𝒙𝟖 ) − 𝒙𝟖 . 𝒇(𝒙𝟗 )


𝒙𝟏𝟎 =
𝒇(𝒙𝟖 ) − 𝒇(𝒙𝟗 )

𝒙𝟏𝟎 =2.4132627.2..
Actividad: Verifique este resultado realizando las operaciones con los valores indicados.

Comparando 𝒙𝟏𝟎  𝒙𝟗 Se repiten seis decimales, el proceso finaliza.


La raíz es: 𝒙𝟏𝟎 =2.4132627.2..
2. Ejemplo tomado de (Burden & Faires, 2003). Utilice el método de la secante para encontrar una
solución de x  cosx.
Aplique el método hasta que cinco cifras decimales se repitan consecutivamente, es decir,
el error debe ser menor de 0.00001.
Métodos Numéricos
sfs

69

Procedimiento

Llevando la ecuación a la forma de una función queda:

f (x)  x  cosx

Al hacer dar valores a x se tiene que hay cambio de signo en x = 0.5 y x = 0.9, se pueden tomar:
𝒙𝟎 =  0.5  𝒙𝟏 =   / 4  0.78539816
Para 𝒏 = 𝟎 se tiene:

(𝒙𝟏 ). 𝒇(𝒙𝟎 ) − 𝒙𝟎 . 𝒇(𝒙𝟏 )


𝒙𝟐 =
𝒇(𝒙𝟎 ) − 𝒇(𝒙𝟏 )

𝒙𝟐  0.7363841

n xn

0 x0  0,5

1 x1  0,78539816

2 x2  0,73638414

3 x3  0,73905814

4 x4  0,73908515

5 x5  0,73908513
Métodos Numéricos
sfs

70
6 x6  0,73908513

7 x7  0,73908513

Se puede observar que después de la quinta iteración se obtiene la respuesta con la exactitud
deseada, por lo tanto, la raíz de esta ecuación es: x5  0.7390851.

1.4.6.1 MÉTODO DE ITERACIÓN DE PUNTO FIJO

 Video: Método del Punto Fijo


 ”Ver Video”: https://www.youtube.com/watch?v=mQogK-O4-o4

 Video: Ejercicio punto fijo


 ”Ver Video”: https://www.youtube.com/watch?v=LlRWk0WXPfo

 Video: Punto Fijo


 ”Ver Video”: https://www.youtube.com/watch?v=SeqMmyCUKDQ

Un punto fijo de una función es un número en el cual se cumple que x  y


Por ejemplo, la función:

y  f x  3x  4

Tiene un punto fijo en x 2 , ya que f  2  3 2 2  6  4  2 .


Según Purcell & Varberg (1993), Supóngase que una ecuación
(x)  0 que nos interesa resolver puede escribirse de la forma x  g(x) . Resolver una ecuación
de la forma
Métodos Numéricos
sfs

71
x  g(x) consiste en encontrar una raíz r que no cambie mediante la función g. Al número r se
le llama punto fijo.

Procedimiento para aplicar el método PUNTO FIJO:

1. Determine intervalos para cada raíz de la ecuación f (x)  0


Recuerde que se sabe que hay una raíz en un intervalo a, b porque:

f (x) Es continua en el intervalo a, b f (b)  f a 0.

𝑎+𝑏
2. Halle la semilla para la primera raíz como: 𝑥1 = 2
3. Lleve la ecuación 𝒇(𝒙) = 𝟎 𝑎 𝑙𝑎 𝑓𝑜𝑟𝑚𝑎 𝒙 = 𝒈(𝒙), es decir despeje 𝒙. Tenga presente que
pueden existir muchas formas de obtener 𝒙 = 𝒈(𝒙).

4. Luego obtenga x2  g x1  , x3  g x2  y en general xn1  g xn , es posible que


xn converja hacia la raíz r cuando n  .

5. La ilustración del método se puede observar en los siguientes ejemplos, es de aclarar


que el método puede fallar, esto depende de la elección de x  g(x).

1.4.7 EJERCICIOS DE APRENDIZAJE


1. Ejemplo propuesto por Purcell & Varberg (1993), Resuelva la ecuación 𝒙 − 𝒆−𝒙 = 𝟎, de tal
manera que los 4 primeros dígitos coincidan dos veces consecutivas. Utilice el método del punto
fijo.
Procedimiento:

Primero se lleva la ecuación a la forma x  g(x) , queda: x  e x con g(x)  e x

Dando valores a x en la ecuación original, es decir, en ecuación tiene una raíz en el intervalo 0, 1, ya
que:

f (x)  x  e x Es continua en todos los reales y f (1)=+0.632 ∧ f (0)=−𝟏

 Cálculo de la semilla:
Métodos Numéricos
sfs

72
𝑎+𝑏 1+0 1
𝑥1 = → 𝑥1 = → 𝑥1 = → 𝑥1 = 0.5
2 2 2
Aplicando el algoritmo, resulta que:

x2  g(x1 )  x2 g0.5 e0.5  0.60653065

x3  g(x2 )  x3 g0.606530659 e0.606530659 0.54523921

𝑥4  g(𝑥3 )  𝑥4  g0.545239211 e0.545239211 0.57970309


𝑥5  g(𝑥4 )  𝑥5  g0.579703094  0.56006462
𝑥6  g(𝑥5 )  𝑥6 = g 0.560064627  e0.560064627 0.57117214
𝑥7 = g(𝑥6 )  𝑥7 = g 0.571172149  e 0.571172149 0.56486294

𝑥8 = g(𝑥7 )  𝑥8 = g 0.564862949  e 0.564862947 0.56843804


𝑥9 = g(𝑥8 )  𝑥9 = g 0.568438047  e 0.568438047 0.56640945
𝑥10 = g(𝑥9 )  𝑥10 = g 0.566409452   e 0.566409452 0.56755963
𝑥11 = g(𝑥10 )  𝑥11 = g0.567559634e0.567559634 0.56690721
𝑥12 = g(𝑥11 )  𝑥12 = g0.566907212 e0.566907212  0.56727719

𝑥13 = g(𝑥12 )  𝑥13 = g 0.567277196)= e0.567277196  0.5670673

𝑥14 = g(𝑥13 )  𝑥14 = g 0.56706735) e0.56706735  0.5671863

𝑥15 = g(𝑥14 )  𝑥15 = g 0.56718636) e0.56718636 0.56711886

raíz  x15  0.567118864 0.5671


2. Ejemplo propuesto por Purcell & Varberg (1993), Resuelva 𝒙 = 𝟐 𝒄𝒐𝒔𝒙 con una
precisión de cuatro cifras decimales.
Procedimiento:

La ecuación tiene una raíz en el intervalo 0.5, 1.5


La semilla para esta raíz es:

𝑎+𝑏 1.5 + 0.5


𝑥1 = → 𝑥1 = → 𝑥1 = 2 → 𝑥1 = 1
2 2
Métodos Numéricos
sfs

73
Luego se lleva la ecuación a la forma x  g(x)
Una posible forma es x  2 cos x pero al calcular los siguientes valores de x se observa una
oscilación, lo que indica que no hay convergencia, por lo tanto con esta forma de:
x  gx  2cosx, no se obtiene una solución a la ecuación:
x  2 cos x  0.
Se debe escribir la ecuación de otra forma:

𝑥+2𝑐𝑜𝑠𝑥
Puede ser: 𝑔(𝑥) = 2

De esta manera se obtiene:

1+2𝑐𝑜𝑠1
x2  g(x1 ) x2𝑔(1) = = 1.040302306
2
1.040302306 +2𝑐𝑜𝑠1.040302306
x3  g(x2 ) x3= = 1.026110667
2
1.026110667+2𝑐𝑜𝑠1.026110667
x4  g(x3 ) x4= =
2
1.031204594+2𝑐𝑜𝑠1.031204594
x5  g(x4 ) x5= =
2
1.029388072+2𝑐𝑜𝑠1.029388072
x6  g(x5 ) x6= =
2
1.03003739+2𝑐𝑜𝑠1.03003739
x7  g(x6 ) x7= =
2
1.029805485+2𝑐𝑜𝑠1
x8  g(x7 ) x8= = 1.02988833.5
2

raíz  x8  1.02988833.5
 CONVERGENCIA DEL MÉTODO

Prueba de convergencia.

La g(x) que usted encontró se puede utilizar si y sólo si g'(x)  1 para todo número en a, b
Actividad: comprobar la convergencia en los dos ejemplos anteriores.

Taller de Entrenamiento
1. Utilizando al menos dos métodos diferentes, con una precisión de 10-4, encuentre todas
Métodos Numéricos
sfs

74
las raíces de las siguientes ecuaciones:
a. 5cosx  7x 2  4  0

b. 4x 4 10x3  5x 2  9x 15  0

c. ln(3x)  5x2  4  0

2. Use el método de punto fijo para encontrar la raíz de x3  3x 1  0 en 1,2 con


una precisión de tres cifras decimales.

3. ANUALIDADES
Al comprar a crédito, en muchos casos el pago se efectúa mediante anualidades; Las
ANUALIDADES, son cuotas de dinero periódicas e iguales que se entregan o se reciben al
comienzo o al final de cada período de tiempo. Para calcular el valor presente de las
anualidades se utiliza la fórmula:

(𝟏 + 𝒊)𝒏 − 𝟏
𝒑 = 𝑨[ ]
𝒊(𝟏 + 𝒊)𝒏
Donde:

A. El valor de la cuota fija periódica o ANUALIDAD.

p : El valor presente, es el dinero que se tendría que pagar de contado, es decir, sin financiación.

n : Es el periodo, es decir el número de anualidades a pagar.

i : Es la tasa de interés.

NOTA: i y n deben tener las mismas unidades de tiempo.

Con una precisión de 10-4 y utilizando el método que desee solucione: Un automóvil de

contado cuesta $ 25’000 000, pero se puede pagar de la siguiente manera: Una cuota

inicial de $5’000 000 y el saldo mediante cuotas mensuales de $

450.000 durante 5 años. ¿Cuál es la tasa de interés mensual pagada?


Métodos Numéricos
sfs

75

4. Utilizando el método que desee, con una precisión de 10-5 solucione:

Encuentre el número 𝒙 garantizado por el teorema del valor medio para la función

𝒙
𝒇(𝒙) = 𝟒 𝒆𝒏 𝒆𝒍 𝒊𝒏𝒕𝒆𝒓𝒗𝒂𝒍𝒐 [𝟎, 𝟏]
𝒙 +𝟏

5. En una fábrica la función de ingreso para cierto producto es:

rq 5000q  0.2q2 Unidades de dinero

La función de costo para el mismo producto es:

cq q3  2q  4 Unidades de dinero


Donde q está en millones de unidades. Encuentre el Punto de Equilibrio.

6. Un apartamento de 133 millones de pesos se compra de la siguiente manera: Una cuota


inicial de 78 millones de pesos, el saldo se financia a 15 años pagando cuotas mensuales
$ 660.000 Determine la tasa de interés de la financiación.

7. A partir de una lámina de cartón de 15 cm por 10 cm se desea construir una caja, para ello
se debe cortar cuadrados idénticos en las cuatro esquinas y se doblan los lados hacia
arriba. La caja debe ser hecha para contener un volumen de 60 cm3, determine las
dimensiones de la caja. Seleccione aquellas dimensiones que permitan obtener la mayor
altura.
Métodos Numéricos

76
2 UNIDAD 2 REGRESIÓN E INTERPOLACIÓN

 Video: Mínimos Cuadrados


 ”Ver Video”: https://www.youtube.com/watch?v=efEsu_9hXmA

 Video: Interpolación polinomial por método diferencias divididas


 ”Ver Video”: https://www.youtube.com/watch?v=GTq_9bq1CsM

2.1.1 OBJETIVO GENERAL


Aplicar métodos para la interpolación numérica de conjuntos de datos y para la regresión
polinómica.

2.1.2 OBJETIVOS ESPECÍFICOS


 Utilizar el método de determinantes para la solución de sistemas lineales de n ecuaciones
con n incógnitas.

 Utilizar el método de los mínimos cuadrados para la regresión polinómica.

 Reconocer algunos métodos para la interpolación numérica de conjuntos de datos

2.2 TEMA 1 SOLUCIÓN DE SISTEMAS DE ECUACIONES LINEALES


UTILIZANDO DETERMINANTES
 Solución de sistemas de ecuaciones lineales utilizando determinantes

Definición de Determinante

Un determinante es un número que se asocia a toda matriz cuadrada.


Métodos Numéricos
sfs

77
Si A es una matriz cuadrada de orden n, su determinante se puede simbolizar como

det A o |𝐀|

Cálculo de determinantes

2.2.1.1 CASO 1: DETERMINANTE DE MATRICES 2X2


Para una matriz 2X2 de la forma:

𝑎 𝑏
𝐴=[ ]
𝑐 𝑑
Recuerde que el determinante se calcula realizando la diferencia entre el producto de los
elementos de la diagonal principal y el producto de los elementos de la diagonal secundaria, esto
es:

𝑑𝑒𝑡 𝐴 = 𝑎. 𝑑 − 𝑏. 𝑐
Ejercicios:

6 −2
1. 𝑑𝑒𝑡 = [ ] = (6). (1) − (−2). (−3) = 6 − 6 = 0
−3 1
2 −3
2. 𝑑𝑒𝑡 = [ ] = (2). (−5) − (−3). (4) = −10 + 12 = 2
4 −5

2.2.1.2 CASO 2 DETERMINANTE DE MATRICES 3X3

 Video: 15 - Determinantes: Definición y Propiedades - Matemática I - Instituto


ISIV
 ”Ver Video”: https://www.youtube.com/watch?v=jEDFLrx9JCw

 Video: Determinantes 3x3 por Sarrus 01


 ”Ver Video”: https://www.youtube.com/watch?v=wQNaEAM5D8A
Métodos Numéricos
sfs

78
Sea la matriz A de orden 3 × 3: Se tiene que para calcular su determinante se repiten las dos
primeras columnas (o las dos primeras filas) y luego se realiza la suma del producto de las
diagonales menos la suma del producto de las diagonales secundarias.

A través de un ejemplo revisaremos este concepto:

𝑎11 𝑎12 𝑎13


𝐴 = [𝑎21 𝑎22 𝑎23 ]
𝑎31 𝑎32 𝑎33
Repitiendo las dos primeras columnas se tiene el siguiente arreglo para calcular el determinante
de una matriz cuadrada de orden 3:

Suma del producto de las diagonales principales:

𝑎11 𝑎12 𝑎13 𝑎11 𝑎12


𝑑𝑒𝑡𝐴 = [𝑎21 𝑎22 𝑎23 ] 𝑎21 𝑎22
𝑎31 𝑎32 𝑎33 𝑎31 𝑎32
𝑎11 . 𝑎22 . 𝑎33 + 𝑎12 . 𝑎23 . 𝑎31 + 𝑎13 . 𝑎21 . 𝑎32
Suma del producto de las diagonales secundarias:

𝑎11 𝑎12 𝑎13 𝑎11 𝑎12


𝑑𝑒𝑡𝐴 = [𝑎21 𝑎22 𝑎23 ] 𝑎21 𝑎22
𝑎31 𝑎32 𝑎33 𝑎31 𝑎32
𝑎13 . 𝑎22 . 𝑎31 + 𝑎11 . 𝑎23 . 𝑎32 + 𝑎12 . 𝑎21 . 𝑎33

Realizando las diferencias indicadas, en forma general se tiene que:


Métodos Numéricos
sfs

79
(𝒂𝟏𝟏 . 𝒂𝟐𝟐 . 𝒂𝟑𝟑 + 𝒂𝟏𝟐 . 𝒂𝟐𝟑 . 𝒂𝟑𝟏 + 𝒂𝟏𝟑 . 𝒂𝟐𝟏 . 𝒂𝟑𝟐 ) − (𝒂𝟏𝟑 . 𝒂𝟐𝟐 . 𝒂𝟑𝟏 + 𝒂𝟏𝟏 . 𝒂𝟐𝟑 . 𝒂𝟑𝟐 +
𝒂𝟏𝟐 . 𝒂𝟐𝟏 . 𝒂𝟑𝟑 )

(Producto Diagonal Principal − (Producto Diagonal Secundaria)

Ejercicio de Aprendizaje

Calcule el determinante de la siguiente matriz:

2 −3 1
𝐴 = [4 0 −2]
3 −1 −3

Se repiten las dos primeras columnas:

2 −3 1 2 −3
𝑑𝑒𝑡𝐴 = [3 0 −2] 4 0
3 −1 −3 3 −1

Se realiza la suma de los productos de la matriz principal:

(2 × 0 × −3) + (−3 × −2 × 3) + (1 × 4 × −1) = 14


Se realiza la suma de los productos de la matriz secundaria:

(𝟏 × 𝟎 × 𝟑) + (𝟐 × −𝟐 × −𝟏) + (−𝟑 × 𝟒 × −𝟑) = 𝟒𝟎

𝑑𝑒𝑡𝐴 = 14 − 40 = −26
Métodos Numéricos
sfs

80
Actividad: calcula el determinante para la siguiente matriz, siguiendo el proceso anterior:

7 −2 5
𝐴 = [9 3 −2]
8 7 5
2.2.1.3 CASO 3: DETERMINANTE DE MATRICES 𝒏 × 𝒏. DESARROLLO POR
COFACTORES
Dada una matriz cuadrada 𝐴 de orden 𝒏 × 𝒏:

La matriz menor resulta de eliminar en la matriz 𝐴 la fila 𝑖 y la columna 𝑗.


Revisemos un ejemplo para una mejor comprensión (te puedes remitir al módulo de Álgebra Lineal ‘para
que realices un repaso más completo).
Dada la matriz:

7 −2 5
𝐴 = [9 3 −2]
8 7 5
Obtenga la matriz menor para cada una de las filas:

Para la primera fila: Se elimina su fila y su columna

7 −2 5
𝐴11 = [9 3 −2]
8 7 5
Quedaría:

3 −2
𝐴11 = [ ]
7 5
Para la segunda fila: Se elimina su fila y su respectiva columna:
Métodos Numéricos
sfs

81
7 −2 5
𝐴21 = [9 3 −2]
8 7 5
7 5
𝐴21 = [ ]
8 5
Para la tercera fila: Se elimina su fila y su respectiva columna:

7 −2 5
𝐴31 = [9 3 −2]
8 7 5
−2 5
𝐴31 = [ ]
3 −2
Nota: Para calcular el determinante por menores, se realiza el siguiente proceso:

 Se determinan los menores por fila o por columna.

 Se realiza el determinante que queda indicado y se multiplica por el respectivo elemento


(pivote)

 Los signos de los productos se colocan alternos empezando por el signo más.

En el ejemplo, se tiene:

3 −2 7 5 −2 5
7[ ] − 9[ ]+8 [ ]=
7 5 8 5 3 −2
𝟕(𝟑 × 𝟓 − (−𝟐) × 𝟕) − 𝟗(𝟕 × 𝟓 − 𝟓 × 𝟖) + 𝟖(−𝟐 × (−𝟐) − 𝟓 × 𝟑) =

𝟕(𝟏𝟓 + 𝟏𝟒) − 𝟗(𝟑𝟓 − 𝟒𝟎) + 𝟖(𝟒 − 𝟏𝟓) =

𝟕 × 𝟐𝟗 − 𝟗(−𝟓) + 𝟖(−𝟏𝟏) = 𝟐𝟎𝟑 + 𝟒𝟓 − 𝟖𝟖 = 𝟏𝟔𝟎


Métodos Numéricos
sfs

82
2.2.2 EJERCICIO DE ENTRENAMIENTO
1. Dada la siguiente matriz, encuentra el determinante por menores de los elementos de
cada una de las columnas y compara el resultado con el ejercicio anterior.

7 −2 5
𝐴 = [9 3 −2]
8 7 5
2. ¿Cómo se calcularía el determinante de una matriz de orden 4 × 4?Realiza un ejemplo.

COFACTOR: Cij

El cofactor es un número que se asocia a cada entrada de una matriz cuadrada 𝐴


CÁLCULO DE COFACTORES:
(http://www.youtube.com/watch?v=xIt4td_4OaY)

𝑪𝒊𝒋 = (−𝟏)𝒊+𝒋 𝑴𝒊𝒋

La única diferencia entre un menor y un cofactor es el factor 1


i j

2.2.3 EJERCICIO DE APRENDIZAJE

1. Para la siguiente matriz encuentre: 𝑪𝟐𝟑 ∧ 𝑪𝟏𝟏

𝟑 𝟓 𝟑
𝑨 = [𝟑 𝟎 𝟖]
𝟐 𝟐 𝟑
Métodos Numéricos
sfs

83
Procedimiento:

𝟑 𝟓
 𝑪𝟐𝟑 = (−𝟏)𝟐+𝟑 [ ] = (−𝟏)𝟓 [𝟑(𝟐) − 𝟐(𝟓)] =
𝟐 𝟐

−𝟏 (𝟔 − 𝟏𝟎) = −𝟏(−𝟒) = 𝟒

---------------------------------------------------------
𝟎 𝟖
 𝑪𝟏𝟏 = (−𝟏)𝟏+𝟏 [ ] = (−𝟏)𝟐 [𝟎(𝟑) − 𝟐(𝟖)] =
𝟐 𝟑

𝟏 (𝟎 − 𝟏𝟔) = 𝟏(−𝟏𝟔) = −𝟏𝟔

2. Para la siguiente matriz encuentre: 𝑪𝟐𝟑 ∧ 𝑪𝟏𝟐

𝟔 −𝟑 𝟓
𝑨 = [−𝟒 −𝟗 𝟖]
𝟓 𝟑 𝟕
Procedimiento:

𝟔 −𝟑
 𝑪𝟐𝟑 = (−𝟏)𝟐+𝟑 [ ] = (−𝟏)𝟓 [𝟔(𝟑) − 𝟓(−𝟑)] =
𝟓 𝟑

−𝟏 (𝟏𝟖 + 𝟏𝟓) = −𝟏(𝟑𝟑) = −𝟑𝟑


---------------------------------------------------------
−𝟒 −𝟖
 𝑪𝟏𝟐 = (−𝟏)−𝟏+𝟐 [ ] = (−𝟏)𝟑 [−𝟒(𝟕) − 𝟓(𝟖)] =
𝟓 𝟕

−𝟏 (−𝟐𝟖 − 𝟒𝟎) = −𝟏(−𝟔𝟖) = 𝟔𝟖


Métodos Numéricos
sfs

84
Para calcular determinantes de matrices n x n se utiliza el método: CÁLCULO DE DETERMINANTES
POR EXPANSIÓN POR COFACTORES que consiste en:

Para encontrar el determinante de cualquier matriz cuadrada A de orden 𝒏 se


selecciona cualquier renglón (o columna) de 𝑨 y se multiplica cada entrada del
renglón (o columna) por su cofactor. La suma de estos productos será el
determinante de A, llamado determinante de orden 𝒏.

2.2.4 EJERCICIO DE APRENDIZAJE


1. Encuentre el determinante de la siguiente matriz:

𝟐 −𝟏 𝟑
𝑨 = [𝟑 𝟎 −𝟓]
𝟐 𝟏 𝟏
Procedimiento

a. Ampliando las dos primeras columnas.

𝟐 −𝟏 𝟑 𝟐 −𝟏
[𝟑 𝟎 −𝟓] 𝟑 𝟎
𝟐 𝟏 𝟏 𝟐 𝟏
𝒅𝒆𝒕 𝑨 = 𝟐(𝟎)(𝟏) + (−𝟏)(−𝟓)(𝟐) + 𝟑(𝟑)(𝟏) −2031 52 131

𝒅𝒆𝒕 𝑨 = 𝟎 + 𝟏𝟎 + 𝟗 − 𝟎 + 𝟏𝟎 + 𝟑 → 𝒅𝒆𝒕 𝑨 = 𝟑𝟐

b. Utilizando cofactores en el primer renglón (fila):

det A  a11C11  a12C12 +a13C13

𝒅𝒆𝒕𝑨 = 𝒂𝟏𝟏 (−𝟏)𝟏+𝟏 𝑴𝟏𝟏 + 𝒂𝟏𝟐 (−𝟏)𝟏+𝟐 𝑴𝟏𝟐 + 𝒂𝟏𝟑 (−𝟏)𝟏+𝟑 𝑴𝟏𝟑
Métodos Numéricos
sfs

85

𝟎 −𝟓 𝟑 −𝟓 𝟑 𝟎
𝒅𝒆𝒕𝑨 = 𝒂𝟏𝟏 (−𝟏)𝟏+𝟏 [ ] + 𝒂𝟏𝟐 (−𝟏)𝟏+𝟐 [ ] + 𝒂𝟏𝟑 (−𝟏)𝟏+𝟑 [ ]
𝟏 𝟏 𝟐 𝟏 𝟐 𝟏

𝒅𝒆𝒕𝑨 = 𝟐(−𝟏)𝟐 × [(𝟎)(𝟏) − (𝟏)(−𝟓)] + (−𝟏)(−𝟏)𝟑 × [(𝟑)(𝟏) − (𝟐)𝟖 − 𝟓)] + 𝟑(−𝟏)𝟒 × [(𝟑)(𝟏) − (𝟐)(𝟎)]

𝒅𝒆𝒕𝑨 = 𝟐(𝟏) × 𝟓 + (−𝟏)(−𝟏) × 𝟏𝟑 + 𝟑(𝟏) × 𝟑

𝒅𝒆𝒕𝑨 = 𝟏𝟎 + 𝟏𝟑 + 𝟗 = 𝟑𝟐

c. Utilizando cofactores en la segunda columna:

det A  a12C12  a22C22 +a32C32

𝒅𝒆𝒕𝑨 = 𝒂𝟏𝟐 (−𝟏)𝟏+𝟏 𝑴𝟏𝟐 + 𝒂𝟐𝟐 (−𝟏)𝟏+𝟐 𝑴𝟐𝟐 + 𝒂𝟑𝟐 (−𝟏)𝟏+𝟑 𝑴𝟑𝟐

𝟑 −𝟓 𝟐 𝟑 𝟐 𝟑
𝒅𝒆𝒕𝑨 = (−𝟏)(−𝟏)𝟑 [ ] + 𝟎(−𝟏)𝟒 [ ] + 𝟏(−𝟏)𝟑+𝟐 [ ]
𝟐 𝟏 𝟐 𝟏 𝟑 −𝟓

𝒅𝒆𝒕𝑨 = (−𝟏)(−𝟏) × [(𝟑)(𝟏) − (𝟐)(−𝟓)] + 𝟎 + 𝟏(−𝟏) × [(𝟐)(−𝟓) − (𝟑)(𝟑)]

𝒅𝒆𝒕𝑨 = (−𝟏)(−𝟏) × [𝟏𝟑] + 𝟎 + 𝟏(−𝟏) × [−𝟏𝟗]

𝒅𝒆𝒕𝑨 = 𝟏𝟑 + 𝟎 + 𝟏𝟗 = 𝟑𝟐

2. Utilizando cofactores a lo largo del segundo renglón (fila) obtenga el determinante de la matriz:

12 −1 3
𝐴=[2 3 2]
0 −1 3
Métodos Numéricos
sfs

86
Det(A)  3(1) 21 1(3)  2*3 1(1) 22 12(3)  (10)3  1123 12* 2  (10)
(1)

det(𝐴) = −1
3. Utilizando cofactores, calcule el determinante de la matriz:

0 1 1
𝐴 = [2 3 2]
0 −1 3
Nota: Como en la primera columna hay dos entradas iguales a cero, vamos a calcular el determinante
utilizando cofactores en la primera columna

det A  2(1)12 1*3  (1)1 8

4. Calcule el determinante de la matriz:


2 0 0 1
𝐵 = [0 1 0 3]
0 0 1 2
1 2 3 0
Vamos a utilizar cofactores en el primer renglón (Fila):

1 0 3 0 0 1
1+1 4
det 𝐵 = 2(−1) × 𝑑𝑒𝑡 [0 1 2] + 1(−1) 𝑑𝑒𝑡 [1 0 3]
2 3 0 0 1 2

det(𝐵) = 2 × [−(3 × 1 × 2 + 2 × 3 × 1]

det(𝐵) = −24 + (−1) = −25


5. Determine el valor de x, de tal manera que el determinante de la matriz A sea igual a -105.
Métodos Numéricos
sfs

87
5 4 −3
𝐴=[ 2 𝑥 0]
−8 6 3
A  det A  5x3 40 8  326  8x 3 605 324 105

15x  36  24x  24  105  9x  50  105  9x  105  50

−45
−9𝑥 = −45 → 𝑥 = →𝑥=5
−9

2.2.5 TALLER DE ENTRENAMIENTO


Dados los siguientes enlaces, realice los ejercicios que allí se plantean, utilizando los métodos
vistos en el desarrollo del tema.

 Video: Determinante de una matriz de 4x4 con cofactores


 ”Ver Video”: https://www.youtube.com/watch?v=vAVcikH4Y5A

 Video: Determinante de una matriz 4x4 (Bachillerato)


 ”Ver Video”: https://www.youtube.com/watch?v=ZO0naBrmgj4

 Solución de sistemas de ecuaciones lineales utilizando determinantes o REGLA DE


CRAMER

Sea un sistema de 𝒏 ecuaciones lineales con 𝒏 incógnitas de la forma:


𝑎11 𝑥1 𝑎12 𝑥2 𝑎1𝑛 𝑥𝑛 = 𝑐1
𝑎21 𝑥1 𝑎22 𝑥2 𝑎2𝑛 𝑥𝑛 = 𝑐2
[ ⋮ ⋮ ⋮ ]
𝑎𝑛1 𝑥𝑛 𝑎𝑛2 𝑥2 𝑎𝑛𝑛 𝑥𝑛 = 𝑐𝑛
Métodos Numéricos
sfs

88

Si el det A, matriz de coeficientes, es diferente de cero, el sistema tiene única solución que se
obtiene como:

𝐷1 𝐷2 𝐷𝑛
𝑥1 = , 𝑥2 = , … 𝑥𝑛 =
𝐷 𝐷 𝐷
Dónde:

D: Es el determinante de la matriz de coeficientes.


D1: Es el determinante de la matriz que resulta de intercambiar en la matriz de coeficientes las
entradas de la primera columna por las entradas de la matriz B ó matriz del término independiente.

D2: Es el determinante de la matriz que resulta de intercambiar en la matriz de coeficientes las


entradas de la segunda columna por las entradas de la matriz B o matriz del término independiente.

D3: Es el determinante de la matriz que resulta de intercambiar en la matriz de coeficientes las


entradas de la tercera columna por las entradas de la matriz B o matriz del término independiente.
Dn: Es el determinante de la matriz que resulta de intercambiar en la matriz de coeficientes las
entradas de la n ésima columna por las entradas de la matriz B ó matriz del término independiente.

2.2.6 EJERCICIO DE APRENDIZAJE


1. Utilizando determinantes, resuelva el sistema:

2x  5 y  4z  13

3x  2 y  3z  30

4x  6 y  5z  4

Procedimiento

a. Para obtener el determinante D utilizamos la matriz de coeficientes:


Métodos Numéricos
sfs

89
2 −5 4
𝐷 = [3 2 3]
4 6 −5
b. Se repiten las dos primeras columnas:

2 −5 4 2 −5
𝐷 = [3 2 3 ]3 2
4 6 −5 4 6

D  22 5  534 436 424 632  53 5

D  151
c. Para calcular el determinante D1 se cambia las entradas de la primera columna en la matriz de
coeficientes por las entradas del término independiente de cada ecuación:

13 −5 4
𝐷1 = [ 30 2 3]
−4 6 −5
Recuerde: Para realizar el determinante se repiten las dos primeras columnas

𝐷1 =132 5  53 4 4306  424 6313  530 5

D1  302

d. Para calcular el determinante D2 se cambia las entradas de la segunda columna en la matriz de


coeficientes por las entradas del término independiente de cada ecuación:

2 13 4
𝑫𝟐 = [3 30 3]
4 −4 −5
Recuerde: Para realizar el determinante se repiten las dos primeras columnas
Métodos Numéricos
sfs

90
𝑫𝟐 =230 5 1334 43 4 4304  432 133 5

D2  453

e. Para calcular el determinante 𝑫𝟑 se cambia las entradas de la tercera columna en la matriz de


coeficientes por las entradas del término independiente de cada ecuación:

2 −5 13
𝐷 = [3 2 30 ]
4 6 −4
Recuerde: Para realizar el determinante se repiten las dos primeras columnas

= D3  22 4  53041336 4213 6302  43 5


D3  906
f. La solución del sistema es:

−302 −453 −906


𝑥= = 2, 𝑦 = = 3, 𝑧 = =6
−151 −151 −151
La solución del sistema es:

x  2, y  3  z  6

g. Prueba:

2x  5 y  4z  13 →22 53 46  13  4 15 24  13  13  13

3x  2 y  3z  30 →32 23 36  30  6  6  18  30  30  30 

4x  6 y  5z  4 →42 63 56  4  8 18 30  4  4  4


Métodos Numéricos
sfs

91
2.3 TALLER DE ENTRENAMIENTO
1. Siguiendo el procedimiento presentado en el ejercicio de aprendizaje, resuelva el
siguiente sistema, dándole la respectiva prueba para verificar los resultados obtenidos:

2x  y  z  0
4x  3 y  2z  2
2x  y  3z  0
Respuesta: x=-1/2, y=2, z=-1
2. Determine el valor de x, de tal manera que el determinante de la matriz sea igual a -1

−6 3 𝑥
𝐴 = [ 7 2 −5] 𝑅𝑒𝑠𝑝𝑢𝑒𝑠𝑡𝑎. 𝑥 = 4
5 9 −𝑥

3. Determine el valor de x, de tal manera que el determinante de la matriz sea igual a -141

𝑥 6 5
𝐴 = [4 𝑥2 2] 𝑅𝑒𝑠𝑝𝑢𝑒𝑠𝑡𝑎. 𝑥 = 3
6 𝑥 5

4. Determine el valor de x, de tal manera que el determinante de la matriz sea igual a 22

𝑥−3 9 3
𝐴=[ 7 8 𝑥 − 2] 𝑅𝑒𝑠𝑝𝑢𝑒𝑠𝑡𝑎. 𝑥 = 8
5 0 4
Métodos Numéricos
sfs

92
2.4 TEMA 2 REGRESIÓN POLINÓMICA
Regresión polinómica.

En esta sección se estudiará el siguiente problema:

Según Asmar (1999), suponga que existe una relación funcional y = f (x) entre dos cantidades
x e y, con f desconocida y se conoce los distintos valores x e y.

Se trata de encontrar la función y =f (x) a partir de un conjunto de valores x e y


Este problema se conoce como problema de ajuste de datos o ajuste de curvas. Interesa el caso en que
f es una función polinómica.
La idea es obtener una función de aproximación que se ajuste a la forma o a la tendencia general
de los datos, sin coincidir necesariamente en todos los puntos.

Suponga que se conocen los datos x0, y0 , x1, y1 ...,xn, yn distintos y se desea encontrar un
polinomio:

𝑃𝑚 (x)  𝑏0  𝑏1 x  𝑏2 x2  ... 𝑏𝑛 xn con m  n


Es decir, el grado del polinomio de regresión no puede ser mayor o igual que el número de parejas x y.
Se debe cumplir que la suma de los cuadrados del error cometido sea mínima.

El grado m del polinomio 𝑷𝒎 (x) puede ser el que cada quien deseé, teniendo en cuenta que
no puede superar a n  1 que es el número de parejas x e y .

A continuación, se explica la forma de encontrar los diferentes polinomios de regresión.

2.4.1.1 REGRESIÓN LINEAL

 Video: MinimosCuadrados
 ”Ver Video”: https://www.youtube.com/watch?v=TiLeQ7rDvSc
Métodos Numéricos
sfs

93
 Video: REGRESION LINEAL
 ”Ver Video”: https://www.youtube.com/watch?v=rF66WNFYAoM

 Video: Problema de regresión lineal


 ”Ver Video”: https://www.youtube.com/watch?v=DdAM9DXfrGI

 Video: REGRESION LINEAL


 ”Ver Video”: https://www.youtube.com/watch?v=zZvOocwYDkA

Consiste en encontrar una fórmula o ecuación lineal que muestre como se relacionan las variables

x e y, permitiendo efectuar estimaciones de y con valores de x propuestos o viceversa. La expresión


matemática para la línea recta es:

𝒚 = 𝒃𝟎 + 𝒃𝟏 𝒙
Donde:

𝒃𝟎 𝒚 𝒃𝟏 𝒔𝒐𝒏 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕𝒆𝒔
𝒃𝟏 : 𝑷𝒆𝒏𝒅𝒊𝒆𝒏𝒕𝒆 𝒅𝒆 𝒍𝒂 𝒍í𝒏𝒆𝒂 𝒓𝒆𝒄𝒕𝒂
𝒃𝟎 : 𝑰𝒏𝒕𝒆𝒓𝒔𝒆𝒄𝒄𝒊ó𝒏 𝒄𝒐𝒏 𝒆𝒍 𝒆𝒋𝒆 𝒚.
Para obtener los valores de b0 y b1 de esta ecuación lineal, y utilizando el método de los mínimos
cuadrados, se tiene que:

𝒆 = 𝒚 − 𝒃 𝟎 − 𝒃𝟏 𝒙
Calculando el error para cada valor, se tiene que:

e0  y0  b0  b1 x0 ,

e1  y1  b0  b1 x1 ,

e2  y2  b0  b1 x2 ,...,

en  yn  b0  b1 xn
Elevando al cuadrado cada lado de la ecuación y sumando, se tiene que:
Métodos Numéricos
sfs

94
𝑛 𝑛

∑(𝑒𝑖 )2 = ∑(𝑦𝑖 − 𝑏0 − 𝑏1 𝑥𝑖 )2
𝑖=1 𝑖=1
Lo que es lo mismo:

𝑺 = ∑(𝒚𝒊 − 𝒃𝟎 − 𝒃𝟏 𝒙𝒊 )𝟐

El objetivo es encontrar los valores de 𝒃𝟎 𝒚 𝒃𝟏


Para determinar los mínimos de esta función se debe derivar de manera parcial con respecto a cada
variable, queda entonces:

𝝏𝑺
= −𝟐 ∑(𝑦𝑖 − 𝑏0 − 𝑏1 𝑥𝑖 )
𝝏𝒃𝟎
𝝏𝑺
= −𝟐 ∑(𝑦𝑖 − 𝑏0 − 𝑏1 𝑥𝑖 )𝒙
𝝏𝒃𝟏
Igualando estas derivadas a cero, se encuentra el S mínimo:

→ ∑(𝑦𝑖 − ∑ 𝑏0 − ∑ 𝑏1 𝑥𝑖 ) = 𝟎

→ ∑ 𝑦𝑖 𝑥𝑖 − ∑ 𝑏0 𝑥𝑖 − ∑ 𝑏1 𝑥𝑖 𝟐 = 𝟎

Plantea un sistema 𝟐 × 𝟐 de la siguiente manera:

∑ 𝑦 = 𝑏0 𝑛 + 𝑏1 ∑ 𝑥


Métodos Numéricos
sfs

95

∑ 𝑥𝑦 = 𝑏0 ∑ 𝑥 + 𝑏1 ∑ 𝑥 2
Cuya solución es:

𝑛 ∑ 𝑥𝑦 − ∑ 𝑥 × ∑ 𝑦
𝑏1 =
𝑛 ∑ 𝑥 2 − (∑ 𝑥)2
𝑏0 = 𝑦̅ − 𝑏1 𝑥̅
EJERCICIOS DE APRENDIZAJE:
1. Este es parte de un ejemplo planteado por Asmar (1999). Para la tabla:

x y
0 -1
2 0
3 2
5 1

Determine la recta regresión.


SOLUCIÓN.

Para hallar las constantes b0 y b1 es conveniente que se llene la siguiente tabla:

x y xy
x2
0 -1 0 0
2 0 4 0
3 2 9 6
5 1 25 5

x 10 y  2 x 2
 38 x  y  11
Con n = 4, Ya que hay cuatro parejas x, y
Cálculo de 𝑥̅ :
Métodos Numéricos
sfs

96
∑ 𝑥 10 5
𝑥̅ = = =
𝑛 4 2
Cálculo de 𝑦̅:
∑𝑦 2 1
𝑦̅ = = =
𝑛 4 2
Cálculo de 𝑏1 :

𝑛 ∑ 𝑥𝑦 − ∑ 𝑥 × ∑ 𝑦
𝑏1 =
𝑛 ∑ 𝑥 2 − (∑ 𝑥)2

Reemplazando se tiene que:

𝑛 ∑ 𝑥𝑦 − ∑ 𝑥 × ∑ 𝑦
𝑏1 =
𝑛 ∑ 𝑥 2 − (∑ 𝑥)2
𝟒 ∗ 𝟏𝟏 − 𝟏𝟎 ∗ 𝟐 𝟐𝟒 𝟔
𝒃𝟏 = = =
𝟒 ∗ 𝟑𝟖 − (𝟏𝟎)𝟐 𝟓𝟐 𝟏𝟑
Cálculo de:

̅ − 𝒃𝟏 𝒙
𝒃𝟎 = 𝒚 ̅
𝟏 𝟔 𝟓 𝟏 𝟑𝟎 𝟏𝟑 − 𝟑𝟎
𝒃𝟎 = − × = − =
𝟐 𝟏𝟑 𝟐 𝟐 𝟐𝟔 𝟐𝟔
𝟏𝟕
𝒃𝟎 = − 𝟐𝟔
La recta regresión queda entonces:

17 6
𝑦= + 𝑥
26 13
Métodos Numéricos
sfs

97
2. Para los datos de la siguiente tabla determine la recta de regresión

x y
-5 15
-2 9
1 3
4 -3
6 -7
PROCEDIMIENTO:

Para encontrar la recta, es conveniente diligenciar la siguiente tabla:

X Y X2 X.Y
-5 15 25 -75
-2 9 4 -18
1 3 1 3
4 -3 16 -12
6 -7 36 -42
x  4  y  17  x 2
 82 x  y  144
Con n = 5, Ya que hay cinco parejas x, y

Cálculo de 𝑥̅ :
∑𝑥 4
𝑥̅ = =
𝑛 5
Cálculo de 𝑦̅:
∑ 𝑦 17
𝑦̅ = =
𝑛 5
Cálculo de 𝑏1 :

𝑛 ∑ 𝑥𝑦 − ∑ 𝑥 × ∑ 𝑦
𝑏1 =
𝑛 ∑ 𝑥 2 − (∑ 𝑥)2
Métodos Numéricos
sfs

98
𝟓 ∗ (−𝟏𝟒𝟒) − 𝟒 ∗ 𝟏𝟕 𝟕𝟖𝟖
𝒃𝟏 = = − = −𝟐
𝟓 ∗ 𝟖𝟐 − (𝟒)𝟐 𝟑𝟗𝟒
Cálculo de:

̅ − 𝒃𝟏 𝒙
𝒃𝟎 = 𝒚 ̅
𝟏𝟕 𝟒 𝟏𝟕 𝟖 𝟐𝟓
𝒃𝟎 = − (−𝟐) × = + = =𝟓
𝟓 𝟓 𝟓 𝟓 𝟓
𝒃𝟎 = 𝟓
La recta regresión queda entonces:

𝑦 = 5 − 2𝑥
REGRESIÓN CUADRÁTICA O REGRESIÓN PARABÓLICA

La parábola que se debe obtener tiene la siguiente forma:

𝒚  𝒃𝟎  𝒃𝟏 x  𝒃𝟐 x 2
Para encontrar los valores de 𝒃𝟎, 𝒃𝟏 , 𝒃𝟐 el método de los mínimos cuadrados plantealas siguientes
ecuaciones:

2
 y  𝒃𝟎 n  𝒃𝟏  x  𝒃𝟐  x
2 3
 xy 𝒃𝟎  x  𝒃𝟏  x  𝒃𝟐  x
0 1 2

2 2 3 4
 x y 𝒃𝟎  x  𝒃𝟏  x  𝒃𝟐  x
Al solucionar este sistema 3X3 se obtienen las constantes b0 , b1  b2
Métodos Numéricos
sfs

99
2.4.2 EJERCICIOS DE APRENDIZAJE
1. Para los datos de la tabla:
x y
-5 113
-2 23
1 5
4 59
6 135
10 383

Determine la parábola de regresión.


Procedimiento:

Para hallar las constantes b0 , b1  b2 es conveniente que llenemos la siguiente tabla:

X Y X2 X.Y x2.y x3 x4
-5 113 25 -565 2825 -125 625
-2 23 4 -46 92 -8 16
1 5 1 5 5 1 1
4 59 16 236 944 64 256
6 135 36 810 4860 216 1296
10 383 100 3830 38300 1000 10000
2 2 3 4
x  14 y  718 x  182 x  y  4270 x  y  47026 x  1148 x  12194

Con n6
Se debe plantear y resolver el siguiente sistema:

718 = 6𝑏0 + 14𝑏1 + 182𝑏2


4270 = 14𝑏0 + 182𝑏1 + 1148𝑏2
47026 = 182𝑏0 + 1148𝑏1 + 12194𝑏2
El sistema dado se resuelve por determinantes:

6 14 182
𝐷 = [ 14 182 1148 ] = 2.840.040
182 1148 −12194
Métodos Numéricos
sfs

100
718 14 182
𝐷1 = [ 4270 182 1148 ] = 8.520.120
47026 1148 −12194
6 718 182
𝐷2 = [ 14 4270 1148 ] = −5.680.080
182 47026 −12194
6 14 718
𝐷3 = [ 14 182 4270 ] = 11.360.160
182 1148 47026

Recuerde que:

𝐷
𝑏𝑛 =
𝐷𝑛
Reemplazando los valores obtenidos, se tiene que:

𝐷1 8.520.120
𝑏0 = = =3
𝐷 2.840.040
𝐷2 −5.680.080
𝑏1 = = = −2
𝐷 2.840.040
𝐷3 11.360.160
𝑏2 = = =4
𝐷 2.840.040
La parábola de regresión está dada por:

𝒚 = 𝟑 − 𝟐𝒙 − 𝟒𝒙𝟐

2. Encuentre la parábola de regresión que mejor aproxime los puntos:

 4,10, 1,2, 3,10, 5,28 8,70


Métodos Numéricos
sfs

101
Procedimiento:

Es conveniente llenar la siguiente tabla:

X Y X2 X.Y x2y x3 x4
-4 10 16 -40 160 -64 256
-1 -2 1 2 -2 -1 1
3 10 9 30 90 27 81
5 28 25 140 700 125 625
8 70 64 560 4480 512 4096
 y  11 6  x x  y  69 2  x  y  542 8 x  599 x
2 3 4
 115  5059
2

x  11

Se debe plantear y resolver el siguiente sistema:

116 = 5𝑏0 + 11𝑏1 + 115𝑏2


692 = 11𝑏0 + 115𝑏1 + 599𝑏2
5428 = 115𝑏0 + 599𝑏1 + 5059𝑏2
El sistema dado se resuelve por determinantes:

5 11 115
𝐷 = [ 11 115 599 ] = 497.376
115 599 5059
116 11 115
𝐷1 = [ 692 115 599 ] = −994.752
5428 599 5059
5 116 115
𝐷2 = [ 11 692 599 ] = 497.376
115 5428 5059
5 11 116
𝐷3 = [ 11 115 692 ] = 497.376
115 599 5248
Métodos Numéricos
sfs

102
Recuerde que:

𝐷
𝑏𝑛 =
𝐷𝑛
Reemplazando los valores obtenidos, se tiene que:

𝐷1 −994.752
𝑏0 = = = −2
𝐷 497.376
𝐷2 497.376
𝑏1 = = =1
𝐷 497.376
𝐷3 497.376
𝑏2 = = =1
𝐷 497.376
La parábola de regresión es:

y  2  x  x2
Polinomio de regresión de grado 3 o polinomio cúbico de regresión

El polinomio de grado 3 o polinomio cúbico que se debe obtener tiene la siguiente forma:

𝒚  𝒃𝟎  𝒃𝟏 x  𝒃𝟐 x2+𝒃𝟑 x3
Para encontrar las constantes b0 , b1 , b2  b3 el método de los mínimos cuadrados plantea las
siguientes ecuaciones:

2 3
 y  𝒃𝟎 n 𝒃𝟏  x  𝒃𝟐  x  𝒃𝟑  x
0 1 2 3

2 3 4
 xy  𝒃𝟎  x 𝒃𝟏  x  𝒃𝟐  x  𝒃𝟑  x
0 1 2 3

2 2 3 4 5
 x y  𝒃𝟎  x  𝒃𝟏  x  𝒃𝟐  x  𝒃𝟑  x
0 1 2 3

3 3 4 5 6
 x y  𝒃𝟎  x  𝒃𝟏  x 𝒃𝟐  x 𝒃𝟑  x
0 1 2 3
Métodos Numéricos
sfs

103
Al solucionar este sistema 4X4 se obtienen las constantes b0 , b1 , b 2  b3 .
Polinomio de regresión de grado n.

El polinomio de grado 3 o polinomio cúbico que se debe obtener tiene la siguiente forma:

𝒚  𝒃𝟎  𝒃𝟏 x  𝒃𝟐 x2+𝒃𝟑 x3+ ⋯ + 𝒃𝒏 xn
Para encontrar las constantes b0 , b1, b2 , b3 .... bn el método de los mínimos cuadrados plantea
las siguientes ecuaciones:

3 n
 y  b0 n  b1  x  b2  x  b3  x ... bn  x
0 1 2 3 n

2 3 4 n1
 xy  b0  x  b1  x  b2  x  b3 x  ... b  x
0 1 2 3 n

2 2 3 4 5 n2
 x y  b0  x  b1  x  b2 x  b3  x  ... bn  x
0 1 2 3 n

3 3 5 6 n3
 x y  b0  x  b1  x  b  x  b3  x  ... bn  x
0 1 2 3 n

En forma general:

n n+ n+2  b n+3  b n+4  ... b nn


 x y  b0  x  b1  x 2 x 3 x n x
0 1 2 3 n

Al solucionar este sistema nXn se obtiene las constantes b0, b1 , b2 , b3 .... bn

2.4.3 TALLER DE ENTRENAMIENTO


1. Encuentre la recta de regresión para el siguiente conjunto de datos:

x -2 3 4 9
y 7 4 -2 -3
Métodos Numéricos
sfs

104
2. Encuentre la parábola de regresión para el siguiente conjunto de datos.

x -3 -1 1 2 4 5
y 9 2 1 3 7 9

3. Encuentre el polinomio de regresión de grado 3 para el siguiente conjunto de datos.

x -9 -6 -4 -1 1 2 3 5 7 10
y -921 -303 -101 7 19 25 39 115 295 865

2.5 TEMA 3 INTERPOLACIÓN POLINOMIAL


Interpolación polinomial

 Este tema consiste en dado un conjunto de valores x e y encontrar una función polinómica
que represente estos valores, que pase exactamente por cada uno de ellos y que permita
estimar valores no observados (es decir, que permita hacer predicciones o interpolaciones o
aproximaciones).
También permite aproximar mediante una función polinómica funciones que no son polinómicas. Se
utiliza la aproximación polinomial ya que los polinomios son continuos en todo su dominio, con esto
se garantizan funciones continuas, que son las que interesa estudiar en el análisis numérico y además
la derivada y la integral de polinomios son muy fáciles de determinar y también son polinómicas.

La forma general de un polinomio es:

𝑷𝒏 (𝒙) = 𝒂𝒏 𝒙𝒏 + 𝒂𝒏−𝟏 𝒙𝒏−𝟏 + ⋯ 𝒂𝟏 𝒙 + 𝒂𝟎


Donde n es un entero no negativo que indica el grado del polinomio y 𝒂𝟎 , 𝒂𝟏 …,𝒂𝒏 Son
constantes reales.
Para determinar el polinomio interpolante, se desarrollarán dos métodos.

 Video: Interpolación 1b
 ”Ver Video”: https://www.youtube.com/watch?v=S3QEyUN3sn0
Métodos Numéricos
sfs

105
Polinomio interpolante de LaGrange

 Video: Interpolación 2
 ”Ver Video”: https://www.youtube.com/watch?v=h5s6lZ9EDbk
 Video: Calculo y métodos Numéricos. Interpolación de Lagrange
 ”Ver Video”: https://www.youtube.com/watch?v=S4KdYhTzaxE

 CASO 1: POLINOMIO DE PRIMER GRADO

La idea es encontrar un polinomio de primer grado que pasa por dos puntos distintos:

x0 , y0   x1 , y1 
O lo que es lo mismo aproximar una función f, para la cual:
f (x0 )  y0 y f (x1 )  y1

Por medio de un polinomio de primer grado que interpole los valores de f en los puntos dados o que
coincida con ellos.

El polinomio interpolante de primer grado es:

P1 (x)  a0  a1 x

Se tiene que:

𝑃1 (𝑥0 ) = 𝑎0 + 𝑎1 𝑥0 = 𝑦0
𝑃1 (𝑥1 ) = 𝑎0 + 𝑎1 𝑥1 = 𝑦1
Se tiene un sistema 2X2, donde las incógnitas son a0  a1
Solución del sistema:

1 𝑥0
𝐷=[ ] = 𝑥1 − 𝑥0
1 𝑥1

𝑦0 𝑥0
𝐷1 = [𝑦 𝑥1 ] = 𝑥1 𝑦0 − 𝑥0 𝑦1
1
Métodos Numéricos
sfs

106
1 𝑦0
𝐷2 = [ ] = 𝑦1 − 𝑦0
1 𝑦1
𝐷1 𝑥1 𝑦0 −𝑥0 𝑦1 𝐷 𝑦 −𝑦
𝑎0 = = , 𝑎1 = 𝐷2 = 𝑥1−𝑥0
𝐷 𝑥1 −𝑥0 1 0

Reemplazando 𝑎0 y 𝑎1 en: P1 (x)  a0  a1 x se tiene:

𝑥1 𝑦0 −𝑥0 𝑦1 𝑦1 −𝑦0
P1 (x)  +( )𝑥
𝑥1 −𝑥0 𝑥1 −𝑥0

Realizando la operación indicada:


𝑥1 𝑦0 −𝑥0 𝑦1 +𝑥𝑦1 −𝑥𝑦0
P1 (x) 
𝑥1 −𝑥0

Factorizando:

𝑦0 (𝑥1 −𝑥)+𝑦1 (𝑥−𝑥0 ) 𝑥1 −𝑥 𝑥−𝑥0


P1 (x)  → P1 (x)= 𝑦0 ( ) + 𝑦1 ( )→
𝑥1 −𝑥0 𝑥1 −𝑥0 𝑥1 −𝑥0

𝒙−𝒙𝟏 𝒙−𝒙𝟎
P1 (x)= ( ) 𝒚𝟎 + ( ) 𝒚𝟏
𝒙𝟎 −𝒙𝟏 𝒙𝟏 −𝒙𝟎

Para ello se definen las funciones:

𝒙−𝒙𝟏 𝒙−𝒙𝟎
L0 (x)= ( ) ∧ 𝑳𝟏 (𝒙) ( )
𝒙𝟎 −𝒙𝟏 𝒙𝟏 −𝒙𝟎

Y se define:

P1 (x)= L0 (x)𝒇(𝒙𝟎 ) + 𝑳𝟏 (𝒙)𝒇(𝒙𝟏 )

Se puede afirmar que P(x) es la única función lineal que pasa por los puntos x0 , y0   x1 , y1  .

2.5.1 TALLER DE APRENDIZAJE


1. Encuentre el polinomio lineal que pasa por los puntos (2,5)  (8,3)
Métodos Numéricos
sfs

107
Procedimiento

Datos del problema:

x0  2, y0  5  x1  8, y1  3



El polinomio o recta de interpolación tiene la forma:

P1 (x)  L0 (x) f (x0 )  L1(x) f (x1 )

𝒙−𝒙𝟏 𝒙−𝟖 𝒙−𝟖 𝑥 8


Cálculo de: L0 (x)= ( )=( ) → L0 (x)= → L0 (x)=−6 − −6
𝒙𝟎 −𝒙𝟏 𝟐−𝟖 −𝟔

𝑥 8 𝟒 𝒙
L0 (x)=− + → L0 (x)= −
6 6 𝟑 𝟔

Cálculo de:

𝒙 − 𝒙𝟎 𝒙−𝟐 𝑥 2
𝑳𝟏 ( 𝒙 ) = ( ) → 𝑳𝟏 ( 𝒙 ) = ( ) → 𝑳𝟏 ( 𝒙 ) = − →
𝒙𝟏 − 𝒙𝟎 𝟖−𝟐 6 6
𝒙 𝟏
L1 (x)=𝟔 − 𝟑
Reemplazando en: P1 (x)  L0 (x) f (x0 )  L1 (x) f (x1 )

𝒙 𝟏𝟕
𝟒 𝒙 𝒙 𝟏
P1 (x) ( − ) 𝟓 + ( − ) 𝟑 →P1 (x)
𝟑 𝟔 𝟔 𝟑
− +
𝟑 𝟑

2. Encuentre el polinomio lineal que pasa por los puntos (1,7)  (6,27)

3. Procedimiento

Datos del problema:


Métodos Numéricos
sfs

108

x0  1, y0  7  x1  6, y1  27



El polinomio o recta de interpolación tiene la forma:

P1 (x)  L0 (x) f (x0 )  L1(x) f (x1 )

𝒙−𝒙𝟏 𝒙−𝟔 𝒙−𝟔 −𝒙+𝟔


Cálculo de: L0 (x)= ( )=( ) → L0 (x)= → L0 (x)=
𝒙𝟎 −𝒙𝟏 𝟏−𝟔 −𝟓 𝟓

Cálculo de:

𝒙 − 𝒙𝟎 𝒙−𝟏 𝒙−𝟏
𝑳𝟏 ( 𝒙 ) = ( ) → 𝑳𝟏 ( 𝒙 ) = ( ) → 𝑳𝟏 ( 𝒙 ) = ( )
𝒙𝟏 − 𝒙𝟎 𝟔−𝟏 𝟓
𝒙 𝟏
L1 (x)=𝟔 − 𝟑
Reemplazando en: P1 (x)  L0 (x) f (x0 )  L1 (x) f (x1 )

𝒙 𝟏𝟕
P1 (x) (
−𝒙+𝟔
𝟓
𝒙−𝟏
) 𝟕 + ( 𝟓 ) 𝟐𝟕 →P1(x) − +
𝟑 𝟑
−7𝑥+42 27𝑥−27
P1 (x) 
5
+
5
→P1 (x)=−7𝑥+42+27𝑥−27
5
=
20𝑥+15
5

Dividiendo por 5, se tiene que: :

20𝑥 15
P1 (x)= + = 4𝑥 + 3
5 5

 CASO 2: POLINOMIO DE GRADO MÁXIMO 2 QUE PASE POR TRES PUNTOS


(http://www.youtube.com/watch?v=Tk_ZVMhHAmY)

Los puntos por los cuales debe pasar el polinomio son: x0 , y0 , x1 , y1 , x2 , y2 
Métodos Numéricos
sfs

109
Recuerde que x0 , x1 , x2 deben ser diferentes.

Se desea encontrar el polinomio interpolante de grado menor o igual que dos, el polinomio se simboliza
como P2 (x) y tiene la forma:

𝑃2 (𝑥 ) = 𝑎0 + 𝑎1 𝑥 + 𝑎2 𝑥 2
Se debe encontrar el valor de las constantes: a0 , a1  a2

Deducción del polinomio de LaGrange de grado 2 (Esta demostración puede ser omitida)

Como los puntos x0 , y0 , x1 , y1 , x2 , y2  pertenecen al polinomio P2 (x) , se debe cumplir que:

𝒚𝟎 = 𝒂𝟎 + 𝒂𝟏 𝒙𝟎 + 𝒂𝟐 𝒙𝟎 𝟐

𝒚𝟏 = 𝒂𝟎 + 𝒂𝟏 𝒙𝟏 + 𝒂𝟐 𝒙𝟏 𝟐

𝒚𝟐 = 𝒂𝟎 + 𝒂𝟏 𝒙𝟐 + 𝒂𝟐 𝒙𝟐 𝟐

Este es un sistema 3 X 3 donde las incógnitas son a0 , a1  a2

SOLUCIÓN DEL SISTEMA

Se soluciona este sistema matricialmente utilizando determinantes.


La matriz de coeficientes A, la matriz de variables X y la matriz de términos independientes son en su
orden:

1 𝑥0 𝑥0 2 𝑎0 𝑦0
𝐴 = [1 𝑥1 𝑥1 2 ] , 𝑥 = [𝑎1 ] ∧ 𝐵 = [𝑦1 ]
1 𝑥2 𝑥2 2 𝑎2 𝑦2
Sea D el det(A)
Para simplificar el cálculo de del determinante D se aplican operaciones elementales en la matriz A:
Métodos Numéricos
sfs

110
R1  R1 , R2  R2  R1 (1)  R3  R3  R1 (1)

1 𝑥0 𝑥0 2
[0 𝑥1 − 𝑥0 𝑥1 2 −𝑥0 2 ]
0 𝑥2 − 𝑥0 𝑥2 2 −𝑥0 2
Se obtiene D por expansión por cofactores a lo largo de la primera columna:

𝑥1 − 𝑥0 𝑥1 2 −𝑥0 2 2 2
𝐷=[ 2 2 ]2 = (𝑥1 − 𝑥0 ) (𝑥2 −𝑥0 ) − (𝑥2 −
𝑥2 − 𝑥0 𝑥2 −𝑥0
𝑥0 )(𝑥1 2 −𝑥0 2 )

D  x1  x0 x2  x0 x2  x0   x1  x0   x1  x0 x2  x0 x2  x0  x1  x0 

𝑫 = (𝒙𝟏 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )


𝑫𝟏 : Es el determinante de la matriz que resulta al intercambiar en la matriz A los elementos de
la primera columna por los elementos del vector:

𝑦0 𝑥0 𝑥0 2
𝑫𝟏 = [𝑦1 𝑥1 𝑥1 2 ]
𝑦2 𝑥2 𝑥2 2
Se obtiene este determinante por expansión a lo largo de la primera columna:

1+1 𝑥1 𝑥1 2 2+1 𝑥0 𝑥0
2
𝑫𝟏 = 𝑦0 (−1) [ ] + 𝑦1 (−1) [ ]
𝑥2 𝑥2 2 𝑥2 𝑥2 2
2
3+1 𝑥0 𝑥0
+ 𝑦2 (−1) [ ]
𝑥1 𝑥1 2
𝑫𝟏 = 𝑦0 (𝑥1 𝑥2 2 − 𝑥1 2 𝑥2 ) − 𝑦1 (𝑥0 𝑥2 2 − 𝑥0 2 𝑥2 ) + 𝑦2 (𝑥0 𝑥1 2 − 𝑥0 2 𝑥1 )

𝑫𝟏 = 𝑦0 𝑥1 𝑥2 (𝑥2 − 𝑥1 ) − 𝑦1 𝑥0 𝑥2 (𝑥2 − 𝑥0 ) + 𝑦2 𝑥0 𝑥1 (𝑥1 − 𝑥0 )

𝑫𝟐 : Es el determinante de la matriz que resulta al intercambiaren la matriz A los elementos de


la segunda columna por los elementos del vector B.
Métodos Numéricos
sfs

111
𝑦0 𝑥0 2
𝟏
𝑫𝟐 = [𝟏 𝑦1 𝑥1 2 ]
𝟏 𝑦2 𝑥2 2
Se soluciona este determinante por expansión a lo largo de la segunda columna:

1+2 1 𝑥1 2 2+2 1 𝑥0
2
3+2 1 𝑥0
2
𝑫𝟐 = 𝑦0 (−1) [ ] + 𝑦1 (−1) [ ] + 𝑦2 (−1) [ ]
1 𝑥2 2 1 𝑥2 2 1 𝑥1 2
𝑫𝟐 = −𝑦0 (𝑥2 2 − 𝑥1 2 ) + 𝑦1 (𝑥2 2 − 𝑥0 2 ) − 𝑦2 (𝑥1 2 − 𝑥0 2 )

𝑫𝟐 = 𝑦0 (𝑥2 − 𝑥1 )(𝑥2 + 𝑥1 ) − 𝑦1 (𝑥2 − 𝑥0 )(𝑥2 + 𝑥0 )


+ 𝑦2 (𝑥1 − 𝑥0 )(𝑥1 + 𝑥0 )
𝑫𝟑 : Es el determinante de la matriz que resulta al intercambiaren la matriz A los elementos de
la tercera columna por los elementos del vector B.

𝟏 𝑥0 𝑦0
𝑫𝟑 = [𝟏 𝑥1 𝑦1 ]
𝟏 𝑥2 𝑦2
Se soluciona este determinante por expansión a lo largo de la segunda columna:

1+3 1 𝑥1 2+3 1 𝑥0 3+3 1 𝑥0


𝑫𝟑 = 𝑦0 (−1) [ ] + 𝑦1 (−1) [ ] + 𝑦2 (−1) [ ]
1 𝑥2 1 𝑥2 1 𝑥1
𝑫𝟑 = 𝑦0 (𝑥2 − 𝑥1 ) + 𝑦1 (𝑥2 − 𝑥0 ) + 𝑦2 (𝑥1 − 𝑥0 )

Se tiene que:

𝐷1 𝐷2 𝐷3
𝑎0 = ; 𝑎1 = ; 𝑎2 =
𝐷 𝐷 𝐷
Reemplazando estas expresiones en el polinomio 𝑷𝟐 (𝒙) = 𝒂𝟎 + 𝒂𝟏 𝒙 + 𝒂𝟐 𝒙𝟐 queda:
Métodos Numéricos
sfs

112
𝑦 𝑥1 𝑥2 (𝑥2 − 𝑥1 ) − 𝑦1 𝑥0 𝑥2 (𝑥2 − 𝑥0 ) + 𝑦2 𝑥0 𝑥1 (𝑥1 − 𝑥0 )
𝑷𝟐 (𝒙) = 0 +
(𝑥2 − 𝑥1 )(𝑥2 − 𝑥0 )(𝑥1 − 𝑥0 )

[−𝑦0 (𝑥2 − 𝑥1 )(𝑥2 + 𝑥1 ) − 𝑦1 (𝑥2 − 𝑥0 )(𝑥2 + 𝑥0 ) + 𝑦2 (𝑥1 − 𝑥0 )(𝑥1 + 𝑥0 )]


𝑥
(𝑥2 − 𝑥1 )(𝑥2 − 𝑥0 )(𝑥1 − 𝑥0 )
+

𝑦0 (𝑥2 − 𝑥1) + 𝑦1 (𝑥2 − 𝑥0) + 𝑦2 (𝑥1 − 𝑥0)


(𝑥2 − 𝑥1 )(𝑥2 − 𝑥0 )(𝑥1 − 𝑥0 )

Factorizando, se obtiene:

𝑦0 (𝑥2 − 𝑥1 )[𝑥1 𝑥2 − (𝑥2 + 𝑥1 )𝑥 + 𝑥2 ]


𝑷𝟐 (𝒙) =
(𝑥2 − 𝑥1 )(𝑥2 − 𝑥0 )(𝑥1 − 𝑥0 )

−𝑦1 (𝑥2 − 𝑥0 )[𝑥0 𝑥2 − (𝑥1 + 𝑥0 )𝑥 + 𝑥2 ]


+
(𝑥2 − 𝑥1 )(𝑥2 − 𝑥0 )(𝑥1 − 𝑥0 )

𝑦2 (𝑥1 − 𝑥0 )[𝑥0 𝑥1 − (𝑥2 + 𝑥0 )𝑥 + 𝑥2 ]


+
(𝑥2 − 𝑥1 )(𝑥2 − 𝑥0 )(𝑥1 − 𝑥0 )

Se tiene que:

𝑦0 (𝑥2 − 𝑥1 )[𝑥1 𝑥2 − (𝑥2 + 𝑥1 )𝑥 + 𝑥 2 ] = 𝑦0 (𝑥2 − 𝑥1 )[𝑥1 𝑥2 − 𝑥2 𝑥 − 𝑥1 𝑥 + 𝑥 2 ] =

𝑦0 (𝑥2 − 𝑥1 )[(𝑥1 𝑥2 − 𝑥2 𝑥) − (𝑥1 𝑥 − 𝑥 2 )] = 𝑦0 (𝑥2 − 𝑥1 )[𝑥2 (𝑥1 − 𝑥) − 𝑥(𝑥1 − 𝑥)] =

𝑦0 (𝑥2 − 𝑥1 )(𝑥1 − 𝑥)(𝑥2 − 𝑥) = 𝑦0 (𝑥2 − 𝑥1 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )

−𝑦1 (𝑥2 − 𝑥0 )[𝑥0 𝑥2 − (𝑥2 + 𝑥0 )𝑥 + 𝑥 2 ] = −𝑦1 (𝑥2 − 𝑥0 )[𝑥0 𝑥2 − 𝑥2 𝑥 − 𝑥0 𝑥 + 𝑥 2 ] =

−𝑦1 (𝑥2 − 𝑥0 )[𝑥 2 − 𝑥0 𝑥 − 𝑥2 𝑥 + 𝑥0 𝑥2 ] = −𝑦1 (𝑥2 − 𝑥0 )[(𝑥 2 − 𝑥0 𝑥) − (𝑥2 𝑥 − 𝑥0 𝑥2 ] =


Métodos Numéricos
sfs

113
−𝑦1 (𝑥2 − 𝑥0 )[𝑥(𝑥 − 𝑥0 ) − 𝑥2 (𝑥 − 𝑥0 )] = −𝑦1 (𝑥2 − 𝑥0 )(𝑥 − 𝑥0 )(𝑥 − 𝑥2 )

𝑦2 (𝑥1 − 𝑥0 )[𝑥0 𝑥1 − (𝑥1 + 𝑥0 )𝑥 + 𝑥 2 ] = 𝑦0 (𝑥1 − 𝑥0 )[𝑥0 𝑥1 − 𝑥1 𝑥 − 𝑥0 𝑥 + 𝑥 2 ] =

𝑦2 (𝑥1 − 𝑥0 )[𝑥 2 − 𝑥0 𝑥 − 𝑥1 𝑥 − 𝑥0 𝑥1 ] = 𝑦2 (𝑥1 − 𝑥0 )[(𝑥2 − 𝑥0 𝑥) − (𝑥1 𝑥 − 𝑥0 𝑥1 )] =

𝑦2 (𝑥1 − 𝑥0 )[𝑥(𝑥 − 𝑥0 ) − 𝑥1 (𝑥 − 𝑥0 ) = 𝑦2 (𝑥1 − 𝑥0 )(𝑥 − 𝑥0 )(𝑥 − 𝑥1 )

De esta manera el polinomio P2 (x) queda:

𝑦0 (𝑥2 − 𝑥1 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )


𝑃2 (𝑥 ) = +
(𝑥2 − 𝑥1 )(𝑥2 − 𝑥0 )(𝑥1 − 𝑥0 )
−𝑦1 (𝑥2 − 𝑥0 )(𝑥 − 𝑥0 )(𝑥 − 𝑥2 )
(𝑥2 − 𝑥1 )(𝑥2 − 𝑥0 )(𝑥1 − 𝑥0 )
𝑦2 (𝑥1 − 𝑥0 )(𝑥 − 𝑥0 )(𝑥 − 𝑥1 )
+
(𝑥2 − 𝑥1 )(𝑥2 − 𝑥0 )(𝑥1 − 𝑥0 )

Simplificando, se tiene:

𝒚𝟎 (𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )
𝑷𝟐 (𝒙) =
(𝒙𝟎 − 𝒙𝟏 )(𝒙𝟎 − 𝒙𝟐 )
𝒚𝟏 (𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )
+
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )
𝒚𝟐 (𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )
+
(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )
Métodos Numéricos
sfs

114
Haciendo el mismo análisis que para el polinomio interpolante lineal, se tiene que:

(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )
𝑳𝟎 (𝒙) =
(𝒙𝟎 − 𝒙𝟏 )(𝒙𝟎 − 𝒙𝟐 )
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )
𝑳𝟏 (𝒙) =
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )
𝑳𝟐 (𝒙) =
(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )
1 2
Y el polinomio interpolante de grado dos sería:

P2 (x)  y0 L0 (x)  y1L1 (x)  y2 L2 (x)

2.5.2 EJERCICIOS DE APRENDIZAJE


1. Encuentre el polinomio de grado 2 o parábola de interpolación que pase por los puntos:

 2, 3, 5, 4 7, 2


Datos:

(x0 = −2, y0 = 3);


(x1 = 5, y1 = 4);
(x2 = 7, y2 = 2)

La parábola tiene la forma:

P2 (x)  y0 L0 (x)  y1L1 (x)  y2 L2 (x)

Cálculo de:
(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )
𝐿0 (𝑥 ) =
(𝒙𝟎 − 𝒙𝟏 )(𝒙𝟎 − 𝒙𝟐 )
Métodos Numéricos
sfs

115
(𝒙 − 𝟓)(𝒙 − 𝟕) 𝒙𝟐 − 𝟕𝒙 − 𝟓𝒙 + 𝟑𝟓
𝐿0 (𝑥 ) = = =
(−𝟐 − 𝟓)(−𝟐 − 𝟕) (−𝟕)(−𝟗)

𝒙𝟐 − 𝟏𝟐𝒙 + 𝟑𝟓
𝐿0 (𝑥 ) =
𝟔𝟑
Cálculo de:

(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )
𝐿1 (𝑥 ) =
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )

(𝒙 − (−𝟐))(𝒙 − 𝟕) (𝒙 + 𝟐))(𝒙 − 𝟕)
𝐿1 (𝑥 ) = = =
(𝟓 − (−𝟐))(𝟓 − 𝟕) (𝟓 + 𝟐)(𝟓 − 𝟕)

𝒙𝟐 − 𝟓𝒙 − 𝟏𝟒
𝐿1 (𝑥 ) =
−𝟏𝟒

Cálculo de:
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )
𝐿2 (𝑥 ) =
(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )

(𝒙−(−𝟐))(𝒙−𝟓) (𝒙+𝟐))(𝒙−𝟓) 𝒙𝟐 −𝟓𝒙+𝟐𝒙−𝟏𝟎


= (𝟕+𝟐)(𝟕−𝟓)
= (𝟗)(𝟐)
=
(𝟕−(−𝟐))(𝟕−𝟓)

𝒙𝟐 − 𝟑𝒙 − 𝟏𝟎
𝐿2 (𝑥 ) =
𝟏𝟖
Cálculo de 𝑃2 (𝑥): Recuerde que:

P2 (x)  y0 L0 (x)  y1L1 (x)  y2 L2 (x)


Reemplazando se tiene que:
Métodos Numéricos
sfs

116

𝒙𝟐 −𝟏𝟐𝒙+𝟑𝟓 𝒙𝟐 −𝟓𝒙−𝟏𝟒 𝒙𝟐 −𝟑𝒙−𝟏𝟎


P2 (x)= 𝟑( ) + 𝟒( ) + 𝟐( )
𝟔𝟑 −𝟏𝟒 𝟏𝟖

Simplificando:

𝒙𝟐 −𝟏𝟐𝒙+𝟑𝟓 𝒙𝟐 −𝟓𝒙−𝟏𝟒 𝒙𝟐 −𝟑𝒙−𝟏𝟎


P2 (x)= 𝟏 ( ) − 𝟐( ) + 𝟏( )
𝟐𝟏 𝟕 𝟗

𝟏 𝟐 𝟏 −𝟏𝟐 𝟏𝟎 𝟑 𝟑𝟓 𝟐𝟖 𝟏𝟎
P2 (x)= ( − + ) 𝒙𝟐 − ( + − )𝒙 + ( + − )
𝟐𝟏 𝟕 𝟗 𝟐𝟏 𝟕 𝟗 𝟐𝟏 𝟕 𝟗
𝟖 𝟏𝟏 𝟒𝟏
P2 (x)=− 𝒙𝟐 + 𝒙+
𝟔𝟑 𝟐𝟏 𝟗

Actividad: Verificar este resultado dando la prueba al ejercicio.


1. Encuentre el polinomio de grado 2 o parábola de interpolación que pase por los puntos:

 4, 85, 10, 499 12, 717


Datos:

x0  4, y0  85, x1  10, y1  499, x2  12, y2  717

La parábola tiene la forma:

P2 (x)  y0 L0 (x)  y1L1 (x)  y2 L2 (x)


Cálculo de:
(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )
𝐿0 (𝑥 ) =
(𝒙𝟎 − 𝒙𝟏 )(𝒙𝟎 − 𝒙𝟐 )
(𝒙 − 𝟏𝟎)(𝒙 − 𝟏𝟐) 𝒙𝟐 − 𝟏𝟐𝒙 − 𝟏𝟎𝒙 + 𝟏𝟐𝟎
𝐿0 (𝑥 ) = = =
(𝟒 − 𝟏𝟎)(𝟒 − 𝟏𝟐) (−𝟔)(−𝟖)

𝒙𝟐 − 𝟐𝟐𝒙 + 𝟏𝟐𝟎
𝐿0 (𝑥 ) =
𝟒𝟖
Métodos Numéricos
sfs

117
Cálculo de:

(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )
𝐿1 (𝑥) =
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )
(𝒙 − (−𝟐))(𝒙 − 𝟕) (𝒙 + 𝟐))(𝒙 − 𝟕)
𝐿1 (𝑥 ) = = =
(𝟓 − (−𝟐))(𝟓 − 𝟕) (𝟓 + 𝟐)(𝟓 − 𝟕)

𝒙𝟐 − 𝟓𝒙 − 𝟏𝟒
𝐿1 (𝑥 ) =
−𝟏𝟒

Cálculo de:

(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )
𝐿2 (𝑥 ) =
(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )

(𝒙−(−𝟐))(𝒙−𝟓) (𝒙+𝟐))(𝒙−𝟓) 𝒙𝟐 −𝟓𝒙+𝟐𝒙−𝟏𝟎


= (𝟕+𝟐)(𝟕−𝟓)
= (𝟗)(𝟐) =
(𝟕−(−𝟐))(𝟕−𝟓)

𝒙𝟐 − 𝟑𝒙 − 𝟏𝟎
𝐿2 (𝑥 ) =
𝟏𝟖
Cálculo de 𝑃2 (𝑥): Recuerde que:

P2 (x)  y0 L0 (x)  y1L1 (x)  y2 L2 (x)


Reemplazando se tiene que:

𝒙𝟐 −𝟏𝟐𝒙+𝟑𝟓 𝒙𝟐 −𝟓𝒙−𝟏𝟒 𝒙𝟐 −𝟑𝒙−𝟏𝟎


P2 (x)= 𝟑( ) + 𝟒( ) + 𝟐( )
𝟔𝟑 −𝟏𝟒 𝟏𝟖

Simplificando:

𝒙𝟐 −𝟏𝟐𝒙+𝟑𝟓 𝒙𝟐 −𝟓𝒙−𝟏𝟒 𝒙𝟐 −𝟑𝒙−𝟏𝟎


P2 (x)= 𝟏 ( ) − 𝟐( ) + 𝟏( )
𝟐𝟏 𝟕 𝟗
Métodos Numéricos
sfs

118
𝟏 𝟐 𝟏 −𝟏𝟐 𝟏𝟎 𝟑 𝟑𝟓 𝟐𝟖 𝟏𝟎
P2 (x)= ( − + ) 𝒙𝟐 − ( + − )𝒙 + ( + − )
𝟐𝟏 𝟕 𝟗 𝟐𝟏 𝟕 𝟗 𝟐𝟏 𝟕 𝟗
𝟖 𝟏𝟏 𝟒𝟏
P2 (x)=− 𝒙𝟐 + 𝒙+
𝟔𝟑 𝟐𝟏 𝟗

Actividad: Verificar este resultado dando la prueba al ejercicio.

1. Encuentre el polinomio de grado 2 o parábola de interpolación que pase por los puntos:

 4, 85, 10, 499 12, 717


Datos:

x0  4, y0  85, x1  10, y1  499, x2  12, y2  717

La parábola tiene la forma:

P2 (x)  y0 L0 (x)  y1L1 (x)  y2 L2 (x)


Cálculo de:
(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )
𝐿0 (𝑥 ) =
(𝒙𝟎 − 𝒙𝟏 )(𝒙𝟎 − 𝒙𝟐 )

(𝒙 − 𝟏𝟎)(𝒙 − 𝟏𝟐) 𝒙𝟐 − 𝟏𝟐𝒙 − 𝟏𝟎𝒙 + 𝟏𝟐𝟎


𝐿0 (𝑥 ) = = =
(𝟒 − 𝟏𝟎)(𝟒 − 𝟏𝟐) (−𝟔)(−𝟖)

𝒙𝟐 − 𝟐𝟐𝒙 + 𝟏𝟐𝟎
𝐿0 (𝑥 ) =
𝟒𝟖

Cálculo de:

(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )
𝐿1 (𝑥 ) =
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )
Métodos Numéricos
sfs

119
(𝒙 − 𝟒)(𝒙 − 𝟏𝟐) 𝒙𝟐 − 𝟏𝟔𝒙 + 𝟒𝟖
𝐿1 (𝑥 ) = = =
(𝟏𝟎 − 𝟒)(𝟏𝟎 − 𝟏𝟐) −𝟏𝟐

𝒙𝟐 − 𝟏𝟔𝒙 + 𝟒𝟖
𝐿1 (𝑥 ) =
−𝟏𝟐

Cálculo de:

(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )
𝐿2 (𝑥 ) =
(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )

(𝒙−𝟒)(𝒙−𝟏𝟎) 𝒙𝟐 −𝟏𝟎𝒙−𝟒𝒙+𝟒𝟎 𝒙𝟐 −𝟏𝟒𝒙+𝟒𝟎


(𝟏𝟐−𝟒)(𝟏𝟐−𝟏𝟎)
== (𝟖)(𝟐)
= =
𝟏𝟔

𝒙𝟐 − 𝟏𝟒 + 𝟒𝟎
𝐿2 (𝑥 ) =
𝟏𝟔
Cálculo: P2 (x)  y0 L0 (x)  y1L1 (x)  y2 L2 (x)

𝒙𝟐 −𝟐𝟐𝒙+𝟏𝟐𝟎 𝒙𝟐 −𝟏𝟔𝒙+𝟒𝟖 𝒙𝟐 −𝟏𝟒+𝟒𝟎


P2 (x) =85( ) + 499( ) + 𝟕𝟏𝟕( )
𝟒𝟖 −𝟏𝟐 𝟏𝟔

85 499 717 −1870 7984 10038 10200 23952


P2 (x) =( −
48 12
+
16
) 𝑥2 + ( 48
+
12

16
)𝑥 +( 48

12
+
28680
16
)

P2 (x) = 𝟓𝒙𝟐 − 𝒙 + 𝟗

Actividad: Verificar este resultado dando la prueba al ejercicio.


Métodos Numéricos
sfs

120
CASO 3: Polinomio de grado máximo tres que pasa por cuatro puntos, donde los valores de x son
todos diferentes.

El polinomio tiene la forma:

P3 (x)  y0 L0 (x)  y1L1 (x)  y2 L2 (x)  y3 L3 (x)

(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )(𝒙 − 𝒙𝟑 )
𝑳𝟎 (𝒙) =
(𝒙𝟎 − 𝒙𝟏 )(𝒙𝟎 − 𝒙𝟐 )(𝒙𝟎 − 𝒙𝟑 )
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )(𝒙 − 𝒙𝟑 )
𝑳𝟏 (𝒙) =
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )(𝒙𝟏 − 𝒙𝟑 )
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟑 )
𝑳𝟐 (𝒙) =
(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )(𝒙𝟐 − 𝒙𝟑 )
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )
𝑳𝟑 (𝒙) =
(𝒙𝟑 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )(𝒙𝟑 − 𝒙𝟐 )

Caso 4: Polinomio de grado máximo cuatro que pasa por cinco puntos con la condición que todos
los valores de x son diferentes.

El polinomio tiene la forma:

P4 (x)  y0 L0 (x)  y1L1 (x)  y2 L2 (x)  y3 L3 (x)  y4 L4 (x)

(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )(𝒙 − 𝒙𝟑 )(𝒙 − 𝒙𝟒 )


𝑳𝟎 (𝒙) =
(𝒙𝟎 − 𝒙𝟏 )(𝒙𝟎 − 𝒙𝟐 )(𝒙𝟎 − 𝒙𝟑 )(𝒙𝟎 − 𝒙𝟒 )
Métodos Numéricos
sfs

121
(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 )(𝒙 − 𝒙𝟑 )(𝒙 − 𝒙𝟒 )
𝑳𝟏 (𝒙) =
(𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 )(𝒙𝟏 − 𝒙𝟑 )(𝒙𝟏 − 𝒙𝟒 )

(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟑 )(𝒙 − 𝒙𝟒 )


𝑳𝟐 (𝒙) =
(𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )(𝒙𝟐 − 𝒙𝟑 )(𝒙𝟐 − 𝒙𝟒 )

(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )(𝒙 − 𝒙𝟒 )


𝑳𝟑 (𝒙) =
(𝒙𝟑 − 𝒙𝟎 )(𝒙𝟐 − 𝒙𝟏 )(𝒙𝟑 − 𝒙𝟐 )(𝒙𝟑 − 𝒙𝟒 )

(𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 )(𝒙 − 𝒙𝟑 )


𝑳𝟒 (𝒙) =
(𝒙𝟒 − 𝒙𝟎 )(𝒙𝟒 − 𝒙𝟏 )(𝒙𝟒 − 𝒙𝟐 )(𝒙𝟒 − 𝒙𝟒𝟑 )

Nota: Lo anterior se puede generalizar para polinomios de grado 5, grado


6 y más.

2.5.3 EJERCICIOS DE APRENDIZAJE


Use los polinomios interpolantes de LaGrange de grado uno, dos y tres, más apropiados, para aproximar:

𝒇(𝟐, 𝟓) 𝒔𝒊:
 𝒇(𝟐) = 𝟎, 𝟓𝟏𝟎𝟑𝟕𝟓
 𝒇(𝟐, 𝟐) = 𝟎, 𝟓𝟐𝟎𝟕𝟖𝟒
 𝒇(𝟐, 𝟒) = 𝟎, 𝟓𝟏𝟎𝟒𝟏𝟒
 𝒇(𝟐, 𝟔) = 𝟎, 𝟒𝟖𝟏𝟑𝟑𝟎
 𝒇(𝟐, 𝟖) = 𝟎, 𝟒𝟑𝟓𝟗𝟏𝟔
Métodos Numéricos
sfs

122
Procedimiento:

Como 𝟐, 𝟓 ∈ [𝟐. 𝟒, 𝟐. 𝟔 ]entonces el polinomio de interpolación de LaGrange de grado 1, más


apropiado es el que se obtiene tomando los nodos 𝒙𝟎 = 𝟐, 𝟒 ∧ 𝒙𝟏 = 𝟐, 𝟔 , lo que nos da:
(𝒙 − 𝒙𝟏 ) (𝒙 − 𝒙𝟎 )
𝑷𝟏 (𝒙) = 𝒇(𝒙𝟎 ) + 𝒇(𝒙𝟏 )
(𝒙𝟎 − 𝒙𝟏) (𝒙𝟏 − 𝒙𝟎 )
Reemplazando los valores dados se tiene que:

(𝟐, 𝟓 − 𝟐, 𝟔) (𝟐, 𝟓 − 𝟐, 𝟒)
𝑷𝟏 (𝟐, 𝟓) = 𝟎, 𝟓𝟏𝟒𝟏𝟒𝟕 + 𝟎, 𝟒𝟖𝟏𝟑𝟑𝟎𝟔
(𝟐, 𝟒 − 𝟐, 𝟔) (𝟐, 𝟔 − 𝟐, 𝟒)

𝑷𝟏 (𝟐, 𝟓) = 𝟎, 𝟒𝟗𝟓𝟖𝟕𝟐𝟕 ≈ 𝒇(𝟐, 𝟓)


Para el caso del grado 2, hay dos polinomios interpelantes igualmente apropiados.
Un primer polinomio se obtiene tomando los nodos: 𝒙𝟎 = 𝟐, 𝟒, 𝒙𝟏 = 𝟐, 𝟒 ∧ 𝒙𝟐 = 𝟐, 𝟔
Lo que nos da:

(𝒙 − 𝒙𝟏 )(𝒙 − 𝒙𝟐 ) (𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟐 ) (𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 )


𝑷𝟐 (𝒙) = 𝒇(𝒙𝟎 ) + 𝒇(𝒙𝟏 ) + 𝒇(𝒙𝟐 )
(𝒙𝟎 − 𝒙𝟏 )(𝒙𝟎 − 𝒙𝟐 ) (𝒙𝟏 − 𝒙𝟎 )(𝒙𝟏 − 𝒙𝟐 ) (𝒙𝟐 − 𝒙𝟎 )(𝒙𝟐 − (𝒙𝟏 )

(𝟐, 𝟓 − 𝟐, 𝟒)(𝟐, 𝟓 − 𝟐, 𝟔)
𝑷𝟐 (𝟐, 𝟓) = 𝟎, 𝟓𝟐𝟎𝟕𝟖𝟒
(𝟐, 𝟐 − 𝟐, 𝟒)(𝟐, 𝟐 − 𝟐, 𝟔)

(𝟐, 𝟓 − 𝟐, 𝟐)(𝟐, 𝟓 − 𝟐, 𝟔)
+𝟎, 𝟓𝟏𝟎𝟒𝟏𝟒
(𝟐, 𝟒 − 𝟐, 𝟐)(𝟐, 𝟒 − 𝟐, 𝟔)

(𝟐, 𝟓 − 𝟐, 𝟐)(𝟐, 𝟓 − 𝟐, 𝟒)
+𝟎, 𝟒𝟖𝟏𝟑𝟑𝟎𝟔
(𝟐, 𝟔 − 𝟐, 𝟐)(𝟐, 𝟔 − 𝟐, 𝟒)
=

𝟎, 𝟒𝟗𝟖𝟐𝟏𝟐𝟎 ≈ 𝒇(𝟐, 𝟓)
El otro polinomio interpolante de grado dos se obtiene tomando los nodos:
Métodos Numéricos
sfs

123
𝒙𝟎 = 𝟐, 𝟒, 𝒙𝟏 = 𝟐, 𝟔 ∧ 𝒙𝟐 = 𝟐, 𝟖
(𝟐, 𝟓 − 𝟐, 𝟔)(𝟐, 𝟓 − 𝟐, 𝟖)
𝑷𝟐 (𝟐, 𝟓) = 𝟎, 𝟓𝟏𝟎𝟒𝟏𝟒
(𝟐, 𝟐 − 𝟐, 𝟒)(𝟐, 𝟐 − 𝟐, 𝟔)

(𝟐, 𝟓 − 𝟐, 𝟐)(𝟐, 𝟓 − 𝟐, 𝟔)
+𝟎, 𝟓𝟏𝟎𝟒𝟏𝟒
(𝟐, 𝟒 − 𝟐, 𝟐)(𝟐, 𝟒 − 𝟐, 𝟔)

(𝟐, 𝟓 − 𝟐, 𝟐)(𝟐, 𝟓 − 𝟐, 𝟒)
+𝟎, 𝟒𝟑𝟓𝟗𝟏𝟔𝟎
(𝟐, 𝟔 − 𝟐, 𝟐)(𝟐, 𝟔 − 𝟐, 𝟒)
==

𝟎, 𝟒𝟗𝟕𝟗𝟏𝟒𝟎 ≈ 𝒇(𝟐, 𝟓)

2.5.4 TALLER DE ENTRENAMIENTO


Para el grado tres el polinomio interpolante más apropiado, se obtiene tomando los nodos:

𝒙𝟎 = 𝟐, 𝟐 𝒙𝟏 = 𝟐, 𝟒 𝒙𝟐 = 𝟐, 𝟔 ∧ 𝒙𝟑 = 𝟐, 𝟖
Ya que 𝟐, 𝟓 ∈ [𝟐. 𝟐, 𝟐. 𝟖] y 𝟐, 𝟐; 𝟐, 𝟒; 𝟐, 𝟔 ∧ 𝟐, 𝟖 son los nodos más cercanos a 2,5

Nota: Resuelve el problema teniendo como modelo los ejercicios realizados, debes tener
presente que es para el grado tres el polinomio interpolante.

Polinomio interpolante de Newton

2.5.4.1 DIFERENCIAS DIVIDIDAS

 Video: Interpolación polinomial por método diferencias divididas


 ”Ver Video”: https://www.youtube.com/watch?v=GTq_9bq1CsM
Métodos Numéricos
sfs

124
Dados n + 1 puntos x0 , y0 , x1 , y1 ,...,xn , yn  con 𝒙𝟎 , 𝒙𝟏 , . . , 𝒙𝒏 números distintos.

Para obtener los coeficientes de los polinomios de interpolación utilizando el método de Newton,
es necesario efectuar cálculos con resultados anteriores, para facilitar esto, se utilizarán las
diferencias divididas:

 Para dos puntos (𝒙𝟎 , 𝒇(𝒙𝟎 )) ∧ (𝒙𝟏 , 𝒇(𝒙𝟏 )) se tiene que las primeras y únicas diferencias
divididas son:

𝒇(𝒙𝟏 ) − 𝒇(𝒙)𝟎
𝒇[𝒙𝟏 , 𝒙𝟎 ] =
𝒙𝟏 − 𝒙𝟎
 Para tres puntos (𝒙𝟎 , 𝒇(𝒙𝟎 )), (𝒙𝟏 , 𝒇(𝒙𝟏 )) ∧ (𝒙𝟐 , 𝒇(𝒙𝟐 )) se tiene que las primeras
diferencias divididas son:

𝒇(𝒙𝟐 )−𝒇(𝒙𝟏 ) 𝒇(𝒙𝟏 )−𝒇(𝒙)𝟎


𝒇[𝒙𝟐 , 𝒙𝟏 ] = ∧ 𝒇[𝒙𝟏 , 𝒙𝟎 ] =
𝒙𝟐 −𝒙𝟏 𝒙𝟏 −𝒙𝟎
Las segundas diferencias divididas son:

𝒇(𝒙𝟐 , 𝒙𝟏 ) − 𝒇(𝒙𝟏 , 𝒙𝟎 )
𝒇[𝒙𝟐 , 𝒙𝟏, 𝒙𝟎 ] =
𝒙𝟏 − 𝒙𝟎

 Para cuatro puntos (𝒙𝟎 , 𝒇(𝒙𝟎 )), (𝒙𝟏 , 𝒇(𝒙𝟏 )), (𝒙𝟐 , 𝒇(𝒙𝟐 )) ∧ (𝒙𝟑 , 𝒇(𝒙𝟑 )) se tiene que las
primeras diferencias

divididas son:

𝒇(𝒙𝟑 ) − 𝒇(𝒙𝟐 )
𝒇[ 𝒙 𝟑 , 𝒙 𝟐 ] = ,
𝒙𝟑 − 𝒙𝟐

𝒇(𝒙𝟐 ) − 𝒇(𝒙𝟏 )
𝒇[ 𝒙 𝟐 , 𝒙 𝟏 ] = ∧
𝒙𝟐 − 𝒙𝟏
Métodos Numéricos
sfs

125
𝒇(𝒙𝟏 ) − 𝒇(𝒙𝟎 )
[ 𝒙𝟏 , 𝒙𝟎 ] =
𝒙𝟏 − 𝒙𝟎
Las segundas diferencias divididas son:

𝒇(𝒙𝟑 , 𝒙𝟐 ) − 𝒇(𝒙𝟐 , 𝒙𝟏 )
𝒇[𝒙𝟑 , 𝒙𝟐, 𝒙𝟏 ] =
𝒙𝟑 − 𝒙𝟏

𝒇(𝒙𝟐 , 𝒙𝟏 ) − 𝒇(𝒙𝟏 , 𝒙𝟎 )
𝒇[𝒙𝟐 , 𝒙𝟏, 𝒙𝟎 ] =
𝒙𝟐 − 𝒙𝟎
 Las Terceras diferencias divididas son:

𝒇(𝒙𝟑 , 𝒙𝟐 , 𝒙𝟏 ) − 𝒇(𝒙𝟐 , 𝒙𝟏 , 𝒙𝟎 )
𝒇[𝒙𝟑 , 𝒙𝟐, 𝒙𝟏 , 𝒙𝟎 ] =
𝒙𝟑 − 𝒙𝟎
Y así sucesivamente.

2.5.5 TALLER DE APRENDIZAJE


Sean los puntos (𝟏, 𝟑), (𝟒, 𝟕) ∧ (𝟖, 𝟐), determine todas las diferencias divididas.

Procedimiento:

𝒙𝟎 = 𝟏, 𝒇(𝒙𝟎 ) = 𝟎
𝒙𝟏 = 𝟒, 𝒇(𝒙𝟏 ) = 𝟏, 𝟑𝟖𝟔𝟑
𝒙𝟐 = 𝟖, 𝒇(𝒙𝟐 ) = 𝟏, 𝟕𝟗𝟏𝟖

 Las primeras diferencias divididas son:


Métodos Numéricos
sfs

126
𝟐−𝟕 𝟓
𝒇[𝒙𝟐 , 𝒙𝟏 ] = =−
𝟖−𝟒 𝟒
𝟕−𝟑 𝟒
[𝒙𝟏 , 𝒙𝟎 ] = =
𝟒−𝟏 𝟑
 Las segundas diferencias divididas son:

𝟓 𝟒 𝟑𝟏
𝒇(𝒙𝟐 , 𝒙𝟏 ) − 𝒇(𝒙𝟏 , 𝒙𝟎 ) − 𝟒 − 𝟑 − 𝟏𝟐 𝟑𝟏
𝒇[𝒙𝟐 , 𝒙𝟏, 𝒙𝟎 ] = = = =−
𝒙𝟐 − 𝒙𝟎 𝟖−𝟏 𝟕 𝟖𝟒

 Forma general interpolante de Newton:

(𝒙𝟎 , 𝒇(𝒙𝟎 )) ∧ (𝒙𝟏 , 𝒇(𝒙𝟏 )), (𝒙𝟐 , 𝒇(𝒙𝟐 )), … , (𝒙𝒏 , 𝒇(𝒙𝒏 ))
Con:

𝒙𝟎 ≠ 𝒙𝟏 ≠ 𝒙𝟐 ≠ ⋯ ≠ 𝒙𝒏
El polinomio que pasa exactamente por todos los puntos es:

𝑷𝒏 = 𝒃𝒐 + 𝒃𝟏 (𝒙 − 𝒙𝟎 ) + 𝒃𝟐 (𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 ) + ⋯ + 𝒃𝒏 (𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 ) … (𝒙 − 𝒙𝒏−𝟏 )

Resumiendo, se puede decir que dados los puntos:

(𝒙𝟎 , 𝒇(𝒙𝟎 )) ∧ (𝒙𝟏 , 𝒇(𝒙𝟏 )), (𝒙𝟐 , 𝒇(𝒙𝟐 )), … , (𝒙𝒏 , 𝒇(𝒙𝒏 ))
Con:

𝒙𝟎 ≠ 𝒙𝟏 ≠ 𝒙𝟐 ≠ ⋯ ≠ 𝒙𝒏
El polinomio de Newton que pasa exactamente por todos los puntos es:

𝑷𝒏 = 𝒃𝒐 + 𝒃𝟏 (𝒙 − 𝒙𝟎 ) + 𝒃𝟐 (𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 ) + ⋯ + 𝒃𝒏 (𝒙 − 𝒙𝟎 )(𝒙 − 𝒙𝟏 ) … (𝒙 − 𝒙𝒏−𝟏 )

Con:

𝒃𝒐 = 𝒇(𝒙𝟎 )
𝒃 𝟏 = 𝒇[ 𝒙 𝟏 , 𝒙 𝟎 ]
Métodos Numéricos
sfs

127
𝒃 𝟐 = 𝒇[ 𝒙 𝟐 , 𝒙 𝟏 , 𝒙 𝟎 ]
𝒃 𝟑 = 𝒇[ 𝒙 𝟑 , 𝒙 𝟐 , 𝒙 𝟏 , 𝒙 𝟎 ]
𝒃 𝟒 = 𝒇[ 𝒙 𝟒 , 𝒙 𝟑 , 𝒙 𝟐 , 𝒙 𝟏 , 𝒙 𝟎 ]
𝒃 𝟓 = 𝒇[ 𝒙 𝟓 , 𝒙 𝟒 , 𝒙 𝟑 , 𝒙 𝟐 , 𝒙 𝟏 , 𝒙 𝟎 ]
𝒃 𝒏 = 𝒇[ 𝒙 𝒏 … , 𝒙 𝟓 , 𝒙 𝟒 , 𝒙 𝟑 , 𝒙 𝟐 , 𝒙 𝟏 , 𝒙 𝟎 ]

EJERCICIO DE APRENDIZAJE:

1. Para  3, 40, 1,  8 5, 40, encuentre el polinomio de Newton de grado máximo que
pasa por todos los puntos.

Procedimiento

Datos:

𝒙𝟎 = −𝟑 ∧ 𝒇(𝒙𝟎 ) = 𝟒𝟎
𝒙𝟏 = 𝟏 ∧ 𝒇(𝒙𝟏 ) = −𝟖
𝒙𝟐 = 𝟓 ∧ 𝒇(𝒙𝟐 ) = 𝟒𝟎

El polinomio de Newton de grado máximo que pasa por todos los puntos es:

P2 (x)  b0  b1 x  x0  b2 x  x0 x  x1 

Con:

𝒃𝟎 = 𝒇(𝒙𝟎 )

𝒃𝟏 = 𝒇[𝒙𝟏 , 𝒙𝟎 ]

𝒃𝟐 = 𝒇[𝒙𝟐 , 𝒙𝟏 , 𝒙𝟎 ]

Determinación de 𝒃𝟎 = 𝒇(𝒙𝟎 ):

𝒃𝟎 = 𝒇(𝒙𝟎 ) = 𝟒𝟎
Métodos Numéricos
sfs

128
Cálculo de todas las diferencias divididas:

Las primeras diferencias divididas son:


𝟒𝟎 − (−𝟖) 𝟒𝟖
𝒇[𝒙𝟐 , 𝒙𝟏 ] = = = 𝟏𝟐
𝟓−𝟏 𝟒

−𝟖 − 𝟒𝟎 −𝟒𝟖
𝒇[𝒙𝟏 , 𝒙𝟎 ] = = = −𝟏𝟐
𝟓−𝟏 𝟒
Las segundas diferencias divididas son:

𝒇[𝒙𝟐 , 𝒙𝟏 ] − 𝒇[𝒙𝟏 , 𝒙𝟎 ] 𝟏𝟐 − (−𝟏𝟐) 𝟐𝟒


𝒇[𝒙𝟐 , 𝒙𝟏 , 𝒙𝟎 ] =] = = =𝟑
𝒙𝟐 − 𝒙𝟎 𝟓 − (−𝟑) 𝟖

Determinación de:

𝒃𝟏 = 𝒇[𝒙𝟏 , 𝒙𝟎 ] = −𝟏𝟐
Determinación de:

𝒃𝟐 = 𝒇[𝒙𝟐 , 𝒙𝟏 , 𝒙𝟎 ] = 𝟑
Determinación de:

P2 (x)  b0  b1 x  x0  b2 x  x0 x  x1 

𝑷𝟐 (𝒙) = 𝟒𝟎 + (−𝟏𝟐)(𝒙 − (−𝟑)) + 𝟑(𝒙 − (−𝟑))(𝒙 − 𝟏)

𝑷𝟐 (𝒙) = 𝟒𝟎 − 𝟏𝟐(𝒙 + 𝟑) + 𝟑(𝒙 + 𝟑)(𝒙 − 𝟏)

𝑷(𝒙) = 𝟒𝟎 − 𝟏𝟐𝒙 − 𝟑𝟔 + 𝟑(𝒙𝟐 + 𝟐𝒙 − 𝟑)

𝑷𝟐 (𝒙) = 𝟒𝟎 − 𝟏𝟐𝒙 − 𝟑𝟔 + 𝟑𝒙𝟐 + 𝟔𝒙 − 𝟗

𝑷𝟐 (𝒙) = 𝟑𝒙𝟐 − 𝟔𝒙 − 𝟓
Métodos Numéricos
sfs

129
2.5.6 ACTIVIDAD: REALIZAR LA PRUEBA DEL EJERCICIO
2. Para  10, 4, 3, 8 5, 10, encuentre el polinomio de Newton de grado máximo que
pasa por todos los puntos.

Procedimiento

Datos:

𝒙𝟎 = −𝟏𝟎 ∧ 𝒇(𝒙𝟎 ) = 𝟒
𝒙𝟏 = 𝟑 ∧ 𝒇(𝒙𝟏 ) = 𝟖
𝒙𝟐 = 𝟓 ∧ 𝒇(𝒙𝟐 ) = 𝟏𝟎
El polinomio de Newton de grado máximo que pasa por todos los puntos es:

P2 (x)  b0  b1 x  x0  b2 x  x0 x  x1 

Con:

𝒃𝟎 = 𝒇(𝒙𝟎 )

𝒃𝟏 = 𝒇[𝒙𝟏 , 𝒙𝟎 ]

𝒃𝟐 = 𝒇[𝒙𝟐 , 𝒙𝟏 , 𝒙𝟎 ]

Determinación de 𝒃𝟎 = 𝒇(𝒙𝟎 ):

𝒃𝟎 = 𝒇(𝒙𝟎 ) = 𝟒
Cálculo de todas las diferencias divididas:

Las primeras diferencias divididas son:


𝟏𝟎 − 𝟖 𝟐
𝒇[𝒙𝟐 , 𝒙𝟏 ] = = =𝟏
𝟓−𝟑 𝟐

𝟖−𝟒 𝟒
𝒇[𝒙𝟏 , 𝒙𝟎 ] = = = −𝟏𝟐
𝟑 − (−𝟏𝟎) 𝟏𝟑

Las segundas diferencias divididas son:


Métodos Numéricos
sfs

130
𝟒 𝟗
𝒇[𝒙𝟐 , 𝒙𝟏 ] − 𝒇[𝒙𝟏 , 𝒙𝟎 ] 𝟏−
𝒇[𝒙𝟐 , 𝒙𝟏 , 𝒙𝟎 ] =] = 𝟏𝟑 = 𝟏𝟑 = 𝟑
𝒙𝟐 − 𝒙𝟎 𝟓 − (−𝟏𝟎) 𝟏𝟓 𝟔𝟓

Determinación de:
𝟒
𝒃𝟏 = 𝒇[𝒙𝟏 , 𝒙𝟎 ] =
𝟏𝟑

Determinación de:

𝟑
𝒃𝟐 = 𝒇[𝒙𝟐 , 𝒙𝟏 , 𝒙𝟎 ] =
𝟔𝟓
Determinación de:

P2 (x)  b0  b1 x  x0  b2 x  x0 x  x1 

𝟒 𝟑
𝑷𝟐 (𝒙) = 𝟒 + (𝒙 − (−𝟏𝟎)) + (𝒙 − (−𝟏𝟎))(𝒙 − 𝟑)
𝟏𝟑 𝟔𝟓
𝟒 𝟑
𝑷𝟐 (𝒙) = 𝟒 + (𝒙 + 𝟏𝟎) + (𝒙 + 𝟏𝟎)(𝒙 − 𝟑)
𝟏𝟑 𝟔𝟓
𝟒 𝟒𝟎 𝟑 𝟐
𝑷(𝒙) = 𝟒 + 𝒙+ + (𝒙 + 𝟕𝒙 − 𝟑𝟎)
𝟏𝟑 𝟏𝟑 𝟔𝟓
𝟒 𝟒𝟎 𝟑 𝟐 𝟐𝟏 𝟗𝟎
𝑷𝟐 (𝒙) = 𝟒 − 𝒙+ + 𝒙 + 𝒙−
𝟏𝟑 𝟏𝟑 𝟔𝟓 𝟔𝟓 𝟔𝟓
𝟑 𝟒𝟏 𝟕𝟒
𝑷𝟐 (𝒙) = 𝒙𝟐 + 𝒙+
𝟔𝟓 𝟔𝟓 𝟏𝟑

2.5.7 ACTIVIDAD: REALIZAR LA PRUEBA DEL EJERCICIO.


3. Para  6, 11, encuentre el polinomio de Newton de grado máximo que pasa por
todos los puntos.

Procedimiento
Métodos Numéricos
sfs

131
Datos:

𝒙𝟎 = 𝟒 ∧ 𝒇(𝒙𝟎 ) = 𝟑
𝒙𝟏 = 𝟔 ∧ 𝒇(𝒙𝟏 ) = 𝟏𝟏

El polinomio de Newton de grado máximo que pasa por todos los puntos es:

P1(x)  b0  b1 x  x0 

Con:

𝒃𝟎 = 𝒇(𝒙𝟎 )

𝒃𝟏 = 𝒇[𝒙𝟏 , 𝒙𝟎 ]

Determinación de 𝒃𝟎 = 𝒇(𝒙𝟎 ):

𝒃𝟎 = 𝒇(𝒙𝟎 ) = 𝟑

Determinación de:

𝟏𝟏 − 𝟑 𝟖
𝒃𝟏 = 𝒇[𝒙𝟏 , 𝒙𝟎 ] = = =𝟒
𝟔−𝟒 𝟐
Determinación de:

P1(x)  b0  b1 x  x0 

𝑷𝟏 (𝒙) = 𝟑 + 𝟒(𝒙 − 𝟒)

𝑷𝟏 (𝒙) = 𝟑 + 𝟒𝒙 − 𝟏𝟔

𝑷𝟏 (𝒙) = −𝟏𝟑 + 𝟒𝒙
Métodos Numéricos
sfs

132
2.5.8 ACTIVIDAD: REALIZAR LA PRUEBA DEL EJERCICIO
2.5.9 TALLER DE ENTRENAMIENTO
Teniendo presente los conceptos vistos y los ejercicios de aprendizaje, resuelva el siguiente taller:

1. Utilizando el polinomio interpolante de LaGrange y el Polinomio interpolante de


Newton, encuentre el polinomio de grado máximo para los datos de la tabla. (Ambos
polinomios deben dar el mismo resultado).

2. x -6 4 10

y -1 10 3

2. Utilizando el polinomio interpolante de LaGrange y el Polinomio interpolante de


Newton, encuentre el polinomio de grado máximo para los datos de la tabla. (Ambos
polinomios deben dar el mismo resultado).

x 10 20

y -3 15

3. Encuentre un polinomio de la forma

y  f x  Ax3  Bx2  Bx  C

Que pase por los puntos:

(−3, 92), (2, 22), (7, 398)

𝒙 -4 -1 4 5 8 9 11 13 20 25

𝒚 10 6 3 4 7 10 9 5 1 -4

Halle el polinomio de regresión de grado 3 que mejor aproxime estos datos.


Métodos Numéricos

133
3 UNIDAD 3 CÁLCULO INFINITESIMAL Y MÉTODOS
NUMÉRICOS

 Video: Aplicación de la Regla de Los Trapecios al cálculo de Integrales


 ”Ver Video”: https://www.youtube.com/watch?v=Zu4AeUAl0R8

 Video: Resolución Ecuaciones Diferenciales 1er orden por Euler, Euler Modificado
y Runge-Kutta 2da Parte
 ”Ver Video”: https://www.youtube.com/watch?v=aypLi3gIepU

3.1.1 OBJETIVO GENERAL


Aplicar algunos métodos numéricos en el cálculo infinitesimal (derivadas, integrales y solución del
problema de valor inicial).

3.1.2 OBJETIVOS ESPECÍFICOS


 Calcular la primera y segunda derivada utilizando algunos métodos numéricos. Calcular
integrales utilizando la regla del trapecio.

 Solucionar ecuaciones diferenciales de primer orden utilizando algunos métodos


numéricos.

3.2 TEMA 1 DIFERENCIACIÓN NUMÉRICA


Cálculo numérico de la primera derivada

FÓRMULA DE DOS PUNTOS: Del cálculo diferencial se sabe que:

Sea la función y  f x, la primera derivada de esta función se obtiene aplicando la siguiente fórmula:
Métodos Numéricos
sfs

134

1(
𝑓 (𝑥 + ℎ) − 𝑓(𝑥)
𝑓 𝑥 ) = lim
ℎ→0 ℎ
Como la derivada se va a hallar numéricamente:

1(
𝑓(𝑥 + ℎ) − 𝑓(𝑥)
𝑓 𝑥 ) = lim −𝐸
ℎ→0 ℎ
Si se quiere evaluar la derivada en un valor 𝑥0 , la derivada se obtiene como:

1(
𝑓(𝑥0 + ℎ) − 𝑓(𝑥0 )
𝑓 𝑥0 ) = lim −𝐸
ℎ→0 ℎ
Donde el error se obtiene como 𝐸 :

ℎ "
𝐸 = 𝑓 (𝑐)
2
Con: 𝑥0 < 𝑐 < 𝑥0 + ℎ
Se toma el valor de c que produce el mayor error, esto se llama encontrar una cota para el error.

Pero este error se puede calcular de manera exacta sólo si se conoce la función f x

3.2.1 EJERCICIO DE APRENDIZAJE

1. Este ejemplo fue tomado de (Burden & Faires, 2003)

Sean 𝒇(𝒙) = 𝐥𝐧 𝒙 ∧ 𝒙𝟎 = 𝟏, 𝟖

a. Obtenga 𝒇(𝟏, 𝟖) de manera exacta.

b. Obtenga 𝒇(𝟏, 𝟖) de manera aproximada para valores de ℎ = 0,1, ℎ = 0,01 ∧ ℎ = 0,001

Obtenga también el error real y una cota para el error.


Procedimiento:

a. Como 𝒇(𝒙) = 𝐥𝐧 𝒙 ∧ 𝒙𝟎 = 𝟏, 𝟖
Métodos Numéricos
sfs

135
1 𝟏
𝑓 1 (𝑥) = ⟹ 𝒇(𝟏, 𝟖) = = 𝟎, 𝟓𝟓𝟓𝟓𝟓𝟓𝟓 …
𝑥 𝟏, 𝟖

b. Para h  0.1

1 𝒇(𝟏, 𝟖 + 𝟎, 𝟏) − 𝒇(𝟏, 𝟖) 𝐥𝐧 𝟏, 𝟗 − 𝐥𝐧 𝟏, 𝟖
𝑓 (1,8) = = = 𝟎, 𝟓𝟒𝟎𝟔𝟕𝟐𝟐𝟏 …
𝟎, 𝟏 𝟎, 𝟏

El Error real es:

𝐸𝑅𝑒𝑎𝑙 = |𝟎, 𝟓𝟓𝟓𝟓𝟓𝟓𝟓 − 𝟎, 𝟓𝟒𝟎𝟔𝟕𝟐𝟐𝟏| = 𝟎, 𝟎𝟏𝟒𝟖𝟖𝟑𝟑𝟒𝟓

El Error Estimado o calculado: Hay que obtener 𝑓 " (𝑥):

1 1
𝑓 1 (𝑥 ) = ⇒ 𝑓 " (𝑥 ) = − 2
𝑥 𝑥
0,1 " 0,1 1 0,1 1
𝐸𝐶𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 = 𝑓 (𝑐) = × |− 2 | ≤ × ≤ 0,015432098
2 2 𝑐 2 1,82

Nota: El error se calculó con c  1.8 y no con c  1.9, ya que con c  1.8 se obtiene
un error mayor y de esta manera se puede garantizar que el error verdadero sea
menor que este error calculado

a. Para h  0.01

1 𝒇(𝟏, 𝟖 + 𝟎, 𝟎𝟏) − 𝒇(𝟏, 𝟖) 𝐥𝐧 𝟏, 𝟖𝟏 − 𝐥𝐧 𝟏, 𝟖


𝑓 (1,8) = = = 𝟎, 𝟓𝟓𝟒𝟎𝟏𝟖𝟎𝟑𝟕 …
𝟎, 𝟎𝟏 𝟎, 𝟎𝟏

El Error real es:

𝐸𝑅𝑒𝑎𝑙 = |𝟎, 𝟓𝟓𝟓𝟓𝟓𝟓𝟓 − 𝟎, 𝟓𝟓𝟒𝟎𝟏𝟖𝟎𝟑𝟕| = 𝟎. 𝟎𝟎𝟏𝟓𝟑𝟕𝟓𝟏𝟖𝟓𝟎, 𝟎𝟏

El Error Estimado o calculado: Hay que obtener 𝑓 " (𝑥):


Métodos Numéricos
sfs

136

1(
1 "( )
1
𝑓 𝑥) = ⇒𝑓 𝑥 =− 2
𝑥 𝑥

0,01 " 0,01 1 0,01 1


𝐸𝐶𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 = 𝑓 (𝑐) = × |− 2 | ≤ × ≤ 0,0015432098
2 2 𝑐 2 1,82

b. Para h  0.01

1 𝒇(𝟏, 𝟖 + 𝟎, 𝟎𝟏) − 𝒇(𝟏, 𝟖) 𝐥𝐧 𝟏, 𝟖𝟏 − 𝐥𝐧 𝟏, 𝟖


𝑓 (1,8) = = = 𝟎, 𝟓𝟓𝟒𝟎𝟏𝟖𝟎𝟑𝟕 …
𝟎, 𝟎𝟏 𝟎, 𝟎𝟏

El Error real es:

𝐸𝑅𝑒𝑎𝑙 = |𝟎, 𝟓𝟓𝟓𝟓𝟓𝟓𝟓 − 𝟎, 𝟓𝟓𝟒𝟎𝟏𝟖𝟎𝟑𝟕| = 𝟎, 𝟎𝟎𝟏𝟓𝟑𝟕𝟓𝟏𝟖𝟓

El Error Estimado o calculado: Hay que obtener 𝑓 " (𝑥):

1(
1 "
1
𝑓 𝑥 ) = ⇒ 𝑓 (𝑥 ) = − 2
𝑥 𝑥

0,01 " 0,01 1 0,01 1


𝐸𝐶𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 = 𝑓 (𝑐) = × |− 2 | ≤ × ≤ 0,001543200987
2 2 𝑐 2 1,82

d. Para h  0.001

1 𝒇(𝟏, 𝟖 + 𝟎, 𝟎𝟎𝟏) − 𝒇(𝟏, 𝟖) 𝐥𝐧 𝟏, 𝟖𝟎𝟏 − 𝐥𝐧 𝟏, 𝟖


𝑓 (1,8) = = = 𝟎, 𝟓𝟓𝟓𝟒𝟎𝟏𝟐𝟗𝟏 …
𝟎, 𝟎𝟎𝟏 𝟎, 𝟎𝟎𝟏
Métodos Numéricos
sfs

137
El Error real es:

𝐸𝑅𝑒𝑎𝑙 = |𝟎, 𝟓𝟓𝟓𝟓𝟓𝟓𝟓 − 𝟎, 𝟓𝟓𝟓𝟒𝟎𝟏𝟐𝟗𝟏| = 𝟎, 𝟎𝟎𝟎𝟏𝟓𝟒𝟐𝟔𝟒𝟓

El Error Estimado o calculado: Hay que obtener 𝑓 " (𝑥):

1(
1 "( )
1
𝑓 𝑥) = ⇒𝑓 𝑥 =− 2
𝑥 𝑥

0,001 " 0,001 1 0,001 1


𝐸𝐶𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 = 𝑓 (𝑐) = × |− 2 | ≤ × ≤ 0,00015432087
2 2 𝑐 2 1,82

___________________________________________
Dada:

𝑓(3+0,1)−𝑓(3)
𝒇(𝒙) = , Halle:
0,1

a) f '3 De manera exacta.

b) f '3 De manera aproximada con h  0.1

c) Error real.
d) Una cota para el error cometido.

Procedimiento:

a. 𝑓 1 (𝑥 ) = 10𝑥 − 5 ⟹ 𝑓(3) = 10(3) − 5 = 25


b. 𝑓 (3) de manera aproximada con ℎ1 = 0,1
𝑓(3+0,1)−𝑓(3) 𝑓(3,1)−𝑓(3) 5(3,1)2 −5(3,1)+8−(3)−[5(3)2 −5(3)+8]
𝑓(3) = = =
0,1 0,1 0,1
Métodos Numéricos
sfs

138

40,55 − 38 2,55
𝑓(3) = = 𝑓(3) = = 25,5
0,1 0,1
c. Error Real:

𝐸𝑅𝑒𝑎𝑙 = |25 − 25,5| = 0,5

d. Una cota para el error cometido:

𝑓 1 (𝑥 ) = 10𝑥 − 5 ⟹ 𝑓"(𝑥) = 10
𝑓"(3) = 10 ∧ 𝑓"(3,1) = 10
0,1 0,1
𝐸𝐶𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 = 𝑓"(𝑐) = × |10|
2 2

 Fórmulas de tres puntos (O formula con error de orden 𝒉𝟐 )


Para obtener diferenciación numérica, existen otras fórmulas o métodos, con los cuales se obtiene una
mayor exactitud. Dichas fórmulas son las siguientes:

𝟏
𝒇𝟏 (𝒙𝟎 ) = [−𝟑𝒇(𝒙𝟎 ) + 𝟒𝒇(𝒙𝟎 + 𝒉) − 𝒇(𝒙𝟎 + 𝟐𝒉)]
𝟐𝒉
Que se conoce como Fórmula de Tres Puntos.

Los tres puntos son:

 𝒙𝟎

 𝒙𝟎 + 𝒉
Métodos Numéricos
sfs

139
 𝒙𝟎 + 𝟐𝒉
El Error está dado por:

𝒉𝟐 𝟏𝟏𝟏
𝑬= 𝒇 (𝒄)
𝟑
Con:

𝒙𝟎 < 𝒄 < 𝒙𝟎 + 𝟐𝒉
El valor de c se toma en este intervalo de tal manera que, con este se, obtenga el mayor error.

𝟏(
𝟏
𝒇 𝒙𝟎 ) = [𝒇(𝒙𝟎 + 𝒉) − 𝒇(𝒙𝟎 − 𝒉)]
𝟐𝒉
También se conoce como fórmula de. Los tres puntos son: x0  h, x0  x0  h aunque
x0 no aparece en la fórmula.

El error es:

𝒉𝟐 𝟏𝟏𝟏
𝑬= 𝒇 (𝒄)
𝟔

Nota: 𝒇𝟏𝟏𝟏 (𝒄): 𝒔𝒆 𝒍𝒆𝒆 𝒄𝒐𝒎𝒐 𝒕𝒆𝒓𝒄𝒆𝒓𝒂 𝒅𝒆𝒓𝒊𝒗𝒂𝒅𝒂 𝒅𝒆 𝒍𝒂 𝒇𝒖𝒏𝒄𝒊ó𝒏.

𝒙𝟎 − 𝒉 < 𝒄 < 𝒙𝟎 + 𝒉
El valor de c se toma en este intervalo de tal manera que, con este, se obtenga el mayor error.

3.2.1.1 CÁLCULO NUMÉRICO DE LA SEGUNDA DERIVADA

𝟏
𝒇" (𝒙𝟎 ) = [𝒇(𝒙𝟎 − 𝒉) − 𝟐𝒇(𝒙𝟎 ) + 𝒇(𝒙𝟎 + 𝒉)
𝒉𝟐
Métodos Numéricos
sfs

140
El Error es:

𝒉𝟐 (𝟒)
𝑬= 𝒇 (𝒄)
𝟏𝟐
Con 𝒙𝟎 − 𝒉 < 𝒄 < 𝒙𝟎 + 𝒉

3.2.2 EJERCICIOS DE APRENDIZAJE

1. Para los datos de la tabla, encuentre f '2 f ''2

x 1 1.5 2 2.5 3
y  f x 6 9.5 15 22.5 32

Procedimiento:

Utilizando para la primera derivada la fórmula:

𝟏(
𝟏
𝒇 𝒙𝟎 )−= [𝒇(𝒙𝟎 + 𝒉) − 𝒇(𝒙𝟎 − 𝒉)]
𝟐𝒉

Con 𝒉 = 𝟎, 𝟓 se tiene que:

𝟏
𝒇𝟏 (𝟐) = [𝒇(𝟐 + 𝟎, 𝟓) − 𝒇(𝟐 − 𝟎, 𝟓)] = 𝟐𝟐, 𝟓 − 𝟗, 𝟓
𝟐(𝟎, 𝟓)
𝒇𝟏 (𝟐) = 𝟏𝟑

Utilizando para la segunda derivada la fórmula:


Métodos Numéricos
sfs

141
𝟏
𝒇" (𝒙𝟎 ) = [𝒇(𝒙𝟎 − 𝒉) − 𝟐𝒇(𝒙𝟎 ) + 𝒇(𝒙𝟎 + 𝒉)
𝒉𝟐
Reemplazando los valores conocidos, se tiene que:

𝟏
𝒇" (𝟐) = [𝒇(𝟐 − 𝟎, 𝟓) − 𝟐𝒇(𝟐) + 𝒇(𝟐 + 𝟎, 𝟓)
(𝟎, 𝟓)𝟐
𝟏
𝒇" (𝟐) = [𝒇(𝟏, 𝟓) − 𝟐𝒇(𝟐) + 𝒇(𝟐, 𝟓)
𝟎, 𝟐𝟓
𝒇" (𝟐) = 𝟒[𝟗, 𝟓 − 𝟐(𝟏𝟓) + 𝟐𝟐, 𝟓] = 𝟖

3.2.3 APLICACIONES EN GEOMETRÍA


Encuentre la ecuación de la recta tangente a la curva en el punto:

𝒙=𝟑
La curva corresponde a los datos que se muestran en la tabla.

x 2 2.5 3 3.5 4 4.5


y 3 7.25 13 20.25 29 39.25

Procedimiento:

Aplicando la fórmula de tres puntos se tiene:

𝟏
𝒎 = 𝒇𝟏 (𝟑) = [𝒇(𝟑, 𝟓) − 𝒇(𝟐, 𝟓)] = 𝟏𝟑
𝟐(𝟎, 𝟓)

Para encontrar la ecuación de una recta se utiliza la fórmula:

𝒚 − 𝒚𝟎 = 𝒎(𝒙 − 𝒙𝟎 )
Métodos Numéricos
sfs

142
Donde:

𝒎 = 𝟏𝟑, 𝒙𝟎 = 𝟑 ∧ 𝒚𝟎 = 𝟏𝟑
Reemplazando en la ecuación de la recta se tiene:

𝒚 − 𝟏𝟑 = 𝟏𝟑(𝒙 − 𝟑) ⟹ 𝒚 − 𝟏𝟑 = 𝟏𝟑𝒙 − 𝟑𝟗 ⟹ 𝒚 = 𝟏𝟑𝒙 − 𝟑𝟗 + 𝟏𝟑 ⟹ 𝒚 = 𝟏𝟑𝒙 − 𝟐𝟔

La recta Tangente está dada por: 𝒚 = 𝟏𝟑𝒙 − 𝟐𝟔


1. Aplicaciones en Física:

La tabla muestra la Posición de un objeto en diferentes tiempos:

Tiempo 1 1.2 1.4 1.6 1.8 2


(en segundos)
Posición 13 21.4 31.4 43 56.2 71
(en metros)

Determine velocidad y aceleración para t  1.4 segundos

Procedimiento

La velocidad se obtiene como la primera derivada de la posición:

𝟏
𝒗 = 𝒇𝟏 (𝒙𝟎 )−= [𝒇(𝒙𝟎 + 𝒉) − 𝒇(𝒙𝟎 − 𝒉)]
𝟐𝒉
Con: 𝒙𝟎 = 𝟏, 𝟒 ∧ 𝒉 = 𝟎, 𝟐

Cálculo de la velocidad:

𝟏 𝟏
𝒗 = 𝒇𝟏 (𝟏, 𝟒)−= [𝒇(𝟏, 𝟒 + 𝟎, 𝟐) − 𝒇(𝟏, 𝟒 − 𝟎, 𝟐)] = [𝒇(𝟏, 𝟔) − 𝒇(𝟏, 𝟐)]
𝟐(𝟎, 𝟐) 𝟎, 𝟒

𝒗 = 𝒇𝟏 (𝟏, 𝟒)−= 𝟐, 𝟓(𝟒𝟑 − 𝟐𝟏, 𝟒) = 𝟐, 𝟓(𝟐𝟏, 𝟔) = 𝟓𝟒


Métodos Numéricos
sfs

143
𝑚
𝑣 = 54
𝑠
La aceleración se obtiene como la segunda derivada de la posición:

𝟏
𝒂 = 𝒇" (𝒙𝟎 ) = [𝒇(𝒙𝟎 − 𝒉) − 𝟐𝒇(𝒙𝟎 ) + 𝒇(𝒙𝟎 + 𝒉)
𝒉𝟐
Cálculo de la aceleración:

𝟏
𝒂 = 𝒇" (𝟏, 𝟒) = [𝒇(𝟏, 𝟒 − 𝟎, 𝟐) − 𝟐𝒇(𝟏, 𝟒) + 𝒇(𝟏, 𝟒 + 𝟎, 𝟐)
𝒉𝟐

𝟏
𝒂 = 𝒇" (𝟏, 𝟒) = [𝒇(𝟏, 𝟐) − 𝟐𝒇(𝟏, 𝟒) + 𝒇(𝟏, 𝟔)]
𝒉𝟐

𝒂 = 𝒇" (𝟏, 𝟒) = 𝟐, 𝟓[𝟐𝟏, 𝟒 − 𝟐(𝟑𝟏, 𝟒) + 𝟒𝟑] = 𝟐𝟓 ∗ 𝟏, 𝟔 = 𝟒𝟎


𝒎
𝒂 = 𝒇" (𝟏, 𝟒) = 𝟒𝟎
𝒔𝟐

3.2.4 TALLER DE ENTRENAMIENTO


De acuerdo a lo desarrollado en el tema y los ejercicios de aprendizaje resueltos:
1. Encuentre la ecuación de la recta tangente en el punto x  1 para la curva cuyos datos se
muestran en la tabla.

x 0.6 1.0 1.4 1.8 2.2 2.6


y 0.6 0.5 0.467 0.45 0.44 0.433

2. La tabla muestra la Posición de un objeto en diferentes tiempos:

Tiempo (en segundos) 0.5 1 1.5 2 2.5


Posición (en metros) 300 320 400 450 500

Determine velocidad y aceleración para t  2.5 segundos


Métodos Numéricos
sfs

144
3.3 TEMA 2 INTEGRACIÓN NUMÉRICA
En algunos casos al intentar resolver una integral por los métodos de integración conocidos es
una tarea difícil y en algunas ocasiones no es posible resolver dicha integral, adicionalmente si se
tiene una serie de datos y es necesario integrarlos es muy práctico aplicar algún método numérico
para resolverla.

A continuación, se explica la integración numérica utilizando la regla del trapecio.

REGLA DEL TRAPECIO PARA RESOLVER UNA INTEGRAL DEFINIDA.

El problema a analizar es el siguiente:

Se desea obtener la solución a la integral

𝒃
∫ 𝒇(𝒙)𝒅𝒙
𝒂
Deducción de la fórmula de la regla del trapecio para integrar.

Suponga que la gráfica de 𝑓(𝑥)𝑒𝑛𝑡𝑟𝑒 𝑎 ∧ 𝑏 es la que se ve en la figura 15.

El problema conduce a encontrar el área bajo la curva f (x) entre x  a  x  b

Para ello se divide el área en áreas más pequeñas, y observe, en la figura 16, que las áreas más
pequeñas corresponden a trapecios.
Métodos Numéricos
sfs

145

Figura 15. Gráfica de una función y= f(x) (Autor Elkin Ceballos Gómez abril 18 de 2011)

Figura 16. División del área en n trapecios iguales. (Autor Elkin Ceballos Gómez abril 18 de 2011)
Métodos Numéricos
sfs

146
Se pude observar que el intervalo a, b está dividido en n sub intervalos de igual longitud. Como
la longitud de a, bes 𝑏 − 𝑎 la longitud de cada sub intervalo es:

𝒃−𝒂
El cual se llamará ℎ.
𝒏

Observe que:

𝒙𝟎 = 𝒂
𝒙𝟏 = 𝒂 + 𝒉
𝒙𝟐 = 𝒂 + 𝟐𝒉
𝒙𝟑 = 𝒂 + 𝟑𝒉

𝒙𝒏 = 𝒂 + 𝒏𝒉 = 𝒃
También se puede ver que en cada uno de los sub intervalos se forma un trapecio, el área A de la
b

región limitada por la curva y el eje x entre x = a y x = b es


𝒃
∫ 𝒇(𝒙)𝒅𝒙
𝒂
La cual puede aproximarse por la suma de las áreas de cada trapecio determinado por los sub intervalos
El área de un trapecio se obtiene como: La suma de sus lados paralelos multiplicado por su base y todo
dividido entre dos.
Para cada trapecio de la figura 16 se tiene que:

𝟏
𝑨𝟏 = 𝒉[𝒇(𝒂) + 𝒇(𝒂 + 𝒉)]
𝟐
El segundo trapecio tiene área igual a:
Métodos Numéricos
sfs

147
𝟏
𝑨𝟐 = 𝒉[𝒇(𝒂 + 𝒉) + 𝒇(𝒂 + 𝟐𝒉)]
𝟐
El tercer trapecio tiene área igual a:

𝟏
𝑨𝟑 = 𝒉[𝒇(𝒂 + 𝟐𝒉) + 𝒇(𝒂 + 𝟑𝒉)]
𝟐
El área total es:

𝒃 𝟏
∫𝒂 𝒇(𝒙)𝒅𝒙 ≈ 𝟐h f (a)  2 f a  h 2 f (a  2h)  2 f (a  3h)  ... 2 f (a  (n 1) h)  f
(b)

Otra forma para el área está dada por:


𝒃 𝟏
∫𝒂 𝒇(𝒙)𝒅𝒙 ≈ 𝟐h f (𝒙𝟎 )  2 f 𝒙𝟏  2 f (𝒙𝟐 )  2 f (𝒙𝟑 )  ... 2 f (𝒙𝒏−𝟏 )  f (𝒙𝒏 )

Con:

𝒃−𝒂
𝒉= 𝒏
El error cometido en la regla trapezoidal está dado por:

ℎ3
𝐸𝑎 = | 𝑁𝑓 " (𝑐)| , 𝑐 ∈ [𝑎, 𝑏]
12
Se escoge el valor de c que produzca el máximo error; esto es lo que se llama: ENCONTRAR UNA COTA
PARA EL ERROR. Para utilizar esta fórmula se debe conocer la función a integrar.

3.3.1 TALLER DE APRENDIZAJE

1. Ejemplo tomado de (Haeussler & Richa1rd, 1997

1 1
Use la regla del trapecio para estimar el valor de ∫0 1+𝑥 2
Métodos Numéricos
sfs

148
Con: 𝒂 = 𝟎, 𝒃 = 𝟏 ∧ 𝑵 = 𝟓
PROCEDIMIENTO

Determinación del valor de h:

𝟏−𝟎
𝒉=
𝟓

Los términos a sumar son:

𝑓(𝑎) = 𝑓(0) = 𝟏

1
2𝑓(𝑎 + ℎ) = 2 × = 𝟏, 𝟗𝟐𝟑𝟏
1 + 0,22
1
2𝑓(𝑎 + 2ℎ) = 2 × = 𝟏, 𝟕𝟐𝟒𝟏
1 + 0,42
1
2𝑓(𝑎 + 3ℎ) = 2 × = 𝟏, 𝟒𝟕𝟎𝟔
1 + 0,62
1
2𝑓(𝑎 + 4ℎ) = 2 × = 𝟏, 𝟐𝟏𝟗𝟓
1 + 0,82

2𝑓(𝑏) = 𝑓(1) = 𝟎, 𝟓

La suma de estos resultados es: 7,8373


Por lo tanto:

1
1 1
∫ 2
≈ × 0,2 × 0,78373 = 0,784
0 1 + 𝑥 2

Error cometido:
Métodos Numéricos
sfs

149
𝟐(𝟑𝒙𝟐 − 𝟏)
𝒇"(𝒙) =
(𝟏 + 𝒙𝟐 )𝟑
Reemplazando los valores respectivos, se tiene:

𝟐(𝟑(𝟎)𝟐 − 𝟏)
𝒇"(𝟎) = = −𝟐
(𝟏 + (𝟎)𝟐 )𝟑

𝟐(𝟑(𝟏)𝟐 − 𝟏) 𝟒 𝟏
𝒇"(𝟏) = = =
(𝟏 + (𝟏)𝟐 )𝟑 𝟖 𝟐

(1 − 0)3 " 𝟏𝟑
𝐸𝑎 = | 𝑓 (0)| < × 𝟐 = 𝟔, 𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟒 × 𝟏𝟎−𝟑 = 𝟎, 𝟎𝟎𝟔𝟔𝟔𝟔 …
12 × 52 12 × 52

Otra forma de calcular el error:

ℎ3 (𝟎, 𝟐)𝟑
𝐸𝑎 = | 𝑁𝑓 " (𝑐)| < × 𝟓 × 𝟐 = 𝟎, 𝟎𝟎𝟔𝟔𝟔𝟔𝟔𝟔𝟔
12 12 ×

El valor verdadero de la integral es:tan−1 1 = 0,78539816


Error verdadero es: 0,007025163

1. Ejemplo tomado de (Purcell & Varberg, 1993)

Use la regla trapezoidal con 𝒏 = 𝟔, para dar el valor aproximado de :

𝟏
𝟐
∫ 𝒆−𝒙 𝒅𝒙
𝟎
Procedimiento

Datos:
Métodos Numéricos
sfs

150
𝟐
𝒇(𝒙) = 𝒆−𝒙 , 𝒂 = 𝟎, 𝒃 = 𝟏, 𝒏 = 𝟔
1−0 1
ℎ= = = 0,1667
6 6
𝟐
𝒇(𝒂) = 𝒇(𝟎) = 𝒆−(𝟎) = 𝟏
𝟐
𝟐 × 𝒇(𝒂 + 𝒉) = 𝒇(𝟎, 𝟏𝟔𝟔𝟕) = 𝟐 × 𝒆−(𝟎,𝟏𝟔𝟔𝟕) = 𝟏, 𝟗𝟒𝟓𝟐
𝟐
𝟐 × 𝒇(𝒂 + 𝟐𝒉) = 𝒇(𝟎, 𝟑𝟑𝟑𝟑) = 𝟐 × 𝒆−(𝟎,𝟑𝟑𝟑𝟑) = 𝟏, 𝟕𝟖𝟗𝟕
𝟐
𝟐 × 𝒇(𝒂 + 𝟑𝒉) = 𝒇(𝟎, 𝟓) = 𝟐 × 𝒆−(𝟎,𝟓) = 𝟏, 𝟓𝟓𝟕𝟔
𝟐
𝟐 × 𝒇(𝒂 + 𝟒𝒉) = 𝒇(𝟎, 𝟔𝟔𝟔𝟕) = 𝟐 × 𝒆−(𝟎,𝟔𝟔𝟔𝟕) = 𝟏, 𝟐𝟖𝟐𝟒
𝟐
𝟐 × 𝒇(𝒂 + 𝟓𝒉) = 𝒇(𝟎, 𝟖𝟑𝟑𝟑) = 𝟐 × 𝒆−𝟎,𝟖𝟑𝟑𝟑) = 𝟎, 𝟗𝟗𝟖𝟖
𝟐
𝒇(𝒃) = 𝒇(𝟏) = 𝒆−(𝟏) = 𝟎, 𝟑𝟔𝟕𝟗
La suma de estos valores es: 𝟖, 𝟗𝟒𝟏𝟔

El valor de la integral es:

𝟏
−𝒙𝟐
𝟏 𝟏
∫ 𝒆 𝒅𝒙 = × × 𝟖, 𝟗𝟒𝟏𝟔 = 𝟎, 𝟕𝟒𝟓𝟏𝟑𝟑𝟑𝟑 ≈ 𝟎, 𝟕𝟒𝟓𝟏
𝟎 𝟐 𝟔
El error cometido en la regla trapezoidal está dado por:

ℎ3
𝐸𝑎 = | 𝑁𝑓 " (𝑐)| , 𝑐 ∈ [𝑎, 𝑏]
12
Estime el máximo error posible para el ejemplo anterior

𝟐 𝟐 𝟐
𝒇(𝒙) = 𝒆−𝒙 , 𝒇𝟏 (𝒙) = −𝟐𝒙𝒆−𝒙 ∧ 𝒇" (𝒙) = 𝒆−𝒙 (𝟒𝒙𝟐 − 𝟐)
Se tiene que:
Métodos Numéricos
sfs

151
" −𝒙 𝟐
|𝒇 (𝒙)| = |𝒆 (𝟒𝒙𝟐 − 𝟐)|

𝟐
|𝒇" (𝟎)| = |𝒆−𝟎 (𝟒(𝟎)𝟐 − 𝟐)| ≤ |𝟏(𝟎) − 𝟐| = |−𝟐| = 𝟐
𝟐
|𝒇" (𝟏)| = |𝒆−𝟏 (𝟒(𝟏)𝟐 − 𝟐)| = |𝟎, 𝟑𝟔𝟕𝟖𝟕𝟗𝟒𝟒𝟏(𝟎 − 𝟐)| = |𝟎, 𝟑𝟔𝟕𝟖𝟕𝟗𝟒𝟒𝟏(−𝟐)|
= 𝟎, 𝟕𝟑𝟓𝟕𝟖𝟖𝟖𝟐

Entonces:

𝟏 𝟑
( 𝟔)
|𝑬| ≤ × 𝟔 × 𝟐 ⇒ |𝑬| = 𝟎, 𝟎𝟎𝟒𝟔𝟐𝟗𝟔𝟗 …
𝟏𝟐

2. Utilizando la regla del trapecio con 𝑛 = 8 , halle el valor de:

𝟏 𝟐 − 𝟏 𝒛𝟐
∫ 𝒆 𝟐 𝒅𝒛
𝟐𝝅 −𝟐
Procedimiento

𝟏 𝟐 − 𝟏 𝒛𝟐
∫ 𝒆 𝟐 𝒅𝒛
𝟐𝝅 −𝟐
𝟐 − (−𝟐) 𝟒
𝒉= = = 𝟎, 𝟓
𝟖 𝟖

𝒙𝟎 = −𝟐

𝒙𝟏 = −𝟐 + 𝟎, 𝟓 = −𝟏, 𝟓

𝒙𝟐 = −𝟏, 𝟓 + 𝟎, 𝟓 = −𝟏

𝒙𝟑 = −𝟏 + 𝟎, 𝟓 = −𝟎, 𝟓

𝒙𝟒 = −𝟎, 𝟓 + 𝟎, 𝟓 = 𝟎
Métodos Numéricos
sfs

152
𝒙𝟓 = 𝟎 + 𝟎, 𝟓 = 𝟎, 𝟓

𝒙𝟔 = 𝟎, 𝟓 + 𝟎, 𝟓 = 𝟏

𝒙𝟕 = 𝟏 + 𝟎, 𝟓 = 𝟏, 𝟓

𝒙𝟑 = 𝟏, 𝟓 + 𝟎, 𝟓 = 𝟐
Métodos Numéricos
sfs

153

3.3.2 TALLER DE ENTRENAMIENTO


Encuentre el área de la siguiente figura (fig.17)

Para ello tome: 𝒂 = 𝟐𝟎𝒄𝒎. , 𝒃 = 𝟑𝟎𝒄𝒎. 𝒚 𝒄 = 𝟏𝟓𝒄𝒎.

FIGURA 17

𝟓 𝟑𝒙
2. Resuelva ∫𝟎 𝒙𝟐 +𝟒
a. De manera exacta, aplicando el método de integración por sustitución.

b. De manera aproximada, utilizando la regla del trapecio con 𝒏 = 𝟏𝟎 𝒕𝒓𝒂𝒑𝒆𝒄𝒊𝒐𝒔.

3.4 TEMA 3 SOLUCIÓN NUMÉRICA DE ECUACIONES


DIFERENCIALES
 Video: METODO DE EULER, EULER MODIFICADO Y RUNGE-KUTTA
 ”Ver Video”: https://www.youtube.com/watch?v=nZ1aYtHusLE
Métodos Numéricos
sfs

154
El problema a resolver en esta sección es que:
Dada una ecuación diferencial de la forma:

𝒅𝒚
= 𝒇(𝒙, 𝒚)
𝒅𝒙
Con la condición inicial:

𝒚(𝒙𝟎 ) = 𝒚𝟎
Halle el valor que da la función f(x) para otro valor dado de x, es decir, halle:

𝑦(𝑥𝑑𝑎𝑑𝑜 )
Para solucionar este problema que se llama P.V.I. (Problema de valor inicial) se analizarán dos
métodos.

3.4.1 MÉTODO DE EULER (O MÉTODO DE LA RECTA TANGENTE)


3.4.2 EJERCICIO DE APRENDIZAJE
Solucione la ED:

𝒅𝒚
= 𝒇(𝒙, 𝒚)
𝒅𝒙
Con la condición inicial:

𝒚(𝒙𝟎 ) = 𝒚𝟎
Se tienen los siguientes datos:

En

𝒙 = 𝒙𝟎 , 𝒚𝟎 = 𝒚(𝒙𝟎 )
Como:
Métodos Numéricos
sfs

155
𝒅𝒚
= 𝒇(𝒙𝟎 , 𝒚𝟎 ) en 𝒙 = 𝒙𝟎 ∧ 𝒚 = 𝒚𝟎
𝒅𝒙

𝒅𝒚 𝒚𝟏 − 𝒚𝟎
𝒎= = 𝒇(𝒙𝟎 , 𝒚𝟎 ) =
𝒅𝒙 𝒙𝟏 − 𝒙𝟎
Tomando la igualdad:

𝒚𝟏 − 𝒚𝟎
𝒇(𝒙𝟎 , 𝒚𝟎 ) =
𝒙𝟏 − 𝒙𝟎
Despejando 𝒚𝟏 , quedaría:
𝒚𝟏 = 𝒚𝟎 + 𝒇(𝒙𝟎 , 𝒚𝟎 )(𝒙𝟏 − 𝒙𝟎 )

Con 𝒚𝟏 determinado se puede calcular 𝒎𝟐 = 𝒇(𝒙𝟏 , 𝒚𝟏 )


𝒚𝟐 − 𝒚𝟏
𝒎 = 𝒇(𝒙𝟏 , 𝒚𝟏 ) =
𝒙𝟐 − 𝒙𝟏

𝒚𝟐 = 𝒚𝟏 + 𝒇(𝒙𝟏 , 𝒚𝟏 )(𝒙𝟐 − 𝒙𝟏 )

En forma general:

𝒚𝒊+𝟏 = 𝒚𝒊 + 𝒇(𝒙𝒊 , 𝒚𝒊 )(𝒙𝒊+𝟏 − 𝒙𝒊 )


Sea: 𝒉 = 𝒙𝒊+𝟏 − 𝒙𝒊 , Valor constante, la ecuación queda:

𝒚𝒊+𝟏 = 𝒚𝒊 + 𝒇(𝒙𝒊 , 𝒚𝒊 )𝒉
Para: 𝒊 = 𝟏, 𝟐, 𝟑, … , 𝒏 − 𝟏
Llamada Fórmula de Euler.
Métodos Numéricos
sfs

156
Enlaces método de Euler.

 Video: Método de euler,kuta


 ”Ver Video”: https://www.youtube.com/watch?v=nZ1aYtHusLE

3.4.3 MÉTODO MEJORADO (O RUG KUTA DE SEGUNDO ORDEN) (O


MÉTODO PREDICTOR – CORRECTOR DE EULER)
 Video: Euler Mejorado
 ”Ver Video”: https://www.youtube.com/watch?v=eV4DLPXXp7w

Se utiliza el método de Euler para obtener un valor predictor:

𝒚∗𝑰+𝟏 = 𝒚𝒊 + 𝒇(𝒙𝒊 , 𝒚𝒊 )𝒉
Este valor se ingresa al corrector, que es la fórmula de EULER MEJORADO.

𝒇(𝒙𝒊 , 𝒚𝒊 ) + 𝒇(𝒙𝒊+𝟏 , 𝒚𝒊+𝟏 )


𝒚𝒊+𝟏 = 𝒚𝒊 + 𝒉
𝟐

3.4.4 EJERCICIOS DE APRENDIZAJE


1. Ejemplo tomado de (Asmar, 1999)

𝒅𝒚 𝟏 𝟐
Para el P.V.I: = (𝒚 + 𝒚) con condición inicial:
𝒅𝒙 𝒙

𝒚(𝟏) = −𝟐 , 𝒉𝒂𝒍𝒍𝒆 𝒚(𝟏, 𝟔), con un tamaño de paso 𝒉 = 𝟎, 𝟐,


utilizando el método de EULER MEJORADO.

Procedimiento
Métodos Numéricos
sfs

157
Con 𝒉 = 𝟎, 𝟐, 𝒙𝟎 = 𝟏 ∧ 𝒚𝟎 = −𝟐

Se tiene que:

𝒙𝟏 = 𝟏, 𝟐; 𝒙𝟐 = 𝟏, 𝟒; 𝒙𝟑 = 𝟏, 𝟔
𝟏
𝒚𝟏 = 𝒚𝟎 + 𝒇(𝒙𝟎 , 𝒚𝟎 )𝒉 ⇒ 𝒚𝟏 = −𝟐 + [ ((−𝟐)𝟐 + (−𝟐))] 𝟎, 𝟐 = 𝟏, 𝟔
𝟏
𝒇(𝒙𝟎 , 𝒚𝟎 ) + 𝒇(𝒙𝟏 , 𝒚𝟏 )
𝒚𝟏 = 𝒚𝟎 + 𝒉⇒
𝟐
𝟏 𝟐 𝟏 𝟐
((−𝟐) + (−𝟐)) + (−𝟏, 𝟔))
𝟏 𝟏, 𝟐 ((−𝟏, 𝟔) +
𝒚𝟏 = −𝟐 + × 𝟎, 𝟐 ⇒
𝟐

𝒚𝟏 = −𝟏, 𝟕𝟐 = 𝒚(𝟏, 𝟐)
𝟏
𝒚∗𝟐 = 𝒚𝟏 + 𝒇(𝒙𝟏 , 𝒚𝟏 )𝒉 ⇒ 𝒚∗𝟏 = −𝟏, 𝟕𝟐 + [ ((−𝟏, 𝟕𝟐)𝟐 + (−𝟏, 𝟕𝟐))] 𝟎, 𝟐 = −𝟏, 𝟓𝟏𝟑𝟔
𝟏, 𝟐

𝒇(𝒙𝟏 , 𝒚𝟏 ) + 𝒇(𝒙𝟐 , 𝒚∗𝟐 )


𝒚𝟐 = 𝒚𝟏 + 𝒉
𝟐
𝟏 𝟐 𝟏 𝟐
(−𝟏, 𝟕𝟐)) + (−𝟏, 𝟓𝟏𝟑𝟔))
𝟏, 𝟐 ((−𝟏, 𝟕𝟐) + 𝟏, 𝟒 ((−𝟏, 𝟓𝟏𝟑𝟔) +
𝒚𝟐 = −𝟏, 𝟕𝟐 + × 𝟎, 𝟐
𝟐

𝒚𝟐 = −𝟏, 𝟓𝟔𝟏𝟐𝟕𝟐𝟓 = 𝒚(𝟏, 𝟒)


𝟏
𝒚∗𝟑 = −𝟏, 𝟓𝟔𝟏𝟐𝟕𝟐𝟓 + [ ((−𝟏, 𝟓𝟔𝟏𝟐𝟕𝟐𝟓)𝟐 + (−𝟏, 𝟓𝟔𝟏𝟐𝟕𝟐𝟓))] 𝟎, 𝟐 ⇒
𝟏, 𝟒

𝒚∗𝟑 = −𝟏, 𝟒𝟑𝟔𝟎𝟖𝟔𝟖𝟖𝟑


𝟏 𝟏
((−𝟏, 𝟓𝟔𝟏𝟐𝟕𝟐𝟓)𝟐 + (−𝟏, 𝟓𝟔𝟏𝟐𝟕𝟐𝟓)) + ((−𝟏, 𝟒𝟑𝟔𝟎𝟖𝟔𝟖𝟖𝟑)𝟐 + (−𝟏, 𝟒𝟑𝟔𝟎𝟖𝟔𝟖𝟖𝟑))
𝟏, 𝟒 𝟏, 𝟔
𝒚𝟑 = −𝟏, 𝟓𝟔𝟏𝟐𝟕𝟐𝟓 + × 𝟎, 𝟐
𝟐

𝒚𝟑 = −𝟏, 𝟒𝟓𝟗𝟓𝟑𝟖𝟓 = 𝒚(𝟏, 𝟔)


Métodos Numéricos
sfs

158
𝒅𝒚
1. Para el P.V.I: = 𝟐𝒙𝒚 + 𝒚 𝒄𝒐𝒏 𝒚(𝟐) = 𝟐
𝒅𝒙

Use la fórmula de Euler mejorado con tamaños de paso de h  0.2 para determinar un valor
aproximado de 𝒚(𝟐, 𝟒)

Procedimiento:

h  0,2; 𝒙𝟎  2  𝒚𝟎  2
𝒙𝟏 = 𝟐 + 𝟎, 𝟐 = 𝟐, 𝟐
𝒙𝟐 = 𝟐, 𝟐 + 𝟎, 𝟐 = 𝟐, 𝟒
Cálculo del Predictor:

𝒚∗𝟏 = 𝒚𝟎 + 𝒇(𝒙𝟎 , 𝒚𝟎 )𝒉

𝒚∗𝟏 = 𝒚𝟎 + 𝒇(𝒙𝟎 , 𝒚𝟎 )𝒉 ⟹ 𝒚𝟏∗ = 𝟐 + [𝟐(𝟐)(𝟐) + (𝟐)]𝟎, 𝟐 = 𝟒

Cálculo del Corrector:

𝒇(𝒙𝟎 , 𝒚𝟎 ) + 𝒇(𝒙𝟏 , 𝒚𝟏 )
𝒚𝟏 = 𝒚𝟎 + 𝒉
𝟐
[𝟐(𝟐)(𝟐) + 𝟐] + [𝟐(𝟐, 𝟐)(𝟒) + 𝟒]
𝒚𝟏 = 𝟐 + 𝟎, 𝟐
𝟐

𝒚𝟏 = 𝟓, 𝟏𝟔 = 𝒚(𝟐, 𝟐)
Cálculo del Predictor:

𝒚∗𝟐 = 𝒚𝟏 + 𝒇(𝒙𝟏 , 𝒚𝟏 )𝒉

𝒚∗𝟐 = 𝟓, 𝟏𝟔 + [𝟐(𝟐, 𝟐)(𝟓, 𝟏𝟔) + (𝟓, 𝟏𝟔)]𝟎, 𝟐 = 𝟏𝟎, 𝟕𝟑𝟐𝟖


Métodos Numéricos
sfs

159
Cálculo del Corrector:

𝒇(𝒙𝟏 , 𝒚𝟏 ) + 𝒇(𝒙𝟐 , 𝒚𝟐 )
𝒚𝟐 = 𝒚𝟏 + 𝒉
𝟐
[𝟐(𝟐, 𝟐)(𝟓, 𝟏𝟔) + 𝟓, 𝟏𝟔] + [𝟐(𝟐, 𝟒)(𝟏𝟎, 𝟕𝟑𝟐𝟖) + 𝟏𝟎, 𝟕𝟑𝟐𝟖]
𝒚𝟐 = 𝟓, 𝟏𝟔 + 𝟎, 𝟐
𝟐

𝟐𝟕, 𝟖𝟔𝟒 + 𝟔𝟐, 𝟐𝟓𝟎𝟐𝟒 𝟗𝟎, 𝟏𝟏𝟒𝟐𝟒


𝒚𝟐 = 𝟓, 𝟏𝟔 + 𝟎, 𝟐 = 𝟓, 𝟏𝟔 + 𝟎, 𝟐
𝟐 𝟐

𝒚𝟐 = 𝟏𝟒, 𝟏𝟕𝟏𝟒𝟐𝟒 = 𝒚(𝟐, 𝟒)


En los siguientes enlaces

 Video: Euler Mejorado


 ”Ver Video”: https://www.youtube.com/watch?v=eV4DLPXXp7w

 Video: METODOS NUMERICOS EULER MODIFICADO


 ”Ver Video”: https://www.youtube.com/watch?v=1j-Um-p06xo

 Video: Euler's Method: Example


 ”Ver Video”: https://www.youtube.com/watch?v=4CqaepeaJHA

3.4.5 TALLER DE ENTRENAMIENTO


Recuerda que para cada uno de los problemas indicados a continuación debes
realizar el procedimiento adecuado, de acuerdo con la temática desarrollada a lo
largo de la unidad.
1. Para el problema de valor inicial:
Métodos Numéricos
sfs

160
𝑑𝑦 𝑥
= 𝑐𝑜𝑛 𝑦(2) = 3
𝑑𝑥 𝑦

Halle: 𝑦(2) 𝑐𝑜𝑛 ℎ = 0,25

2. Para el problema de valor inicial:

𝑑𝑦 𝑥
= 𝑐𝑜𝑛 𝑦(2) = 4
𝑑𝑥 𝑦

Halle: 𝑦(2,5) 𝑐𝑜𝑛 ℎ = 0,1

Las preguntas 3 y 4 se resuelven con la siguiente integral:

𝟔
𝟐𝒙
∫ 𝟐 𝒅𝒙
𝟎 𝒙 + 𝟏
3. Aproximando a cuatro cifras decimales, el resultado exacto de la integral es:

a. 3.5835

b. 3.6109

c. 3.4174

d. 6.8348

4. Aplicando la regla del trapecio con 𝑛 = 6 y dando la respuesta final aproximada con
cuatro decimales, el resultado obtenido es:

a. 3.5835

b. 2.6333

c. 1.7918

d. 3.4174

5. Los datos de la tabla corresponden al tiempo en horas y a la distancia en kilómetros


recorrida por un objeto.
Métodos Numéricos
sfs

161
X Y
(Tiempo) (distancia)
2 -46
2.5 -1,875
3 62
3.5 149,375

Utilizando la diferenciación numérica encuentre la velocidad y la aceleración para un tiempo de 2.5


horas.
6. Para el P.V.I: y' ex, y (0)  0, h  0.5, el valor y (2) aproximando el último cálculo a cinco
2

cifras decimales, el resultado es:

a. y (2)  1.36218

b. y (2)  1.13520

c. y (2)  1.25996

d. y (2)  1.08250

e. y (2)  1.12604
Métodos Numéricos

162
4 PISTAS DE APRENDIZAJE

No olvide Que para encontrar la raíz de una ecuación no lineal f x  0 , Se sabe que hay una
raíz en un intervalo a, bporque: La función f x es continua en a, by porque f a f b 0

Traer a la memoria Que para hallar un polinomio de interpolación a un conjunto de parejas:

x0 , y0 , x1 , y1 , x2 , y2 ,...,xn , yn 

Los valores de x deben ser todos diferentes.

Traer a la memoria La división entre cero no está definida en ningún campo numérico. Cuando
en el numerador hay un número diferente de cero y en el denominador está el cero se dice que
el resultado no existe; si en el numerador y en el denominador está el cero, se dice que el
resultado es indefinido.

Tener en cuenta Cuando se desea interpolar un conjunto de puntos de la forma:

x0 , y0 , x1 , y1 , x2 , y2 ,...,xn , yn 

El grado máximo del polinomio de interpolación debe ser menor que el número de puntos.

Tener en cuenta Que cuando se realiza la interpolación por el método de LaGrange y por el
método de Newton, el polinomio obtenido por ambos métodos debe ser el mismo.
Tenga Presente Que, al utilizar el método del Punto Fijo para solucionar una ecuación no lineal,
para que la

función gx, encontrada por usted, funcione, se debe cumplir que |𝑔1 (𝑥)| < 1 en el intervalo
donde encuentra la raíz

Traer a la memoria Que en física si f xes la función de posición de un objeto, la velocidad instantánea del
objeto se halla como v  f 'x y función de aceleración instantánea se obtiene como a  f ''x.
Métodos Numéricos
sfs

163
Traer a la memoria Para encontrar la ecuación de una recta conocidos; un punto de coordenadas
x0 , y0 y la pendiente m de la recta, se utiliza la ecuación punto pendiente de la línea recta
cuya forma es:

y  y0  mx  x0 

Traer a la memoria Si f xes una función, la pendiente de la recta tangente a la curva en un


punto de coordenadas x0 , y0  , se obtiene como:

m  f ' x0 
Tenga presente El orden en que se efectúa operaciones es: Primero potencias o raíces, luego
multiplicaciones o divisiones y por último sumas y restas.

Tenga presente Al utilizar la calculadora para resolver una ecuación lineal que involucra funciones
trigonométricas, la calculadora se debe llevar al modo radianes
Métodos Numéricos

164
5 GLOSARIO
 Métodos numéricos: Son técnicas que permiten solucionar problemas utilizando operaciones
repetitivas de suma, resta, multiplicación y división.

 Interpolación: Dado un conjunto de valores 𝑥, 𝑦 encontrar una función que represente estos
valores, que pase exactamente por cada uno de ellos y que permita estimar valores no observados
(es decir, que permita hacer predicciones o interpolaciones o aproximaciones).

 Regresión: Dado un conjunto de parejas 𝑥, 𝑦 encontrar un una función que represente el


conjunto de datos, en la regresión la función encontrada, necesariamente no pasa por todos los
puntos.
 Dominio: El dominio de una función 𝑦 = 𝑓(𝑥) está formado por el conjunto de todos los
números
 reales que puede tomar la variable x.
 Rango: El rango de una función 𝑦 = 𝑓(𝑥) está formado por el conjunto de todos los números
reales que puede tomar la variable y.
 Continuidad: Intuitivamente una función es continua cuando al observar o recorrer su gráfica,
esta no presenta huecos, espacios, separaciones; si algo de esto sucede, se dice que la función es
discontinua.

 Semilla: Es una estimación inicial de la raíz de una ecuación de la forma f x  0. Es el valor
inicial con el cual se empieza a operar el método numérico que se va a aplicar.

 Punto fijo: Dada una función y=f xun punto fijo es un punto en el cual 𝑦 = 𝑥
 Convergencia: Se dice que un método converge cuando se llega al valor buscado, esto se puede
identificar por qué se puede observar que cada cálculo efectuado lleva a un valor cada vez más
parecido al anterior.

 Divergencia: Un método diverge cuando no se llega al valor buscado, esto se identifica, porque
cada cálculo efectuado lleva a un valor cada vez diferente.
 Asíntota vertical: Es un valor de x en el cual el denominador de una función se hace igual a cero.
Métodos Numéricos

165
6 BIBLIOGRAFÍA
 Asmar, I. (1999). Métodos numéricos un 
primer curso. Medellín. http://books.google.com.co/books?id=65J
SdxLG4M4C&printsec=frontcover&dq=m%
 Burden, & Faires, D. (2003). Análisis C3%A9todos
numérico. México: International
 +num%C3%A9ricos&source=bl&ots=oe_U
Thomson Editores, S.A. 1iXJDU&sig=rfk15RFvg6S3O1Vx8wvgZXFC
8sE&hl=es&ei=
 Chapra, S., & Canale, R. (2003). Métodos a_qLTeuuKsaCtgfahsm1DQ&sa=X&oi=boo
numéricos para ingenieros. . México D.F.: k_result&ct=result&resnum=10&ved=0CG
Mc Graw-Hill Interamericana. AQ6AEwCQ#v
 =onepage&q&f=false
 Haeussler, E., & Richard, P. S. (1997).
 Mendoza, F. Matemáticas para ingenieros.
Matemáticas para Administración,
Recuperado el 24 de marzo de 2011 de,
Economía, Ciencias sociales y de la vida. http://matematicas5.galeon.com/index.htm
México: Prentice hall. l
 Métodos numéricos. Recuperado
 Purcell, E., & Varberg, D. (1993). Cálculo
el 24 de
con geometría analítica. México: Prentice
marzo de 2011,
Hall.
de
 Asmar,I. (2000). Métodos numéricos. http://www.mitecnologico.com/Main/Met
Recuperado el 24 odosNumericos
de marzo de 2011,
 Quintana, P, Villalobos, E “y” Cornejo, M.
de
http://www.unalmed.edu.co/~ifasmar/int (2005). Métodos numéricos con
ro.html aplicaciones en Excel. Recuperado el 24 de
marzo de 2011 de,
 Díaz,W. (1998). Métodos numéricos.
http://books.google.com.co/books?id=wX
Recuperado el 24
de marzo de 2011, P3VzSa-
de http://www.uv.es/diaz/mn/fmn.html HIC&printsec=frontcover&dq=m%C3%A9to
dos+num%C3%A9ricos&source=bl&ots=j9J
4hx8T_Q&si
 Huerta, A, Sarrate, J “y” Rodríguez, A. g=BnEOwJ1CgzurjEudVbSP2qAdTJk&hl=es
(2001). Métodos numéricos introducción,
&ei=a_qLTeuuKsaCtgfahsm1DQ&sa=X&oi
aplicaciones y programación.
Recuperado el =book_result&
24 de ct=result&resnum=9&ved=0CF0Q6AEwCA#
marzo de v=onepage&q&f=false
2011 de,

También podría gustarte