function S = e(f,x0,N)

S = x0; %Startvärde (den första approximationen)
x = x0;
h = 0.001;

for k=1:N
   df = (feval(f,x+h)-feval(f,x))./h; % Derivatan av f beräknas i x
    x = x - feval(f,x)./df;  % Ny approximation beräknas enligt Newton-Raphsons metod
    S = [S;x]; % Vi samlar upp alla approximationer
end

% Testa på följande exempel;
% fun = inline('x.^2-3.*x+2')
% x = e(fun,1.4,5)
