clear, clc

fh = @(x) x^2-2;
a0 = [0 -2 2];
b0 = [2  0 0];
prec = 0.1;

for i=1:length(a0)
    [x pit] = puleni(fh,a0(i),b0(i))
end