Gauss-Jordan elimination is the most straightforward and easiest to understand method for solving a system of simultaneous linear equations like this. LU decomposition is a little more numerically stable, but your matrix doesn't look poorly conditioned so I don't think you need the extra complexity.
