clear all, close all

% solution of a system of linear equations using the "backslash" operator
N = 5;
A = rand(N)
b = (1:N)'

x = A \ b


% LU decomposition with and withour row permutations
disp('========================================')

[ L, U ] = lu(A)

disp('========================================')

[ L, U, P ] = lu(A)

disp('========================================')
