  The theory of two-dimensional linear quaternion-valued differential equations
(QDEs) was recently established (see Kou and Xia, SAPM). Some profound
differences between QDEs and ODEs were observed.
  Also, an algorithm to evaluate the fundamental matrix by employing the
eigenvalues and eigenvectors was presented in [Kou and Xia, SAPM]. However, the
fundamental matrix can be constructed providing that the eigenvalues are
simple. If the linear system has multiple eigenvalues, how to construct the
fundamental matrix? In particular, if the number of independent eigenvectors
might be less than the dimension of the system. That is, the numbers of the
eigenvectors is not enough to construct a fundamental matrix. How to find the
"missing solutions"? The main purpose of this paper is to answer this question.
Furthermore, Caley determinant for Quaternion-valued matrix was adopted to
proceed the theory of QDEs in [Kou and Xia, SAPM]. One big disadvantage of
Caley determinant is that it can be expanded along the different rows and
columns. This may lead to different results due to non-commutativity of the
quaternions. This approach is not convenient to be used. The novel definition
of determinant for Quaternion-valued matrix based on permutation is introduced
to analyze the theory. This newly definition of determinant has great advantage
compare to Calay determinant.
