% Sitzung S. 170

startintlab; format long e; intvalinit('DisplayInfsup')

E=inline('1/AGM(sqrt(1-4*(pEW+pNS)^2),sqrt(1-4*(pEW-pNS)^2))',...
    'pEW','pNS');
f=inline('E(sqrt((0.25+epsilon)*(0.25-epsilon)),0.25)-2',...
    'epsilon','E');
epsilon=IntervalBisection(f,infsup(0.06,0.07),1e-15,E)