The method of numerical analysis of the singular integral equations for crack analysis by using the method of continuously distributed dislocations model is presented. This method is based on the concept of boundary collocation and reduces singular integral equation into linear algebraic equations. Presented method made it possible to obtain the solution of the singular integral equations in symmetrical or asymmetrical problems including two or more unknown functions. The accuracy of stress intensity factors obtained by this method was verified by comparing these values with the exact solution or the reliable numerical solution obtained by the other researchers.