--- a/scripts/statistics/distributions/laplace_inv.m
+++ b/scripts/statistics/distributions/laplace_inv.m
@@ -32,19 +32,16 @@
     usage ("laplace_inv (x)");
   endif
 
-  [r, c] = size (x);
-  s = r * c;
-  x = reshape (x, 1, s);
-  inv = (-Inf) * ones (1, s);
+  inv = (-Inf) * ones (size (x));
 
   k = find (isnan (x) | (x < 0) | (x > 1));
   if (any (k))
-    inv(k) = NaN * ones (1, length (k));
+    inv(k) = NaN;
   endif
 
   k = find (x == 1);
   if (any (k))
-    inv(k) = Inf * ones (1, length (k));
+    inv(k) = Inf;
   endif
 
   k = find ((x > 0) & (x < 1));
@@ -53,6 +50,4 @@
 	      - (x(k) > 1/2) .* log (2 * (1 - x(k))));
   endif
 
-  inv = reshape (inv, r, c);
-
-endfunction
\ No newline at end of file
+endfunction