This work uses a new boundary integral equation (BIE) and finite element method (FEM) to analyze an infinite anisotropic plate containing two elliptic/circular holes subjected to remote bending or twisting moments. The foundation of the boundary integral equation is the classical plate theory with Cauchy integral formula. The BIE is used to calculate the curvatures and moments on the boundaries directly. Numerical examples are given for orthotropic and isotropic plates with circular or elliptic holes under uniform bending and twisting moments. Comparison of the numerical results with the analytic solution for one hole shows that in general BIE can achieve higher accuracies in evaluating moments while BIEs and FEM have comparable accuracies for computing deflections.