% Session p. 179

startintlab, format long e

f = inline('theta(exp(-pi^2*t),k)-pi/2/sqrt(5)','t','k'); kt = [];
for k=1:4, kt = [kt; k IntervalNewton(f,infsup(0.4,0.5),k)]; end
kt