Metode Newton Raphson

function fsNewtonRaphson(f,a,n,tol)

% ----------- FUNGSI Newton Raphson ---------------

% variabel f :fungsi f(x)

%          a :nilai awal yang didefinisikan

%          n :banyaknaya iterasi

%          toll: nilai toleransi (epsiolon) yang membatasi iterasi

syms x

N=[];

x0 = a;

tur_f = diff(f,x);

for i=1:n,

    fx=subs(f,x,x0);

    tur_fx=subs(tur_f,x,x0);

    xi=x0-(fx/tur_fx);

    N=[N; i, x0, xi,fx];

    if abs(fx)<tol,

        break;

    end

    x0=xi;

end

N

% ----------------------------------------------------------------------

% Untuk memanggil di Command  Window

% Pilih nilai awal untuk mencari akar a=2, n=10 dan toll=10^(-5)

% syms x;

% f=x^2-1

% f =

% x^2 - 1

% fsNewtonRaphson(f,2,10,1e-5)

N =

    1.0000    2.0000    1.2500    3.0000

    2.0000    1.2500    1.0250    0.5625

    3.0000    1.0250    1.0003    0.0506

    4.0000    1.0003    1.0000    0.0006

    5.0000    1.0000    1.0000    0.0000