a=2*ones(5,1);
b=-ones(5,1); b(1)=0;
c=-ones(5,1); c(5)=0;
f=[1 0 0 0 1]';
a,b,c,f

disp(' ')
disp('push return to continue')
pause;

[alpha,beta]=trilu(a,b,c);
alpha,beta

disp(' ')
disp('push return to continue')
pause;

disp(' ');
disp('     alpha     beta')
disp('    ----------------')
disp([alpha beta])


x=trilu_solve(alpha,beta,c,f);
x

disp(' ')
disp('push return to continue')
pause;

A=diag(2*ones(5,1))-(diag(ones(4,1),1)+diag(ones(4,1),-1));
A

disp(' ')
disp('push return to continue')
pause;

lu(A)
