Suppose we have matrix of indices A=[1,3;2,1] and array/cell array/struct/table of values

I have attempted to implement the algorithm given in my ‘Numerical Analysis’ textbook (Burden and Faire, 9th edition, page 311.)