function [a,b] = iterpi(err)
%   Berechnung eines einschliessenden Intervalls [a,b] fuer die 
%   Kreiszahl pi mit abs(b-a) < err

a = 3/2*sqrt(3);
b = 2*sqrt(3);
while abs(b-a) >= err
    a = sqrt(a*b);
    b = 2*a*b/(a+b);
end