Как должен выглядеть алгоритм для метода половинного деления на языке программирования?

Question

[Аня Фролова]

Алгоритм метода половинного деления:

Вводим отрезок [a,b]

Вводим eps (достаточно маленькая величина)

Вычисляем f(a) и f(b)

Цикл пока abs(a-b)>eps

c=(a+b)/2

Если f(a)f(b)<0

b=c

f(b)=f(c)

Иначе

a=c

f(a)=f(c)

Всеесли

Всецикл

x=c

Вывод х

А как он будет выглядеть на языке программирования в максимально коротком виде?
Любой язык, интересует возможность уменьшения количества кода при сохранении его работоспособности.

Comments

[Василий Банников]

> языке программирования
А на каком именно?

[Kovalsky]

"Любой язык, интересует возможность уменьшения количества кода при сохранении его работоспособности" - это вам на code golf,

[Василий Банников]

Kovalsky, оу. Что-то я в глаза долблюсь и не увидел эту строку.

[Василий Банников]

> Любой язык
Думаю, в вольфраме или каком-нибудь другом математическом языке уже есть готовая функция, которая это делает.

PS: Да BisectionMethodFindRoot

Answer 1

[Василий Банников]

Не уверен, на каком языке вы хотите, по тому опишу на F#

open System

let epsilon = Double.Epsilon
let abs: float -> float = Math.Abs
let equals a b = abs (a - b) <= epsilon

let rec findRoot f segment =
    let bisect f (a, b) =
        let avg = (a + b) / 2.0
        if ((f a) * (f b)) < 0.0 then // знак f(a) != знак f(b)
            a, avg
        else
            avg, b

    let a, b = bisect f segment
    if equals a b then a else findZero f (a, b)

let segment = -10.0, 10.0 // Отрезок [-10,10]
let f x = x + 5.0 // Пусть f(x) = x + 5
let result = findRoot f segment
printfn "%f" result

Или ещё короче на Wolfram language

ResourceFunction["BisectionMethodFindRoot"][x + 5, {x, -10, 10}, 5, 100]

https://resources.wolframcloud.com/FunctionReposit...