  In this work, we consider the Dirichlet boundary value problem for nonlinear
triharmonic equation. Due to the reduction of the nonlinear boundary value
problem to operator equation for the nonlinear term and the unknown second
normal derivative we design an iterative method at both continuous and discrete
level for numerical solution of the problem. Some examples demonstrate that the
numerical method is of fourth order convergence.
