catch case n <= 0

Expander · Sep 30, 2023 · 9fda0d0 · 9fda0d0
1 parent c5d6adf
commit 9fda0d0
Showing 1 changed file with 8 additions and 4 deletions.
diff --git a/src/Harmonic.jl b/src/Harmonic.jl
@@ -22,9 +22,13 @@ function harmonic_f64(n::Integer)
 end
 
 function harmonic_big(n::Integer)
-    sum = one(BigFloat)
-    for k in 2:n
-        sum += inv(big(k))
+    if n <= 0
+        throw(DomainError(n, "harmonic not implemented for n <= 0"))
+    else
+        sum = one(BigFloat)
+        for k in 2:n
+            sum += inv(big(k))
+        end
+        sum
     end
-    sum
 end