function newtonraphson

  clc;

  clear;

  disp(' Metode Newton Raphson');

  disp('=============================');

 

  x0=input('Masukkan Nilai X awal :');

  E=input('Masukkan Nilai Toleransi :');

  i=0;

    M=input('Masukkan Nilai Iterasi Maksimum :');

 

  xb=0;

  disp('_______________________________________________');

  disp(' i           xi          f(xi)     epsilon');

  disp('_______________________________________________');

  while (E<M)

    fx=x0^3-x0^2+4*x0-1;

    gx=3*x0^2-2*x0+4;

    xb=x0-(fx/gx);

    M= abs(x0-xb);

    x0=xb;

    i=i+1;

    fprintf('%3.0f %12.6f %12.6f %12.6f\n',i,xb,fx,M);

  end

  disp('------------------------------')

  fprintf('Akarnya Adalah = %10.8f\n',xb);

end