%Lab 3
%Written by Jeannette Spuhler, 2011-10-04

%Exercise 3
%Newtons metod
%Gerd p74
clear all; clc;
home
format long;
disp('>> Titta först på grafen för att få startgissningar för rötterna');
disp('>> Tryck mellanslag för att fortsätta');
x=[-1:0.05:8];
y=x-4*sin(2*x)-(243/80);

plot(x,y);
title('Lab1 -Exercise Newton metod'); 


grid on;
pause
x = input('>> Ange en startgissning: ');
tol = input('>> Ange en tolerans som stopkristerium (t.ex. 10^-3): ');
iter = input('>> Ange max antal iterationer: ');
dx = 1;
it = 0;
X = []; %för att spara alla iterationsteg
while (abs(dx/x)>tol && it<iter)
    f = x-4*sin(2*x)-(243/80);
    fder = 1-8*cos(2*x);
    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);
else
    disp('Ej konvergens');
end
