TECNOLOGICO DE ESTUDIAOS SUPERIORES DE
HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
“Métodos de solución de ecuaciones. “
Alumno: Herrera Galindo Brenda Berenice
Matricula: 19090013
Profesor: Ing. Víctor Esteban Santiago Trejo
Huixquilucan Edo. De México a 13 de marzo de 2021
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
Introducción
La resolución de un sistema de ecuaciones lineales aparece frecuentemente en la
rama de las matemáticas e ingeniería siendo un subproblema de un problema más
complicado del análisis numérico, por ejemplo, cuando se resuelve iterativamente
un sistema de ecuaciones no lineales por el método de Newton-Raphson, donde en
cada etapa de ese proceso iterativo se requiere resolver un sistema de ecuaciones
lineales, o en procesos de optimización tanto lineales como no lineales. Los
sistemas de ecuaciones tienen una estructura muy especial que puede ser utilizados
para la aplicación de un problema en particular. En la presente investigación se
muestra una pequeña clasificación de la solución de los problemas de ecuaciones
lineales siendo métodos directos y métodos iterativos al mismo tiempo presentan
los métodos de cada clasificación.
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
Objetivo general
Uno de los propósitos en general de la presente investigación es forzar el estudio e
investigación, permitiéndonos entender su teoría y aplicación; es decir, resolver el
sistema. El objetivo de este tema es la resolución de un sistema de 𝑛 ecuaciones
lineales con 𝑛 incógnitas donde son conocidos la matriz de coeficientes y el vector
de términos independientes puesto en notación matricial, el sistema de ecuaciones
lineales se escribe: 𝐴𝑥 = 𝑏 y se puede resolver por dos tipos de métodos siendo
los directos e iterativos permitiéndonos encontrar el resultado con distintos métodos.
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
Desarrollo
Antecedentes Históricos
Una muestra de las matemáticas, más específicamente de la aritmética y del álgebra,
que se desarrollaron en Egipto se encuentran en el Papiro de Rhind, escrito por el escriba
Ahmés, hacia 1650 a.C. Este documento que es una copia de otra más antigua (2000 –
1800 a.C) arroja evidencia del uso de ecuaciones lineales y más aún de sistemas de
ecuaciones simultáneas, no necesariamente lineales por otro lado, los babilonios
llamaban a las incógnitas con palabras tales como longitud, anchura, área, o volumen,
sin que tuvieran relación con problemas de medida. Los griegos también resolvían
algunos sistemas de ecuaciones, pero utilizando métodos geométricos. Thymaridas (400
a. de C.) había encontrado una fórmula para resolver un determinado sistema de 𝑛
ecuaciones con 𝑛 incógnitas. Diophante resuelve también problemas en los que
aparecían sistemas de ecuaciones, pero transformándolos en una ecuación lineal. Los
sistemas de ecuaciones aparecen también en los documentos indios. No obstante, no
llegan a obtener métodos generales de resolución, sino que resuelven tipos especiales
de ecuaciones. En el libro El arte matemático, de un autor chino desconocido (siglo III a.
de C.), contiene algunos problemas donde se resuelven ecuaciones. En ellos
encontramos un esbozo del método de las matrices para resolver sistemas de ecuaciones
lineales. Uno de dichos problemas equivale a resolver un sistema de tres ecuaciones
lineales por dicho método matricial.
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
Antecedentes Matemáticos
La primera fase, que comprende el periodo de 1700 a. de C. a 1700 d. de C., se
caracterizó por la invención gradual de símbolos y la resolución de ecuaciones. Dentro
de esta fase encontramos un álgebra desarrollada por los griegos (300 a. de C.), llamada
álgebra geométrica, rica en métodos geométricos para resolver ecuaciones algebraicas.
La introducción de la notación simbólica asociada a Viéte (1540-1603), marca el inicio de
una nueva etapa en la cual Descartes (1596-1650) contribuye de forma importante al
desarrollo de dicha notación. En este momento, el álgebra se convierte en la ciencia de
los cálculos simbólicos y de las ecuaciones. Posteriormente, Euler (1707-1783) la define
como la teoría de los "cálculos con cantidades de distintas clases" (cálculos con números
racionales enteros, fracciones ordinarias, raíces cuadradas y cúbicas, progresiones y
todo tipo de ecuaciones).
Para llegar al actual proceso de resolución de la ecuación 𝑎𝑥 + 𝑏 = 𝑐 han pasado más
de 3.000 años.
Aceptando el simbolismo actual x 2 , se plantearon problemas que implicaban la solución
de sistemas como
x2 +y2 =100
3
𝑦= 𝑥.
4
𝑥
En tiempos modernos el problema pediría resolver la ecuación 𝑥 + 7 = 19. Para este
problema, se toma 7 como si fuera la solución, quizás porque le permite realizar de
manera más sencilla los cálculos, y a partir de este supuesto encuentra el valor real.
7
Calcula 7 + = 8 ≠ 19 , la cual es una falsa posición. Posteriormente, dice que hay
7
que encontrar un número de tal forma que al multiplicarlo por el resultado de aplicar el
19
valor remplazado resulte 19, es decir se necesitaría obtener .
8
Los sistemas de ecuaciones lineales.
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
Un sistema de m ecuaciones lineales con n incógnitas es un conjunto de 𝑚 relaciones
lineales entre 𝑛 incógnitas 𝑥1, 𝑥2, … , 𝑥𝑛, es decir:
donde los elementos 𝑎𝑖𝑗 , 𝑐𝑜𝑛 𝑖 = 1, … 𝑚 𝑦 𝑗 = 1, … , 𝑛. reciben el nombre de coeficientes
del sistema de ecuaciones lineales. y 𝑏i el de términos independientes del sistema de
ecuaciones lineales. Tanto los coeficientes como los términos independientes son en
principio escalares pertenecientes a un cuerpo determinado. En esta asignatura nos
restringiremos al caso en el que tanto unos como otros son números reales (y, por tanto,
las incógnitas serán también números reales, en este caso desconocidos en principio).
Definición: Una solución del S.E.L. es toda colección de números reales: 𝛼1, 𝛼2, . . . , 𝛼𝑛
tal que al ser sustituidos en el sistema: 𝑥1 = 𝛼1, 𝑥2 = 𝛼2, … , 𝑥𝑛 = 𝛼𝑛, lo reducen a 𝑚
identidades.
Utilidad de los sistemas de ecuaciones lineales
La utilidad radica en analizar los problemas ficticios y de la vida real en que las
ecuaciones lineales nos proporcionan métodos donde buscamos despejar una incógnita.
Un ejemplo fácil en la vida diaria podemos usar las ecuaciones lineales para determinar
la cantidad de integrantes que deben pasar a un banco, dada una relación con otros
servicios, mientras hacen operaciones en cualquier cajero automático, para determinar
el tiempo de espera que hicimos, entre otros nos permite estudiar costos, ofertas,
demandas, gastos, consumos como el consumo y costo de la electricidad se puede medir
con una ecuación lineal la cual nos puede ayudar a predecir los pagos que realizaremos
o cómo disminuir los gastos por el consumo.
Sistema de Ecuaciones Lineales
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
Un sistema de ecuaciones lineales es compatible si posee alguna solución.
Dentro de los sistemas compatibles, podemos distinguir:
• Sistema de ecuaciones lineales compatible determinado: Se trata del caso en el
que existe una única solución (𝑥1, . . . , 𝑥𝑛) del sistema.
• Sistema de ecuaciones lineales compatible indeterminado: Se trata de los
sistemas que poseen más de una solución.
I. Un sistema de ecuaciones lineales es incompatible si no tiene soluciones.
II. Un sistema de ecuaciones lineales es homogéneo si todos sus términos
independientes son nulos, es decir:
III. En caso contrario se dice que el sistema es no-homogéneo.
Propiedades básicas
• Un sistema de ecuaciones lineales homogéneo siempre admite la solución: 𝑥1 =
𝑥2 = … = 𝑥𝑛 = 0, que recibe habitualmente el nombre de solución trivial. Por
lo tanto, los sistemas homogéneos son siempre compatibles.
• El conjunto de soluciones de un sistema homogéneo forma un subespacio
vectorial del especio ℝn
Por esta razón, si un sistema homogéneo es compatible indeterminado, admitirá
necesariamente infinitas soluciones. Si es compatible determinado, obviamente admitirá
sólo la solución trivial.
Esta proposición se convierte en obvia si interpretamos el sistema como una aplicación
lineal 𝑓: ℝn → ℝm , 𝑐𝑜𝑛𝑓(𝑥⃗) = ⃗0⃗. Las soluciones son entonces el núcleo de 𝑓, Ker 𝑓, es
decir un subespacio vectorial de ℝn .
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
• 𝑆𝑖𝛼⃗ = (𝛼1, … , 𝛼𝑛) es una solución de un sistema de ecuaciones lineales no
homogéneo, 𝐴 · 𝑋 = 𝐵, entonces cualquier vector de la forma ⃗𝛼 + 𝛽 ⃗, siendo 𝛽⃗
una solución del sistema homogéneo asociado: 𝐴 · 𝑋 = 0, también será solución
del sistema 𝐴 · 𝑋 = 𝐵. Además, es fácil demostrar que si ⃗𝛼 𝑦 ⃗𝛾 son dos
soluciones del sistema de ecuaciones lineales:
𝐴 · 𝑋 = 𝐵, entonces ⃗𝛼 − ⃗𝛾 es una solución del sistema homogéneo asociado.
Ambos resultados nos llevan a concluir que si se conoce una solución concreta de
un sistema 𝐴 · 𝑋 = 𝐵, entonces todas y cada una de las soluciones del sistema
se obtienen sumando a dicha solución las soluciones del sistema homogéneo
asociado. En particular, queda demostrado que un sistema compatible
indeterminado posee necesariamente infinitas soluciones.
𝑆𝑖 ⃗𝛼 = (𝛼1, . . . , 𝛼𝑛) es una solución de un sistema de ecuaciones lineales. no
homogéneo 𝐴 · 𝑋 = 𝐵1 𝑦 𝛽⃗ = (𝛽1, . . . , 𝛽𝑛) es una solución del sistema de
ecuaciones lineales no ℎ𝑜𝑚𝑜𝑔é𝑛𝑒𝑜 𝐴 · 𝑋 = 𝐵2 entonces ⃗𝛼 + 𝛽⃗ es solución del
S.E.L.: 𝐴 · 𝑋 = 𝐵1 + 𝐵2.
En resumen, podría decirse de la siguiente forma:
Un sistema de n ecuaciones lineales y n incógnitas generalmente se representa de la
siguiente forma:
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
O también en forma matricial:
Figura1: Interpretación de la forma matricial 𝐴𝑥 = 𝑏
Donde los 𝑥i son desconocidos y son lo que se desea encontrar. Para encontrar los 𝑥i
se pueden usar dos métodos principalmente: Métodos directos. Métodos iterativos.
Métodos Directos
Son aquellos en los cuales se obtendrá en un número finito de pasos la respuesta exacta
(suponiendo que se trabaja con infinitos dígitos). Estos métodos directos, estrictamente
hablando, pertenecen más al álgebra lineal que a los métodos numéricos; pero, debido a
su importancia y a la necesidad de hacer algunas observaciones importantes cuando se
desarrolla en los computadores, se verá solamente el de Gauss-Jordán, ya que otros
métodos, tal como el de factorización triangular, se pueden considerar una modificación
de aquel. En resumen, son métodos que, en un número finito de operaciones, obtienen
la solución exacta.
Métodos Gauss-Jordán
Es uno de los métodos más importantes para resolver sistemas de ecuaciones. Se tiene
el sistema matricial 𝐴𝑥 = 𝑏.
El método consiste en poner 1 en la diagonal de la matriz A y 0 en el resto, usando las
siguientes propiedades: Una fila se puede multiplicar por una constante. A una fila se le
puede sumar otra fila.
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
Para obtener la matriz final:
Donde:
𝑣i = 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒 𝑥i
𝐼 = 𝑚𝑎𝑡𝑟𝑖𝑧 𝑖𝑑𝑒𝑛𝑡𝑖𝑑𝑎𝑑
El elemento en el cual se pone 1 se llama elemento pivote app (𝑝 = 𝑝𝑖𝑣𝑜𝑡𝑒), de donde el
elemento pivote estará en la diagonal. La fila asociada al elemento pivote se llama fila
pivote (𝑓𝑖𝑙𝑎p ). La columna asociada al elemento pivote se llama columna pivote
(𝑐𝑜𝑙𝑢𝑚𝑛𝑎p ).
• Poner un 1 en el elemento pivote:
En el siguiente sistema de ecuaciones supongamos que el elemento pivote es el
a2,2, por lo que la fila pivote (𝑓𝑖𝑙𝑎p ) será la segunda, al igual que la columna pivote
(𝑐𝑜𝑙𝑢𝑚𝑛𝑎p ).
Ejemplo de aplicación:
Para poner un 1 en el elemento pivote = app = a2,2 basta con dividir la fila 2 por a2,2 = 4.
Quedando:
Generalizando: Para poner un 1 en la 𝑓𝑖𝑙𝑎p : 𝑓𝑖𝑙𝑎p ← 𝑓𝑖𝑙𝑎p /𝑎p,p.
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
Poner un 0 en la 𝒇𝒊𝒍𝒂k
En el mismo sistema de ecuaciones anterior se desea poner un 0 en la fila 4. Seguir
considerando que la fila pivote (𝑓𝑖𝑙𝑎p) es la segunda, al igual que la columna pivote
(𝑐𝑜𝑙𝑢𝑚𝑛𝑎p ).
Primero hay que hacer que aparezca un 7 en el elemento pivote 𝑎p,p. = 𝑎2,2. .
Para lo anterior, basta con dividir la fila2 (𝑓𝑖𝑙𝑎p) por cuatro (𝑎p,p. ) y multiplicarla por siete
(𝑎k,p. ), quedando el elemento 𝑎2,2. = 4/4 * 7 = 7. Es decir, la 𝑓𝑖𝑙𝑎k se multiplica por:
multiplicadores de la eliminación
Quedando la fila2:
Para poner un 0 en la 𝑓𝑖𝑙𝑎k = 𝑓𝑖𝑙𝑎4, simplemente restar la fila obtenida, de la 𝑓𝑖𝑙𝑎k . El
resultado es la nueva 𝑓𝑖𝑙𝑎k . Quedando:
De donde: Para poner un 0 en la 𝑓𝑖𝑙𝑎k .
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
Como ya se dijo:
multiplicadores de la eliminación
Recordemos que para poner un 1 en la 𝑓𝑖𝑙𝑎p esta fue dividida por ap,p . Por tal razón, si
usamos esta filap, en la cual ya fue puesto un 1, la ecuación anterior se transforma en:
Para poner un 0 en la 𝑓𝑖𝑙𝑎k : 𝑓𝑖𝑙𝑎k ← 𝑓𝑖𝑙𝑎k - ak,p. * 𝑓𝑖𝑙𝑎p. Resumiendo:
Poner un 1 en la 𝑓𝑖𝑙𝑎p.
Poner un 0 en la 𝑓𝑖𝑙𝑎𝑝
Procedimiento
Poner un 1 en el elemento pivote ap,p. Poner 0 en todos los elementos de la 𝑐𝑜𝑙𝑢𝑚𝑛𝑎p
(excepto en ap,p .). Hacerlo para 𝑝 = 1, después 𝑝 = 2, etc. Al programar este método,
tenga los siguientes cuidados:
• Asigne ap,p a otra variable, ej.: aux1.
• Asigne ak,p . a otra variable, ej.: aux2.
• Use estas variables en las fórmulas:filap ← ffilap /aux1, etc.
Determinación de la inversa y del determinante
Es conveniente recodar que usando el método de Gauss-Jordán se pueden obtener la
inversa de la matriz y su determinante.
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
Procedimiento para calcular la inversa de la matriz A
Se amplía la matriz 𝐴 en la identidad (𝐼).
[𝐴 | 𝐼 ] = [𝑏]
O también: [𝐴 | 𝐼 |𝑏]
Al finalizar el procedimiento, donde estaba 𝐴 queda la matriz identidad 𝐼, lo cual es
equivalente a haber multiplicado por A-1
Es decir;
A-1 ∗ 𝑨 = 𝑰
A-1 ∗ 𝑰 = A-1
A-1 ∗ 𝒃 = 𝑿
𝐸𝑠 𝑑𝑒𝑐𝑖𝑟, 𝑞𝑢𝑒𝑑𝑎: [𝐼 |𝐴– 1] = [𝑋].
𝐷𝑜𝑛𝑑𝑒 𝑒𝑠𝑡𝑎𝑏𝑎 𝐴 𝑞𝑢𝑒𝑑𝑎 𝐼.
𝐷𝑜𝑛𝑑𝑒 𝑒𝑠𝑡𝑎𝑏𝑎 𝐼 𝑞𝑢𝑒𝑑𝑎 𝐴– 1.
𝐷𝑜𝑛𝑑𝑒 𝑒𝑠𝑡𝑎𝑏𝑎 𝑏 𝑞𝑢𝑒𝑑𝑎 𝑋.
Procedimiento para calcular el determinante de la matriz A Hay que tener en cuenta que:
Cada vez que se divide una fila por una constante el determinante queda dividido por ese
número. Cada vez que se intercambian filas el signo del determinante cambia. El
determinante de la matriz I = 1. Por lo tanto, el valor del determinante es igual a la
multiplicación de todos los AP,P por lo que se dividió cada 𝑓𝑖𝑙𝑎P para poder obtener un 1
en el elemento pivote.
Signo del determinante:
• Negativo: si el número de intercambios de filas es
impar.
• Positivo: si el número de intercambios de filas es par.
Determinante R = # de intercambios de la fila que hubo.
Pivoteo parcial Consiste en escoger como fila pivote aquella en la que se encuentre el
elemento de mayor valor absoluto dentro de la columna pivote (desde la diagonal
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
principal hacia abajo). Se recomienda fuertemente su uso para reducir la propagación de
errores.
Método de Jacobi
De 𝐴𝑥 = 𝑏 se tiene que (𝐷 − 𝐿 − 𝑈)𝑥 = 𝑏 y, por tanto, 𝐷𝑥 = (𝐿 + 𝑈)𝑥 + 𝑏. Luego:
𝑥 (𝑚+1) = 𝐵j 𝑥 (𝑚) + 𝑐j .
Siendo
Por tanto, el valor de cada incógnita en cada paso del método 𝑚 es
Teniendo en cuenta la anterior igualdad se deduce el siguiente método iterativo.
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
Método de descomposición LU
Si 𝐴 = 𝐿𝑈, con 𝐿 una matriz triangular inferior (𝑙𝑜𝑤) y 𝑈 una matriz triangular superior
(𝑢𝑝𝑝𝑒𝑟), entonces el sistema de ecuaciones lineales 𝐴𝑥 = 𝑏 se puede resolver en dos
fases:
𝑎) 𝐿𝑦 = 𝑏,
𝑏) 𝑈𝑥 = 𝑦,
con un coste de 2n 2 operaciones.
Ejemplo:
Aplicando el algoritmo de resolución de un sistema factorizado en 𝐿𝑈, obtenemos
Figura2: Representación de la Factorización LU
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
4𝑛 3 + 2𝑛 2𝑛3
Coste general de la factorización: = 𝑂( ), (a falta de añadir 𝑛2 + 𝑛2
6 3
).
Algunas variantes:
• Variante de Doolittle: L tiene diagonal de 1.
• Variante de Crout: U tiene diagonal de 1.
Figura3: Representación de la Factorización LU para matrices tridiagonales.
Coste: 4𝑛 − 3 para factorizar y 2(3𝑛 − 2) para resolver.
Método de Gauss-Seidel
De 𝐴𝑥 = 𝑏 se tiene que (𝐷 − 𝐿 − 𝑈)𝑥 = 𝑏 y, por tanto, (𝐷 − 𝐿)𝑥 = 𝑈𝑥 + 𝑏. Luego
𝑥 = (𝐷 − 𝐿)-1 𝑈𝑥 + (𝐷 − 𝐿)-1 𝑏
Teniendo en cuenta la anterior igualdad se deduce el siguiente método iterativo
(𝐷 − 𝐿)𝑥 (𝑚+1) = 𝑈𝑥 (𝑚) + 𝑏
que se trata de un sistema triangular inferior en cada paso, que se resuelve por sustitución
progresiva. Por tanto, el valor de cada incógnita en cada paso del método 𝑚 es
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
Método de Cramer
Un sistema de ecuaciones lineales recibe el nombre de sistema de Cramer cuando se
cumplen las dos condiciones siguientes: El número de ecuaciones es igual al número
de incógnitas. El determinante de la matriz de los coeficientes (matriz del sistema) es
distinto de cero ( 𝑑𝑒𝑡 ( 𝐴 ) # 0 ) Un sistema de Cramer es, por definición, compatible
determinado, puesto que se cumple que rango (𝐴) = 𝑟𝑎𝑛𝑔𝑜 (𝐴 ∗) =
𝑛 (𝑛º 𝑑𝑒 𝑖𝑛𝑐ó𝑔𝑛𝑖𝑡𝑎𝑠). Consideremos un sistema de Cramer, es decir, un sistema de n
ecuaciones lineales con 𝑛 incógnitas, cuya expresión general es la siguiente:
Sean A la matriz del sistema (matriz de los coeficientes), entonces det (A) # 0.
Llamaremos matriz asociada a la incógnita xi y la designaremos por Ai a la matriz que se
obtiene al sustituir en la matriz del sistema la columna i por la matriz columna de los
términos independientes. Es decir:
Todos los sistemas de Cramer son compatibles determinados. El valor de cada
incógnita se obtiene dividiendo el determinante de la matriz asociada a dicha incógnita
por la matriz del sistema (matriz de los coeficientes de las incógnitas).
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
1
TECNOLOGICO DE ESTUDIAOS SUPERIORES DE HUIXQUILUCAN
Departamento de Ingeniería en Sistema computacionales
Métodos Numéricos
Bibliografía
Chapra, S. (2011). Métodos Numéricos para Ingenieros. México: Mc Graw Hill.
Luthe, O. (1996). Métodos Numéricos. México: Limusa.
Nieves Hurtado, A. (2015). Métodos numéricos: aplicados a la ingeniería. Grupo Editorial
Patria. [Link]
Saver, T. (2013). Análisis Numérico. México: Pearson.
Vázquez, L. (2009). Métodos numéricos para la física y la ingeniería. McGraw-Hill España.
[Link]
