%possible solution to exercise 2.4.4 in Bradie p. 104

clear all; clc;
%plotta funktionen
x=[1:0.01:2];
y=x.^7-3;
plot(x,y);
title('x.^7-3'); 

x = 1; %startgissning
tol = 10^-5;%tolerans
iter = 100; % antal iterationer
dx = 1;
it = 0;
X = []; %för att spara alla iterationsteg
while (abs(dx/x)>tol && it<iter)
    f = x.^7-3; 
    fder = 7*x.^6;
    dx =  -f/fder;
    x= x+dx;
    X = [X; x dx];
    it = it+1;
end
if(it<iter)
    disp(['>> Man anropade fixpunkiterationen ', num2str(it),' gånger.']);
    disp(['>> Fixpunkten är ', num2str(x), '.']);
    disp('>> Iterationssteg:');
    disp('       x    x_n-x_n+1');
    disp(X);
    disp(['Given lösning:', num2str(3^(1/7))])
else
    disp('Ej konvergens');
end
