2013 09 14 Sturm

2013
Sistema de polinomios de Sturm
RAICES DE ECUACIONES
Yamil Armando Cerquera Rojas

Universidad Surcolombiana
15/09/2013
RAICES DE ECUACIONES
Teorema De Sturm
Ing Yamil Armando Cerquera Rojas – yacerque@gmail.com
Especialista en Sistemas Universidad Nacional
Especialista en Administración de la Informática Educativa
Neiva - Huila
Contenido
Preámbulo ....................................................................................................................................... 4
Ecuaciones de resolución ......................................................................................................... 4
Definición de Ecuación ............................................................................................................. 5
Desarrollo del Concepto ....................................................................................................... 5
Clasificación de las ecuaciones con una incógnita: ........................................................ 6
Raíces de ecuaciones: introducción ........................................................................................... 6
Historia de las ecuaciones lineales. ....................................................................................... 7
Historia de los sistemas de ecuaciones lineales. ............................................................... 12
Teoría de ecuaciones .................................................................................................................. 13
Polinomios ................................................................................................................................. 13
Grado de un polinomio ........................................................................................................... 14
Ecuaciones con una sola incógnita ....................................................................................... 14
Ecuaciones con varias incógnitas .......................................................................................... 15
Raíces de un polinomio .......................................................................................................... 16
Factorización de un polinomio .............................................................................................. 17
Representación gráfica de las raíces de un polinomio ..................................................... 18
Raíces Únicas y Múltiples: ...................................................................................................... 20
Teorema fundamental del Álgebra....................................................................................... 23
Todo polinomio de grado n tiene n raíces. ......................................................................... 23
Regla de los signos de Descartes .......................................................................................... 24
Conjunto de posibles raíces................................................................................................... 25
¿Qué hacer cuando se tenga una raíz? ................................................................................. 26
Deflación polinomial............................................................................................................... 27
TEOREMA DE STURM .................................................................................................................... 29
Charles Sturm: Biografía ........................................................................................................ 29
Introducción ............................................................................................................................. 29
Sistema de polinomios de Sturm........................................................................................... 30
EJEMPLO 1: Calcule el sistema de polinomios de Sturm para
f ( x) = x 3 − 9 x 2 + 24 x − 36 . .................................................................................................. 31
EJEMPLO 2: Sistema de polinomios de Sturm ................................................................ 32
El Teorema de Sturm .............................................................................................................. 34
TEOREMA 1: Teorema de Sturm ....................................................................................... 34
EJEMPLO 3: Teorema de Sturm ........................................................................................ 34
Cota superior para los ceros de un polinomio .................................................................... 36
TEOREMA 2: Cota superior de los ceros de un polinomio. .......................................... 36
EJEMPLO 4: Búsqueda de cotas f ( x) = x 3 − 9 x 2 + 24 x − 36 = 0 ..................................... 36
TEOREMA 3: Búsqueda de cotas........................................................................................ 37
EJEMPLO 5: Búsqueda de cotas f ( x) = x 3 − 9 x 2 + 24 x − 36 = 0 ..................................... 37
EJEMPLO 6: Ejemplo completo para búsqueda de ceros .................................................. 37
Ejemplo 7: p(x) = x 4 − x 3 − x 2 + x − 1 .................................................................................... 40
Reflexiones para un programa computacional ................................................................... 42
Programa para ejecutar en MatLab: .................................................................................... 50
Apéndice ................................................................................................................................... 53
Métodos para aproximar los ceros de un polinomio ..................................................... 53
Método de bisección ........................................................................................................... 53
Método de bisección acelerada ........................................................................................ 54
Método de Newton .............................................................................................................. 54
Corolario del Teorema de Rolle y la Regla de Descartes ................................................. 54
TEOREMA 4 Corolario del Teorema de Rolle ............................................................... 54
TEOREMA 5 La Regla de Descartes................................................................................. 54
Recursos Bibliográficos ........................................................................................................... 56
Bibliografia OnLine: ................................................................................................................ 57
Preámbulo
Muchas veces al desarrollar un tema en clase, no se consideran ciertos aspectos que
son importantes para el estudiante, como la parte histórica detrás de un concepto o
teorema, o la explicación de cómo se aplican estos conceptos en nuestra realidad. La
Matemática es una ciencia que no está desligada del mundo, todo lo contrario, si se
observa atentamente el paisaje que le rodea, encontrará la aplicación de algún
concepto matemático.
Por esta razón, el papel del estudiante es ahondar en los temas desarrollados y
relacionarlos con su entorno inmediato.
Cuando se encuentra ante una situación problemática, a veces se plantean

ecuaciones, y es el Álgebra que en uno de sus aspectos más elementales, brinda el
conocimiento necesario para resolverla. No hay que olvidar, además, que los
algoritmos de resolución de ecuaciones y de sistemas de ecuaciones han ocupado a
muchos matemáticos a lo largo de la historia. Por ejemplo, se conoce la existencia de
problemas resueltos por procedimientos algebraicos que datan del año 3000 AC.
Los babilonios y egipcios fueron los más entusiastas cultivadores del Álgebra. A partir
de ellos, los griegos, y luego los árabes, elaboraron los cimientos de la Matemática
actual, así como el contexto de todos los axiomas y teoremas que se fueron
desarrollando.
El presente libro está centrado en mostrar cómo las culturas milenarias abordaron el
tema de resolución de ecuaciones y cómo en la actualidad, de algún modo muy sutil,
se están relacionando problemas ancestrales con situaciones reales y actuales.
Resolver una ecuación o un sistema de ecuaciones son problemas que se plantearon en

la antigüedad y que, salvo las diferencias en los actuales métodos más complicados de
resolución, lograban. Determinar la solución adecuada a cada situación.
Ecuaciones de resolución
En la antigüedad, algunas culturas representaron los números mediante letras. Los
griegos usaban su alfabeto para representar los números, al igual que la numeración
romana.
En realidad, el Álgebra comienza cuando los matemáticos se empiezan a interesar por

operaciones que se pueden realizar con cualquier número más que por los mismos
números, lo que los llevó a generalizar un número cualquiera a través de una letra
llamada variable.
Si bien en un principio no existían las variables, los problemas se plantearon mediante
palabras, por lo que se le llamó Álgebra Retórica1, y la variable era llamada cosa, de
ahí que el Álgebra fuera conocida como “la regla de la cosa”. A partir del siglo XII los
árabes introducen el Álgebra Simbólica, la cual asigna símbolos a la variable buscada
en el problema.
Tras muchos milenios, las ecuaciones matemáticas de la forma + = 0 se resuelven

actualmente por diversos métodos, incluso, mediante el uso de computadoras. Existe
software que resuelven diversos tipos de ecuaciones.
Definición de Ecuación
Una ecuación es una igualdad que contiene una o más incógnitas.
Desarrollo del Concepto

En una ecuación existen cantidades desconocidas (incógnitas), que en general se
designan por letras minúsculas de la parte final del alfabeto: x, y, z y cantidades
conocidas (coeficientes), que pueden designarse por letras minúsculas iniciales del
alfabeto: a, b, c. Lo anterior lo introdujo el matemático René Descartes en 1637.
En la ecuación: + = , a, b y c son coeficientes, es la incógnita.
En la ecuación 5 − 4 = 16, los coeficientes son los enteros 5, -4, y 16 y la incógnita es

.
Se llama raíces o soluciones de la ecuación a los valores de las incógnitas que cumplen
la igualdad.
Ejemplos:
1. Si voy a la panadería con $5000 y deseo comprar 3 panes de $1000 ¿qué vueltas
recibiré? Si v representa el valor de las vueltas, éste tiene que cumplir:
5000 = 3 ∗ 1000 +
En la ecuación anterior v es la incógnita y el valor = 2000 es la solución.
2. Un norteamericano llega a Neiva – Huila - Colombia y manifiesta que la

temperatura ambiental es de 85° Farenheit ¿cuál es la temperatura ambiental en
grados Celsius? Si C representa la temperatura en grados Celsius se cumple la
relación:
1
Se denomina así, en el campo de la educación matemática, a la primera fase en el desarrollo histórico
del álgebra, debido a que los problemas y sus soluciones se describían mediante lenguaje natural, sin
incluir ningún símbolo, ni siquiera de las operaciones. Es el álgebra de la edad clásica. Desde los
egipcios y babilonios hasta la obra de Diofanto (Siglo III)
85 = 9 / 5o C + 32 , C es la incógnita y el valor = 265/9 es la solución.
En: ax + 8 = 3 , el valor = −5/ es la solución.
En 4 – 3 + 2 = 2 los valores =0y = 3/4 son soluciones.
Clasificación de las ecuaciones con una incógnita:

Las ecuaciones se catalogan según el exponente o potencia más alto que tenga la
incógnita.
Así,
6 + 34 = 5, Es una ecuación de primer grado.

8 + 7 + 45 = 3, Es una ecuación de segundo grado.
4 + 35 – 3 + 2 = 7, Es una ecuación de tercer grado.
Raíces de ecuaciones: introducción2

Constantemente, al plantear en términos matemáticos problemas de distintas áreas
(economía, física, ingeniería, biología, etc.), aparece la siguiente cuestión: tratar de
determinar los ceros de ciertas funciones, es decir valores para los cuales la función se
anula.
Resolver una ecuación es encontrar el ó los valores de la incógnita que satisface la

igualdad.
Por ejemplo la ecuación: 500 = 450 + , se satisface para = 50
Después de las funciones lineales, las funciones polinomiales (en 1 variable) son las
más simples.
Estudiar los ceros (raíces) de funciones polinomiales tiene un gran interés por lo
menos por las dos razones siguientes:
– No se puede pretender poder resolver el problema para funciones más

generales si no se logra resolverlo en este caso más sencillo.
– La mayoría de veces es posible traducir de alguna manera el problema original

de hallar ceros de una función cualquiera, al de calcular las raíces de ciertos
polinomios (que “aproximan” a la función original).
2
Polinomios y Raíces, Teresa Krick., Departamento de Matemática. Facultad de Ciencias Exactas y
Naturales. Universidad de Buenos Aires. -1428- Buenos Aires. ARGENTINA.
Generalmente, para las aplicaciones, se trabaja con funciones reales, y se trata de
encontrar ceros reales. Más aún, debido a la estructura de los números con los cuales
trabajan las computadoras, las funciones suelen tener coeficientes racionales y los
ceros que se calculan serán números racionales que aproximan suficientemente una
verdadera solución al problema.
En situaciones reales la solución de la ecuación debe tener sentido en el contexto en

que se trabaja.
Nótese los siguientes casos:
a) Pertinencia de la solución: Se quiere repartir equitativamente 24 dulces a 5

niños. Sea la cantidad de dulces que corresponde a cada niño, x debe ser un
número natural que satisfaga la ecuación: 5 = 24
La ecuación anterior no tiene solución en los naturales (N).
b) Existencia de la solución: La ecuación 4 x + 5 = 2, no tiene solución en los

naturales (N) ni en los enteros (Z) sino que en los racionales y en los reales.
La ecuación 4 * x * x = −7 , no tiene solución en los reales (R) ya que no existe ningún

número real que la satisfaga.
c) Infinitas soluciones: La ecuación 2 + + = 2( + 1), es una ecuación que es

satisfecha por cualquier valor que tome x, luego tiene infinitas soluciones.
En este libro, se trata de profundizar sobre raíces de polinomios con coeficientes en Q

(el cuerpo de los números racionales), ℜ (el cuerpo de los números reales) y C (el
cuerpo de los números complejos). Se verá en qué medida la teoría sobre Q está más
resuelta, y la de C, más teórica, permite aclarar la estructura de los polinomios en
R[X] .
Además presupone cierta familiaridad con la teoría general básica de polinomios en

una variable, aunque todos los resultados que se usan están. También se requieren
ciertas nociones de Análisis elemental, para el caso de funciones reales continuas y
derivables (como por ejemplo los Teoremas de Bolzano y de Rolle).
Historia de las ecuaciones lineales.

