report/pyp_clustering/acl09-short/code/crppred_geom.m


1
2
3
4
5
6
7
8
9
10
11
12

function output = crppred_geom(input,lengths,b)


output = zeros(length(input),1);

p0=(1/52).^lengths;
a=b*p0;
for i = 1:length(input)
  output(i) = a(i)*sum(1./((1:input(i))+a(i)-1));
end