/Applications/Julia-1.5.app/Contents/Resources/julia/bin/julia using LinearAlgebra A = [1 3 4 -4; 1 1 2 2; -1 -5 -6 10] lu(A) L,U,p = lu(A, check=false) L U p L*U == A A[[1,3,2],:] A = rand(5, 5) lu(A) A = rand(50, 50) @time lu(A) A = rand(100,100) @time lu(A); A = rand(1000,1000); @time lu(A); A = rand(2000,2000); @time lu(A); A = rand(5000,5000); @time lu(A); A = rand(10000,10000); @time lu(A); A = rand(20000,20000); @time lu(A); n = 5000 n = 10000 n = 20000 A = rand(n,n); b = rand(n,1); @time x = A\b; A*x == b norm(A*x-b,1) norm(A*x-b,2)