Page 219 - Chapra y Canale. Metodos Numericos para Ingenieros 5edición_Neat

7.7 LOCALIZACIÓN DE RAÍCES CON BIBLIOTECAS Y PAQUETES DE SOFTWARE 195

TABLA 7.2 Rutinas de IMSL para localizar raíces.

Categoría Rutina Capacidad
Raíces de una función
ZREAL Encuentra los ceros reales de una función real
con el método de Müller.
ZBREN Encuentra un cero de una función real que
cambia de signo en un intervalo dado.
ZANLY Encuentra los ceros de una función compleja
univariada usando el método de Müller.
Raíz de un sistema de ecuaciones
NEQNF Resuelve un sistema de ecuaciones no lineales
usando un algoritmo híbrido de Powell
modifi cado (una variación del método de
Newton) y una aproximación en
diferencias fi nitas del Jacobiano.
NEQNJ Resuelve un sistema de ecuaciones no lineales
usando un algoritmo híbrido de Powell
modifi cado (una variación del método de
Newton) con el Jacobiano propuesto por
el usuario.
NEQBF Resuelve un sistema de ecuaciones no lineales
usando la actualización de la secante
factorizada y una aproximación en
diferencias fi nitas del Jacobiano.
NEQBJ Resuelve un sistema de ecuaciones no lineales
usando la actualización de la secante
factorizada con el Jacobiano propuesto
por el usuario.
Raíces de polinomios
ZPORC Encuentra los ceros de polinomios con
coefi cientes reales con el algoritmo de
Jenkins-Traub.
ZPLRC Encuentra los ceros de polinomios con
coefi cientes reales con el método de
Laguerre.
ZPOCC Encuentra los ceros de polinomios con
coefi cientes complejos con el algoritmo
de Jenkins-Traub.

Donde

F = Una función definida por el usuario para la cual van a encontrarse las raíces
ERABS = Primer criterio de terminación, termina si |ƒ(x )| < ERABS. (Entrada)
i
ERREL = Segundo criterio de terminación, termina si |(x – x )/x |< ERREL. (Entrada)
i–1
i
i
EPS = Véase ETA. (Entrada)
ETA = Criterio de extensión para raíces múltiples. (Entrada)
Si la raíz x se ha calculado y |x – x | < EPS, donde x es una raíz previamen-
j
j
i
i
te calculada, se reinicia el cálculo con un nuevo valor inicial de x + ETA.
i
NR = Número de raíces a ser encontradas. (Entrada)
IMAX = Máximo número permitido de iteraciones por raíz. (Entrada)
6/12/06 13:51:29
Chapra-07.indd 195
Chapra-07.indd 195 6/12/06 13:51:29