It is shown how boundary conditions are numerically imposed to obtain three-dimensional magnetic field solutions in electrical devices. The effect of changing the size of the solution volume and type of boundary condition imposed on the numerically obtained parameters of the device is examined in detail.