La primera fase, que comprende el periodo de 1700 AC a 1700 DC., se
caracterizó por la invención gradual de símbolos y la resolución de
ecuaciones. Dentro de esta fase se encuentra un álgebra desarrollada por
los griegos (300 AC), 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ète3 (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.
Los egipcios nos dejaron en sus papiros4 (sobre todo en el de Rhid-1.650 AC y el de

Moscú-1.850 AC) multitud de problemas matemáticos resueltos. La mayoría de ellos
son de tipo aritmético y respondían a situaciones concretas de la vida diaria; sin
embargo, encontramos algunos que podemos clasificar como algebraicos, pues no se
refiere a ningún objeto concreto. En éstos, de una forma retórica, obtenían una
solución realizando operaciones con los datos de forma análoga a como hoy se
resuelven dichas ecuaciones.
Las ecuaciones más utilizadas por los egipcios eran de la forma:
x + ax = b y x + ax + bx = 0
Donde a, b y c eran números conocidos y x la incógnita que ellos denominaban aha o

montón.
Una ecuación lineal que aparece en el papiro de Rhid responde al problema siguiente:
"Un montón y un séptimo del mismo es igual a 24".
En notación moderna, la ecuación sería: x + 1 x = 24

7
La solución la obtenía por un método que hoy se conoce con el nombre de "método de
la falsa posición" o "regula falsi". Consiste en tomar un valor concreto para la
3
Fue un matemático francés (Fontenay-le-Comte, 1540 - París, 1603). Se le considera uno de los
principales precursores del álgebra. Fue el primero en representar los parámetros de una ecuación con
letras. François Viète también fue conocido en su época como súbdito del rey fiel y competente. Fue
consejero privado de los reyes de Francia Enrique III y de Enrique IV.
4
Papiro es el nombre que recibe el soporte de escritura elaborado a partir de una planta acuática,
también denominada papiro, muy común en el río Nilo, en Egipto, y en algunos lugares de la cuenca
mediterránea, una hierba palustre de la familia de las ciperáceas, el Cyperus papyrus.
Etimológicamente, la palabra papiro proviene del término griego πάπυρος papiros que en latín es
papyrus, cuyo plural es papyri, tomada del egipcio antiguo per-peraâ, que significa “flor del rey”.
También es el origen de la palabra papel.
incógnita, se prueba con él y si se verifica la igualdad ya se tiene la solución, si no,
mediante cálculos se obtiene la solución exacta.
Suponga que fuera 7 la solución, al sustituir en la x daría: 7 + 1 7 = 8 , y como nuestra

7
solución es 24, es decir, 8*3, la solución es 21 = 3 ∗ 7, ya que 3 * (7 + 1 / 7 * 7) = 24 .
Generalmente, el cálculo de la solución correcta no era tan fácil como en este caso e
implicaba numerosas operaciones con fracciones unitarias (fracciones con numerador
la unidad), cuyo uso dominaban los egipcios. En cuanto el simbolismo, solamente en
algunas ocasiones utilizaban el dibujo de un par de piernas andando en dirección de la
escritura o invertidas, para representar la suma y resta, respectivamente.
Los babilonios (el mayor número de documentos corresponde al periodo 600 AC. a 300
DC.) casi no le prestaron atención a las ecuaciones lineales, quizás por considerarlas
demasiado elementales, y trabajaron más los sistemas de ecuaciones lineales y las
ecuaciones de segundo grado.
Entre las pocas que aparecen, Se tiene la ecuación 5x = 8 . En las tablas en base
sexagesimal hallaban el recíproco de cinco que era 12/60 y en la tabla de multiplicar
por 8, se encuentra que 8 ∗ 12/60 = 1 ∗ 36/60.
Los matemáticos griegos no tuvieron problemas con las ecuaciones lineales y,
exceptuando a Diophante5 (250 DC), no se dedicaron mucho al álgebra, pues su
preocupación era como se ha visto, mayor por la geometría. Sobre la vida de
Diophante aparece en los siglos V o VI un epigrama algebraico que constituye una
ecuación lineal y dice:
"Transeúnte, ésta es la tumba de Diophante: es él quien con esta sorprendente

distribución te dice el número de años que vivió. Su juventud ocupó su sexta
parte, después durante la doceava parte su mejilla se cubrió con el primer
vello. Pasó aún una séptima parte de su vida antes de tomar esposa y, cinco
años después, tuvo un precioso niño que, una vez alcanzada la mitad de la
edad de su padre, pereció de una muerte desgraciada. Su padre tuvo que
sobrevivirle, llorándole durante cuatro años. De todo esto, deduce su edad. "
Los primeros documentos matemáticos que existen (datan del siglo III d. de C.) son los
Sulvasütras6, donde se recogen todos los conocimientos necesarios para construir los
templos. En éstos aparece el siguiente problema:
5
Nació alrededor del 200 d.C. y Murió: alrededor del 284 DC, Diofanto, a menudo conocido como el
'padre del álgebra', es mejor conocido por su Aritmética, un trabajo sobre la solución de ecuaciones
algebraicas y sobre la teoría de los números. Sin embargo, esencialmente no se sabe nada de su vida y
ha existido mucho debate respecto de la fecha en que vivió.
http://ciencia.astroseti.org/matematicas/articulo.php?num=3629
6
Los Sulvasutras son apéndices de los Vedas en los que se dan normas para la construcción de altares.
Para que el sacrificio ritual tuviera éxito, el altar tenía que ajustarse a medidas muy precisas. La gente
hacía sacrificios a sus dioses para que estuvieran contentos y le diesen abundante comida, buena
suerte, salud, una vida larga y muchos otros beneficios materiales. Para tener contentos a los dioses
“Hallar el lado de un rectángulo, conociendo el otro lado y sabiendo que su área
es igual al área de un cuadrado dado "
Esto es:
Es decir, ∗ = .
Lo resolvían utilizando el método de la falsa posición, como los egipcios.
Posteriormente, Brahmagupta7 (siglo VII) expresa, ya de forma sincopada, cómo

resolver ecuaciones lineales. La incógnita la representaba por la abreviatura ya, y las
operaciones con la primera sílaba de las palabras.
Dada la ecuación ax + b = cx + d , la solución vendrá dada dividiendo la diferencia de

los términos conocidos entre la diferencia de los coeficientes de los desconocidos,
esto es,
d −b
x=
a−c
Estos métodos pasaron a los árabes que los extendieron por Europa. Al algebrista
“Abu-Kamil”8 (siglo IX y X) se le atribuye una obra donde trata la solución de
ecuaciones lineales por simple y doble falsa posición.
todo se tenía que realizar siguiendo una fórmula muy precisa, por lo que la exactitud matemática se
veía como algo de importancia capital. Todo lo que se sabe de las matemáticas védicas se halla en los
Sulbasutras.
7
Fue un matemático y astrónomo indio. Su padre fue Jisnugupta. Nació en el año 598, posiblemente en
Ujjain, donde vivió. En esta ciudad de la zona central de la India se encontraba el más famoso y antiguo
observatorio de astronomía del que Brahmagupta era el director. Está considerado el más grande de los
matemáticos de esta época. Murió en el año 670. Es posible que Brahmagupta haya sido el idealizador
del concepto del "cero" ya que en su obra Brahmasphutasiddhanta del año 628 aparece por primera vez
esta idea. La obra trataba también sobre aritmética y números negativos en términos muy parecidos a
los de la moderna matemática
8
La muerte de al-Jwarizmi coincide aproximadamente con el nacimiento en Egipto de Abu Kamil ibn
Aslam ibn Mohammed, llamado el calculista egipcio. Vivió ochenta años y nos dejó numerosas obras
matemáticas. Entre ellas un tratado de álgebra, cuyo original árabe se ha perdido, pero del que nos han
llegado dos traducciones, una latina y otra hebrea. Las ecuaciones de segundo grado las resuelve
geométricamente, como su predecesor de Bagdad, pero se apoya más directamente en los Elementos.
Demuestra una proposición cuyo equivalente algebraico es la célebre fórmula de suma por diferencia
El método de la doble falsa posición es el siguiente:
Sea la ecuación ax + b = 0 y suponga dos valores para la x:

Ecuación 1
x=m am + b = p
Ecuación 2
x=n an + b = q Ec. 2
Restando, a ( m − n ) = p − q
Por otra parte, eliminando a del sistema de ecuaciones anteriores se tiene.
Ecuación 1*n → amn + bn = pn

Ecuación 2*m → amn + bm = qm
Que restando, b( n − m) = pn − qm
Y dividiendo ambos resultados, − a/b = (p − q)/(pn − qm)
O también − b /a = (pn − qm)/(p − q)
Siendo esto último el valor de x.
Vea un ejemplo. Sea la ecuación 5 x − 10 = 0 , si se toma como valor de x: x = 3 y x = 4

, y sustituyendo,
5 * 4 − 10 = p
5 * 3 − 10 = q
Se tiene que: x = (10 * 3 − 5 * 4) /(10 − 5) = (30 − 20) / 5 = 10 / 5 = 2
Este principio fue posteriormente presentado en una forma ligeramente modificada

por el método de las escalas. El nombre proviene de un diagrama que permitía escribir
la solución rápidamente:
igual diferencia de sus cuadrados. Con este resultado, la ecuación x + bx = c puede se escrita de este
2
modo.
Las dos líneas de la izquierda representan p y q y las de la derecha m y n y la cruz del
centro indica que hay que multiplicar.
El método puede ser sintetizado como sigue:
1. Consideran dos valores cualesquiera de la incógnita m, n.

2. Calculan los errores correspondientes a ellos p, q.
3. Hallan el valor de la incógnita en función de los valores dados y sus errores.
En el ejemplo,
A partir de aquí se dedican al estudio de ecuaciones de grado superior.
Historia de los sistemas de ecuaciones lineales.

Los sistemas de ecuaciones lineales fueron ya resueltos por los babilonios, los cuales
llamaban a las incógnitas con palabras tales como longitud, anchura, área, o volumen,
sin que tuvieran relación con problemas de medida.
Un ejemplo tomado de una tablilla babilónica9 plantea la resolución de un sistema de

ecuaciones en los siguientes términos:
1/4 ! ℎ#$ + %&!'()#* = 7 + !&,

%&!'()#* + ! ℎ#$ = 10 + !&,
Para resolverlo comienzan asignando el valor 5 a una mano y observaban que la

solución podía ser: anchura=20, longitud=30. Para comprobarlo utilizaban un método
parecido al de eliminación. La notación, sería:
- + 4 = 28 y - + = 10
Restando la segunda de la primera, se obtiene 3 = 18, es decir = 6 e - = 4.
También resolvían sistemas de ecuaciones, donde alguna de ellas era cuadrática.
9
Los Babilonios vivieron en Mesopotamia, en unos claros de tierras fértiles entre los ríos Tigris y
Éufrates, hacia finales del milenio IV antes de Cristo. Desarrollaron una forma abstracta de escritura
basada en símbolos cuneiformes. Sus símbolos fueron escritos en tablas de arcilla mojadas cocidas al sol
Miles de estas tablillas han sobrevivido hasta nuestros días. Gracias a ello, se ha podido conocer, entre
otras cosas, gran parte de las matemáticas babilónicas. El uso de una arcilla blanda condujo a la
utilización de símbolos cuneiformes sin líneas curvas porque no podían ser dibujadas.
Los griegos también resolvían algunos sistemas de ecuaciones, pero uti1izando
métodos geométricos. Thymaridas10 (400 AC) había encontrado una fórmula para
resolver un determinado sistema de n ecuaciones con n incógnitas.
Diophante resuelve también problemas en los que aparecían sistemas de ecuaciones,

pero transformándolos en una ecuación lineal.
Diophante sólo aceptaba las soluciones positivas, pues lo que buscaba era resolver
problemas y no ecuaciones. Utilizó ya un álgebra sincopada como se ha señalado
anteriormente. Sin embargo, unas de las dificultades que se encuentra en la
resolución de ecuaciones por Diophante es que carece de un método general y utiliza
en cada problema métodos a veces excesivamente ingeniosos.
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.
El libro El arte matemático, de autor chino desconocido (siglo III a. de C.), contiene
algunos problemas donde se resuelven ecuaciones. En ellos se encuentra 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.
Teoría de ecuaciones
Rama de las matemáticas que estudia la naturaleza de las raíces de ecuaciones
polinómicas y los métodos de búsqueda de dichas raíces. La teoría de las ecuaciones
tiene aplicaciones en todas las ramas de las matemáticas y de las ciencias.
Polinomios
Un polinomio es una expresión algebraica que se obtiene al expresar cualquier suma
de monomios no semejantes. Una expresión de variable real de la forma
n −1
a n x + a n −1 x + ... + a 2 x + a1 x + a 0 , siendo n un entero no negativo a n , a n −1 ,..., a 2 , a1 , a 0
n 2
dados como números reales o complejos se llama: Polinomio de coeficientes en ℜ (o

en C) y con indeterminada x. Si se recuerda la suma de monomios, cuando estos no
son semejantes no se pueden sumar. En este caso lo que se obtiene es por tanto un
polinomio. A los polinomios suele denotárseles por: .( ), /( ), 0( ), Así, se hablará del
polinomio: P( x) = a n x n + a n−1 x n −1 + ... + a 2 x 2 + a1 x + a0 .
a 0 : Se llama término independiente del polinomio /( ).
10
Thymaridas was a follower of Pythagoras who wrote on prime numbers and gave methods for solving
simultaneous linear equations
a n : Es el coeficiente principal. En particular, cuando 1 = 1, el polinomio se llama
mónico.
En resumen, un polinomio es una suma de términos llamados monomios. Un monomio

es el producto de un coeficiente (un número real), una variable (casi siempre x o y)
elevada a un exponente (entero positivo).
Existen polinomios con uno, dos o más términos, por ejemplo:
Monomio (un término): F ( x ) = 5 x 2 Es un monomio pues tiene un solo

coeficiente no nulo. En este caso el coeficiente es 5, la variable es x el
exponente 2. También se puede decir que 2( ) es un polinomio de un solo
término.
Binomio (dos términos): 6 x − 2
7
Trinomio (tres términos): 3 x 5 + 4 x 3 − x 2
En el caso de los siguientes ejemplos:
1. 4a 2 x 4 y 3 + x 3 by 2 + 3ab 2 xy 3
2. 4 x 4 − 2 x 3 − 3x 2 + 5 x − 6
Se puede decir que en el caso de la primera expresión se trata de un polinomio de 3

términos y en el caso de la segunda expresión se trata de un polinomio de 5 términos.
Se diferencia el segundo del primero en que este último tiene un término
independiente y en este caso es -6.
Grado de un polinomio
Ecuaciones con una sola incógnita

Dado un polinomio / en una cierta variable , su grado es el máximo de los
exponentes de x en los distintos monomios del polinomio. Se suele denotar como
'$(/( )), y se puede omitir la variable si no hay posibilidad de confusión.
Dado el polinomio P( x) = a n x n + a n−1 x n −1 + ... + a 2 x 2 + a1 x + a0 y siendo a n ≠ 0 se dice
que el grado de P(x) es n, siendo n un entero no negativo.
gr (P( x) ) = n
Se debe tener en cuenta que:
Si P ( x) = a 0 , su gr (P( x)) = 0
Si P ( x) = 0 , el polinomio no tiene grado
Ejemplo:
El polinomio F ( x ) = 5 x 2 − 6 x + 3 tiene como gr (F ( x) ) = 2
En resumen se puede decir que el grado de un polinomio es igual al exponente mayor

de la variable. Por ejemplo:
5x 2 Es un monomio de grado 2
6x7 − 2 Es un binomio de grado 7

3x 5 + 4 x 3 − x 2 Es un polinomio de grado 5
3x15 − 2 x 4 + x ¿De qué grado es?
6 x 4 − 3x 2 + x ¿De qué grado es?
3x12 + x10 + x 5 ¿De qué grado es?
13 ¿De qué grado es?
Al polinomio: g ( x ) = 0 + 0 x + 0 x 2 + .. + 0 x n se le llama polinomio nulo
Si a 0 ≠ 0 , al polinomio g ( x) = a0 + 0 x + 0 x 2 + .. + 0 x n −1 + 0 x n se le denomina polinomio

constante.
A un polinomio, cuyo grado es un entero positivo, se le llama polinomio no

constante.
Ecuaciones con varias incógnitas

Cuando se tiene una ecuación algebraica con varias incógnitas, se estudia el grado de
distinta manera. Un monomio es un producto de incógnitas, multiplicadas a su vez por
números. Por ejemplo, - es un monomio, porque sería la multiplicación de las
incógnitas x e y, y a su vez está multiplicado todo por 1 (que no se coloca porque
multiplicar por 1 es como no hacer nada). Otro ejemplo de monomio sería − 4 x 3 y 2 z 4 .
Aquí las incógnitas son x, y, z, se multiplican así: la x se multiplica tres veces a sí
misma (porque x 3 = x * x * x ), la y se multiplica dos veces a sí misma, la z se multiplica
cuatro veces a sí misma, y los tres resultados se multiplican entre sí. Finalmente se
multiplica todo por el valor numérico -4.
Para calcular el grado de una ecuación con varias incógnitas, antes hemos de calcular
los grados de cada uno de los monomios que aparecen en la ecuación. El grado de un
monomio se calcula sumando los exponentes de las incógnitas que aparecen en el
monomio. Por ejemplo, el grado del monomio - es 2, porque es la suma del
exponente de x (que es 1, porque x = x 1 ) y del exponente de y (que también es 1). El
grado del monomio − 4 x 3 y 2 z 4 es 9, que es la suma de 3(exponente de x), 2(exponente
de y) y 4(exponente de z). Nótese que el grado del monomio 6x 2 sería 2, o sea, sería
el exponente de la incógnita, y que siempre se puede considerar que en un monomio
aparecen todas las incógnitas que hay en la ecuación, con sólo considerar que están
elevadas al exponente 0. Por ejemplo, en la ecuación - − 13- = 4 los monomios son
- (aparecen las dos incógnitas de la ecuación, y su grado es 2), −13- (aparece sólo
la incógnita y, pero se puede considerar que aparece también x con exponente 0,
puesto que 3 = 1) y 4 (no aparecen ni ni -, pero se puede considerar que aparecen
como 3 - 3 ). Así, se puede ver la ecuación como - − 13 3 - = 4 3 - 3 . Esto no cambia
el grado de ninguno de los monomios. El monomio 4 tiene entonces grado 0.
Ahora se esta en condiciones de calcular el grado de una ecuación de varias

incógnitas.
Este es el mayor de los grados de todos los monomios que aparecen en la ecuación.
Por ejemplo, en la ecuación - − 13- = 4 el grado es 3, que el grado más grande
entre los grados de todos los monomios de la ecuación (que son 2, 3 y 0).
Es fácil ver que el grado de una ecuación con una incógnita no es otra cosa que un
caso particular del grado de una ecuación con varias incógnitas.
Nota: ¿Cómo se deben escribir los polinomios?. Se deben escribir en orden decreciente
con respecto al grado de cada término.
En este libro se utilizarán polinomios con coeficientes enteros y potencias enteras

positivas.
Raíces de un polinomio
En muchos campos de las matemáticas es necesario hallar las raíces de un polinomio,
por ejemplo, para calcular la integral de una función racional, en la transformada de
Laplace, etc. Solamente existen fórmulas si el polinomio tiene un grado igual o
inferior a cuatro. Excepto para los polinomios de primer y segundo grado, las fórmulas
son complicadas, por lo que se emplean procesos de aproximación numérica. Entre los
numerosos métodos numéricos que existen, el más conocido es quizá el método de
Newton. Sin embargo, se describirán otros métodos, que proporcionan exactitud en las
raíces de un polinomio.
La raíz de un polinomio es un valor numérico tal que hace que el polinomio al ser
evaluado en dicho valor, valga cero. Es decir que, cuando se iguala un polinomio a
cero, las soluciones son las raíces del polinomio.
Se puede decir entonces, que un número a es una raíz del polinomio P (x ) si el valor
numérico de /( ), para x = a es cero. Es decir, P ( a ) = 0 . Las raíces enteras de un
polinomio con coeficientes enteros son divisores del término independiente.
Por ejemplo el polinomio: f ( x) = x 2 + x − 12 .
Cuando se iguala a cero y se resuelve se tiene:
x 2 + x − 12 = 0 Igualando a cero.
( x + 4)( x − 3) = 0 Factorizando.
x = −4 Raíz 1
x=3 Raíz 2
Puesto que x1 = −4 y x2 = 3 , son soluciones de f (x) entonces f (−4) = 0 y f (3) = 0 . Se

dice entonces que x1 = −4 y x2 = 3 , son raíces del polinomio f ( x) = x 2 + x − 12
Las raíces de f ( x) = x 3 − 4 x 2 + x + 6 son = −4, =5y = 6 ¿Por qué?
Factorización de un polinomio
Teorema del factor (factorización):
Si 7 es una raíz del polinomio P (x ) , entonces − 7 es un factor de P (x ) .
Por el contrario, si − 7 es un factor de P (x ) , entonces 7 es una raíz de P (x ) .
El número de factores en que se puede descomponer un polinomio es igual al grado

del polinomio.
Para poder factorizar un polinomio es necesario encontrar sus raíces. Cuando se

tengan estas, los factores correspondientes a cada raíz son de la forma ( − 7) donde r
es una de las raíces.
Esto es, si $8 , $ , . . . , $1 son raíces del polinomio f (x) entonces la factorización de f (x)
es:
f ( x) = ( x − r1 )( x − r2 )...( x − rn )
Por ejemplo, si
1. f ( x) = x 3 − 4 x 2 + x + 6 ,
Como sus raíces son x = −1 , x = 2 y x = 3 entonces f (x) se ha factorizado como:

f ( x) = ( x − ( −1))( x − 2)( x − 3) = ( x + 1)( x − 2)( x − 3)
2. :( ) = + − 12,
Como sus raíces son x = −4 y x = 3 entonces f (x) se ha factorizado como
f ( x) = ( x − ( −4))( x − 3) = ( x + 4)( x − 3)
Si un factor − 7 ocurre k veces en la factorización completa de un polinomio P (x ) , entonces r

es una raíz de /( ) = 0 con multiplicidad k.
Ejemplos:
En el polinomio P ( x ) = x 2 − 10 x + 25 es un polinomio con raíz 5 de multiplicidad 2.

Observa que, x 2 − 10 x + 25 = ( x − 5)( x − 5) = ( x − 5) 2 .
Un polinomio /( ) de menor grado, con coeficiente principal 1 que tiene las siguientes
raíces: -7 de multiplicidad 3 y 5 de multiplicidad 4 queda expresado de la forma
factorizada como: P( x) = ( x + 1) 3 ( x − 5) 4 .
El polinomio P ( x) = 6( x − 3) 3 ( x + 1) 2 ( x − i )( x + j ) es de grado siete y por lo tanto

tiene siete raíces, no todas diferentes. Observa que 3 es una raíz de multiplicidad 3; -
1 es una raíz de multiplicidad 2; i y -j es de multiplicidad 1.
Así que este polinomio de grado siete tiene exactamente siete raíces tomando en
cuenta al 3 y al -1 con sus respectivas multiplicidades.
Representación gráfica de las raíces de un polinomio

Como las raíces de un polinomio, son valores que al evaluar el polinomio en dichos
valores, hacen que éste valga cero, en un plano cartesiano esto se identifica como las
intersecciones de la gráfica del polinomio con el eje de las X (abscisas), si y solo si
son raíces únicas o repetidas impar veces.
Esto es, los puntos en donde cruza la gráfica al eje horizontal tienen como abscisa la
raíz del polinomio graficado.
A continuación se presentan algunas funciones con sus raíces, factores y gráficas:
Descripción Gráfica
Función f ( x) = x 2 + x − 12
Raíces -4y3
Factorización f ( x) = ( x + 4)( x − 3)
Figura 1
Función f ( x) = x 3 − 4 x 2 + x + 6
Raíces - 1, 2 y 3
Factorización f ( x) = ( x + 1)( x − 2)( x − 3)
Figura 2
Función f ( x) = x 4 − 5 x 2 + 4
Raíces - 2, - 1, 1 y 2
Factorización f ( x) = ( x + 1)( x + 2)( x − 1)( x − 2)
Figura 3
Función f ( x) = x 3 + 4 x 2 + 3 x
Raíces ¿Cuáles son?
Factorización f(x) =
Figura 4
Función f ( x) = x 3 − 2 x 2 − 5 x + 6
Raíces 1, - 2 y 3
Factorización f ( x) = ( x − 1)( x + 2)( x − 3)
Figura 5
Raíces Únicas y Múltiples:
Los polinomios pueden tener raíces únicas en un punto determinado del eje x o raíces
que se repiten en un número par o impar veces, es decir una raíz sobre el eje x puede
ser la misma par o impar veces. Dicho de otra manera una raíz por ejemplo dos (2)
puede ser la misma raíz pero repetida tres (3) veces (impar) o repetida 2 veces (par).
La tabla siguiente muestra la función f ( x) = x 2 − 4 , y en la gráfica se observa como

esta corta el eje x tanto en el punto -2 y 2 de un lado a otro siguiendo la forma de la
figura.
Se puede decir que la figura pasa de un lado al otro el eje x en el punto de corte o
raíz. Por lo tanto la raíz es única en el valor de menos dos (-2) y única en el valor de
dos (2).
Es un polinomio de orden dos por lo tanto solo tendrá dos (2) raíces.
En la parte donde se muestra el polinomio factorizado f ( x) = ( x − 2)( x + 2) se puede

observar que si la variable x toma el valor de 2 o toma el valor de -2 la función tomará
el valor de cero.
Función f ( x) = x 2 − 4
Raíces - 2, 2
Factorización f ( x) = ( x − 2)( x + 2)
Figura 6
Para el caso del siguiente ejemplo, la función f ( x) = x 2 − 2 x + 1 tiene dos (2) raíces, y
si observa la gráfica, esta no corta el eje x en ningún sector. Ahora si observa el valor
de uno (1) en el eje x, es un punto donde la función se vuelve cero (0).
Se debe considerar al valor de uno (1) como raíz de la función.
Lo que pasa es que dicha raíz se repite par veces (para el caso del ejemplo 2 veces), por
esta razón la gráfica no corta el eje x, sino que lo toca tangencialmente en el punto
raíz y cambia su pendiente. Se puede decir matemáticamente que en el punto raíz, la
derivada de la función es igual a cero (0) ó dicho de otra manera en este punto la
tangente es igual a cero (0).
En la factorización que se hace de la función, se puede observar que el único valor

que hace que la función tome el valor de cero (0) es el punto uno (1) sobre el eje x, y
puede ser en el término de la izquierda o en el término de la derecha, es decir dos
veces.
Función f ( x) = x 2 − 2 x + 1
Raíces 1, 1
Factorización f ( x ) = ( x − 1)( x − 1)
Figura 7
Para el caso del siguiente ejemplo la función f ( x) = x 3 − 6 x 2 + 12 x + 8 tiene 3 raíces, y

si observa la gráfica, esta corta el eje x aparentemente en varios puntos cercanos a 2.
Ahora si observa el valor de 2 en el eje x, es un punto donde la función se vuelve cero.

Se debe considerar al valor de 2 como raíz de la función.
Lo que pasa es que dicha raíz se repite impar veces (para el caso del ejemplo 3 veces), por
esta razón la gráfica corta el eje x de la forma como se observa en la figura.
En el punto de corte sobre el eje x, este y la gráfica son paralelos superpuestos.
Se puede decir matemáticamente que en el punto raíz la derivada de la función es

igual a cero.

que hace que la función tome el valor de cero es 2, y puede ser en el término de la
izquierda, en el término del centro o en el término de la derecha, es decir tres veces.
Función f ( x) = x 3 − 6 x 2 + 12 x + 8
Raíces 2, 2, 2 Figura 8
Factorización f ( x) = ( x − 2)( x − 2)( x − 2)
Para el caso del siguiente ejemplo la función f ( x) = x 3 − 3x 2 + 4 tiene 3 raíces, y si

observa la gráfica, esta corta el eje x aparentemente en un punto igual a menos uno (-
1) y toca tangencialmente dicho eje en un valor igual a uno (1).
Ahora si observa el valor de -1 en el eje x, es un punto donde la función cruza el eje
de las x con cierta pendiente, esto indica que ese punto de corte es una raíz única.
Ahora en le punto 1 sobre el eje de las x la curva o grafica de la función toca
tangencialmente el eje y cambia de pendiente. Esto debe asumirse como una raíz que
se repite par veces. Como el polinomio de es orden 3 y ya se sabe de una raíz única se
puede decir que dicha raíz es par veces repetida.
Se puede decir matemáticamente que en el punto 1 considerado como raíz repetida
par veces, la derivada de la función es igual a cero.
que hace que la función tome el valor de cero es -1, en el término de la izquierda o 1
en los dos términos de la derecha.
En el caso de que la raíz 1 se repitiera 4 veces diferenciaría la forma de la gráfica en
que la pendiente de esta es mayor o menor al acercarse al eje.
Función f ( x) = x 3 − 3 x 2 + 4
Raíces - 1, 2, 2
Factorización f ( x) = ( x + 1)( x − 2)( x − 2)
En el siguiente ejemplo muestra la combinación de las tres formas que toma la gráfica
dependiendo si sus raíces se repiten par o impar veces o son raíces únicas.
f ( x) = x 6 − 17 x 5 + 102 x 4 − ....
Función
− 248 x 3 + 160 x 2 + 240 x − 288
Raíces - 1, 2, 2, 2, 6, 6
Factorización f ( x) = ( x + 1)( x − 2) 3 ( x − 6) 2
Que análisis de acuerdo con lo mostrado anteriormente le puede realizar a las

siguientes funciones.
f ( x) = ( x + 1)( x − 1) 2 ( x − 1) f ( x) = ( x + 1)( x − 2) 5 ( x − 6) 4
f ( x) = ( x + 2)( x − 2) 3 ( x) f ( x) = ( x − 1) 2 ( x + 1) + 1
f ( x) = ( x + 1)( x − 2) 5 ( x − 6) 4 + 1 f ( x) = ( x − 1) 2 − 1
f ( x) = ( x − 1) 2 + ( x − 1) f ( x) = (( x − 1) 2 + 1) 2 ( x − 1)
Teorema fundamental del Álgebra

Carl Friedrch Gauss ha sido uno de los matemáticos más grandes de todos los tiempos.
Contribuyó a muchas ramas de las matemáticas. En 1798, ¡a los 20 años de edad!,
Gauss demostró el teorema fundamental del Álgebra que dice lo siguiente:
Todo polinomio de grado n tiene n raíces.

Si se toma una ecuación en términos generales, tal como la ecuación siguiente:
a n x n + a n −1 x n −1 + a n − 2 x n − 2 + a n −3 x n −3 + ... + a3 x 3 + a 2 x 2 + a1 x 1 + a 0 = 0 .
Se puede decir entornes que es una ecuación de orden n y por tanto tiene n
soluciones.
Recuerde que en es este apartado sólo se tiene polinomios con coeficientes enteros.
Observa la tabla anterior, donde se da la función, las raíces y la gráfica y verifica que
efectivamente para cada polinomio de grado n hay n raíces.
Una forma en la que se puede interpretar este teorema es como sigue, ya que se
puede factorizar un polinomio, dadas las raíces y hay n raíces para todo polinomio de
este grado, entonces si:
f ( x) = a n x n + a n−1 x n −1 + a n− 2 x n −2 + a n−3 x n−3 + ... + a3 x 3 + a 2 x 2 + a1 x1 + a0 ,
Se puede decir que:
f ( x) = ( x − r1 )( x − r2 )...( x − rn )
Donde r1, r2, ... , rn son las raíces de f (x) .

La demostración de este teorema queda lejos del objetivo de esta página sin embargo
daremos algunas herramientas para encontrar las n raíces.
Regla de los signos de Descartes

Rene Descartes encontró un método para indicar el número de raíces positivas en un
polinomio.
Esta regla dice lo siguiente: “El número de raíces reales positivas (+) de un polinomio
f (x) es igual al número de cambios de signo de término a término de f (x) "
Hay que recordar que los polinomios se deben escribir en orden decreciente conforme
al grado de cada término.
Por ejemplo los polinomios:
f ( x ) = x 2 + x − 12 Tiene un cambio de signo, del segundo al tercer término, por

lo tanto tiene una raíz positiva.
g ( x ) = + x 3 − 4 x 2 + x + 6 Tiene dos cambios de signo, tiene dos raíces positivas
h( x ) = + x 4 − 5 x 2 + 4 Tiene dos raíces positivas
i ( x) = x 3 + 4 x 2 + 3 x No tiene cambios de signo, por lo tanto no tiene raíces
reales positivas.
j ( x ) = x 3 − 2 x 2 − 5 x + 6 ¿Cuántas raíces positivas tiene?
También puede evaluar la expresión en los valores 1 y -1 teniendo en cuenta solo el

signo de cada término de la expresión. En caso de que un término no exista se toma
como 0 (positivo).
En el ejemplo de la siguiente tabla para el caso de la primera función al evaluar la

función en 1, el primer término de la función toma un valor positivo y el segundo toma
un valor negativo, por tanto se dice que tiene una raíz positiva en razón a un solo
cambio de signo (de positivo a negativo). Como se trata de una ecuación de una recta
pues tan solo tiene una raíz.
Pero por probar se ha evaluado la función en -1, resultando en ambos términos un

signo -, o sea que no hay cambio de signo, indicando con esto que no hay raíces
negativas, en la función f ( x ) = x − 1 .
Nro Ecuación Signo Rai_Pos Rai_Neg

f (1) = + −
1 f ( x) = x − 1 1 0
f (−1) = − −
f (1) = + + −
2 f ( x) = x 2 + x − 12 1 1
f (−1) = + − −
f (1) = + − − +
3 f ( x) = x 3 − 2 x 2 − 5 x + 6 2 1
f ( −1) = − − + +
f (1) = + + + +
4 f ( x) = x 3 + 4 x 2 + 3x 0 3
f (−1) = − + − +
f (1) = + + − + +
5 f ( x) = x 4 − 5 x 2 + 4 2 2
f (−1) = + + − + +
f (1) = + − + +
6 f ( x) = x 3 − 4 x 2 + x + 6 2 1
f (−1) = − − − +
Sea f (x) un Polinomio que tiene Coeficientes Reales y un término Independiente

diferente de Cero,
• El número de Raíces Positivas de f ( x ) = 0 es igual al Número de Variaciones de

signo de f (x) o es Menor que este número en una Cantidad Par.
• El número de Raíces Negativas de f ( x) = 0 es igual al Número de Variaciones de
signo de f (− x) o es Menor que este número en una Cantidad Par.
Veremos un ejemplo para aclarar dudas: f ( x) = 2 x 5 + 3 x 4 − 2 x 3 + 4 x 2 − 2 x − 5
Como se puede observar entre el Segundo y el Tercero, entre el Tercero y el Cuarto y

entre el Cuarto y el Quinto Término hay Variaciones de Signo, por lo tanto tiene 3
Variaciones,
Ahora si se calcula f ( − x) se tiene: f (− x) = −2 x 5 + 3 x 4 + 2 x 3 + 4 x 2 − 2 x − 5
Como se puede observar entre el Primero y el Segundo y entre el Cuarto y el Quinto

Término hay Variaciones de Signo, por lo tanto tiene 2 Variaciones, por lo tanto la
Regla de los Signos de Descartes nos quedaría:
Las Raíces pueden ser:
Positivas Negativas Imaginarias

3 2 0
1 2 2
3 0 2
1 0 4
Conjunto de posibles raíces

Existe un método para encontrar un conjunto de números, los cuales pueden ser raíces
de un polinomio. La regla que mencionaremos aquí es aplicable sólo para polinomios
con el coeficiente de la potencia mayor de x igual a 1. Es decir, si
f ( x) = a n x n + a n −1 x n −1 + a n − 2 x n − 2 + a n −3 x n −3 + ... + a3 x 3 + a 2 x 2 + a1 x1 + a 0 , se toma a a0 = 1 .
Esto es que sólo se trabaja con polinomios de la siguiente forma:
f ( x ) = x n + a n −1 x n −1 + a n − 2 x n − 2 + a n −3 x n −3 + ... + a3 x 3 + a 2 x 2 + a1 x1 + a 0
El conjunto de posibles raíces de f (x) se forma con los divisores de a0 (del término
independiente), hay que considerar estos divisores tanto con signo positivo como con
negativo.
La forma en que se puede usar esta información del término independiente es la

siguiente, puesto que cualquier elemento de este conjunto puede ser raíz de f (x) hay
que evaluar a f (x) en algún valor de este conjunto y si el resultado de la evaluación
es cero, entonces ese valor escogido es raíz de f (x) .
En la siguiente tabla se muestran varios polinomios, los divisores del término

independiente y las raíces de los polinomios:
Función Divisores del término independiente Raíces

1, 2, 3, 4, 6, 12,
f ( x) = x 2 + x − 12 -1, -2, -3, -4, -6, -12
-4y3
1, 2, 3, 6,
f(x)= x3 - 4 x2 + x + 6 - 1, 2 y 3
-1, -2, -3, -6
1, 2, 4,
f(x)= x4 - 5 x2 + 4 - 2, - 1, 1 y 2
-1, -2, -4
1, 2, 3, 6,
f(x)= x3 - 2 x2 - 5 x + 6 1, - 2 y 3
-1, -2, -3, -6
¿Qué hacer cuando se tenga una raíz?

Con lo visto en los apartados anteriores se tiene las herramientas necesarias para
encontrar n raíces de un polinomio. Recuerde que para encontrar una raíz es
necesario saber los divisores del término independiente y evaluar el polinomio en con
el valor escogido.
Además de haber encontrado una raíz usando el método anterior se ha hallado un

factor del polinomio. Se puede estar seguro de que si r es una raíz de f (x) entonces
al dividir f ( x) /( x − r ) tendrá como resultado un polinomio de un grado menor a f (x) y
como residuo cero.
Así se ha reducido el problema de encontrar n raíces en otro problema, el encontrar

sólo n-1 raíces.
Deflación polinomial.
Suponga que conoce, una de las raíces de un polinomio, se puede realizar la división
de este polinomio para obtener un polinomio de grado menor. Así por ejemplo si se
tiene:
f n ( x) = a n x n + a n −1 x n −1 + ... + a 2 x 2 + a1 x + a 0
Si una de las raíces es s, entonces se puede escribir
f n ( x ) = ( x − s ) * ( a ' n −1 x n −1 + ... + a ' 2 x 2 + a '1 x + a ' 0 )
En este caso el residuo de la división es cero y se puede calcular un polinomio de

grado n-1. Para un polinomio de orden 3 se tiene que
a3 x 2 (a 2 − a3 s ) x a1 − (a 2 − a3 s ) s
x+s a3x3 +a2x2 + a1x +a0
-a3x2 - a3 s x2
0 (a2-a3s) x2 + a1x +a0
-(a2-a3s) x2 -(a2-a3s)sx
0 (a1-(a2-a3s)s)x +a0
-(a1-(a2-a3s)s)x -(a1-(a2-a3s)s)s
a0-(a1-(a2-a3s)s)s
La división sintética, es una manera de hacer lo mismo, pero de forma compacta. Así
se tiene:
x+s a3 a2 a1 a0
- a3 s -(a2- a3 s)s -(a1-(a2-s a3)s)s
a3 a2- a3 s a1-(a2- a3 a0-(a1-(a2- a3
s)s s)s)s
Ejemplo.
Dado el polinomio f 5 ( x) = x 5 − 7 x 4 − 3x 3 + 79 x 2 − 46 x − 120 encontrar la división

sintética con el monomio (x-4).
x-4 1 -7 -3 79 -46 -120
4 -12 -60 76 120
1 -3 -15 19 30 0
El polinomio resultante, en este caso, es f 4 ( x) = x 4 − 3 x 3 − 15 x 2 + 19 x + 30 . Note que el

residuo es cero.
TEOREMA DE STURM
Charles Sturm: Biografía
(Nació en Ginebra-Suiza en 1803 y murió en París – Francia hacia
1855) Matemático suizo, nacionalizado francés. En colaboración con
J.D. Colladon, midió la velocidad del sonido en el agua del lago
Leman (1827). Enunció (1829) un famoso teorema sobre las raíces de
las ecuaciones algebraicas. Escribió diversas obras de geometría,
análisis y álgebra.
Introducción
Uno de los objetivos principales de la teoría de ecuaciones es la búsqueda de ceros de

un polinomio. Existen varios métodos relacionados con este tema, algunos tratan sobre
la búsqueda de las cotas de los ceros (los valores entre los cuales se encuentran todos
los ceros), otros sobre la aproximación de los ceros y hay otros, como el teorema de
Sturm, que sirven para saber cuántos ceros tiene la función entre dos valores de " x'' y
así poder aislarlos.
El conocimiento de estos teoremas vuelve más sencilla la búsqueda de los ceros del
polinomio y permiten realizar, por ejemplo, algoritmos para hacer programas de
computación que realicen todo el trabajo.
Para la aproximación de los ceros existen varios métodos, algunos de ellos son la
bisección, la bisección acelerada, el método de Newton o el de punto fijo, el
problema que se presenta es que se deben tener valores iniciales entre los cuales se
esté seguro que se encuentra un cero, estos valores iniciales son esenciales para tener
éxito en el proceso y encontrarlos no es nada sencillo, para lograrlo es que se utiliza
el teorema de Sturm.
Por ejemplo, si se tiene la seguridad que un polinomio tiene un cero entre 0 y 1, se

puede aplicar alguno de los métodos con esos valores y así encontrar ese cero; pero
cómo saber que el polinomio tiene un cero entre estos valores? y cómo se puede estar
seguro que es el único cero en ese intervalo?, ya que si hubieran más ceros podrían
haber problemas y se tendrían que utilizar límites más cercanos, como 0 y 0.5, hasta
tener un solo cero aislado.
En realidad hay varios métodos para lograr aislar los ceros, una primer forma sería
usar la gráfica del polinomio, pero para esto se debe tener una buena gráfica que
permita observar todos los ceros y ésta puede cambiar para cada polinomio, además
que no serviría para realizar un programa computacional, la computadora no puede
ver donde están los ceros en una gráfica.
El método basado en un corolario del Teorema de Rolle se vuelve poco práctico

cuando el grado del polinomio es mayor que 3 y aunque la regla de Descartes es muy
sencilla, es muy común que falle al calcular el número total de ceros reales.
Teniendo esto en cuenta, el método más efectivo para aislar los ceros es el Teorema
de Sturm, lo primero que se debe ver es el sistema de polinomios de Sturm, éste
sistema lo necesita para hacer evaluaciones en él, ya que el número de cambios de
signo entre las distintas evaluaciones dará la cantidad de ceros entre los valores.
Para hacer más sencilla la búsqueda de los ceros, también hay una sección que trata
sobre el cálculo de las cotas, estos valores permiten saber entre que números se
encuentran todas las raíces, esto es muy importante porque se sabría exactamente
donde buscar.
Al final se muestran algunas ideas para poder hacer un programa computacional, no se

muestra el código paso a paso sino que se explican las funciones necesarias para
hacerlo, lo que recibe cada función, lo que devuelve y una explicación general de
cómo se hacen.
Sistema de polinomios de Sturm
Suponga que se van a buscar los ceros de un polinomio f (x) , lo primero que se tiene
que buscar al utilizar el Teorema de Sturm es un sistema de polinomios especial para
poder evaluar en él, vea cómo calcular este sistema.
Se dirá que una sucesión de funciones f 0(x), f 1(x), . . . , f m (x ) definidas y continuas en [ a, b ] es

una sucesión de Sturm en [ a, b ] si:
• f 0 es diferenciable,
• f m (x) no se anula en (a,b),
• Si f k (α ) = 0 con k = 1,..., m − 1 entonces f k-1(α)f k +1(α) < 0
• Si f 0 (α ) = 0 entonces f1 (α ) f ' (α ) > 0 .
Se denomina método del Sturm, al procedimiento para calcular todas las raíces reales de un
polinomio p n (x ) y que utiliza una sucesión de Sturm para separar las raíces reales y una
combinación del método de bisección y Newton para obtener con precisión su valor.
El primer polinomio es el propio f (x) ; el segundo polinomio es f ' ( x ) , es decir, la

derivada del polinomio f (x) , a este polinomio se lo denota como f1(x); para el tercer
polinomio se debe realizar la división de f (x) entre f1(x), el residuo que se obtiene
con signo contrario será f2(x); este proceso se continúa, es decir, f3(x) será el residuo
con signo contrario de la división de f1(x) entre f2(x); el proceso acaba cuando se
obtiene una constante.
Si los cálculos se van hacer a mano, existe una forma de hacer que no aparezcan
fracciones en la división, haciéndola más sencilla, si el dividendo es
P ( x ) = a n x n + ... + a 0 y el divisor es Q ( x ) = bm x m + ... + b0 , entonces multiplique P (x ) por
bmn −m +1 . Si la división la hará la computadora, no hace falta multiplicar el dividendo ya
que esto no afecta el resultado del Teorema de Sturm y a la computadora le da lo
mismo hacer los cálculo con fracciones o sin ellas.
Vea un ejemplo,
EJEMPLO 1:
Calcule el sistema de polinomios de Sturm para
f ( x) = x − 9 x + 24 x − 36 .
3 2
Solución
f ( x) = x 3 − 9 x 2 + 24 x − 36
f 1 ( x) = f ' ( x) = 3 x 2 − 18 x + 24 ≡ ax 2 + bx + c, con a = 3, b = −18, c = 24
Al realizar la división de a 2 f ( x) ≡ 3 2 f ( x) entre f 1 ( x) y esquemáticamente

2
representado por a f ( x) f 1 ( x ) se obtiene:
Primer paso de la división Segundo paso de la división

9 x − 81x 2 + 216 x − 324 3 x 2 − 18 x + 24
3
9 x 3 − 81x 2 + 216 x − 324 3 x 2 − 18 x + 24 − 9 x 3 + 54 x 2 − 72 x 3x − 9

− 9 x + 54 x − 72 x
3 2
3x → −27 x 2 + 144 x − 324
−27 x 2 + 144 x − 324 27 x 2 − 162 x + 216
−18 x − 108
Tomando el residuo con signo contrario se tiene que f 2 ( x) = 18x + 108 ≡ a' x + b' ,
con a ' = 18, b ' = 108 .
Resumiendo se tendría hasta el momento el siguiente sistema de ecuaciones de Sturm:
f ( x) = x 3 − 9 x 2 + 24 x − 36
f 1 ( x) = 3x 2 − 18 x + 24
f 2 ( x) = 18x + 108
Ahora se realiza la división de a' 2 f 1 ( x) ≡ 18 2 f 1 ( x) entre f 2 ( x) , realizando

2
a' f 1 ( x) f 2 ( x) , al igual como se realizó en el paso anterior y obteniendo,
972 x 2 − 5832 x + 7776 18 x + 108

− 972 x 2 − 5832 x 54 x − 648
−11664 x + 7776
11664 x + 69984
77760
Tomando el residuo con signo contrario se tiene a f 3 ( x) = −77760
Así, el sistema de polinomios de Sturm que se obtiene para la función

f ( x) = x 3 − 9 x 2 + 24 x − 36 es:
f ( x) = x 3 − 9 x 2 + 24 x − 36
f 1 ( x) = f ' ( x) = 3 x 2 − 18 x + 24
f 2 ( x) = 18x + 108
f 3 ( x) = −77760
Vea otro ejemplo para dejar más claro el proceso,
EJEMPLO 2: Sistema de polinomios de Sturm
Busque el sistema de Sturm para f ( x) = x 3 + 4 x 2 − 7

Solución
f ( x) = x 3 + 4 x 2 − 7
f 1 ( x ) = f ' ( x) = 3 x 2 + 8 x
Al realizar la división de 3 2 f ( x) entre f1 ( x) se obtiene
9 x 3 + 36 x 2 + 0 x − 63 3 x 2 + 8 x
− 9 x 3 − 24 x 2 3x + 4
12 x + 0 x − 63
2
−12 x 2 − 32 x
−32 x − 63
Así f 2 ( x) = 32 x + 63
Ahora se realiza la división de 32 2 f 1 ( x) entre f 2 ( x)
3072 x 2 + 8192 x + 0 32 x + 63
− 3072 x 2 − 6048 x 96 x + 67
+2144 x + 0
−2144 x − 4221
−4221
Por lo tanto f 3 ( x) = 4221
El sistema que se obtiene para f ( x) = x 3 + 4 x 2 − 7 es:
f ( x) = x 3 + 4 x 2 − 7
f1 ( x) = f ' ( x) = 3x 2 + 8 x
f 2 ( x) = 32 x + 63
f 3 ( x) = 4221
El Teorema de Sturm
TEOREMA 1: Teorema de Sturm
Sea f (x) un polinomio de coeficientes reales tal que f ( x) = 0 no tiene raíces

múltiples. Construya el sistema de Sturm para f (x) . Sean a y b números reales, a < b a
y ninguno de los dos es raíz de f ( x) = 0 , entonces la cantidad de raíces reales de
f ( x) = 0 entre a y b es la diferencia entre el número de variaciones de signo del
sistema de Sturm.
f (x), f 1(x), f 2(x),..., f k - 1(x), f k (x)
Para x = b y el número de variaciones del sistema para x = a , los términos que den
cero deben ser descartados antes de contar los cambios de signo.
Observe que el teorema pide que la ecuación no tenga raíces múltiples, en este
escrito no se tomará en cuenta esta limitación, primero porque la probabilidad que
una persona evalúe un polinomio así es bajo, además que el teorema funciona bien
para muchos polinomios con raíces múltiples, es decir, son muy pocos los polinomios
de raíces múltiples que fallan.
Ahora se ilustra el teorema con un ejemplo.
EJEMPLO 3: Teorema de Sturm
Aísle los ceros de f (x) = x 3 − 9 x 2 + 24 x − 36.
Solución: El sistema de Sturm para este polinomio se calculó anteriormente.
f (x) = x 3 − 9 x 2 + 24 x − 36
f1(x) = 3x2 - 18x + 24
f2(x) = 18x + 108
f3(x) = - 77760
Lo primero para aislar los ceros es revisar cuántos hay en total, así que se evalúa el
sistema en − ∞ y + ∞ , también se evalúa en x = 0 para verificar cuántos ceros
positivos y cuántos negativos hay, vea
Signos al evaluar − ∞ 0 + ∞
f (x) - - +
f1(x) + + +
f2(x) - + +
f3(x) - - -
Cambios de signo 2 2 1
Como en − ∞ hay 2 cambios de signo y en + ∞ solo hay uno entonces de − ∞ a + ∞ hay

2-1=1 ceros, es decir, el polinomio solo tiene un cero en ℜ .
Ahora, de − ∞ a 0 no hay ceros pues ambos tienen igual número de cambios de signo,
así que el cero se encuentra de 0 a + ∞ pues al realizar la resta de cambios de signo
se obtiene 2 − 1 = 1 . Ahora se puede ir evaluando de uno en uno empezando desde 0
hasta encontrar en qué intervalo se encuentra.
Signos 1 2 3 4 5 6
f (x) - - - - - 0
f1(x) + 0 - 0 + +
f2(x) + + + + + +
f3(x) - - - - - -
Cambios 2 2 2 2 2 1
De aquí se obtiene que el cero esta entre 5 y 6. De hecho, en este caso se encontró
directamente que f (6) = 0 , es decir, el único cero que tiene el polinomio es x = 6 . Si
el cero no se hubiera encontrado directamente, se hubiera tenido que utilizar algún
método de aproximación en el intervalo ]5,6 [ .
Ahora, si el cero se encuentra en x = 100 , se tiene que evaluar 100 veces para
encontrarlo?, no hay alguna forma para tener un límite superior, algún valor del que
se pueda estar seguro que el cero se encuentra antes?
Cota superior para los ceros de un polinomio
Existen algunas formas para encontrar la cota superior para hallar la raíz de un
polinomio, aquí se presentan dos, no se puede decir cual de las dos sea mejor,
simplemente un método da la mejor cota para algunas funciones y el otro método es
el mejor para otras.
Se debe notar que estos teoremas son para encontrar la cota superior, sin embargo,
para encontrar la cota inferior simplemente se debe aplicar estos teoremas al
polinomio f (− x) (esta evaluación invierte el polinomio, convirtiendo los ceros
positivos en negativos y viceversa), vea los dos teoremas
TEOREMA 2: Cota superior de los ceros de un polinomio.
Si en una ecuación polinomial de coeficientes reales
f ( x ) = a n x n + a n −1 x n −1 + ... + a 0 = 0 (a n > 0)
El primer coeficiente negativo es precedido por k coeficientes positivos o cero y si G

denota el valor más grande de los valores absolutos de los coeficientes negativos,
G
entonces cada cero real es menor de 1 + k
an
Una condición del teorema es que a n > 0 , es decir, que el coeficiente del primer
término debe ser positivo, pero si esto no sucede, tan solo se debe multiplicar el
polinomio por -1 ya que − f ( x) = 0 tiene los mismos ceros que f ( x) = 0 .
EJEMPLO 4: Búsqueda de cotas f ( x) = x 3 − 9 x 2 + 24 x − 36 = 0
Encontrar las cotas de los ceros de la ecuación f ( x) = x 3 − 9 x 2 + 24 x − 36 = 0
Solución
Este es el ejemplo que se utilizó para aislar los ceros con el teorema de Sturm. El
primer coeficiente negativo es -9, este está precedido por un solo término, es decir
que k = 1 , además G es el mayor de los coeficientes de los valores negativos, en este
36
caso G = 36 , así que la cota es 1 + 1 = 1 + 36 = 37.
1
Para buscar la cota de los ceros negativos, se calcula f (− x) = − x 3 − 9 x 2 − 24 x − 36 ,

como el primer coeficiente es negativo, entonces − f (− x) = x 3 + 9 x 2 + 24 x + 36 , este
polinomio tiene todos sus coeficientes positivos, cuando esto sucede, el polinomio no
tiene ceros negativos, es decir que la cota es 0.
Por los dos hechos anteriores se puede asegurar que los ceros del polinomio
f ( x) = x 3 − 9 x 2 + 24 x − 36 = 0 están en el intervalo ]0,37 [ .
Para este ejemplo es mejor usar el segundo teorema ya que da una cota menor, vea:
TEOREMA 3: Búsqueda de cotas
Si en una ecuación polinomial de orden n, con sus coeficientes reales

f ( x) = a n x n + a n −1 x n −1 + ... + a 0 = 0 ( a n > 0) , el valor numérico de cada coeficiente
negativo se divide por la suma de todos los coeficientes positivos que lo preceden, la
cota superior será el mayor de estos valores aumentado en 1, observe mejor en el
siguiente ejemplo:
EJEMPLO 5: Búsqueda de cotas f ( x) = x 3 − 9 x 2 + 24 x − 36 = 0
Busque las cotas de los ceros de la ecuación f ( x) = x 3 − 9 x 2 + 24 x − 36 = 0
Solución
El mismo ejemplo que se utilizó con el teorema anterior, calcule las posibles valores a
9 36
tomar, c1 = = 9, c 2 = = 1.44 . Como se debe tomar el mayor aumentado en uno,
1 1 + 24
la cota es 9 + 1 = 10.
Para la negativa se toma − f (− x) = x 3 + 9 x 2 + 24 x + 36 , que como ya observó no tiene

soluciones negativas. Así, los ceros están en el intervalo ]0,10 [ , vea que esta cota es
mucho mejor que con el teorema anterior.
Observe que una vez que se tiene el intervalo en donde se encuentra el cero, se puede
ir partiendo ese intervalo a la mitad y así ir aislando cada solución. Por ejemplo, se
puede verificar con el teorema de Sturm en x = 0, 5, 10 , así se eliminan 5 valores de
una vez y quedaría el intervalo ]5,10 [ , luego se repite el proceso varias veces hasta
tener bien aislado el cero (se puede ver que este proceso se asemeja mucho al
proceso de bisección).
Ahora un ejemplo completo que haga uso de todos los teoremas.
EJEMPLO 6: Ejemplo completo para búsqueda de ceros
Aísle todos los ceros de la ecuación f ( x) = x 3 + 4 x 2 − 7 = 0

Solución
Se busca primero las cotas de la función, inicie con el primer teorema.

7
La cota superior es 1 + 3 ≈ 3
1
4
Para la cota inferior se usa el polinomio − f (− x) = x 3 − 4 x 2 + 7 . La cota es 1 + 1 =5
1
Por lo tanto los ceros se encuentran en el intervalo ]−5,3 [ .
Con el segundo teorema se obtiene
7
Para la cota superior, como solo hay un coeficiente negativo, la cota es 1 + = 2.4 .
1+ 4
Se toma 3 como cota superior.
Para la inferior, se utiliza el polinomio − f (− x) = x 3 − 4 x 2 + 7 y se obtiene como cota

4
1 + = 5 . Por lo que la cota inferior es -5.
1
Así, los ceros están en el intervalo ] − 5,3[ , este intervalo es igual al que se obtuvo
anteriormente, por lo que es el que se utilizará al aplicar el teorema de Sturm.
El sistema de Sturm ya se había calculado anteriormente, obteniendo que
f ( x) = x 3 + 4 x 2 − 7
f1 ( x ) = 3 x 2 + 8 x
f 2 ( x) = 32 x + 63
f 3 ( x) = 4221
Primero se evalúa la función en los valores − ∞ , 0 y + ∞ para ver cuántas raíces o

ceros hay en total y cuántos son positivos y cuántos negativos.
Signos −∞ 0 +∞
f (x) - - +
f1(x) + 0 +
f2(x) - + +
f3(x) + + +
Cambios 3 1 0
Por lo tanto, la ecuación al ser de orden 3, tiene un total de 3 ceros o raíces, hay 2
raíces negativas (en el intervalo ]−5,0 [ pues -5 es la cota inferior) y 1 raíz positiva (en
]0,3[ pues 3 es la cota superior), inicie aislando cada uno de ellos, evalúe en otros
valores centrales para ir reduciendo el intervalo, inicie con la positiva
Signos f (x) f1(x) f2(x) f3(x) Cambios

x=1 - + + + 1
Por lo que el primer cero está entre 1 y 3, taL como se puede observar en la grafica de
la función mostrada anteriormente.
Ahora aísle las negativas

x=-2 + - - + 2
Así que hay un cero en ]−5,−2 [ y otro en ]−2,0 [ , ahora

x=-1 - - + + 1
Por lo que una solución está en ]−2,−1 [ .
Para la última solución se tiene

x=-3 + + - + 2
Así que el último cero se encuentra en ]−5,−3[ .
Por último, si se quiere verificar estos resultados (tal vez los resultados no están muy
claros), vea una tabla completa
Signos − ∞ -4 -3 -2 -1 0 1 2 + ∞
f (x) - - + + - - - + +
f1(x) + + + - - 0 + + +
f2(x) - - - - + + + + +
f3(x) + + + + + + + + +
Cambios 3 3 2 2 1 1 1 0 0
Ahora ya se ve claramente que una solución está en ]−4,−3 [ , otra en ]−2,−1[ y la

última en ]1,2 [ . Para terminar de encontrar los ceros sólo falta aproximarlos aplicando
algún método para ello, en el apéndice se pueden ver algunos métodos para hacerlo.
Ejemplo 7: p(x) = x 4 − x 3 − x 2 + x − 1
Dado el polinomio p(x) = x 4 − x 3 − x 2 + x − 1
Determinar:
a) Determine una sucesión de Sturm asociada a este polinomio.

