% Session p. 151

session_p148

x = zeros(n,1);
perm = symrcm(A);        % reverse Cuthill-McKee ordering
R = chol(A(perm,perm));
x(perm) = R\(R'\b(perm));
x1 = x(1)