Метод половинного деления

[serhan]

Дана функция, дано а и b и точность e. Нужно решить её методом половинного деления.

program mdp;
        function f(x: real): real;
        begin
             f:=4*(sqr(x))-16*x-9;
        end;
var  a,b,e,c,x:real;
begin
     writeln ('a,b,e');
     read(a,b,e);
     c:=(a+b)/2;
     while (abs(b-a)<=2*e) do
     begin
          if (f(c)=0)  then writeln(c) else
          begin
               if (f(a)*f(c)<0) then
               b:=c
               else
               a:=c;
          end;
     c:=(a+b)/2;
     writeln ('x=',c)  ;
     readln;
     end;

end.

Что не так? Подскажите.

[AnKor94]

Сообщение от serhan

Дана функция, дано а и b и точность e. Нужно решить её методом половинного деления.

Кого решить? Функцию?? Как можно решить функцию?

Сообщение от serhan

Что не так? Подскажите.

А что-то не так? Компилируется вроде. Запускается, пишет непонятные буквы, чего-то спрашивает, считает, выдаёт. Что не так-то, действительно?

[9i.]

(abs(b-a)<=2*e)

здесь должно быть больше или равно, а не меньше или равно

[serhan]

Все равно чтото зацикливается...

[9i.]

в каком смысле зацикливается?

     writeln ('x=',c)  ;
                readln;

это и должно быть в цикле?
или вам нужен только конечный результат?

program mdp;
        function f(x: real): real;
        begin
           f:=4*(sqr(x))-16*x-9;
        end;
var  a,b,e,c,x:real;
begin
     writeln ('a,b,e');
     read(a,b,e);
     c:=(a+b)/2;
     while (abs(f(c))>e) do
     begin
          if f(c)=0  then writeln(c)
          else
               if f(c)>0 then
               	b:=c
               else
               	a:=c;

     c:=(a+b)/2;
     end;
     writeln ('x=',c:5:3);
     readln;
     readln
end.

так вроде работает

[serhan]

У меня вообще чтото так не работает. Ввожу а,в,е ентер и все - виснет

[_Engine_]

попробуйте посмотреть под отладчиком, может тогда все встанет на свои места?

[Alex11223]

А вводите как? Через пробел?
Я бы советовал избегать ввода с помощью readln(a,b,c,еще 10050 переменных), лучше по одной.

[9i.]

а какие данные вы вводите??

[serhan]

Дано система уравнений.
4*(х^2)-16*x-9=0
(e^x)-e^(-3*x)-4=0
(x^5)-x-1=0
Вот её нужно решить методом половинного деления.
Можете под сказать какой код сюда нужно?
И ещо такой вопрос, вот есть три графика функций, как узнать где а и где b?

Вот код, на котором я остановился.

program mdp;
function f(x: real): real;
begin
f:=4*x*x-16*x-9;
{f:=exp(x)-exp(-3*x)-4};
{f:=x*x*x*x*x-x-1};
end;
var
a, b, e, c, x: real;
begin
write('a=');
read(a);
write('b=');
read(b);
write ('e=');
read(e);
c:=(a+b)/2;
while abs(b-a)>e do
begin
if f(a)*f(c)<0 then
b:=c
else
a:=c;
c:=(a+b)/2;
end;
x:=(a+b)/2;
writeln ('x=',x:3:3,' f(x)=',f(x):4:4);
readln;
end.