BisectionMethod-Python/BisectionMethod.py at master · josgard94/BisectionMethod-Python · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
#
#   Autor: Edgard Díaz
#
#   Esta clase contiene la lógica del método de bisección para aproximación de  raices a partir de un intervalo [a,b]
#
#

import math as m

class BisectionMethod:

    def __init__(self,a,b): #constructor de la clase
        self.a = a;
        self.b = b;

    def punto_medio(self, a, b): #calcula el punto medio del intervalo [a, b]
        return (a+b)/2

    def funcion(self,x): #Función matematica a la que se le aproxima la raiz
        return 10 - (m.pow(x,2))
        #return m.pow(x,3) + 4 * pow(x,2) - 10;
        #return pow(x,3) + 4 * pow(x,2)

    def Aproximar_Raiz(self): #Lógica del metodo de Bisección
        n = 1
        iteraciones = 1000
        print("n \t\ta \t\t\tb \t\t\tc \t\t  f(a) \t\t  f(b) \t\t  f(c)")
        while (True):
            c = self.punto_medio( self.a, self.b); #calcula el punto medio del intervalo dado
            #print(c, self.funcion(self.a), self.funcion(self.b), self.funcion(c))
            print(str(n),"\t","{0:.5f}".format(self.a),"\t","{0:.5f}".format(self.b),"\t","{0:.5f}".format(c), "\t{0:.5f}".format(self.funcion(self.a)),"\t{0:.5f}".format(self.funcion(self.b)),"\t{0:.5f}".format(self.funcion(c)))
            if n > iteraciones: #Criterio de paro en caso de que el algoritmo no llege a converger con el intervalo dado
                print("El programa no converge con el intervalo dado. . .")
                break
            else:
                if m.fabs(self.funcion(c)) <= 0.01: #Criterio de paro el error absoluto es muy cercano a un valor definido

                    print("\nResultado: \nRaiz aproximada: {0: .5f}".format(c)," Iteraciones realizadas: ", str(n), " Error: {0: .5f}".format(m.fabs(self.funcion(c))))
                    break

                else: #criterio de asignación para el nuevo intervalo [a, b]
                    if self.funcion(c) * self.funcion(self.a) > 0: #Asignar el punto medio al intervalo a
                        self.a = c
                    elif self.funcion(c) * self.funcion(self.b) > 0: #Asignar el punto medio al intervalo b
                        self.b = c

            n = n + 1 #incrementar el contador de iteraciones