Newton-Raphson 법

function root=newton2(func1,func2,x0, threshold)

%func1: input function

%func2: derivative of input

%func2=diff('func1',1);

fprintf('| iteration  |     x     |    error  |\n')

j=0;

while(1)

    j=j+1;

    x=x0 - feval(func1, x0)/feval(func2, x0);

    rel_error=abs((x-x0)/x);

    x0=x;

    fprintf('|      %d     |% 10.5f | %9.5f |\n',j,x0,rel_error)

    if(rel_error < threshold)

        break

    end

end

root=x0;