Then, we divide the entire region into many hexahedron elements, and interpolate conductivity and electric potential with a trilinear function in each element.