UNIVERSIDAD TECNOLÓGICA DE PANAMÁ CAMPUS

DR. VÍCTOR LEVI SASSO FACULTAD DE INGENIERÍA

MECÁNICA MÉTODOS NÚMERICOS (7895)

ASIGNACIÓN # 3:
“METODOS CERRADOS Y ABIERTOS”

PROFESOR:

Dr. EDSON GALAGARZA

INTEGRANTES:
JOACIM GÓMEZ: 8-955-1073
(joacim.gomez@utp.ac.pa)

ANDRÉS OLARTE: E-8-157021

(andres.olarte@utp.ac.pa)

ROGELIO SINCLAIR: 8-950-1836

(rogelio.sinclair@utp.ac.pa)

ENRIQUE LAU: 8-950-1682

(enrique.lau@utp.ac.pa)

EDWIN MARÍN: 8-942-2461

(edwin.marin1@utp.ac.pa)

PERIODO:
SEGUNDO SEMESTRE DEL 2021

GRUPO: 1IM-122

FECHA DE ENTREGA: 3 DE NOVIEMBRE DE 2021

Universidad Tecnológica de Panamá Grupo: 1IM-122

INDICACIONES GENERALES
La asignación esta dividida en 2 partes, métodos cerrados y métodos abierto y cada parte
tiene 2 algoritmos. Por lo tanto, la tarea expira en 2 semanas luego de asignada.
• Debe ´presentar una Grafica del polinomio y, en base a ella, debe escoger los valores
iniciales (Cada grupo escoge los valores a discreción).
• Todas las raíces calculadas se deben evaluar en la función original para comprobar
su exactitud.
• El estudiante tiene la libertad de mostrar los datos y resultados en el formato que
desee, siempre que cumpla con los requisitos solicitados y presente todos los datos
necesarios.
• Ambos códigos deben alcanzar un error de aproximación menor a 0.01% para cada
raíz.

REQUISITOS DE LOS CÓDIGOS

Cada uno de los códigos a presentar debe cumplir con los siguientes requisitos:
1. Permitir introducir la función a evaluar manualmente.
i. En el caso del método de Newton-Raphson, el programa debe permitir
también introducir la derivada de esta función manualmente.

2. Permitir introducir los parámetros del cálculo

i. Limites de los intervalos (métodos cerrados, xl y Xu) o los valores iniciales
(métodos abiertos).
ii. Error máximo aceptable(εs)
iii. Máximo número de iteraciones

3. Calcular e imprimir una raíz de la función a la vez.

4. Imprimir el error de aproximación y el número de iteraciones que se usaron.
5. Realizar una comprobación de la raíz calculada en la función original.

Problemas

Parte 1- Métodos Cerrados

1-Método de Bisección
1. Desarrolle un código en SCILAB que implemente el método de la bisección para
encontrar las raíces. El programa debe estar basado en el pseudocódigo de la figura
5.11 del libro (Pagina 105) que se nuestra a continuación.
Universidad Tecnológica de Panamá Grupo: 1IM-122

Figura N°1. Pseudocódigo de la figura 5.11 del libro (Página 105).

2-Método de la Falsa Posición Modificado

1. Desarrolle un código en SCILAB que implemente el método de la falsa posición
modificado encontrar raíces. El programa debe estar basado en el pseudocódigo de la
figura 5.15 del libro (página 109) que se muestra a continuación:

Figura N°2. Pseudocódigo de la figura 5.15 del libro (página 109)

Universidad Tecnológica de Panamá Grupo: 1IM-122

Parte 2- Métodos Abiertos

Para los problemas a continuación, los algoritmos a construir utilizan como base el
pseudocódigo de la Figura 6.4 (Pseudocódigo para la iteración de punto fijo) del libro
(Chapra & Canale, 7ed), pagina 117, que también se muestra a continuación:

Figura N°3. Pseudocódigo de la Figura 6.4 (Pseudocódigo para la iteración de punto fijo) del libro
(Chapra & Canale, 7ed), página 117

NOTA: Este código es una forma base para el algoritmo de iteración de punto fijo, usted
debe modificarlo para que funcione con los 2 métodos de abajo solicitados.

1-Método de Newton-Raphson
Desarrolle un código en SCILAB que implemente el método de Newton Raphson para
encontrar raíces.

2-Método de la Secante

Desarrolle un código en SCILAB que implemente el método de la Secante para encontrar

