Numerical study on tunneling splitting in biaxial spin systems is done by performing diagonalization of the Hamilton operator.It is found that the calculated energy splitting agrees quantitatively with theoretical prediction of instanton method.Our result shows that both the instanton method and the large spin limit work well for the total spin around 10.By including the fourth-order term in Hamiltonian,experimental observation can be re-covered quantitatively.