digamma function file

tinybike · Apr 19, 2013 · d6c1a98 · d6c1a98
1 parent 77fdea8
commit d6c1a98
Showing 1 changed file with 17 additions and 0 deletions.
diff --git a/digamma.m b/digamma.m
@@ -0,0 +1,17 @@
+function y = digamma(x,s)
+
+x_max = max(x);
+u = zeros(x_max,1);
+
+u(1) = 1;
+for jj = 2:x_max
+    u(jj) = 1/(1+s*(jj-1));
+end
+
+if min(x)
+    w = cumsum(u);
+    y = w(x);
+else
+    w = [0;cumsum(u)];
+    y = w(x+1);
+end