raíces.
Universidad Tecnológica de Panamá Grupo: 1IM-122

RESULTADOS
Parte 1- Métodos Cerrados

1-Método de Bisección

Figura N°4. Código para el método de bisección elaborado en Scilab para la resolución del
problema.
Universidad Tecnológica de Panamá Grupo: 1IM-122

Figura N°5. Resultados obtenidos para el método de bisección.

2. Método de la Falsa Posición Modificado

Figura N°6. Código para el método de falsa posición elaborado en Scilab para la resolución del
problema.
Universidad Tecnológica de Panamá Grupo: 1IM-122

Figura N°7. Resultados obtenidos para el método de falsa posición.

Parte 2- Métodos Abiertos

1-Método de Newton-Raphson

Figura N°8. Código para el método de Newton-Raphson elaborado en Scilab para la resolución del
problema.
Universidad Tecnológica de Panamá Grupo: 1IM-122

Figura N°9. Resultados obtenidos para el método de Newton-Raphson para X0=4.52.

Figura N°10. Resultados obtenidos para el método de Newton-Raphson para X0=4.54.

Universidad Tecnológica de Panamá Grupo: 1IM-122

2-Método de la Secante

Figura N°11. Código para el método de la secante elaborado en Scilab para la resolución del
problema.

Figura N°12. Resultados obtenidos para el método de la secante para X0=4.52.

Universidad Tecnológica de Panamá Grupo: 1IM-122

Figura N°13. Resultados obtenidos para el método de la secante para X0=4.54.

EXPLICACIÓN DE LOS CÓDIGOS

Parte 1- Métodos Cerrados
1-Método de Bisección
En este codigo, lo primero que procedimos a hacer fue definir la funcion de la biseccion. Una vez
definida esta funcion procedimos a declarar las variables que vamos a utilizar, luego de declarar
las variables entramos a un ciclo while en donde si iter es mayor que imax ea sera igual a 0; luego
del ciclo while vamos escribiendo las fórmulas del método de bisección y luego entramos a un
ciclo for en donde si xr es menos o mayor a 0 se utiliza esa formula, luego se entra nuevamente a
un ciclo if en donde si test es menir a 0 xu es igual a xr y si test es mayor a 0 x1 es igual a xr y f1
es igual a fr, si no es ninguno de los casos entonces ea es igual a 0.

En la consola podemos ver los siguientes datos:

Los valores correspondientes a los parámetro de entrada (de izquierda a derecha) son:
1. Nombre de la función que se va a evaluar.
2. Valor incial del rango a evaluar.
3. Valor final del rango a evaluar.
4. Número máximo de iteraciones permitidas.
5. Error máximo aceptable.

Los valores correspondientes a los parámetro de salida (de izquierda a derecha) son:
1. Valor de la raíz hallada.
2. Comprobación de la raiz calculada en la función original.
3. Error absoluto calculado.
4. Número de iteraciones realizadas.
Universidad Tecnológica de Panamá Grupo: 1IM-122

2. Método de la Falsa Posición Modificado

En este codigo, procedimos primero a declarar la funcion de la falsa posicion, luego declaramos
las variables que vamos a utilizar junto con algunas formulas, despues de esto entramos a un ciclo
while mientras que ea sea mayor que es y si iter es mayor que imax entonces ea es igual a 0. Luego
procedimos a escribir la formula para ese metodo y entramos a un ciclo if en donde mientras que
xr sea menor o mayor a 0 ea sera igual a esa formula usando los valores ya declarados previamente,
luego entramos a un un ciclo if en donde si test es menor a 0 xu es igual a xr, fu es igual a f por
xu, iu es igual a 0 y il es igual a il mas 1, luego si test es mayor a 0 x1 es igual a xr, f1 es igual f
por x1, il es igual a 0, iu es igual a iu mas 1 y por ultimo si iu es mayor a igual a 2 fu es igual a fu
entre 2 sino, ea es igual a 0

En la consola podemos ver los siguientes datos:

Los valores correspondientes a los parámetro de entrada (de izquierda a derecha) son:
1. Nombre de la función que se va a evaluar.
2. Valor incial del rango a evaluar.
3. Valor final del rango a evaluar.
4. Número máximo de iteraciones permitidas.
5. Error máximo aceptable.

Los valores correspondientes a los parámetro de salida (de izquierda a derecha) son:
1. Valor de la raíz hallada.
2. Comprobación de la raiz calculada en la función original.
3. Error absoluto calculado.
4. Número de iteraciones realizadas.

