asd
%cálculo de F(x0)
x=x0;
Fx=eval(funcao);
s=sprintf(' *  x( 0) * %10.6f * %10.6f *',x,Fx);
disp(s)

while (abs(Fx)>eps)
    i=i+1;
    x=eval(fnewton);
    Fx=eval(funcao);
    s=sprintf(' *  x(%2.0f) * %10.6f * %10.6f *',i,x,Fx);
    disp(s)
end
zero=x;