A new easy method has been presented to calculate the variable intervals corresponding to the stable univariant curves and to discriminate the stabilities of invariant points. This method and the one reported previously constitute a simple and universal theory for the computer-plotting of the equilibrium phase diagrams of a multisystem sign function matrix (SFM) discrimination method. Its main steps are: determining the stable univariant scheme according to the derivative (or difference) of r G m ; grouping the univariant curves by comparisons of the mutual relations among them; determining the existing intervals of the variables for the stable curves by comparisons of coordinate values of the curves about the invariant point; determining the stabilities of invariant points by comparisons of relations between the common curves and the invariant points. This method is suitable for any kind of phase diagram of closed or open systems in a phase diagram "space" with either 2 or more than 2 dimensions. Keywords: multisystem, computer-plotting of a phase diagram, stability of invariant point, sign function, discrimination matrix.The auto-discrimination of the stabilities of invariant points, the stable univariant scheme about an invariant point, and the existing intervals of the variables for the stable univariant curves are the key and knotty problems in computer-plotting of the stable equilibrium phase diagrams of a multisystem. To solve these problems, some authors [1 3] presented an algorithm on the basis of the theorem of minimum Gibbs free energy. At the same time, they developed a software package for computer plotting of the p-T-X and p-T-a phase diagrams. Although this algorithm is general, it is theoretically complex and difficult to understand. Previously, Vielzeuf and Boivin [4] developed a simple-easy method to deal with the problems involved in the construction of the petrogenetic grids (p-T phase diagrams) by means of sign function matrix. Following this method, Cheng, et al. [5] proposed a new approach to determining the equilibrium p-T bundle of the (n+2)-phase assemblage with sign function matrixes VF and SF. Based on a new and more general Gibbs free energy equation, Yin [6] extended the applicable range of this method from qualitative description to quantitative plotting. To simplify the discrimination of the metastable curves, Yin [6] also gave a theoretical analysis of the correlation of the topological properties of the curves with the elements in VF, SF and the slopes of the curves. Yin and Han [7] put forward a method to judge the stabilities of invariant points in the light of Guo's theorem [8]