Parte 2- Métodos Abiertos

1-Método de Newton-Raphson
En este código, inicialmente se declaran la función de Newton Raphson y se declaran las variables
que se utilizarán. Podemos ver que declaramos que xr = x0 por lo que las funciones se estarán
evaluando en el punto x0. Luego, iniciamos un ciclo while en el cual, si iter es mayor a imax,
entonces ea será 0, igual que en los programas anteriores. Cuando se cumpla la condición de que
ea es 0. Se realiza el primer cálculo utilizando la fórmula iterativa de Newton en el cual se resta el
valor anterior obtenido menos el f en el valor anterior entre la derivada de la función en el valor
anterior. En esta variable xr se guardarán todas las aproximaciones de la variable. El error se
obtiene a partir de la fórmula del valor absoluto de la resta entre la aproximación actual y la anterior
entre la actual por cien.

En la consola podemos ver los siguientes datos:

Los valores correspondientes a los parámetros de entrada (de izquierda a derecha) son:
1. Nombre de la función que se va a evaluar.
2. Representa la derivada.
3. Límite inferior del intervalo.
4. Número máximo de iteraciones permitidas.
5. Error máximo aceptable.
Universidad Tecnológica de Panamá Grupo: 1IM-122

Los valores correspondientes a los parámetros de salida (de izquierda a derecha) son:
1. Valor de la raíz hallada.
2. Comprobación de la raiz calculada en la función original.
3. Error absoluto calculado.
4. Número de iteraciones realizadas.

2-Método de la Secante
Para este código, inicialmente se declaran la función de Newton Raphson y se declaran las
variables que se utilizarán. En el programa asignaremos valores iniciales a nuestras variables y
mediante un ciclo while implementamos el método de la secante en el cual se estarán almacenando
nuevos valores en cada iteración para xr, mediante la fórmula que se puede ver en el programa
correspondiente, como el programa anterior, el error se obtiene del mismo modo.

En la consola podemos ver los siguientes datos:

Los valores correspondientes a los parámetros de entrada (de izquierda a derecha) son:
1. Nombre de la función que se va a evaluar.
2. Límite inferior del intervalo.
3. Número máximo de iteraciones permitidas.
4. Error máximo aceptable.

Los valores correspondientes a los parámetros de salida (de izquierda a derecha) son:
1. Aproximación de la raíz actual.
2. Comprobación de la raíz calculada en la función original.
3. Error absoluto calculado.
4. Número de iteraciones realizadas.
Universidad Tecnológica de Panamá Grupo: 1IM-122

CONCLUSIONES
Mediante esta asignación, la cual es una continuación de los métodos empleados para
encontrar raíces, como lo son, los métodos cerrados y abiertos, los cuales en esta asignación
nos pudimos percatar de las diferencias que existen entre estos métodos, en las cuales algunas
ya se habían mencionado en clases anteriores, pero mediante la realización de códigos en
Scilab para emplear los métodos que corresponden tanto a abiertos como cerrados, pudimos
darnos de cuenta en realidad de sus diferencias principales. Entre estas diferencias podemos
mencionar que el empleo de los métodos cerrados se necesitan conocer dos valores iniciales
que encierren a la raíz, mientras que los métodos abiertos solo necesitan un solo valor inicial.
Otra diferencia que hay entre estos métodos cerrados necesitan de muchas iteraciones hasta
llegar a la precisión deseada, por otro lado, los métodos abiertos convergen de forma más
rápida que las sucesiones generadas por los métodos cerrados.

Respecto a los códigos y los resultados obtenidos se pueden destacar los inconvenientes a la
hora de compilar el código e introducir los datos de los diferentes problemas, ya que se nos
presentaron problemas de declaración de variables así como problemas de poder correr
correctamente las iteraciones necesarias para llegar al resultado final del problema. En cuanto
a los resultados obtenidos podemos observar que se llega a la misma respuesta utilizando
ambos métodos abiertos, sin embargo se requieren muchas menos iteraciones en el método
de bisección (8 iteraciones) que en el método de falsa posición (48 iteraciones), pero el
método de falsa posición lanzaba menor error; mismo escenario se presenta en los métodos
cerrados, ya que ambos métodos también llegan al mismo resultado, pero el método de la
secante requirió de menos iteraciones para ambas X0 que el método de Newton-Raphson,
sin embargo el método de Newton-Raphson arrojó menor error.