[Pascal] repeat-until и for-do

[Richi106]

Задание:
Установить, на каком из 10 интервалов отрезка [1, 4]  функция f(x)=ln(x)-x+1,8  имеет корень.
Разделив отрезок на 10 интервалов, найти значение функции на интервале путем сравнения правой и левой координаты интервала.
есть наработки, получается неправильное значение:

Код: [Выделить]


program prim4;
 var
  f,h,x,a,b,g,xl,xr:real;
 begin
   writeln('введите границы отрезка');
   readln(a,b);
   h:=(b-a)/10;
   xl:=a;
   xr:=a+h;
    repeat
      f:=ln(xl)-xl+1.8;
      g:=ln(xr)-xr+1.8;
      xl:=xl+h;
      xr:=xr+h;
      if (f*g)<0 then writeln('функция имеет значение на [',xl:4:2,';',xr:4:2,']');
    until xr=b;
  end.

просьба найти ошибку. вроде правильное значение:2,845

[Данила]

с точки зрения кода ошибок не вижу. единственное - не проверяете вводимые границы, вдруг кто то введет a больше b... будет печально.
с математической точки зрения уверены, что реализация такая?

[testtest]

с точки зрения кода ошибку вижу

Код: [Выделить]

      f:=ln(xl)-xl+1.8;
      g:=ln(xr)-xr+1.8;
      xl:=xl+h;
      xr:=xr+h;
-->      if (f*g)<0 then writeln('функция имеет значение на [',xl:4:2,';',xr:4:2,']');
здесь xl, xr принимают уже следующие значения. нужно так:

Код: [Выделить]

      f:=ln(xl)-xl+1.8;
      g:=ln(xr)-xr+1.8;
      if (f*g)<0 then writeln('функция имеет значение на [',xl:4:2,';',xr:4:2,']');
      xl:=xl+h;
      xr:=xr+h;
ну и условие выхода

Код: [Выделить]

until (xr > b + h);

а во-вторых, логарифм при x < 1 быстро убывает, а в окрестности x = 0 вообще бешено несется вниз, поэтому a и b нужно задавать или 0-1, или 1-10, т.к. корня на отрезке, который ты вероятнее всего задаешь, два.

[Richi106]

с точки зрения кода ошибок не вижу. единственное - не проверяете вводимые границы, вдруг кто то введет a больше b... будет печально.
с математической точки зрения уверены, что реализация такая?

ну да.

с точки зрения кода ошибку вижу

Код: [Выделить]

      f:=ln(xl)-xl+1.8;
      g:=ln(xr)-xr+1.8;
      xl:=xl+h;
      xr:=xr+h;
-->      if (f*g)<0 then writeln('функция имеет значение на [',xl:4:2,';',xr:4:2,']');
здесь xl, xr принимают уже следующие значения. нужно так:

Код: [Выделить]

      f:=ln(xl)-xl+1.8;
      g:=ln(xr)-xr+1.8;
      if (f*g)<0 then writeln('функция имеет значение на [',xl:4:2,';',xr:4:2,']');
      xl:=xl+h;
      xr:=xr+h;
ну и условие выхода

Код: [Выделить]

until (xr > b + h);

а во-вторых, логарифм при x < 1 быстро убывает, а в окрестности x = 0 вообще бешено несется вниз, поэтому a и b нужно задавать или 0-1, или 1-10, т.к. корня на отрезке, который ты вероятнее всего задаешь, два.

значение одно в том то и дело.

[Richi106]

с перестановкой "if" значение совпадает с искомым. спасибо за помощь.