secent - 할선법

function root=secent(func, x1, x2, threshold)

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

j=0;

while(1)

    j=j+1;

    xn=x2-feval(func, x2)*(x2-x1)/(feval(func, x2)-feval(func,x1));

    if xn ~=0

        rel_error=abs((xn-x1)/xn)*100;

    end

    x1=x2;

    x2=xn;

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

    if(rel_error < threshold)

        break

    end

end

root=x2;