The goal of this paper is to show that under some assumptions, for a d-dimensional fractional Brownian motion with Hurst parameter H > 1/2, the density of the solution of the stochastic differential equationUnder the framework of this present work, the Laplace method can be obtained in general hypoelliptic case and without imposing the structure equations on vector fields in Theorem 1.1. These two assumptions are imposed to obtain the correct Riemannian distance in the kernel expansion.Remark 1.3. When H > 1/2, to obtain a short-time asymptotic formula for the density of solution to equation (1.1) but with drift, one need to work on a version of Laplace method with fractional powers of ε, which will be very heavy and tedious in computation.Remark 1.4. When the present work was almost completed, we noticed that a proof for the Laplace method for stochastic differential equation driven by fractional Brownian motion with Hurst parameter 1/3 < H <