Non-equilibrium fission has been described by diffusion model. In order to describe the diffusion process analytically, the analytical solution of Smoluchowski equation in harmonic oscillator potential is obtained. This analytical solution is able to describe the probability distribution and the diffusive current with the variable x and t. The results indicate that the probability distribution and the diffusive current are relevant to the initial distribution shape, initial position, and the nuclear temperature T; the time to reach the quasi-stationary state is proportional to friction coefficient β, but is independent of the initial distribution status and the nuclear temperature T. The prerequisites of negative diffusive current are justified. This method provides an approach to describe the diffusion process for fissile process in complicated potentials analytically.