Octave / fun1

function y = fun1(f, a, b, e) #метод бисекции
  flag = false;
  if (e <= 0 || a > b) #проверка на правильность числа е и построения интервала
      flag = true;
      disp('Блок 1')
  else
      fa = f(a);
      fb = f(b);
      if (fa * fb < 0)
        while (true)
          fa = f(a);
          c = (a + b) / 2;
          fc = f(c);
          if (abs(fc) < e)
            break;
          else
            if (fa * fc < 0)
              b = c;
            else
              a = c;
            endif
          endif
        endwhile
      else
        flag = true;
      endif
  endif
  if (flag)
    disp(flag)
    y = 'Проверьте корректность данных';
  else
    y = c;
  endif
endfunction