Решить уравнение методом половинного деления - QBasic

Формулировка задачи:

решить уравнение x3 + x2 + x + 1 = 0; на отрезке [–2,1] методом половинного деления на QBasic помогите составить

Решение задачи: «Решить уравнение методом половинного деления»

Листинг программы

REM
REM OTBET: x = -1
REM
 
DECLARE FUNCTION f! (x!)
DECLARE FUNCTION D! (a!, b!)
 
CLS
CONST epsilon = .000001
 
PRINT D(-2, 1)
END
 
FUNCTION D (a, b)
  
   c = (a + b) / 2
 
   IF f(a) * f(c) <= 0 THEN b = c ELSE a = c
   IF ABS(a - b) > epsilon THEN D = D(a, b) ELSE D = c
END FUNCTION
 
FUNCTION f (x)
   f = x ^ 3 + x ^ 2 + x + 1
END FUNCTION

Объяснение кода листинга программы

Задача программы - решить уравнение методом половинного деления.
Уравнение, которое нужно решить: x^3 + x^2 + x + 1 = 0.
Используются две функции: D(a, b) и f(x).
Функция D(a, b) используется для деления отрезка [a, b] пополам и нахождения точки, в которой происходит пересечение графика функции f(x) с осью x.
Функция f(x) вычисляет значение функции x^3 + x^2 + x + 1.
В начале программы выводится -2 и 1, это начальные приближения для решения уравнения.
Далее в цикле происходит деление отрезка [a, b] пополам и проверка, в какой половине находится решение уравнения.
Если решение находится в левой половине, то значение переменной b остается неизменным, а значение переменной a уменьшается до тех пор, пока разница между a и b не станет меньше заданной точности epsilon.
Если решение находится в правой половине, то значение переменной a остается неизменным, а значение переменной b уменьшается до тех пор, пока разница между a и b не станет меньше заданной точности epsilon.
Значение переменной c равно среднему арифметическому a и b.
Значение переменной f_a и f_c вычисляются как значения функции f(x) в точках a и c соответственно.
Если f_a * f_c <= 0, то решение находится в левой половине, иначе в правой.
Если |a - b| > epsilon, то цикл продолжается, иначе решение уравнения равно c и выводится на экран.
Значение переменной f вычисляется как значение функции f(x) в точке x.
В конце программы выводится значение переменной f(-2, 1).