O que há de errado?

luccasrodge
(usa Ubuntu)

Enviado em 14/04/2015 - 20:59h

function [Raiz,Iter,Erro] = newtonR(funcao,x0,Toler,IterMax);

x=x0;

Fx = eval(funcao, x0);

DFx = eval(diff(x0));

Iter = 0;

disp('      Iter        x     DFx      Fx        DeltaX');

disp([Iter           x        DFx      Fx     ]);

while (1)

    DeltaX = -Fx/DFx;

    x = x+DeltaX;

    Fx = eval(funcao);

    DFx = eval(diff(x));

    Iter = Iter+1;

    disp([Iter     x     DFx      Fx     DeltaX]);      

    if(abs(DeltaX)<Toler && abs(Fx)<Toler)||abs(DFx)==0 || Iter>= IterMax

        break;

    end

end

Raiz = x;

if abs(Fx)<Toler

    Erro = 0;

else

    Erro =1;

end

end