Efficient, robust and accurate algorithms are proposed for solving the multiple scattering problem by M circular obstacles for the whole spectrum of frequency. The representation of the solution is based on an integral equation formulation next solved by using Fourier basis. Numerical examples are provided to show that the approaches are efficient