In this paper a computational algorithm for nonlinear balanced realization and model reduction based on Taylor series expansion is proposed. This algorithm requires recursive computations with respect to the order of the Taylor series in which we need to solve linear equations with unknown parameters in each step. Furthermore, the proposed method is applied to a double pendulum system. Some numerical simulations demonstrate the effectiveness of the proposed algorithm.