Abstract: A new implicit finite volume frequency domain (FVFD) method for solving the Maxwell equations was proposed to accurately simulate the spatial distribution and scattering characteristic of electromagnetic fields for continuous incident wave. The FVFD is a natural extension of finite volume time domain (FVTD) method, in which the computational domain is transformed from time to frequency, and the variables are changed from real type to three-dimensional complex type. In the FVFD method, there are two semi-discrete iteration processes, stable pseudo time marching and LU-SGS implicit iteration of flux residual. Thus, acceleration techniques can be easily implemented into it. In the calculation, large Courant-Friedrichs-Lewyor (CFL) numbers can be adopted, which confirms its unconditional stable performance. The FVFD was used to simulate a series of canonical 2D and 3D perfect conductor, composite dielectric and complex shaped objects. The results agreed well with those simulated by method of moment (MOM) and Mie series, which demonstrated its high precision and extensive application scenarios of FVFD method.