A numerical scheme is proposed to obtain the individual stress intensity factors in an axisymmetric crack and in a three dimensional mixed mode crack. The method is based on the path independence of J and M integral and mutual or two-state conservation integral , which involves two elastic fields. Some numerical example are presented to investigate the effectiveness and applicability of the method for and axisymmetric crack and a three dimensional penny shaped crack problem under mixed mode.