A hybrid method coupling photoelasticity to a numerical-computational method which implements the William s (modes I and II) and Creager s (mode I) equations has been developed to determine the mixed-mode stress intensity factors in sharp notches and blunt cracks. The equations take into account the presence or not of a radius in the tip of the notch. To the proposed equations was added a complete polynomial, which represents the non-singular stress field. Three computational programs were developed for both formulations (William s and Creager): one to determine KI, KII and the non-singular terms, the second to evaluate the error between the actual situation and the results obtained and the third to regenerate the isochromatic fringes. The method has been tested in bars with cracks and notches (taking or not into account the existence of the radius at the depth of the notch) and the results were compared to experimental and analytical data found in the literature. Factors which have influence on the determination of KI, KII, were discussed: the number of data points and sector angle where those points are collected as well as the influence of the number of non-singular terms and some aspects of the numeric method. Two applications were studied: Charpy type test specimen under tension and bending and a beam with deep simetrycal grooves. With the values of KI and KII obtained by the Creager s equations one can determine the valeu of sigma x and sigma y in the analyzed situation, as well as the Kt value.