As a continuation to the work reported by Yin in 1992, a new procedure is presented for computer plotting of the stable equilibrium phase diagram of an n-component system composed of (n + k) stoichiometric phases (or fluid species) where 2≤k≤4. The main points of the procedure are: (i) using the technique of sequential-absence of phases (SAP) to determine the possible invariant and univariant assemblages in a given multisystem; (ii) using the matrix inverse technique to generate and balance the univariant reactions from the corresponding univariant assemblages; (iii) comparing the phase assemblage at each invariant point with that of each univariant reaction to select the univariant curves about the corresponding invariant point; (iv) locating the invariant points with the technique of finding common equilibrium relation (CER); (v) using the sign function matrix (SFM) technique to discriminate between the stable portion of a univariant curve and its metastable extension about the corresponding invaria As a continuation to the work reported by Yin in 1992, a new procedure is presented for computer plotting of the stable equilibrium phase diagram of an n-component system composed of (n + k) stoichiometric phases (or fluid species) where 2 ≦ k ≤4. The main points of the procedure are: (i) using the technique of sequential-absence of phases (SAP) to determine the possible invariant and univariant assemblages in a given multisystem; (ii) using the matrix inverse technique to generate and balance the univariant reactions from the corresponding univariant assemblages; (iii) comparing to the phase assemblage at each invariant point with that of each univariant reaction to select the univariant curves about the corresponding invariant point; (iv) locating the invariant points with the technique of finding common equilibrium relation (CER); (v) using the sign function matrix (SFM) technique to discriminate between the stable portion of a univariant curve and its metastable extension about the c orresponding invaria