%Joonis alglähendite leidmiseks
fplot(@(x)x.^3-6*x.^2+3*atan(x)+3,[-1,7])
grid on
legend()
title('Joonis alglähendite leidmiseks')

%Jooniselt on näha 3 nullkohta, valime need alglähenditeks.
%Meile sobivad -0.5 ja 1 ja 5.6.

%Defineerime võrrandi vasaku poole funktsioonina
f=@(x)x^3-6*x^2+3*atan(x)+3;
%ja selle tuletise
ftuletis=@(x)3*x^2-12*x+3/(1+x^2);

%Lahendame võrrandi Newtoni meetodiga
%kasutame x-i vektorina, alglähend on vektori esimene element
x(1)=-0.5; %alglähend jooniselt, vektori esimene element
i=1; %vektori elementide arv
while abs(f(x(i)))>10^-6
    x(i+1)=x(i)-f(x(i))/ftuletis(x(i)); %Newtoni meetodi eeskiri
    i=i+1; %vektori elementide arv suurenes
end
vastus1=x(i-1)
sammud1=i-1 %alustasime loendust 1
vektor1=x

%teine lahend
%kasutame x-i vektorina, alglähend on vektori esimene element
x=1; %alglähend jooniselt
i=0; %vektori elementide arv arv
while abs(f(x))>10^-6
    x=x-f(x)/ftuletis(x); %Newtoni meetodi eeskiri
    i=i+1; %vektori elementide arv suurenes
end
vastus2=x
sammud2=i

%kolmas lahend
%kasutame x-i vektorina, alglähend on vektori esimene element
x=5.6; %alglähend jooniselt
i=0; %vektori elementide arv arv
while abs(f(x))>10^-6
    x=x-f(x)/ftuletis(x); %Newtoni meetodi eeskiri
    i=i+1; %vektori elementide arv suurenes
end
vastus3=x
sammud3=i