Methode newton rapshon

%metode Newton Raphson
%persaman yang diketahui f(x) = exp(x) - 4x
clc
clear all
x0 = 0; tol = 0.00001;
fprintf('tol = %10.8f, x0 = %8.6f\n', tol, x0);
%f = exp(x) - 4*x;
%f1 = exp(x) - 4;Turunan dari f
[f, f1] = fs(x0); % memanggil function fs.m
fprintf('f(x0) = %10.8f, f1(x0)= %10.8f,   ',f,f1);
i=1;x = x0 - f/f1;
fprintf('maka x%g = %10.8f\n',i,x);
while abs((x - x0)/x) > tol
   x0 = x;
   [f, f1] = fs(x0); % memanggil function fs.m
   fprintf('%3g %10.8f %10.8f %10.8f\n',i,f,f1,x);
   i=i+1;
   x = x0 - f/f1;
   %fprintf('maka x%g = %10.8f\n',i,x);
end;
fprintf('|(x - x0)/x| = %10.8f <= tol = %10.8f\n', abs((x - x0)/x),tol);
fprintf('Akarnya = %8.6f, banyak iterasi = %g \n',x,i);

untuk fungsinya

function [f,f1]=fs(x)
f=exp(x)-4*x;
f1=exp(x)-4;