The dynamic calculations of slender marine risers, such as Finite Element Method (FEM) or Modal Expansion Solution Method (MESM), are mainly for the slender structures with their both ends hinged to the surface and bottom. However, for the re-entry operation, risers held by vessels are in vertical free hanging state, so the displacement and velocity of lower joint would not be zero. For the model of free hanging flexible marine risers, the paper proposed a Finite Difference Approximation (FDA) method for its dynamic calculation. The riser is divided into a reasonable number of rigid discrete segments. And the dynamic model is established based on simple Euler-Bernoulli Beam Theory concerning tension, shear forces and bending moments at each node along the cylindrical structures, which is extendible for different boundary conditions. The governing equations with specific boundary conditions for riser’s free hanging state are simplified by Keller-box method and solved with Newton iteration algorithm for a stable dynamic solution. The calculation starts when the riser is vertical and still in calm water, and its behavior is obtained along time responding to the lateral forward motion at the top. The dynamic behavior in response to the lateral parametric excitation at the top is also proposed and discussed in this paper. The dynamic calculations of slender marine risers, such as Finite Element Method (FEM) or Modal Expansion Solution Method (MESM), are mainly for the slender structures with their both ends hinged to the surface and bottom. However, for the re-entry operation , risers held by vessels are in vertical free hanging state, so the displacement and velocity of lower joints would not be zero. of the model of free hanging flexible marine risers, the paper proposed a Finite Difference Approximation (FDA) method for its dynamic calculation The riser is divided into a reasonable number of rigid discrete segments. And the dynamic model is established based on simple Euler-Bernoulli Beam Theory concerning tension, shear forces and bending moments at each node along the cylindrical structures, which is extendible for different boundary conditions. The governing equations with specific boundary conditions for riser’s free hanging state are simplified by Keller-box method and solved with Newton iterati on algorithm for a stable dynamic solution. The calculation starts when the riser is vertical and still in calm water, and its behavior is obtained along time to to the lateral forward motion at the top. The dynamic behavior in response to the lateral parametric excitation at the top is also proposed and discussed in this paper.