A numerical method is presented for obtaining individual stress intensity factor in arc-shaped interfacial crack, where the stress intensity factor is calculated using the L-integral and the mutual parameter $M_1$ proposed by Chen and Shield. This method is useful when there is a large contact zone alog the interfacial crack. Some numerical examples are presented here for various combination of elastic constants of matrix and inclusion by varying load angle.