Traditional temporal relations propagating is based on Allen's Interval Algebra. This paper proposes an alternative method to propagate temporal relations among intervals, in which 5 5 matrices are used to represent temporal relations of intervals. Hence, the propagation of temporal relations is transformed into a numerical computation. For efficiency, we use the special values of the thirteen matrices to determine the possible temporal relations between two given intervals by using only the final resultant matrix so as to optimize the propagation. To evaluate the utility of the proposed technique, we have implemented the matrix representation in Java. The experimental results demonstrate that the approach is efficient and promising.