b) Obtenga la tabla de signos asociada a dicha sucesión de polinomios.
c) ¿Cuantas raíces tiene dicho polinomio?
d) ¿Cuántas raíces reales tiene dicho polinomio?
e) ¿Cuántas raíces positivas? Escriba un intervalo de longitud unidad en el que se
encuentren cada una de ellas.
f) ¿Cuántas raíces negativas? Escriba un intervalo de longitud unidad en el que se
encuentren cada una de ellas.
g) Según el criterio de los signos de Descartes, ¿cuántas raíces positivas hay?
h) Según el criterio de los signos de Descartes, ¿cuántas raíces negativas hay?
i) Aplique el criterio de Cauchy y determine una región en el plano complejo en la
que se encuentren todas las raíces del polinomio.
j) ¿Cuántas raíces complejas hay?
k) Determine sendos intervalos para la parte imaginaria de cada una de las raíces
complejas.
l) El algoritmo de Bairstow consiste en dividir el polinomio
p(z) = a n z n + a n −1 z n −1 + · · · + a0 , entre el factor cuadrático c(u , v) = z 2 − uz − v .
Determine el cociente y el resto de dicha división, especificando sus
coeficientes b j como una relación de recurrencia en función de u y v.
m) Que ecuaciones hay que resolver mediante el método de Newton como parte
del método de Bairstow.
n) Aplique el método de Newton y obtenga las relaciones de recurrencia para las
derivadas parciales de los coeficientes b j entre u y v.
o) Escriba el método de Newton en formulación delta a partir de las expresiones
del apartado anterior.
Solución. Para resolver este problema se utiliza el método de las sucesiones de Sturm
(que tiene como único defecto que sólo determina el número de raíces distintas). La
sucesión de Sturm más simple es. Tenga en cuenta que para el presente ejercicio se
toma como p1 ( x) el valor de la función original, es decir p ( x) = x 4 − x 3 − x 2 + x − 1 .
p1(x) = p(x) = x 4 − x 3 − x 2 + x − 1 , p 2(x) = − p' (x) = −4 x 3 + 3x 2 + 2 x − 1

