With the LU decomposition it is rather simple to solve a system of linear equations a11x1+a12x2+a13x3+a14x4=w1a21x1+a22x2+a23x3+a24x4=w2a31x1+a32x2+a33x3+a34x4=w3a41x1+a42x2+a43x3+a44x4=w4.
This can be written in matrix form as Ax=w.
where A and w are known and we have to solve for x. Using the LU dcomposition we write Ax≡LUx=w.