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

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

%Lahendame Newtoni meetodiga, vaatame x kui vektorit
x(1)=1; %alglähend jooniselt on vektori esimene element
i=1; %näitab vektori liikmete arvu
while abs(f(x(i)))>10^-6
    x(i+1)=x(i)-f(x(i))/ftuletis(x(i));
    i=i+1;
end
vektor=x
lahendNM=x(i-1)
sammudNM=i-1

%Lahendame lõikajate meetodiga
y(1)=x(1); %alglähend jooniselt
y(2)=x(2); %NM leitud lähend
i=2; %vektoris y on juba 2 elementi
while abs(f(y(i)))>10^-6
    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 %alguses oli 2 elementi olemas
vektorLM=y

