clear %vahemälu puhtaks
%Joonis alglähendite leidmiseks

fplot(@(x)x.^5-x.^2-2,[-2,2])
grid on
legend()
title('Joonis alglähendite leidmiseks')
%Jooniselt on näha, et võrrandil on 1 lahend (üks nullkoht).
%Valime alglähendiks 1.1.

%Defineerime võrrandi vasaku poole funktsioonina
f=@(x)x^5-x^2-2;
%ja selle tuletise
ftuletis=@(x)5*x^4-2*x;

%Newtoni meetod
x(1)=1.1; 
%alglähend jooniselt on lähendite vektori esimene element
i=1; %vektoris on 1 element
while abs(f(x(i)))>10^-5
    x(i+1)=x(i)-f(x(i))/ftuletis(x(i));
    i=i+1;
end
lahendNM=x(i-1)
sammudNM=i-1

%Lahendame lõikajate meetodiga, tähistame lahendit y
y(1)=x(1); %alglähend, kasutame NM tähistusi
y(2)=x(2); %NM abil leitud esimene element
i=2; %vektoris y on juba 2 elementi
while abs(f(y(i)))>10^-5
    y(i+1)=y(i)-(y(i)-y(i-1))/(f(y(i))-f(y(i-1)))*f(y(i));
    i=i+1;
end
lahendLM=y(i-1)
sammudLM=i-2 %2 esimest elementi olid juba defineeritud