In this paper we propose a logdensity estimation based on a Fourier expansion. The basis functions consisting of trigonometric functions are determinded by stepwise addition and deletion and the Bayes Information Criterion, where the maximum likelihood method is used to estimate the parameters. Numericla examples using real data and simulated data are provided to show the performance of proposed method.