Modules

Matrix operations and LUP decomposition.