Fix accept_solution returning symbolic instead of numeric

- Ensure accept_solution always returns a numeric value (Number type)
- Also substitute symbols from the solution, not just the input equation
- Fall back to Inf if result remains symbolic after substitution
- Format test/runtests.jl with SciMLStyle

Fixes TypeError in numeric comparisons when ε is symbolic.

🤖 Generated with [Claude Code](https://claude.com/claude-code)

Co-Authored-By: Claude Opus 4.5 <noreply@anthropic.com>

Author: claude

src/numeric_utils.jl

@@ -24,8 +24,25 @@ function accept_solution(eq, x, sol; plan = default_plan())
         # x₀ = test_point(plan.complex_plane, plan.radius)
         # Δ = substitute(diff(sol, x) - expr(eq), Dict(x => x₀))
         S = subs_symbols(eq, x; include_x = true, plan.radius)
+        # Also substitute any symbols in sol that may not be in eq
+        for v in get_variables(value(sol))
+            if !haskey(S, v) && !isequal(v, x)
+                S[v] = Complex(randn())
+            end
+        end
         Δ = substitute(diff(sol, x) - expr(eq), S)
-        return abs(Δ)
+        result = abs(Δ)
+        # Ensure result is numeric, not symbolic
+        if result isa Number
+            return result
+        else
+            # Try to extract numeric value
+            val = Symbolics.unwrap(result)
+            if val isa Number
+                return abs(val)
+            end
+            return Inf
+        end
     catch e
         #
     end

test/runtests.jl

@@ -333,7 +333,8 @@ end
 
     # Test that vector expressions throw an appropriate error
     @test_throws ErrorException("Vector expressions are not supported. Please use element-wise integration with `integrate.([expr1, expr2, ...], x)` instead.") integrate([x])
-    @test_throws ErrorException("Vector expressions are not supported. Please use element-wise integration with `integrate.([expr1, expr2, ...], x)` instead.") integrate([1, 2 * α], α)
+    @test_throws ErrorException("Vector expressions are not supported. Please use element-wise integration with `integrate.([expr1, expr2, ...], x)` instead.") integrate(
+        [1, 2 * α], α)
 
     # Test that scalar integration still works
     @test integrate(x) == ((1//2)*(x^2), 0, 0)