fplot(@(x)x.^3-6*x.^2+3*atan(x)+3,[-2,6])
grid on

%jooniselt on näha, et võrrandil on 3 lahendit
%valime alglähenditeks -0,5 ja 1 ja 6

%kasutame Newtoni meetodit
x=-0.5; %alglähend
while abs(x^3-6*x^2+3*atan(x)+3)>10^-6
    x=(2*x^3+2*x^5-9*x^2-6*x^4+3*x-3*atan(x)-3*x^2*atan(x)-3)/(3*x^2+3*x^4-12*x-12*x^3+3);
end
vastus1=x

x=1; %alglähend
while abs(x^3-6*x^2+3*atan(x)+3)>10^-6
    x=(2*x^3+2*x^5-9*x^2-6*x^4+3*x-3*atan(x)-3*x^2*atan(x)-3)/(3*x^2+3*x^4-12*x-12*x^3+3);
end
vastus2=x

x=6; %alglähend
while abs(x^3-6*x^2+3*atan(x)+3)>10^-6
    x=(2*x^3+2*x^5-9*x^2-6*x^4+3*x-3*atan(x)-3*x^2*atan(x)-3)/(3*x^2+3*x^4-12*x-12*x^3+3);
end
vastus3=x