In this work, equilibrium molecular dynamics (MD) simulations in a microcanonical ensemble are performed to evaluate the friction coefficient of a Brownian particle (BP) in a Lennard-Jones (LJ) solvent. The friction coefficients are determined from the time dependent friction coefficients and the momentum autocorrelation functions of the BP with its infinite mass at various ratios of LJ size parameters of the BP and solvent, ${sigma}_B/{sigma}_s$. The determination of the friction coefficients from the decay rates of the momentum autocorrelation functions and from the slopes of the time dependent friction coefficients is difficult due to the fast decay rates of the correlation functions in the momentum-conserved MD simulation and due to the scaling of the slope as 1/N (N: the number of the solvent particle), respectively. On the other hand, the friction coefficient can be determined correctly from the time dependent friction coefficient by measuring the extrapolation of its long time decay to t=0 and also from the decay rate of the momentum autocorrelation function, which is obtained by time integration of the time dependent friction coefficient. It is found that while the friction coefficient increases quadratically with the ratio of ${sigma}_B/{sigma}_s$ for all ${sigma}_B$, for a given ${sigma}_s$ the friction coefficient increases linearly with ${sigma}_B$.