In this paper we establish the existence results and construct fourth order numerical method for solving fully third order nonlinear differential equation with integral boundary conditions. The method is based on the discretization of an iterative method on continuous level with the use of the trapezoidal quadrature formulas with corrections. Some examples demonstrate the applicability of the theoretical results of existence and uniqueness of solution and the convergence order 4 of the proposed numerical method. The approach used for the third order nonlinear differential equation with integral boundary conditions can be applied to differential equations of any orders.