metoda bisekcji

Question

Module Module1

    Sub Main()
        Dim a, b, s, eps As Decimal
        Dim fa, fb, fs As Decimal

        Console.WriteLine("Podaj początek i koniec przedziału:")
        a = Console.ReadLine()
        b = Console.ReadLine()
        eps = 0.1
        fa = a / 2 - 5
        fb = b / 2 - 5

        If fa * fb < 0 Then
            Do
                s = (a + b) / 2
                fs = s / 2 - 5
                
                If fa * fs < 0 Then
                    b = s
                Else
                    a = s
                End If

            Loop While (a - b) > eps

            Console.WriteLine("Miejsce zerowe fukncji =" + CStr(s))

        Else
            Console.WriteLine("Funkcja nie ma miejsc zerowych w podanym przedziale")
        End If

        Console.ReadKey()

    End Sub


End Module

 

Witam mam problem z metodom bisekcji. Funkcja może być dowolna ale przyjąłem f(x)=x/2-5, niestety mój kod nie zawsze zwraca prawidłową wartość i nie wiem co robię źle.

Answer 1

Problem może tkwić w tym:

            Loop While (a - b) > eps

ponieważ gdy a < 0 i b > 0 to (a - b) jest zawsze mniejsze od eps więc ja bym dał tam wartość bezwzględną abs(a - b) bo np. dla a = -10 i b = 1 warunek (a - b) > eps nie jest spełniony a spełniony powinien być bo -10 - 1 = -11 więc wartość bezwzględna z tego jest większa od narzuconej tolerancji eps