%*** general settings
clc;clear all;close all;

%*** initialize A and b
A = [ 2,1,5,3,6;...
     -2,3,1,4,1;...
      0,3,0,4,1;...
      9,4,3,0,2;...
      8,4,6,0,7];
    
b = [1;3;2;5;4];
%*** compute LR-decomposition
LR = lr(A);
%*** solve sle
x1 = solve(LR,b); 

%*** solve sle with \-operator
x2 = A\b;

[x1,x2]
