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

%Exercise 3
%Sekantmetod
%Prohl p.81
clear all; clc;
home
f= inline('x-4*sin(2*x)-(243/80)');

disp('>> Ange den första och andra startgissningen genom att kicka i grafen.');

x=[-1:0.05:8];
y=feval(f, x);
plot(x,y); hold on;
title('Lab1 -Exercise Sekantmetod'); 
grid on;

x1 =ginput(1);
plot(x1(1), x1(2),'*');
x2 =ginput(1);
plot(x2(1), x2(2),'*');

x1 = x1(1);
x2 = x2(2);

f1=feval(f, x1);
h=1;
delta = 10^-6;
it =0;
%sekantmetoden
while abs(h)>abs(x1)*delta && it<50
    f2=feval(f, x2);
    h=((x2-x1)/(f2-f1))*f2;
    x1=x2;
    f1=f2;
    x2=x2-h;
    it=it+1;
end

if(it == 50)
    disp('>> Avslutat pga iterationsantal');
else
    disp(['>> Den beräknade roten är ', num2str(x2)]);
end

hold off;