A novel unified 2-D program with the boundary element method(BEM) for MZ-JTE is developed.The distributions of electric field and potential and the effect of-surface charge on high-voltage MZ-JTE are simulated.The BEM results are in agreement with ones obtained by the finite difference method (FDM) and the finite element method (FEM).The presentation of BEM and the successful development of the program have opened a new way for analyses and designs of junction terminations.