A 3-dimensional numerical method of determining the vector fields within a closed lossless electromagnetic system is described. The method is used to find the dominant-mode resonant frequencies for two cases of a rectangular cavity partially loaded with dielectric.