%script solver
A = round(10*rand(8));
b = round(10*rand(8,1));
[L,D,U,count1] = nopivot(A);
[y,count2] = Lower(L,b);
[z,count3] = Diagonal(D,y);
[x,count4] = Upper(U,z);
x1 = A\b;
count1
count2
count3
count4
count1+count2+count3+count4
x-x1