Bounded-input bounded-output stability issues for fractional-order linear time invariant (LTI) system with double noncommensurate orders for the matrix case have been established in this paper. Sufficient and necessary condition of stability is given, and a simple algorithm to test the stability for this kind of fractional-order systems is presented. Based on the numerical inverse Laplace transform technique, time-domain responses for fractional-order system with double noncommensurate orders are shown in numerical examples to illustrate the proposed results.MSC 2010 : Primary 26A33; Secondary 33E12, 34A08, 34K37, 35R11, 60G22Key Words and Phrases: fractional-order systems, double noncommensurate orders, bounded-input bounded-output stability, numerical inverse Laplace transform technique