Implement the LU-decomposition of a real square matrix.

Implement also the corresponding functions to solve linear equations, calculate the determinant, and calculate the inverse matrix.

Is your LU faster than your QR?