In this paper the boundary collocation method is presented for computing the stress intensity factors for an internal crack in a finite anisotropic plate. The stress functions are assumed such that they can represent the stress singularity at the crack tips, satisfying not only the governing equations of the anisotropic plate theory in the domain, but also the stress-free conditions on the crack surfaces. Therefore, only the boundary conditions of the plate need to be considered, and they can be satisfied approximately by the Boundary Collocation Method. Numerical examples demonstrated that the proposed method gives satisfactory results compared with the existing solutions.