% 1. Первая фаза программы нахождения действительного корня
% трансендентного уравнения методом дихотомии 
eps=0.000001;
xb1=0;
xe1=10;
f1='x^2-4';
for n=1:64
  x=(xe1+xb1)/2;
  y=eval(f1);
  if abs(y)<=eps
    disp('Корень найден');
    x, n
    break;
  end;
  if y>0
    xe1=x;
  else
    xb1=x;
  end;
end;


% 2. Вторая фаза программы нахождения действительного корня
% трансендентного уравнения методом дихотомии
clear all; 
eps=0.000001; xb1=0; xe1=0; f='';
xb1=input('Введите начальную границу поиска корня xb : ');
xe1=input('Введите конечную границу поиска корня xe : ');
f=input('Введите левую часть уравнения f(x)=0. (Например : x^2-4) : ','s');
fplot(f,[xb1 xe1]);
disp('Начат поиск корня');
for n=1:128
  x=(xe1+xb1)/2;
  [y]=eval(f);
  if abs(y)<=eps
    s1='Корень найден x = ';
    s2=num2str(x);
    disp(strcat(s1,s2));
    break;
  end;
  if y>0
    xe1=x;
  else
    xb1=x;
  end;
end;
s1='Число шагов n = ';
s2=num2str(n);
disp(strcat(s1,s2));