Пожалуйста, помогите проверить программу

[Марина357]

Помогите пожалуйста проверить правильно ли записана программа к заданию на языке Паскаль
Вычислить корни функционального уравнения.
Составить программу для решения уравнения f(x)=0 на отрезке [A,B] с точностью Е методом деления отрезка пополам. Метод состоит в последовательном приближении к корню за счет уменьшения отрезка, на котором находится корень. Каждое новое положение Х находится как середина текущего отрезка. Концы текущего отрезка выбираются из условия противоположности знака f(х) на его концах. Вычисления корня заканчивается, когда длина отрезка станет меньше Е.
Для проверки программы задать: А≈0, В≈2, Е≈10-4.

program lab17;
var C,A,B,E:real;
begin
writeln('Введите A,B,E');
readln(A,B,E);
while(B-A)>E do 
begin
C:=(A+B)/2;
if (1/(sqrt(5)+sin(0.1*c)+ln(1+c))-c)*(1/(sqrt(5)+sin(0.1*b)+ln(1+b))-b)>0
then b:=c
else A:=C;
end;
writeln((A+B)/2);
end.

Решение.

Введите A,B,E
0
2
0.0001
0.384613037109375

[ViktorR]

Опишите функцию:

function MyF(x: real): real;
begin
   MyF := 1/(sqrt(5)+sin(0.1*x)+ln(1+x))-x);
end;

Тогда оператор:

if (1/(sqrt(5)+sin(0.1*c)+ln(1+c))-c)*(1/(sqrt(5)+sin(0.1*b)+ln(1+b))-b)>0
then b:=c
else A:=C;

можно будет записать так:

if (MyF(c)*MyF(b) > 0) then
   b:=c
else
   a:=c;

И вот ещё. В этом коде предполагается, что на левом конце интервала значение функции меньше нуля, а на правом - больше.
В общем случае это следует проверять, и перед циклом выполнить соответствующую перестановку границ.


Как-то так, ...

[Марина357]

А не подскажете, как сделать эту перестановку границ?

[ViktorR]

Например, переменным a и b присвоены некоторые значения.

a := <что-то>;
b := <что-то>;
if (MyF(a) * MyF(b) > 0) then
   writeln("Уточните границы)
else
   if (MyF(a) > MyF(b)) then
    begin
      p := a;
      a:= b;
      b := p;
   end;

Как-то так, ...