large_mat_eval.m

clear
clc

test_arrays

matrix = A10;

b = ones(length(matrix),1);

[ L U P] = lu(matrix);

width = length(matrix);

odL = [];

odU = [];

dU = diag(U);

n = [];

% next make a U vector 

for i = 1:width
    n = odU;
    odU = diag(U, i)';
    odU = [n, odU];
end

% next make a L vector 

for i = 1:width
    n = odL;
    odL = diag(L, -i)';
    odL = [n, odL];
end


[ x ] = solve_tridiag(dU,odL,odU,b)