New matrix module
New test program: test_matrix