Abstract: In this paper the predictor-corrector difference schemes of three-dimensional Poisson. equation with Dirichlet or/and Neumann boundary value problems are studied. It contains a complete set of schemes of the difference nodes which arc in inner region or in boundary planes, boundary lines and boundary corners. The new finite difference schemes can attain four-order accuracy with the least computational error by optimizing the selected factor.