Tensors have a wide application in control theory, data mining, chemistry, information sciences, documents analysis and medical engineering. The material here is motivated by the development of the efficient numerical methods for solving the coupled tensor equations{ A1 ∗M X ∗N B1 + C1 ∗M Y ∗N D1 = E1, A2 ∗M X ∗N B2 + C2 ∗M Y ∗N D2 = E2, with Einstein product. We propose the tensor form of the LSQR methods to find the solutions X and Y of the coupled tensor equations. Finally we give some numerical examples to illustrate that our proposed methods are able to accurately and efficiently find the solutions of tensor equations with Einstein product