function [T,L,U,P,Tinv] = demo_symtridiag_2(n,initialization)
% set T to a sparse random tridiagonal matrix and spy the L, U, P, Tinv matrices
d = randn(n,1);
e = randn(n-1,1);
T = symtridiag(d,e,initialization); % 2 means set T to a sparse matrix
[L,U,P] = lu(T);
Tinv = inv(T);
