%Joonis alglähendite leidmiseks
fplot(@(x)[x.^3-8*x+2,0])
grid on
legend()
title('Joonis alglähendite leidmiseks')
%Jooniselt on näha, et funktsioonil on kolm nullkohta =>
%võrrandil on kolm lahendit. 
%Valime alglähenditeks -3 ja 0.2 ja 2.7.

%Kontrollime koonduvust.
%Kui g(x)=1/8*(x^3+2), siis g'(x)=3/8*x^2
x=-3; %alglähend jooniselt
tuletis=abs(3/8*x^2) %ei sobi, sest on suurem kui 1
x=0.2;
tuletis=abs(3/8*x^2) %SOBIB
x=2.7;
tuletis=abs(3/8*x^2) %ei sobi
%Kui g(x)=nthroot(8*x-2,3), siis g'(x)=8/(3*nthroot((8*x-2)^2,3)))
x=-3;
tuletis=8/(3*nthroot((8*x-2)^2,3)) %SOBIB
x=2.7;
tuletis=8/(3*nthroot((8*x-2)^2,3)) %SOBIB

%Lahendame HIM kasutades
f=@(x)x.^3-8*x+2; %võrrandi vasak pool funktsioonina
x=-3; %alglähend jooniselt
while abs(f(x))>10^-6 %10^-6 on meie valitud täpsus
    x=nthroot(8*x-2,3);
end
vastus1=x

x=0.2; %alglähend jooniselt
while abs(f(x))>10^-6 %10^-6 on meie valitud täpsus
    x=1/8*(x^3+2);
end
vastus2=x

x=2.7; %alglähend jooniselt
while abs(f(x))>10^-6 %10^-6 on meie valitud täpsus
    x=nthroot(8*x-2,3);
end
vastus3=x