% Harjutustund 4, ülesanne 4
y=@(x) x^3-6*x^2+3*atan(x)+3;
fplot(y,[-5,6])
yline(0)
grid on
% alglähenditeks võime võtta -0.5, 1 ja 5.8
% 1) lahendame võrrandi x^3-6x^2+3atanx+3=0 Newtoni meetodiga
ytuletis=@(x) 3*x^2-12*x+3/(1+x^2);

x=-0.5;
sammud=0;
while abs(y(x))>10^(-5)
  x=x-y(x)/ytuletis(x) % eeskiri
  sammud=sammud+1;
end
lahend1newton=x
sammud1newton=sammud

x=1;
sammud=0;
while abs(y(x))>10^(-5)
  x=x-y(x)/ytuletis(x) % eeskiri
  sammud=sammud+1;
end
lahend2newton=x
sammud2newton=sammud

x=5.8;
sammud=0;
while abs(y(x))>10^(-5)
  x=x-y(x)/ytuletis(x) % eeskiri
  sammud=sammud+1;
end
lahend3newton=x
sammud3newton=sammud

% 2) lahendame modifitseeritud Newtoni meetodiga
x=-0.5;
sammud=0;
while abs(y(x))>10^(-5)
  x=x-y(x)/ytuletis(-0.5) % eeskiri
  sammud=sammud+1;
end
lahend1newtonmod=x
sammud1newtonmod=sammud

x=1;
sammud=0;
while abs(y(x))>10^(-5)
  x=x-y(x)/ytuletis(1) % eeskiri
  sammud=sammud+1;
end
lahend2newtonmod=x
sammud2newtonmod=sammud

x=5.8;
sammud=0;
while abs(y(x))>10^(-5)
  x=x-y(x)/ytuletis(5.8) % eeskiri
  sammud=sammud+1;
end
lahend3newtonmod=x
sammud3newtonmod=sammud
[lahend1newton,lahend2newton,lahend3newton]
[sammud1newton,sammud2newton,sammud3newton]
[sammud1newtonmod,sammud2newtonmod,sammud3newtonmod]