We describe a multipole method for calculating the modes of microstructured optical fibers. The method uses a multipole expansion centered on each hole to enforce boundary conditions accurately and matches expansions with different origins by use of addition theorems. We also validate the method and give representative results.