We developed a MatLab c code that implements the multicomponent diffusion approaches outlined by Morral and Thompson [1,2] and by Glicksman and Lupulescu [3]. Zero flux planes (ZFP's), first introduced by Dayananda and Kim [4], occur commonly in ternary and higher-order multicomponent alloys. ZFP's for a component manifest as stationary or symmetrically moving planes disposed about the Matano interface. In contrast to the normal global mixing that occurs in a binary diffusion couple, the presence of a stationary ZFP requires that mixing of the blocked component occurs unilaterally through the release of pairs of coupled depletion-repletion waves that increase the average concentration in the component-poor alloy, but decrease it in the component-rich alloy [3]. We previously analyzed the kinetics of multicomponent diffusion near ZFP's for single-phase couples based on the behavior of diffusion couples between the fixed end-member alloy Cr-10 at.%, Al-10 at.%, Ni-80 at.%, with a series of end-member alloys having a concentration difference vector with a magnitude of 1 at.%. It was demonstrated that the spreading rates of the ternary diffusion zones could be reduced by choosing the orientation of the composition vectors in composition space to be near the Euler angle for the minor eigenvector of the diffusivity matrix. Moreover, a critical Euler angle, ψ * ≈ 29 • , was found in ternary Cr-Al-Ni alloys, near which diffusion couples selected with their composition vectors within ±5 • of ψ * exhibit minimal rates of interdiffusion [3]. The reduced atomic transport predicted in these alloys was associated with the proximity of their Euler angles to that of the ZFP for the major elemental component. Similar analyses performed on ternary couples with the fixed end member alloy Ni-43.5 at.%, Zn-25 at.%, Cu-31.5 at.% and unit composition vectors, show that diffusive spreading is reduced for couples located in composition space close to the ZFP of the minor component.