JuliaでLU分解を利用して連立一次方程式を数値計算するとき,枢軸選択 (pivoting) に注意する必要があります.みるべき情報はここですが,備忘録も兼ねてまとめておきます.
julia> using LinearAlgebra julia> A = [0 0 2; 1 0 1 ; 2 -1 1] 3×3 Array{Int64,2}: 0 0 2 1 0 1 2 -1 1 julia> L, U = lu(A) # 素直にLU分解 LU{Float64,Array{Float64,2}} L factor: 3×3 Array{Float64,2}: 1.