The exact diagonalization method using a spin-adapted basis is employed to calculate the electronic structure of a multi-electron quantum dot. By isolating spin and orbital angular momentum eigenstates, we have significantly reduced the size of the matrices required in comparison with the standard configuration interaction method. A novel approach to the simplification of the interaction integrals that arise in the calculation is also presented, which allows exact evaluation of the Hamiltonian matrix required in the calculations. This Mathematica package permits accurate calculation of energy levels and wave functions for both ground and excited states of multiple electrons confined in a circular quantum dot. (c) 2005 Wiley Periodicals, Inc.