--- a/scripts/specfun/betaincinv.m	Sat Aug 21 16:20:40 2021 +0200
+++ b/scripts/specfun/betaincinv.m	Sat Aug 21 21:54:34 2021 +0200
@@ -272,15 +272,15 @@
 
 ## Extreme values for y, a, b that really test the algorithm
 %!assert (betaincinv ([0, 1], 1, 3), [0, 1])
-%!assert <*60528> (betaincinv (1e-6, 1, 3), 3.3333344444450657e-07, 2*eps)
-%!assert <*60528> (betaincinv (1-1e-6, 3, 1), 0.999999666666556, 2*eps)
+%!assert <*60528> (betaincinv (1e-6, 1, 3), 3.3333344444450617e-7, 2*eps)
+%!assert <*60528> (betaincinv (1-1e-6, 3, 1), 0.9999996666665555, 2*eps)
 %!assert (betainc (betaincinv (0.9, 1e-3, 1), 1e-3, 1), 0.9, 2*eps)
 %!assert (betainc (betaincinv (.01, 1, 1e-3), 1, 1e-3), .01, 6*eps)
 %!assert (betainc (betaincinv (0.5, 100, 1), 100, 1), 0.5, 8*eps)
 %!assert (betainc (betaincinv (0.5, 1, 100), 1, 100), 0.5, 22*eps)
 %!assert (betaincinv ([0, 1], 1, 3, "upper"), [1, 0])
 %!assert <*60528> (betaincinv (1e-6, 1, 3, "upper"), 0.99, 2*eps)
-%!assert <*60528> (betaincinv (1-1e-6, 3, 1,"upper"), .01000000000003977, 2*eps)
+%!assert <*60528> (betaincinv (1-1e-6, 3, 1,"upper"), .01, 250*eps)
 %!assert (betainc (betaincinv (0.1, 1e-3, 1, "upper"), 1e-3, 1, "upper"),
 %!        0.1, 2*eps)
 %!assert (betainc (betaincinv (.99, 1, 1e-3, "upper"), 1, 1e-3, "upper"),