The investigation of the heating mechanism of a material sample in an energized electromagnetic cavity requires the understanding of the interaction of the electromagnetic fields with the material sample in a cavity. The key step for this understanding is to quantify the distribution of the induced electric field by solving an Electric Field Integral Equation (EFIE) or a Magnetic Field Integral Equation (MFIE) inside the material sample placed in the cavity. In this study, a complete set of vector wave functions which include both solenoidal and irrotational functions are employed. EFIE and MFIE are derived based on the expansion of these vector wave functions. Due to the slow convergence of the dyadic Green's function, some schemes were developed to improve the convergence rate of the numerical calculation in order to save computation time. Numerical results of some material samples agree well with theoretical estimations.