The aim of the present work is to develop a pair of Lommel matrix functions of two variables suggested by the Bessel matrix functions and some of their properties are studied to be special cases of our results.