Generalizando pi = mod pi −2(x), pi −1(x),
Que se obtiene fácilmente utilizando el algoritmo de Euclides de división de

polinomios, dando
16 p1(x) = ( 1 − 4 x)p 2(x) − 11x 2 + 10 x − 15 Donde p 3 ( x) = −11x 2 + 10 x − 15

121 p 2 ( x) = (44 x + 7) p 3 ( x) + 832 x − 16 Donde p 4 ( x) = 832x − 16
43264 p3 ( x) = (509 − 572 x) p 2 ( x) − 640816 Donde p 5 ( x) = −640816
x p1 p2 p3 p4 p5 total
¡Error! No se pueden crear objetos modificando códigos de + + - - - 1
campo.
+∞ + - - + - 3
0 - - - - - 0
-1 - + - - - 2
-2 + + - - - 1
1 - - - + - 2
2 + - - + - 3
Que indica que existen 2 raíces reales distintas, una de ellas (negativa) en [-1,-2] y la
otra (positiva) en [1, 2].
Aplicando el criterio de los signos de Descartes, se puede observar que hay tres
cambios de signo v p y por tanto el número de raíces positivas n p cumple que
v p − n p ∈ {0,2,4} , luego o hay 1 raíz positiva o hay tres. En cuanto al numero de raíces
negativas, aplicando Descartes a q ( p ) = p ( − x ) , se obtiene v q = 1 y por tanto
v p − n q ∈ {0,2,4} hay exactamente una raíz negativa.
Aplicando el criterio de Cauchy que dice que las raíces (compleja o real) están
incluidas en el disco de radio
−1
ρ = 1 + an max a k
0≤ k ≤ n
Para p(x) obtenemos ρ = 2 , luego las raíces tienen la cota superior xi < 2 . Para
p (1 / x ) se obtiene ρ = 2 , luego las raíces tienen la cota inferior xi > 1 / 2 . Como se
sabe por el criterio de Descartes que hay una raíz negativa, esta estaría en el intervalo
[− 2,−1/ 2] .
La raíz, o las tres, reales positivas del polinomio estarán en el intervalo [1 / 2,2] .
a) Hay una raíz (real) positiva (Sturm).

