The edge function method for anisotropic elasticity is based on the complex variable formulation and on the superposition of analytical solutions to the field equations. Herein the method is extended to the study of fracture interaction in anisotropic materials. Numerical results indicate that echelon crack formations are preferred. The effect of variation in Young's modulus is noted, while the variation of Poisson's ratio or shear modulus is not significant. Studies of multiple cracks suggest that the method is well suited to the efficient analysis of such problems, either on its own or in combination with other conventional schemes.