# HG changeset patch
# User jwe
# Date 1171613725 0
# Node ID bae85c1e0e2ad2f1dc5befef1550f210259d0060
# Parent  8c67f8be341d7deb03fb066532ac9b9e987dad2d
[project @ 2007-02-16 08:15:25 by jwe]

diff -r 8c67f8be341d -r bae85c1e0e2a scripts/specfun/nchoosek.m
--- a/scripts/specfun/nchoosek.m	Fri Feb 16 08:10:54 2007 +0000
+++ b/scripts/specfun/nchoosek.m	Fri Feb 16 08:15:25 2007 +0000
@@ -61,22 +61,23 @@
 function A = nchoosek (v, k)
 
   if (nargin == 2)
-    n = length (v);
-    if (n == 1)
-       A = round (exp (sum (log (k+1:v)) - sum (log (2:v-k))));
-    elseif (k == 0)
-      A = [];
-    elseif (k == 1)
-      A = v(:);
-    elseif (k == n)
-       A = v(:).';
-    else
-      m = round (exp (sum (log (k:n-1)) - sum (log (2:n-k))));
-      A = [v(1)*ones(m,1), nchoosek(v(2:n),k-1);
-	   nchoosek(v(2:n),k)];
-    endif
-  else
     print_usage ();
   endif
 
+  n = length (v);
+
+  if (n == 1)
+     A = round (exp (sum (log (k+1:v)) - sum (log (2:v-k))));
+  elseif (k == 0)
+    A = [];
+  elseif (k == 1)
+    A = v(:);
+  elseif (k == n)
+     A = v(:).';
+  else
+    m = round (exp (sum (log (k:n-1)) - sum (log (2:n-k))));
+    A = [v(1)*ones(m,1), nchoosek(v(2:n),k-1);
+	 nchoosek(v(2:n),k)];
+  endif
+
 endfunction