b) Hay una raíz (real) negativa (Descartes o Sturm).
c) La raíz positiva está en [1, 2] (Sturm) y la negativa en [-2,-1](Sturm).
d) Hay un par de raíces complejas conjugadas de parte real en [1/2, 2] y de parte
imaginaria en [0, 2] y [-2, 0], respectivamente (Cauchy).
Calcule el resto de preguntas que se realizaron.
Por último, se podría hacer algunas reflexiones sobre cómo realizar un programa de
computación con estas ideas.
Reflexiones para un programa computacional
Si lo que se quiere es realizar un programa computacional que utilice el Teorema de

Sturm para aproximar los ceros de un polinomio se deben seguir una serie de pasos
que trataremos de enunciar a continuación, nosotros realizamos el programa en Java,
sin embargo se puede programar en cualquier otro lenguaje, aquí la idea no es realizar
un programa completo donde se muestre paso a paso como realizarlo, sino mostrar las
ideas principales sobre cómo se debe realizar, es decir, que no se va a mostrar el
código del programa completo, sino solo algunas de las funciones que son necesarias,
el programa lo puede bajar al final de este artículo para analizarlo y tenerlo
completo.
Un primer problema que se nos puede presentar es cómo leer el polinomio, leer un
polinomio como una tira de caracteres y transformarlo a algo que la computadora
entienda no es algo trivial, a esto se le llama parsear una expresión y esto será algo de
lo que escribiré en otra ocasión, ahora se debe buscar una forma sencilla, se puede
poner una caja de texto para cada coeficiente que se quiera leer y hacer el polinomio
de un tamaño fijo (otra forma es preguntar el tamaño primero y después construir las
cajas de texto que se requieran).
La mejor forma de guardar el polinomio es en un arreglo de números, en la posición

