def D(n, j): pi = pi_float() a = sqrt(Float(2)/3) * pi / j b = Float(n) - Float(1)/24 c = sqrt(b) expa = exp(a*c) iexpa = Float(1)/expa ch = (expa + iexpa)*0.5 sh = (expa - iexpa)*0.5 return sqrt(j) / (2*sqrt(2)*b*pi) * (a*ch-sh/c)
def test_import(): assert Float(3) + Float(6) == 9 assert isinstance(exp(3), Float) assert isinstance(log(3), Float)