We aimed to evaluate the ground-state and excite-state energy eigenvalue (En), wave function, and the time-independent correlation function of the atomic density fluctuation of a particle under the harmonics oscillator Cosine asymmetric potential (Saad et al. 2013). Instead of using the 6-point kernel of 4 Green's function (Cherroret and Skipetrov, 2008), averaged over disorder, we use the numerical shooting method (NSM) to solve the Schrödinger equation of quantum mechanics system with Cosine asymmetric potential. Since our approach does not use complicated formulas, it requires much less computational effort when compared to the Green functions techniques (Cherroret and Skipetrov, 2008). We show that the idea of the program of evaluating time-independent correlation function of atomic density is underdamped motion for the Cosine asymmetric potential from the numerical shooting method of this problem. Comparison of the time-independent correlation function obtained from numerical shooting method by and correlation function experiment by Kasprzak et al. (2008). We show the intensity of atomic density fluctuation ( ( ) =̃( ) −̃( )) in harmonics oscillator Cosine asymmetric potential by numerical shooting method.