The concept of D-stability is significant for matrices of any order, especially when they appear in ordinary differential systems modeling physical problems. The problem of the characterization of a D-stable matrix was solved for low order matrices only, because of the computational complexity required to check the conditions assuring D-stability. Here we present an approach based on numerical linear algebra theorems, which provide conditions easily checked by computer algebra.