Sqrt[b^2+1] EllipticE[b^2/(b^2+1)]=pi