n = 200;
a = 10*speye (n) + sparse (1:n,ceil ([1:n]/2),1,n,n) + ...
    sparse (ceil ([1:n]/2),1:n,1,n,n);
r1 = chol (a);
spy (r1);
axis ("ij");
