# HG changeset patch # User Jordi GutiƩrrez Hermoso # Date 1316936155 18000 # Node ID 544304a09e42200b1a989537508c8a7438fb1e45 # Parent e81b93284605132535fe15e2df2747a6ef94373c Fix offbyones and typos in sprandsym.m diff --git a/scripts/sparse/sprandsym.m b/scripts/sparse/sprandsym.m --- a/scripts/sparse/sprandsym.m +++ b/scripts/sparse/sprandsym.m @@ -66,7 +66,7 @@ offdiag = randperm (n*(n - 1)/2, m); ## Row index - i = lookup (cumsum (1:n), offdiag) + 2; + i = lookup (cumsum (0:n), offdiag - 1) + 1; ## Column index j = offdiag - (i - 1).*(i - 2)/2; @@ -125,7 +125,7 @@ d = [mod(k,2):2:min(n,k)-2]; m = (k - d)/2; q = (n - d).*(n - d - 1).*m ./ (d + 2)./(d + 1)./(a - m + 1); - c = [1 cumprod (q)]; + c = [1 cumprod(q)]; s = sum (c); p = c/s;