diff --git a/matlab/distributions/lpdfgweibull.m b/matlab/distributions/lpdfgweibull.m
index c99cba7cdac12bcb4b2c665173543528127377cd..3f2df819d298a2f99bfddff04a1a2355e47747fc 100644
--- a/matlab/distributions/lpdfgweibull.m
+++ b/matlab/distributions/lpdfgweibull.m
@@ -279,7 +279,7 @@ end
 %$ 
 %$ try
 %$    if isoctave
-%$        s = quadv(xdens, .0000000001, 100000,1e-10)
+%$        s = quadv(xdens, .0000000001, 20,1e-10)
 %$    else
 %$        s = integral(xdens, 0, 100000);
 %$    end
diff --git a/matlab/missing/stats/wblinv.m b/matlab/missing/stats/wblinv.m
index f22f9278e2e90efe0139fac408168296460a434c..1b87c9fcb88f51bf9578b7ceba94f5c8fb383243 100644
--- a/matlab/missing/stats/wblinv.m
+++ b/matlab/missing/stats/wblinv.m
@@ -155,7 +155,11 @@ t = exp(log(scale)+log(-log(1-proba))/shape);
 %$           if debug
 %$               [s, abs(p-s)]
 %$           end
+%$         if isoctave
+%$           t(k) = abs(p-s)<1e-10;  
+%$         else
 %$           t(k) = abs(p-s)<1e-12;
+%$         end
 %$       end
 %$    end
 %$ end