In this paper, we combine the Muskhelishvili's complex variable method and boundary collocation method, and choose a set of new stress function based on the stress boundary condition of crack surface, the higher precision and less computation are reached. This method is applied to calculating the stress intensity factor for a finite plate with an inclined crack. The influence of θ (the obliquity of crack) on the stress intensity factors, as well as the number of summation terms on the stress intensity factor are studied and graphically represented.