function ex05b()
% Lösung von Aufgabe 5b
% Ergebnis wird gleich geplottet
f = @(x) sin(1./x);
a = 0.107;
b = 0.28;  
epsRel = 1e-5;

x = -0.1:0.0001:0.3;
y = f(x);
set(0,"defaultfigurecolor","white");
plot(x, y, "b-", x, 0*x, "k-", [a b], f([a b]), "k*")
hold("on")

error = Inf;
while (error > epsRel)
  c = sekantenSchritt(f, a, b);
  a = b;
  b = c;
  fprintf("[%10.6f %10.6f]\n", a, b);
  plot(b, f(b), "r*")
  error = 0.5*abs(b-a)/abs(b);
end
hold("off")

F = getframe(gcf);
imwrite(F.cdata, "bild82.png");
