In this work a method for two-dimensional dynamic analysis of a cylindrical buoy under surfaced waves and moored to the seafloor by a flexible riser is presented. The finite element model used to describe the dynamics of the maritime line was developed by Lustosa (2002). This is a two node beam element based on the Euler-Bernoulli theory, with the element kinematics referred to a co-rotacional coordinate system attached to the element. The model considers the representation of geometric non-linearities due to the line motion large displacements resulting in axial-bending couplings in the element total stiffness representation. The element displacements in the line are due to the following external loadings: self weight, buoyancy force, hydrodynamic forces due to maritime currents, surface buoy imposed displacements, the forces of inertia, and, if present, the action of floaters. The buoy two displacements and one rotation degrees-of-freedom are measured at its center of gravity and are obtained from the solution of the non-linear algebraic system resulting from the Newmark temporal integration of the differential equations of motion, which solution convergence is granted by a Newton- Rapshon iterative method. The linear wave theory - Airy Theory -, associated to the Wheeler empirical modification, is used in this work to represent the monochromatic bi-dimensional surface waves. To obtain the wave forces acting on the surface buoy while in motion, the numerical integration of the equation presented by Chitrapu et al. (1998) is employed. This is the result of a two part equation: the first is evaluated as the Morison viscous drag force and, in a second part, is the sum of the - Froude-Krylov - and the Morison inertia forces. The numerical integration of the resulting equation is evaluated over the instantaneous wet surface of the buoy, in each time-step, using the adaptive recursive Simpson´s method. Finally the temporal integration of the buoy equations of motion are evaluated using a fifth-order Runge-Kutta method, with an adaptive control of the temporal time step. This procedure has been implemented in a computer program and the results obtained from some numerical analysis are compared to the ones obtained from independent simulations, so as to verify the ability of this proposed method to describe the dynamic of floating buoys.