The Riccati equation method is used for study the oscillatory and non oscillatory behavior of solutions of systems of two first order linear two by two dimensional matrix differential equations. An integral and an interval oscillatory criteria are obtained. Two non oscillatory criteria are obtained as well. On an example one of the obtained oscillatory criteria is compared with some well known results.