  The "Iterative Laplace Transform Method" is used to solve the Fokker-Planck
equation for finding the time evolution of the heavy quarks distribution
functions such as charm and bottom in quark gluon plasma. These solutions will
lead us to calculation of nuclear suppression factor RAA. The results have good
agreement with available experiment data from the PHENIX collaboration.
