The determination of the stress intensity factors is investigated by using the surface integral defined around the crack tip of the structure. In this work, the integral method is derived naturally from the standard path integral J. But the use of the surface integral is also extended to the case where body forces act. Computer program for obtaining the stress intensity factors $K_I$ and $K_{II}$ is developed, which prepares input variables from the result of the conventional finite element analysis. This paper provides a parabolic smooth curve function. By the use of the function and conventional element meshes in which the aspect ratio (element length at the crack tip/crack length) is about 25 percent, relatively accurate $K_I$ and K_{II}$ values can be obtained for the outer integral radius ranging from 1/3 to 1 of the crack length and for inner one zero.