Quaternion-valued differential equations (QDEs) is a new kind of differential equations. In this paper, an algorithm was presented for solving linear nonhomogeneous quaternionic-valued differential equations. The variation of constants formula was established for the nonhomogeneous quaternionic-valued differential equations. Moreover, several examples showed the feasibility of our algorithm. Finally, some open problems end this paper.