diff --git a/matlab/distributions/lpdfgweibull.m b/matlab/distributions/lpdfgweibull.m
index a046c1dd2928e5c7f89d47e5be796a577b50ef5e..8e882f5e17f7659bbec70703ec63fcc7d4bccb25 100644
--- a/matlab/distributions/lpdfgweibull.m
+++ b/matlab/distributions/lpdfgweibull.m
@@ -206,7 +206,7 @@ end
 %$ 
 %$ try
 %$    if isoctave
-%$        s = quadv(density, .0000000001, 100000);
+%$        s = quadv(density, .0000000001, 100000,1e-10);
 %$    else
 %$        s = integral(density, 0, 100000);
 %$    end
@@ -229,7 +229,7 @@ end
 %$ 
 %$ try
 %$    if isoctave
-%$        s = quadv(density, .0000000001, 100000);
+%$        s = quadv(density, .0000000001, 100000,1e-10);
 %$    else
 %$        s = integral(density, 0, 100000);
 %$    end
@@ -252,7 +252,7 @@ end
 %$ 
 %$ try
 %$    if isoctave
-%$        s = quadv(density, .0000000001, 100000);
+%$        s = quadv(density, .0000000001, 100000,1e-10)
 %$    else
 %$        s = integral(density, 0, 100000);
 %$    end
@@ -275,7 +275,7 @@ end
 %$ 
 %$ try
 %$    if isoctave
-%$        s = quadv(xdens, .0000000001, 100000);
+%$        s = quadv(xdens, .0000000001, 100000,1e-10)
 %$    else
 %$        s = integral(xdens, 0, 100000);
 %$    end
@@ -298,7 +298,7 @@ end
 %$ 
 %$ try
 %$    if isoctave
-%$        s = quadv(xdens, .0000000001, 100000);
+%$        s = quadv(xdens, .0000000001, 100000,1e-10)
 %$    else
 %$        s = integral(xdens, 0, 100000);
 %$    end
@@ -321,7 +321,7 @@ end
 %$ 
 %$ try
 %$    if isoctave
-%$        s = quadv(xdens, .0000000001, 100000);
+%$        s = quadv(xdens, .0000000001, 100000,1e-10)
 %$    else
 %$        s = integral(xdens, 0, 100000);
 %$    end
@@ -347,7 +347,7 @@ end
 %$    s = NaN(n, 1);
 %$    for i=1:n
 %$        if isoctave()
-%$            s(i) = quadv(density, .0000000001, .1*i);
+%$            s(i) = quadv(density, .0000000001, .1*i,1e-10)
 %$        else
 %$            s(i) = integral(density, 0, .1*i);
 %$        end
@@ -378,7 +378,7 @@ end
 %$    s = NaN(n, 1);
 %$    for i=1:n
 %$        if isoctave()
-%$            s(i) = quadv(density, .0000000001, .1*i);
+%$            s(i) = quadv(density, .0000000001, .1*i,1e-10)
 %$        else
 %$            s(i) = integral(density, 0, .1*i);
 %$        end
@@ -409,7 +409,7 @@ end
 %$    s = NaN(n, 1);
 %$    for i=1:n
 %$        if isoctave()
-%$            s(i) = quadv(density, .0000000001, .1*i);
+%$            s(i) = quadv(density, .0000000001, .1*i,1e-10)
 %$        else
 %$            s(i) = integral(density, 0, .1*i);
 %$        end
diff --git a/matlab/missing/stats/wblinv.m b/matlab/missing/stats/wblinv.m
index 93ac84ffaab20499bd6aff6b56bea68348675567..f22f9278e2e90efe0139fac408168296460a434c 100644
--- a/matlab/missing/stats/wblinv.m
+++ b/matlab/missing/stats/wblinv.m
@@ -148,7 +148,7 @@ t = exp(log(scale)+log(-log(1-proba))/shape);
 %$               [shape, scale, x(k-1)]
 %$           end
 %$           if isoctave
-%$               s = quadv(density, 0, x(k-1));
+%$               s = quadv(density, 0, x(k-1),1e-10);
 %$           else
 %$               s = integral(density, 0, x(k-1));
 %$           end