# HG changeset patch
# User Rik <octave@nomad.inbox5.com>
# Date 1308283510 25200
# Node ID e7cc2d4a6db32713930224b83a5c3fb0455508db
# Parent  901d466ee55ac902a875ec0ade6f1eccef0841dc
Fix range of sigma in normal distribution to exclude 0.

* normcdf.m, normpdf.m: Correct 's >= 0' to 's > 0'.

diff -r 901d466ee55a -r e7cc2d4a6db3 scripts/statistics/distributions/normcdf.m
--- a/scripts/statistics/distributions/normcdf.m	Wed Jun 15 10:35:37 2011 -0400
+++ b/scripts/statistics/distributions/normcdf.m	Thu Jun 16 21:05:10 2011 -0700
@@ -50,18 +50,18 @@
   cdf = zeros (sz);
 
   if (isscalar (m) && isscalar(s))
-    if (find (isinf (m) | isnan (m) | !(s >= 0) | !(s < Inf)))
+    if (find (isinf (m) | isnan (m) | !(s > 0) | !(s < Inf)))
       cdf = NaN (sz);
     else
       cdf =  stdnormal_cdf ((x - m) ./ s);
     endif
   else
-    k = find (isinf (m) | isnan (m) | !(s >= 0) | !(s < Inf));
+    k = find (isinf (m) | isnan (m) | !(s > 0) | !(s < Inf));
     if (any (k))
       cdf(k) = NaN;
     endif
 
-    k = find (!isinf (m) & !isnan (m) & (s >= 0) & (s < Inf));
+    k = find (!isinf (m) & !isnan (m) & (s > 0) & (s < Inf));
     if (any (k))
       cdf(k) = stdnormal_cdf ((x(k) - m(k)) ./ s(k));
     endif
diff -r 901d466ee55a -r e7cc2d4a6db3 scripts/statistics/distributions/normpdf.m
--- a/scripts/statistics/distributions/normpdf.m	Wed Jun 15 10:35:37 2011 -0400
+++ b/scripts/statistics/distributions/normpdf.m	Thu Jun 16 21:05:10 2011 -0700
@@ -50,18 +50,18 @@
   pdf = zeros (sz);
 
   if (isscalar (m) && isscalar (s))
-    if (find (isinf (m) | isnan (m) | !(s >= 0) | !(s < Inf)))
+    if (find (isinf (m) | isnan (m) | !(s > 0) | !(s < Inf)))
       pdf = NaN (sz);
     else
       pdf = stdnormal_pdf ((x - m) ./ s) ./ s;
     endif
   else
-    k = find (isinf (m) | isnan (m) | !(s >= 0) | !(s < Inf));
+    k = find (isinf (m) | isnan (m) | !(s > 0) | !(s < Inf));
     if (any (k))
       pdf(k) = NaN;
     endif
 
-    k = find (!isinf (m) & !isnan (m) & (s >= 0) & (s < Inf));
+    k = find (!isinf (m) & !isnan (m) & (s > 0) & (s < Inf));
     if (any (k))
       pdf(k) = stdnormal_pdf ((x(k) - m(k)) ./ s(k)) ./ s(k);
     endif