This article describes the use of the SACI package-a package for calculating the energy levels and wavefunctions of a multi-electron quantum dot modelled as a 2D harmonic well with electrons interacting through a Coulomb potential and under the influence of a perpendicular magnetic field. ‡ Introduction Quantum dots are artificially fabricated atoms, in which charge carriers are confined in all three dimensions just like electrons in real atoms. Consequently, they exhibit properties normally associated with real atoms such as quantised energy levels and shell structures. These properties are described by the electron wavefunctions whose evolution is governed by the Schrödinger equation and the Pauli exclusion principle.There are many methods available to solve the Schrödinger equation for multiple electrons. They roughly fall into the categories of the diagonalisation method, mean-field density-functional theory, and the self-consistent field approach. One of the first theoretical studies of quantum dots was by Pfannkuche et al. [1], who compared the results of Hartree-Fock self-consistent calculations and exact diagonalisation of the Hamiltonian for two electrons in a circularly symmetric parabolic potential. They found good agreement between the two methods for the triplet state but marked differences for the singlet state, indicating important spin correlations were not included properly in their HartreeFock model. This suggests that the proper treatment of electron spins is crucial for correctly obtaining the electronic structures in quantum dots.Examples of self-consistent field approaches in the literature include Yannouleas and Landman [2,3], who studied circularly symmetric quantum dots using an unrestricted spin-space Hartree-Fock approach, and McCarthy et al. [4], who developed a Hartree-Fock Mathematica package. Macucci et al. [5] studied quantum dots with up to 24 electrons using a mean-field local-density-functionalThe Mathematica Journal 10:2