Abstract: The Boundary Element Method with 8-nodes quatratic elements is used for the computation of stress intensity factors of surface crack problems. The 8-nodes singlar elements are used along the crack front. In order to transmit easily from densed network to sparse network the 4 to 8 nodes element are used on the other boundary. The integral method of equal precision and other numerical techniques used in the paper improve the effectiveness of the solution in some way. The comparsion with Newman's high degree finite element solution shows that a satisfactory result has been obtained by much less nodes.