Fix error check in Elliptic3 function

cpmech · Apr 3, 2024 · 683911d · 683911d
1 parent 1abd5bf
commit 683911d
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/fun/ellipticintegral.go b/fun/ellipticintegral.go
@@ -121,7 +121,7 @@ func Elliptic3(n, φ, k float64) float64 {
 	if math.Abs(k*s-1.0) < 1e-15 {
 		return math.Inf(1)
 	}
-	if math.Abs(n*s-1.0) < 1e-15 {
+	if math.Abs(n*s*s-1.0) < 1e-15 {
 		return math.Inf(1)
 	}
 	t := N * s * s