Matrix Liouville equation