Метод деления отрезка пополам

[Витяака]

Проблемма такая, что есть программа, но мне её надо немного доработать как - не понимаю( программа взята с вашего форума, и отредактирована на отетах),можно ли как то её исправить и сделать уравнение "у" поменьше?(подставлял "sqr(x)-2*x+1" что то не работало. заранее благодарен.

program rabot; 

{$APPTYPE CONSOLE} 

uses SysUtils; 

Var a,b,e,x:real; 

function y(x:real): real; 
begin 
y:=-cos(exp(0.49*ln(x))+sqrt(30/7))+exp(1/5*ln(x))/x-x; 
end; 

Begin 
write('a='); readln(a); 
write('b='); readln(b); 
write('E='); readln(E); 

if (a>b) or (y(a)*y(b)>0) then write ('HET KOPHEU') {кстати в этой строчке может быть ошибка в этой части y(a)*y(b)>0} 
else 
begin 
while abs(a-b)>=E do 
begin 
x:=(a+b)/2; 
if y(a)*y(x)<0 then b:=x {и здесь тоже, кстати по той же самой причине нужна отдельная переменная для умножения, а после переменную сравнивать с 0} 
else a:=x; 
end; 
writeln('x=', x); 
end; 

readln; 
end.

[ArtGrek]

отредактироват можно по разному, что твоя програма должна делат?

[Витяака]

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