cero se guarda el coeficiente que está solo (es decir que va multiplicado por x0), en la
posición 1 se guarda el coeficiente que acompaña a x1 y así consecutivamente hasta
guardar todo el polinomio.
Una vez que se ha guardado el polinomio deseado, se necesitan varias funciones para
trabajar con él. La primera función será la que evalúe un valor de x en el polinomio,
esta función debe recibir el arreglo del polinomio y devolver la evaluación, para esta
simplemente hay que sumar los término multiplicando cada uno por xn, donde n es la
posición en el arreglo (hay otras formas mejores para hacerlo, pero la idea aquí no es
mostrar la mejor forma, sino la más simple). Veamos el código:
La función que sigue es la que evalúa el polinomio
public double f(Vector fun, double x)

{ double total = (double)((Double)fun.elementAt(0)).doubleValue();
for(int i = 1; i < fun.size(); i++)
{ total+=((double)((Double)fun.elementAt(i)).doubleValue())*Math.pow(x,i*1.0);
}
return total;
} //Fin de la función f
Una segunda función calcula la derivada de un polinomio, esta también recibe un

arreglo con el polinomio y devuelve otro arreglo con su derivada, es muy sencillo pues
cada valor se debe multiplicar por su posición en el arreglo y se elimina el primer
término (el que está multiplicado por x0).
Función que deriva un polinomio
public Vector derivePoli(Vector fun)

{ Vector derivado = new Vector();
derivado.addElement(new double((double)((Double)fun.elementAt(i)).doubleValue()*i));
derivado.trimToSize();
return derivado;
} //derivePoli
La tercera función es la que calcula la división de dos polinomios, esta no es tan

sencilla, se deben recibir dos arreglos que contengan el dividendo y el divisor (son
polinomios) y se puede hacer que solo devuelva el cociente, el residuo o ambos, en
este caso sólo se necesita el residuo, así que se declara un arreglo para guardarlo, se
realiza la división entre los últimos coeficientes del cociente y el divisor (recuerde que
los guardamos empezando con el que acompaña a x0, luego x1,etc), este valor se
multiplica por cada elemento del divisor y se le suma a su correspondiente elemento
del dividendo, estos valores se van guardando en el arreglo declarado.
Función que divide dos polinomios y devuelve el residuo, Si no se puede dividir

devuelve el vector vacío
public Vector dividaPolis(Vector dividendo, Vector divisor)

{ Vector dividendoTemp=new Vector();
Vector cocienteTemp=new Vector();
Vector residuoTemp=new Vector();
dividendoTemp=(Vector)dividendo.clone();
//el grado del dividendo debe ser mayor o igual que el grado del divisor
if(dividendo.size()>=divisor.size()){
//Cantidad de veces que se debe hacer el algoritmo
for(int i=0; i<=(dividendo.size()-divisor.size()); i++)
{ if(residuoTemp.size()>0)
residuoTemp.removeAllElements();
cocienteTemp.addElement(new Double((double)((Double)dividendoTemp.elementAt(
dividendoTemp.size()-1)).doubleValue()/(double)((Double)divisor.elementAt(
divisor.size()-1)).doubleValue()));
//Hace la división de los más grandes para ir llenando el cociente
for(int j=2; j<=dividendoTemp.size(); j++)
{ if (j<=divisor.size())
residuoTemp.addElement(new Double((double)((Double)dividendoTemp.elementAt(
dividendoTemp.size()-j)).doubleValue()-((double)((Double)divisor.elementAt(
divisor.size()-j)).doubleValue()*(double)((Double)cocienteTemp.elementAt(i)
).doubleValue())));
else
residuoTemp.addElement(new Double((double)((Double)dividendoTemp.elementAt(
dividendoTemp.size()-j)).doubleValue()));
} //for
residuoTemp=volverVector(residuoTemp);
dividendoTemp=(Vector)residuoTemp.clone();
} //for
cocienteTemp=volverVector(cocienteTemp);
residuoTemp=quitaCeros(residuoTemp);
if(residuoTemp.size() == 0)
residuoTemp.addElement(new Double(0.0));
} //if
residuoTemp.trimToSize();
cocienteTemp.trimToSize();
return residuoTemp;
//Aquí se puede retornar el cociente o el residuo cambiando la palabra
} //dividaPolis
Se podría hacer otra función para simplificar algunos cálculos, esta recibe un
polinomio y le quita los ceros que se encuentran al final, es decir, que el coeficiente
que acompaña a la variable de mayor grado sea diferente de cero.
Procedimiento que elimina el número de grado mayor si es 0
public Vector quitaCeros(Vector fun)

