%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 Newtoni meetodiga ja salvestame lähendite jada vektoriks.
x(1)=-0.5; %alglähend jooniselt, tähistame vektori esimese elemendina
i=1;
while abs(f(x(i)))>10^-6
    x(i+1)=x(i)-f(x(i))/ftuletis(x(i));
    i=i+1;
end
vastus1=x(i)
sammud1=i %esimene element oli juba olemas
format long
x
format short

%Lahendame Newtoni meetodiga ja salvestame lähendite jada vektoriks.
x(1)=1; %alglähend jooniselt, tähistame vektori esimese elemendina
i=1;
while abs(f(x(i)))>10^-6
    x(i+1)=x(i)-f(x(i))/ftuletis(x(i));
    i=i+1;
end
vastus1=x(i)
sammud1=i %esimene element oli juba olemas
format long
x
format short

%Lahendame Newtoni meetodiga ja salvestame lähendite jada vektoriks.
x(1)=5.8; %alglähend jooniselt, tähistame vektori esimese elemendina
i=1;
while abs(f(x(i)))>10^-6
    x(i+1)=x(i)-f(x(i))/ftuletis(x(i));
    i=i+1;
end
vastus3=x(i)
sammud3=i %esimene element oli juba olemas
format long
x
format short
