%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
x=-0.5; %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
vastus1=x
sammud1=i 

%teine lahend
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
x=5.6; %alglähend jooniselt
i=1; %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