Numerical techniques were discussed for its numerical solution, to deal with the problems of singularity and oscillation in the numerical integration in spectral doman, etc.