diff --git a/test/NumberTypeTest.jl b/test/NumberTypeTest.jl
index 15c769ff4..71c6ad95e 100644
--- a/test/NumberTypeTest.jl
+++ b/test/NumberTypeTest.jl
@@ -55,21 +55,19 @@ include(joinpath(@__DIR__, "testutils.jl"))
     end
 
     @testset "BigFloat relative tolerance bug test Issue #557" begin
-		if false # disabled due to Julia issue https://github.com/JuliaLang/julia/issues/25853
-	        x = Fun(BigFloat(0)..BigFloat(20_000));
-	        ν = BigFloat(1568)
-	        f = x^(ν/2-1) * exp(-x/2)
-	        # TODO: JacobiWeight should support general types to avoid warning here
-	        @test f(10.0) ≈ 10.0^(ν/2-1) * exp(-10.0/2) rtol=eps()
-	        ex_mathematica = big"8.6494114955713766301430915207861966687081153e778"
-	        @test (cumsum(f)(10.0) - ex_mathematica)/ex_mathematica ≤ eps()
-
-	        x = Fun(BigFloat(0)..BigFloat(Inf))
-	        ν = BigFloat(1568)
-            f = x^(ν/2-1) * exp(-x/2)
-            @test f(10.0) ≈ 10.0^(ν/2-1) * exp(-10.0/2) rtol=eps()
-            @test_skip (cumsum(f)(10.0) - ex_mathematica)/ex_mathematica ≤ eps()
-        end
+        x = Fun(BigFloat(0)..BigFloat(20_000));
+        ν = BigFloat(1568)
+        f = x^(ν/2-1) * exp(-x/2)
+        # TODO: JacobiWeight should support general types to avoid warning here
+        @test f(10.0) ≈ 10.0^(ν/2-1) * exp(-10.0/2) rtol=eps()
+        ex_mathematica = big"8.6494114955713766301430915207861966687081153e778"
+        @test (cumsum(f)(10.0) - ex_mathematica)/ex_mathematica ≤ eps()
+
+        x = Fun(Ray(BigFloat(0)..BigFloat(Inf)))
+        ν = BigFloat(1568)
+        f = x^(ν/2-1) * exp(-x/2)
+        @test f(10.0) ≈ 10.0^(ν/2-1) * exp(-10.0/2) rtol=eps()
+        @test_skip (cumsum(f)(10.0) - ex_mathematica)/ex_mathematica ≤ eps()
     end
 
     @testset "#641" begin