{ Vector sinCeros=new Vector();
boolean salir=false; //para ver si encontró un valor de cierto grado o si no poner 0
int i=fun.size()-1;
sinCeros=(Vector)fun.clone();
while(i>=0 && salir==false)
{ if((double)((Double)sinCeros.elementAt(i)).doubleValue() == 0.0)
sinCeros.removeElementAt(i);
else
salir = true;
i--;
}
sinCeros.trimToSize();
return sinCeros;
}//quitaCeros
Otra función que puede ser necesaria es la que calcula - f (x), es decir, una función
que le cambie el signo a todos los coeficientes del polinomio. También es necesaria la
que calcula f (- x), esta función debe cambiarle el signo a los campos impares (1, 3, 5,
7, ...).
Función que le cambia el signo a todas las entradas del polinomio
public Vector cambieSignoPoli(Vector fun)

{ Vector signoCambiado=new Vector();
for(int i=0; i<fun.size();i++)
signoCambiado.addElement(new Double(-
1*(double)((Double)fun.elementAt(i)).doubleValue()));
signoCambiado.trimToSize();
return signoCambiado;
}//cambieSignoPoli
Función que le cambia el signo a las potencias impares del polinomio
public Vector cambiaSignoImpares(Vector fun)

{ Vector nuevafun = new Vector();
{ if(i%2 == 0)
nuevafun.addElement((Double)fun.elementAt(i));
else
nuevafun.addElement(new Double(-
1*((double)((Double)fun.elementAt(i)).doubleValue())));
}
return nuevafun;
} //cambiaSignoImpares
Para el teorema de Sturm puede ser útil realizar una función que devuelva el signo de
algún valor que recibe, si es negativo puede devolver -1, si es positivo devuelve 1 y si
es cero entonces devuelve ese valor, es decir 0.
Se deben hacer las funciones que reciben un polinomio y devuelven la cota superior,
esta también es sencilla pues tan solo es de buscar los valores negativos y seguir las
fórmulas.
Busca la cota superior de lo ceros positivos de una función
public double cota(Vector fun)

{ double valorCota=0,mayorCoef;
mayorCoef = Math.abs((double)(((Double)fun.elementAt(0)).doubleValue()));
for(int i=1; i<fun.size()-1;i++)
mayorCoef = Math.max(mayorCoef,
Math.abs((double)(((Double)fun.elementAt(i)).doubleValue())));
if(Math.abs((double)(((Double)fun.elementAt(fun.size()-1)).doubleValue()))==0)
valorCota=0;
else
valorCota=1+(mayorCoef/Math.abs((double)(((Double)fun.elementAt(fun.size()-
1)).doubleValue())));
return valorCota;
} //cota
Falta también la función (esta ya usa las anteriores) que calcula el sistema de Sturm
para un polinomio dado, esta es bastante simple si ya tenemos la funciones anteriores.
Otra función es la que evalúa un valor en el sistema de Sturm y devuelve la cantidad
de cambios de signo que hubo.
La función que sigue es la que hace el teorema de sturm primero deriva el polinomio,
divide el polinomio entre la derivada y recoge el residuo, luego divide la derivada
entre el residuo y nuevamente recoge el residuo y sigue dividiendo los dos últimos
residuos para obtener otro residuo hasta que ya no se pueda. Luego evalúa las
funciones en -infinito y +infinito y con la resta de los cambios de signo se sabe cuantos
ceros hay...
public Vector sturm(Vector fun)

{ boolean salir=false;
int cerosPositivos=0, cerosNegativos=0;
Vector funciones=new Vector();
Vector raices=new Vector();
Vector cotas=new Vector();
if(fun.size()>1)
{ fun = quitaFallo(fun);
if(fun.size()>1)
{ //Se buscan las funciones del sistema de sturm
funciones.addElement(fun);
funciones.addElement(derivePoli(fun));
int i = 0;
while(!salir && i < (fun.size()-2)){
funciones.addElement(cambieSignoPoli(dividaPolis((Vector)funciones.elementAt(i),(Vector)fun
ciones.elementAt(i+1))));
if(((Vector)funciones.elementAt(funciones.size()-1)).size() == 1)
{ if(((double)(((Double)(((Vector)funciones.elementAt(funciones.size()-1))).
elementAt(0))).doubleValue()) == 0){
salir = true;
}
}
i++;
}
//Se buscan los intervalos con el sistema de sturm
cerosPositivos = numeroCerosPositivos(funciones);
cerosNegativos = numeroCerosNegativos(funciones);
cotas = busqueCotas((Vector)funciones.elementAt(0));
raices = unaVectores(busqueRaices(funciones, cerosNegativos, (double)(((Double)cotas.
elementAt(0)).doubleValue()), (double)(((Double)cotas.elementAt(1)).doubleValue())),
busqueRaices(funciones, cerosPositivos,
(double)(((Double)cotas.elementAt(2)).doubleValue()),
(double)(((Double)cotas.elementAt(3)).doubleValue())));
if(falloSturm)
raices.addElement(new Double(0.0));
}
else
raices.addElement(new Double(0.0));
}
raices.trimToSize();
return raices;
}//sturm
En este programa se presentan algunas funciones de las que no hemos hablado, por
ejemplo, el número de ceros positivos y negativos, estas funciones simplemente se
evalúa el sistema de Sturm de cero a infinito para saber los ceros positivos y de menos
infinito a cero para los negativos. busqueRaices es la función que se encarga de ir
acotando las raíces y buscarlas. La función quitaFallo lo que hace es eliminar la
multiplicidad del cero.
Función que busca los intervalos donde hay ceros
public Vector busqueRaices(Vector fun,int numCeros,double minimo,double maximo)

{ Vector raiz=new Vector();
int cerosEncontrados=0;
double ci=minimo, cs=maximo,cotaux;
cotaux=(ci+cs)/2;
if(numCeros==1)
raiz.addElement(new Double(miRedondeo(buscarCero((Vector)fun.elementAt(0),ci,cs))));
else
{ while(cerosEncontrados<numCeros)
{ if(Math.abs(cambiosDeSigno(fun,cotaux)-cambiosDeSigno(fun,ci))==0)
{ ci = cotaux;
cotaux=(cs+cotaux)/2;
}
else if(Math.abs(cambiosDeSigno(fun, ci)-cambiosDeSigno(fun, cotaux)) > 1)
cotaux = (ci+cotaux)/2;
else if(Math.abs(cambiosDeSigno(fun, ci)-cambiosDeSigno(fun, cotaux)) == 1){
raiz.addElement(new Double(miRedondeo(buscarCero((Vector)fun.elementAt(0), ci,
cotaux))));
ci = cotaux;
cotaux = (cs+ci)/2;
cerosEncontrados++;
}
} //while
} //if
raiz.trimToSize();
return raiz;
} //busqueraices
Función que quita el fallo de Sturm, específicamente, la multiplicidad del cero
public Vector quitaFallo(Vector fun)

{ Vector nuevaFun=new Vector();
falloSturm=false;
int i = 0;
while((double)(((Double)fun.elementAt(i)).doubleValue()) == 0.0)
{ i++;
}
if(i > 0)
{ falloSturm=true;
while(i < fun.size())
{ nuevaFun.addElement((Double)fun.elementAt(i));
i++;
}
}
else
nuevaFun = fun;
return nuevaFun;
} //quitaFallo
Por último, se debe hacer el programa principal que maneje todas las funciones
anteriores, es decir que debe leer el polinomio, acomodarlo, calcular el sistema de
sturm e imprimir los ceros.
public boolean action(Event evt,Object arg)

{ if (evt.target instanceof Button)
{ if (evt.target == panel.b1)
{ panel.lista.removeAll();
leaParametros();
funcion = leePoli(formulaPoli);
panel.info.setText("");
if(funcion.size() > 0)
{ funcion = acomodaPoli(funcion);
raicesPol = sturm(funcion);
imprimaCeros(raicesPol);
}
lienzo.repaint();
}
}//if del botón
return true;
}//action
Hay que hacer la función que reciba el polinomio y los valores iniciales y aproxime el
cero, es decir, la función que realice el método de bisección, bisección acelerada o el
de Newton.
Buscar cero es la función que aproxima el cero, utiliza el método de bisección
public double buscarCero(Vector fun, double a, double b)

{ double xmax=b, xmin=a;
double xn = a;
boolean cero = false;
if(f(fun, xmin) == 0)
xn = xmin;
else if(f(fun, xmax) == 0)
xn = xmax;
else{
int i = 0;
while(!cero && i < 1000){
xn = (xmin+xmax)/2;
if(Math.abs(f(fun, xn)) == 0)
cero = true;
else if(f(fun, xn)*f(fun, xmin) < 0)
xmax = xn;
else
xmin = xn;
i++;
}//while
}
//Si este método no funciona se hace la llamada al método newton
if(Math.abs(f(fun, xn)) > 0.000001)
xn = newton(fun, a, b);
return xn;
}//buscarCero
Buscar cero es la función que aproxima el cero, utiliza el método de Newton
public double newton(Vector fun, double a, double b){

double xn = (a+b)/2;
Vector derivada = derivePoli(fun);
boolean cero = false;
int i = 0;
while(!cero && i < 500){
xn = xn-(f(fun, xn)/f(derivada, xn));
if(Math.abs(f(fun, xn)) == 0)
cero = true;
i++;
}//while
return xn;
}//buscarCero
Ya para finalizar, se deben mostrar en pantalla los ceros obtenidos.
//Imprime los ceros en una lista
public void imprimaCeros(Vector raices){
if(raices.size() == 0)
panel.lista.addItem("No hay ceros reales.");
else
for(int i = 0; i < raices.size(); i++)
panel.lista.addItem("x = "+(double)((Double)raices.elementAt(i)).doubleValue());
}//imprimaCeros
Espero que estas reflexiones hayan servido para que puedan realizar un programa
computacional eficaz y eficiente, cualquier duda se pueden comunicar a mi correo y
con mucho gusto les puedo contestar y ayudar con algún problema que tengan. Más
adelante voy a escribir como hacer un parseador simple para leer un polinomio y otro
más complicado sobre como leer cualquier expresión metemática y evaluar un valor
en ella. El programa que se puede bajar tiene mucho más que solo la búsqueda de
ceros, pues grafica el polinomio y marca los ceros en la gráfica.
Programa para ejecutar en MatLab:
Acotación de raíces en polinomios empleando el método de las sucesiones de STURM.

La salida se imprime en pantalla. Pruebe con un ejemplo y saque sus propias
conclusiones
function sturm(p)
n = size(p);
n = n(1,2);
s = sucsturm(p);
a1 = CSRP(p);
if mod(n,2) == 0 l = 1;
else l = 2;
end
p2 = p;
for k = [l:2:n-1]
p2(1,k) = p(1,k)*(-1);
end
a2 = - CSRP(p2);
t = [a2:a1];
f = size(s);
f = f(1,1);
tabla1 = [];
for m = [1:f]
pol = s(m,:);
tabla1(m,:) = sign(polyval(pol,t));
end
[p,q] = size(tabla1);
for i = [1:p]
for j = [1:q]
if tabla1(i,j)==0
tabla1(i,j)= 1;
end
end
end
% menos infinito
neginf=[];
for m = [1:f];
neginf(m)= sign(polyval(s(m,:),-inf));
end
neginf = neginf';
% infinito
posinf=[];
for m = [1:f];
posinf(m)= sign(polyval(s(m,:),inf));
end
posinf = posinf';
% TABLAS
T = [-inf,t,inf]
TABLA = [neginf,tabla1,posinf]
c = contc(TABLA)
%
m = size(c);
m = m(1,2);
for k = [1:m-1]
a = c(k);
b = c(k+1);
if a ~= b
fprintf('\nExiste una raiz en el intervalo: [ %d , %d ] \n',T(k),T(k+1));
end
end
Función que genera las sucesiones de STURM
function s = sucsturm(p)
n = size(p);
n = n(1,2);
f = p;
f1 = polyder(p);
s = [];
s(1,:) = f;
s(2,:) = rellena(n,f1);
q = deconv(f,f1);
r = -f + conv(q,f1);
r = roots(r);
r = poly(r);
t = size(r);
t = t(1,2);
k = 2;
if t == 1 break;
end
while t > 1
k = k+1;
f = f1;
f1 = r;
s(k,:)=rellena(n,f1);
q = deconv(f,f1);
r = -f +conv(q,f1);
r = roots(r);
r = poly(r);
s(k+1,:)= rellena(n,r);
t = size(r);
t = t(1,2);
k = k+1;
end
s;
Función que calcula la Cota Superior de la Raíces negativas del polinomio
function c = CSRP(p)
if sign(p(1,1)) == -1
p = p* (-1);
end
k = 1;
p1 = deconv(p,[1 -k]);
res = polyval(p,k);
if (res > 0) & (todospos(p1)==0) c = k;
break
end
while or((todospos(p1)~= 0),(res < 0))
k = k +1;
p1 = deconv(p,[1 -k]);
res = polyval(p,k);
if (res > 0) & (todospos(p1)==0)
break
end
end
c = k;
Funciones Accesorias
function f = rellena(n,p)
N = size(p);
N = N(1,2);
if N > n fprintf('\n error en llamada a función\n');
pause;
end
f = zeros(1,n);
for k = [n:-1:n-N+1]
f(1,k) = p(1,(k-(n-N)));
end
function c = contc(a)
[fil,col] = size(a);
c = zeros(1,col);
for n = [1:col]
s1 = a(1,n);
for m = [2:fil]
s2 = a(m,n);
if s1 ~= s2
c(n) = c(n)+1;
end
s1= s2;
end
end
function c = todospos(p)
n = size(p);
n = n(1,2);
t = 0;
for k = [1:n]
if sign(p(1,k))== sign(-1) t = 1;
end
end
c = t;
Apéndice
Métodos para aproximar los ceros de un polinomio
Aunque existen muchos métodos para aproximar los ceros de un polinomio, aquí solo
vamos a mencionar tres muy simples y conocidos: el método de la bisección, el
método de bisección acelerada y el método de Newton. Si se quiere investigar más
sobre estos métodos o buscar más se puede investigar en algún libro de Métodos
Numéricos.
Método de bisección
Para poder utilizar este método se necesitan tener dos valores iniciales a y b tales que
f (a) y f (b) tienen signo contrario. Una vez que se cumple esto, entonces tome
c = (a + b ) / 2 y evalúe f (c), si f (c) = 0 entonces ya se encontró que el cero es x = c,
sino si f (a) y f (c) tienen signo contrario, entonces tome a = a y b = c y repita el
proceso, sino tome a = c y b = b. El proceso se repite hasta que f (c) sea menor que un
valor pequeño predeterminado, puede ser f (c) < 0.00001.
Método de bisección acelerada
Este método es similar al anterior, se necesitan dos valores iniciales a y b tales que
f(a) y f(b) tengan diferente signo. La diferencia está en que ahora se toma
f (a )(b − a )
c=a− . Si f(c)=0 entonces ya se encontró que el cero es x=c, sino si f (a) y
f (b) − f (a )
f (c) tienen signo contrario, entonces tome a=a y b=c y repita el proceso, sino tome
a=c y b=b. El proceso se termina cuando f (c) es menor que un valor pequeño
predeterminado, tal vez f (c) = 0.00001.
Método de Newton
Para este método solo se necesita un valor inicial a, el siguiente valor que se toma es
c = a − f ( a ) / f ' ( a ) . Si f (c) = 0 entonces ya encontramos el cero, es x=c, sino se
continúa el proceso hasta que f(c) sea menor que un valor pequeño predeterminado,
puede ser f(c)<0.00001.
El problema que presentan estos métodos es que los valores iniciales deben ser muy
buenos, el primer método es muy lento, el segundo es un poco más rápido y el tercero
es el más sencillo y puede ser el más rápido, sin embargo, puede pasar que se pierda y
nunca encuentre el cero. Si se realizan en un programa de computadora es
conveniente que se ponga un número máximo de iteraciones por si nunca se llega a un
resultado satisfactorio.
Corolario del Teorema de Rolle y la Regla de Descartes

