Решить уравнение методом деления отрезка пополам - Pascal ABC (16567)

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

Решить уравнение методом деления отрезка пополам. Каждый этап решения выводить в файл: sin(1-x)=x на отрезке

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

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

uses crt;
function F(x:real):real;
begin
F:=sin(1-x)-x
end;
 
var a,b,c,x,e:real;
    k:integer;
    t:text;
begin
a:=0;
b:=pi/2;
e:=0.0001;
writeln('Решение уравнения sin(1-x)-x=0');
writeln('на интервале [0;pi/2] с погрешностью 0.0001');
assign(t,'output.txt');
rewrite(t);
writeln(t,'Решение уравнения sin(1-x)-x=0');
writeln(t,'на интервале [0;pi/2] с погрешностью 0.0001');
k:=0;
repeat
 c:=(a+b)/2;
 k:=k+1;
 if F(a)*F(c)<=0 then b:=c
 else a:=c;
 writeln(t,'Приближение ',k:2,'  x=',c:7:5);
until abs(b-a)<e;
x:=(a+b)/2;
writeln(t,'С точностью 0.0001  x=',x:7:5);
close(t);
write('Результат записан в файл output.txt');
end.

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

В начале кода подключается библиотека crt, которая используется для работы с числами с плавающей точкой в языке Pascal ABC.
Создается функция F, которая принимает аргумент x типа real и возвращает значение типа real. Внутри функции происходит вычисление синуса значения x минус x.
Создаются переменные a, b, c, x и e типа real. Переменная k является целочисленной.
Выводится сообщение о решении уравнения sin(1-x)-x=0 на интервале [0;pi/2] с погрешностью 0.0001.
Устанавливается файл вывода output.txt.
Запускается цикл repeat, который выполняется до тех пор, пока разница между текущим и предыдущим приближением меньше заданной погрешности e.
Внутри цикла вычисляется среднее значение a и b, и обновляется значение переменной c.
Проверяется условие F(a)*F(c)<=0. Если оно выполняется, то обновляется значение переменной b, иначе обновляется значение переменной a.
Выводится текущее приближение и его точность.
Значение переменной x вычисляется как среднее значение a и b.
Записывается точное значение x в файл output.txt.
Конец программы.