By increasing the number of Gaussian integrating points,the calculation precision of conductivity matrix of compound elements was improved and the free surface oscillation during the iteration was overcome. Adopting this method, program MESeep was coded and cases were carried out. The results showed that the method and program are feasible.