Integrable ODEs with matrix variables. Bi-Hamiltonian approach