
% From Higham&Higham, MATLAB Guide, Sec.15.2
A=gallery('wathen',8,8);
subplot(1,2,1), spy(A), subplot(1,2,2), spy(chol(A))
figure
r=symrcm(A);
subplot(1,2,1), spy(A(r,r)), subplot(1,2,2), spy(chol(A(r,r)))
figure
m=symmmd(A);
subplot(1,2,1), spy(A(m,m)), subplot(1,2,2), spy(chol(A(m,m)))


