Método de bisección para calcular raíces de una función
El método de bisección para calcular una raíz de una función en el intervalo [a,b] se basa en el teorema de Bolzano:
Si f(x) es una función continua en el intervalo \([a, b]\), y si, además, en los extremos del intervalo la función \(f\) toma valores de signo opuesto \((f(a)f(b) < 0)\), entonces existe al menos un valor \(c\) en \((a, b)\) para el que \(f(c) = 0\)».
El método para calcular una raíz de la función \(f\) en el intervalo \([a,b]\) con un error menor que \(e\) consiste en tomar el punto medio del intervalo \(c = \frac{a+b}{2}\) y considerar los siguientes casos:
- Si \(|f(c)| < e\), hemos encontrado una aproximación del punto que anula \(f\) en el intervalo con un error aceptable.
- Si \(f(c)\) tiene signo distinto de \(f(a)\), repetir el proceso en el intervalo \([a,c]\).
- Si no, repetir el proceso en el intervalo \([c,b]\).
Definir la función
|
1 |
biseccion :: (Double -> Double) -> Double -> Double -> Double -> Double |
tal que biseccion f a b e es una aproximación del punto del intervalo [a,b] en el que se anula la función f, con un error menor que e, calculada mediante el método de la bisección. Por ejemplo,
|
1 2 3 4 |
biseccion (\x -> x^2 - 3) 0 5 0.01 == 1.7333984375 biseccion (\x -> x^3 - x - 2) 0 4 0.01 == 1.521484375 biseccion cos 0 2 0.01 == 1.5625 biseccion (\x -> log (50-x) - 4) (-10) 3 0.01 == -5.125 |
Soluciones
A continuación se muestran las soluciones en Haskell y las soluciones en Python.
|
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 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 |
module Metodo_de_biseccion_para_calcular_ceros_de_una_funcion where import Test.Hspec (Spec, hspec, it, shouldBe) -- 1ª solución -- =========== biseccion :: (Double -> Double) -> Double -> Double -> Double -> Double biseccion f a b e | abs (f c) < e = c | f a * f c < 0 = biseccion f a c e | otherwise = biseccion f c b e where c = (a+b)/2 -- 2ª solución -- =========== biseccion2 :: (Double -> Double) -> Double -> Double -> Double -> Double biseccion2 f a b e = aux a b where aux a' b' | abs (f c) < e = c | f a' * f c < 0 = aux a' c | otherwise = aux c b' where c = (a'+b')/2 -- Verificación -- ============ verifica :: IO () verifica = hspec spec spec :: Spec spec = do it "e1" $ biseccion (\x -> x^2 - 3) 0 5 0.01 `shouldBe` 1.7333984375 it "e2" $ biseccion (\x -> x^3 - x - 2) 0 4 0.01 `shouldBe` 1.521484375 it "e3" $ biseccion cos 0 2 0.01 `shouldBe` 1.5625 it "e4" $ biseccion (\x -> log (50-x) - 4) (-10) 3 0.01 `shouldBe` -5.125 it "e5" $ biseccion2 (\x -> x^2 - 3) 0 5 0.01 `shouldBe` 1.7333984375 it "e6" $ biseccion2 (\x -> x^3 - x - 2) 0 4 0.01 `shouldBe` 1.521484375 it "e7" $ biseccion2 cos 0 2 0.01 `shouldBe` 1.5625 it "e8" $ biseccion2 (\x -> log (50-x) - 4) (-10) 3 0.01 `shouldBe` -5.125 -- La verificación es -- λ> verifica -- -- e1 -- e2 -- e3 -- e4 -- e5 -- e6 -- e7 -- e8 -- -- Finished in 0.0008 seconds -- 8 examples, 0 failures |
|
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 48 49 50 51 |
from math import cos, log from typing import Callable # 1ª solución # =========== def biseccion(f: Callable[[float], float], a: float, b: float, e: float) -> float: c = (a+b)/2 if abs(f(c)) < e: return c if f(a) * f(c) < 0: return biseccion(f, a, c, e) return biseccion(f, c, b, e) # 2ª solución # =========== def biseccion2(f: Callable[[float], float], a: float, b: float, e: float) -> float: def aux(a1: float, b1: float) -> float: c = (a1+b1)/2 if abs(f(c)) < e: return c if f(a1) * f(c) < 0: return aux(a1, c) return aux(c, b1) return aux(a, b) # Verificación # ============ def test_biseccion() -> None: assert biseccion(lambda x : x**2 - 3, 0, 5, 0.01) == 1.7333984375 assert biseccion(lambda x : x**3 - x - 2, 0, 4, 0.01) == 1.521484375 assert biseccion(cos, 0, 2, 0.01) == 1.5625 assert biseccion(lambda x : log(50-x) - 4, -10, 3, 0.01) == -5.125 assert biseccion2(lambda x : x**2 - 3, 0, 5, 0.01) == 1.7333984375 assert biseccion2(lambda x : x**3 - x - 2, 0, 4, 0.01) == 1.521484375 assert biseccion2(cos, 0, 2, 0.01) == 1.5625 assert biseccion2(lambda x : log(50-x) - 4, -10, 3, 0.01) == -5.125 print("Verificado") # La comprobación es # >>> test_biseccion() # Verificado |