Diagonal matrix method for fast computing positive region