Both experimental and estimated multicomponent tray efficiencies are known to become unbounded and exhibit bizarre variation from tray to tray in a distillation column even for ideal systems. Besides lacking the physical significance, they pose numerical difficulties in computation of the component tray efficiencies from the point efficiencies whenever the latter assume negative values (Lockett, M. J. (1986). Distillation tray fundamentals. Cambridge: Cambridge University Press). The use of the Stefan-Maxwell approach for estimation of point efficiency leads to a matrix form. Recently, the tray-efficiency matrix for well-mixed, dispersed and plug flow of liquid has been obtained in terms of point-efficiency matrix similar to binary systems. However, there are no method to use directly the tray-efficiency matrix in simulation. Here, we present a method for direct incorporation of the efficiency matrix in the well-known Naphtali-Sandholm method. We have studied the convergence characteristics of the modified and the original methods for distillation, absorption and extraction columns. The results show the modified one is superior to the original method.