This study proposed a technique for inverse problem, linear approximation of contact position and loading in single and double meshing of spiral bevel gear , using 2-dimension model considered near the tooth by root stress. Determine root stress is carried out far the gear tooth by finite element method and boundary element method. Boundary element discretization near contact point is carefully performed to keep high computational accuracy. And from those estimated results, the comparing estimate value with boundary element method value was discussed.