Abstract: The direct solution of mixed boundary value problem plays an important role in the accuracy improvement of Earth's geopotential model. In this paper the authors discuss the computation method of wigner 3-j symbol and present a direct solution of mixed boundary value problem. Simulated computation results shows high accuracy of this method.