Решение нелинейных уравнений на Фортране. Метод половинного деления.

Постановка задачи: Уточнить корни методом половинного деления с точностью 10^-6

program lab_one

real a,b,eps

integer N

eps=0.000001

a= !значение а

b= !значение b          

N=0

do while (abs(a-b)>eps.or.N<1000)

c=(a+b)/2

    if (f(a)*f(c)<0)then 

        b=c

    else

        a=c

    end if

N=N+1

end do

c=(a+b)/2

print*,c,N !вывод корня и количество делений

pause

end

real function f(z)

    real z

 f= !Ваша функция

end function f