% Eksempel paa bruk av Newtons metode
clear all; close all; clc
format long

f = @(x) x.*cos(x) + sin(x);        % funksjonen vi vil finne nullpunkt til
fderiv = @(x) 2*cos(x)-x*sin(x);    % eksakt derivert. 

xn=1.5;                             % startverdi
xverdier=xn;                        % vektor med iterasjonsverdier
nmaks=25; epsilon=0.0001;           % maks antall iterasjoner og feilmargin

n=0;
while abs(f(xn))>epsilon && n<=nmaks
    n = n+1;
    xn = xn - f(xn)/fderiv(xn);     % <---- Newton
    xverdier=[xverdier, xn];
end

% Skriv ut resultater: (I egen lokke for ikke aa forvanske lokken over)
fprintf('%s %15s %18s\n\n', 'Iterasjon', 'xn', 'f(xn)')
for i = 1:length(xverdier)
    fprintf('%5d %25.10f %15.4e\n', i-1, xverdier(i), f(xverdier(i)))
end