Matlab - busqueda de raices

clear all; clc

syms Y(x)

Y=input('ingrese la función: ');

a=input('ingrese el límite inferior: ' );

b=input('ingrese el límite superior: ');

t=input('ingrese Tolerancia');

c=1;

n=0;

while c

    n=n+1;

    x=(a+b)/2;

    if eval(Y)>0

        a=x;

    else

        b=x;

    end

    c=abs(eval(Y))>t;

    D(n,1:4)=[a b x abs(eval(Y))];

end

uitable('Data',D,'ColumnName',{'a' 'b' 'valor asumido' 'error'});