TEOREMA 4 Corolario del Teorema de Rolle
Entre dos raíces reales consecutivas a y b de f (x) = 0, hay un número impar de raíces
reales de f'(x) = 0, una raíz de multiplicidad m es contada m veces.
TEOREMA 5 La Regla de Descartes
El número de raíces reales positivas de una ecuación real es igual al número v de sus
variaciones de signo o es menor que v por un número entero positivo par. Una raíz de
multiplicidad m se cuentan m veces.
Dickson, L. New First Course in the Theory of Equations. Jhon Wiley and Sons, New
York, 1962.
Recursos Bibliográficos
MATHEUS. John H. Fink Kurtis D. Métodos Numéricos con MATLAB. Editorial
Prentice Hall
ALTZ, Franz L. Electronic. Digital. computers: Their use in science and
Engineering. 1958 Academic Press inc. New York.
BURDEN Richard L., J. Douglas Faires; Análisis numérico. tr. Efrén Alatorre
Miguel; Revisión Técnica. Ildefonso. 1998 (Biblioteca USCO. Nro Topográfico:
515 / B949a.)
CHAPRA Steven C., CANALE Raymond P, Numerical Methods for engineers.
McGraw Hill, Inc. 1988. 839p. ISBN 0-07-909944-0.
CHAPRA Steven C., CANALE Raymond P. Métodos numéricos para ingenieros: con
aplicaciones en computadoras personales. 1988 (Biblioteca USCO Nro
Topográfico: 519.5 / C467m)
CONDE S. D, Carl de Boor. Análisis numérico elemental: Un enfoque
algorítmico. Mc. Graw-Hill 1972, (Biblioteca USCO Nro Topográfico: 511.8 /
C761 Biblioteca).
CORMICK MC., John M. and SALVADOR M.C. Numerical Methods in FORTRAN.
1964. Prentice-Hall Inc Englewood Cliffs N:J.
CURTIS, F. Gerald, WHEATLEY, O. Patrick. Análisis numérico con aplicaciones.
Tr. Hugo Villagomez Vasquez. 6 Ed. Pearson Educación. 2000, 698p. ISBN 968-
444-393-5
FADDEEVA, V.N. Computacional methods of linear algebra, Dover Publications.
1969, New York.
GASTINEL Noél; Análisis numérico lineal. tr. Javier Ruiz Fernández de Pinedo.
1975. (Biblioteca USCO Nro Topográfico: 511.7 / G255).
GREENSPAN, D. Theory and solutions of Ordinary Differencial Equations. 1960
The. Mc Millan Co. New York.
KINCAID David y Ward Cheney; Análisis numérico: Las matemáticas del cálculo
científico. tr. Rafael. 1994 (Biblioteca USCO Nro Topográfico: 515 / K51a).
LUTHE. Rodolfo, OLIVERA Antonio, SCHUTZ Fernando, Métodos numéricos. 1986
(Biblioteca USCO Nro Topográfico: 511.7 / L973m).
McCRACKEN, Daniel D., Métodos numéricos y programación fortran: con
aplicaciones en ingeniería y ciencias. 1986. Editorial Limusa. México.
(Biblioteca USCO Nro. Topográfico: 001.6424 / M117).
NAKAMURA Shoichiro; Análisis numérico y visualización gráfica con MATLAB. tr.
Roberto Escalona García. 1998 (Biblioteca USCO N ro Topográfico: 515.1 /
N163a).
NAKAMURA Shoichiro; Métodos numéricos aplicados con software. tr. Oscar
Alfredo Palmas Velasco. Prentice Hall Hispanoamericana S.A. 1995. 570p.
(Biblioteca USCO. Nro. Topográfico: 511.8 / N163m) ISBN 968-880-263-8
NIETO RAMIREZ José A., Métodos numéricos en computadoras digitales. Editorial
Limusa 1980. (Biblioteca USCO Nro Topográfico: 001.64042 / N677).
RALSTON Anthony; Introducción al análisis numérico. tr. Carlos E. Cervantes de
Gortari. Editorial Limusa. Mexico. 1978. 629p. (Biblioteca USCO Nro
Topográfico: 511.7 / R164.)
SCARBOROUGH, J.B Numerical mathematics analysis
SIERRA ROMERO, Alberto. Manual de Métodos Numéricos. Universidad
Tecnológica de Pereira.
SMITH, W. Allen; Análisis numérico. tr. Francisco Javier Sánchez Bernabe; Rev.
Téc. José Luis Turriza Pinto. Prentice Hall Hispanoamericana S.A. 1988. 608p.
(Biblioteca USCO Nro Topográfico: 515 / S664a) ISBN 968-880-119-4.
STANTON, Ralp G. Numerical Methods for Science and Engineering. 1967.
Prentice-Hall Inc. Englewood Cliffs N.J
Bibliografia OnLine:
http://sai.uam.mx/apoyodidactico/mn/
http://www.unalmed.edu.co/~ifasmar/libro.shtml
http://www.iesrodeira.com/metodos_numericos/index-2.htm
http://uprhmate01.upr.clu.edu/~pnm/notas4061/index.htm
http://mailweb.udlap.mx/~ccastane/Analisis_Numerico_html/Lindley.html#RegresaGral1
http://www.ii.uam.es/~pedro/ccii/teoria/
http://webdiee.cem.itesm.mx/web/servicios/archivo/tutoriales/metodos/algoritmos/index.html
http://www.ucsc.cl/~kdt/numerico/index.htm
http://luda.azc.uam.mx/curso2/cp2indic.html
http://www.geocities.com/SiliconValley/Pines/7894/metodos/
http://www.uv.es/~diaz/mn/fmn.html
http://proton.ucting.udg.mx/posgrado/cursos/metodos/temario.html
http://webdiee.cem.itesm.mx/web/servicios/archivo/tutoriales/metodos/algoritmos/index.html
http://www.damtp.cam.ac.uk/user/fdl/people/sd/lectures/nummeth98/contents.htm
http://anamat1.csi.ull.es/anamat_p/Titulaciones/matematicas.htm
http://www.cnice.mecd.es/Descartes/
http://mathworld.wolfram.com/
http://www.elprisma.com/
http://www-gap.dcs.st-and.ac.uk/~history/index.html
http://www.ciencia-hoy.retina.ar/indice.htm
http://omega.ilce.edu.mx:3000/sites/ciencia/html/fisica.htm
http://books.pdox.net/
http://arxiv.org/
http://pessoal.sercomtel.com.br/matematica/resumos.htm
http://www.itlp.edu.mx/publica/tutors.htm
http://www.rinconmatematico.com/libros.htm

2013 09 14 Sturm

Cargado por

Copyright:

Formatos disponibles

2013 09 14 Sturm

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

2013 09 14 Sturm

Cargado por

Copyright:

Formatos disponibles

2013

Sistema de polinomios de Sturm

Yamil Armando Cerquera Rojas

Cuando se encuentra ante una situación problemática, a veces se plantean

Resolver una ecuación o un sistema de ecuaciones son problemas que se plantearon en

En realidad, el Álgebra comienza cuando los matemáticos se empiezan a interesar por

Tras muchos milenios, las ecuaciones matemáticas de la forma + = 0 se resuelven

Desarrollo del Concepto

En la ecuación: + = , a, b y c son coeficientes, es la incógnita.

En la ecuación 5 − 4 = 16, los coeficientes son los enteros 5, -4, y 16 y la incógnita es

En la ecuación anterior v es la incógnita y el valor = 2000 es la solución.

2. Un norteamericano llega a Neiva – Huila - Colombia y manifiesta que la

En: ax + 8 = 3 , el valor = −5/ es la solución.

En 4 – 3 + 2 = 2 los valores =0y = 3/4 son soluciones.

Clasificación de las ecuaciones con una incógnita:

6 + 34 = 5, Es una ecuación de primer grado.

Raíces de ecuaciones: introducción2

Resolver una ecuación es encontrar el ó los valores de la incógnita que satisface la

Por ejemplo la ecuación: 500 = 450 + , se satisface para = 50

– No se puede pretender poder resolver el problema para funciones más

– La mayoría de veces es posible traducir de alguna manera el problema original

En situaciones reales la solución de la ecuación debe tener sentido en el contexto en

Nótese los siguientes casos:

a) Pertinencia de la solución: Se quiere repartir equitativamente 24 dulces a 5

La ecuación anterior no tiene solución en los naturales (N).

b) Existencia de la solución: La ecuación 4 x + 5 = 2, no tiene solución en los

La ecuación 4 * x * x = −7 , no tiene solución en los reales (R) ya que no existe ningún

c) Infinitas soluciones: La ecuación 2 + + = 2( + 1), es una ecuación que es

En este libro, se trata de profundizar sobre raíces de polinomios con coeficientes en Q

Además presupone cierta familiaridad con la teoría general básica de polinomios en

Historia de las ecuaciones lineales.

Para llegar al actual proceso de resolución de la ecuación + = , han pasado más

Los egipcios nos dejaron en sus papiros4 (sobre todo en el de Rhid-1.650 AC y el de

Las ecuaciones más utilizadas por los egipcios eran de la forma:

Donde a, b y c eran números conocidos y x la incógnita que ellos denominaban aha o

"Un montón y un séptimo del mismo es igual a 24".

En notación moderna, la ecuación sería: x + 1 x = 24

Suponga que fuera 7 la solución, al sustituir en la x daría: 7 + 1 7 = 8 , y como nuestra

"Transeúnte, ésta es la tumba de Diophante: es él quien con esta sorprendente

Lo resolvían utilizando el método de la falsa posición, como los egipcios.

Posteriormente, Brahmagupta7 (siglo VII) expresa, ya de forma sincopada, cómo

Dada la ecuación ax + b = cx + d , la solución vendrá dada dividiendo la diferencia de

Sea la ecuación ax + b = 0 y suponga dos valores para la x:

Por otra parte, eliminando a del sistema de ecuaciones anteriores se tiene.

Ecuación 1*n → amn + bn = pn

Y dividiendo ambos resultados, − a/b = (p − q)/(pn − qm)

O también − b /a = (pn − qm)/(p − q)

Siendo esto último el valor de x.

Vea un ejemplo. Sea la ecuación 5 x − 10 = 0 , si se toma como valor de x: x = 3 y x = 4

Se tiene que: x = (10 * 3 − 5 * 4) /(10 − 5) = (30 − 20) / 5 = 10 / 5 = 2

Este principio fue posteriormente presentado en una forma ligeramente modificada

El método puede ser sintetizado como sigue:

1. Consideran dos valores cualesquiera de la incógnita m, n.

A partir de aquí se dedican al estudio de ecuaciones de grado superior.

Historia de los sistemas de ecuaciones lineales.

Un ejemplo tomado de una tablilla babilónica9 plantea la resolución de un sistema de

1/4 ! ℎ#$ + %&!'()#* = 7 + !&,

Para resolverlo comienzan asignando el valor 5 a una mano y observaban que la

Restando la segunda de la primera, se obtiene 3 = 18, es decir = 6 e - = 4.