Abstract In this paper, an interface is built to compute the normal modes splitting parameters, based on the MINEOS software, ode45 toolbox of MATLAB, and interpolation function. The hydrochloric coefficient of J2 and J4 are used to compute the first and second order term of Earth’s ellipticity. The splitting parameters and singlets frequencies of 0S2 are generated based on the PREM and 1066A. It is proven to be correct when compared with the results of Dahlen. Our research shows that the second order ellipticity has a little effect on the splitting parameters, but almost no influence on the singlet frequencies. Finally, the coupling matrix and coupling strength between 0S2 and selected modes are given. This lays a good base to use the splitting function to inverse the interior Earth structure.