In this research, numerical solution of integral algebraic - equations are investigated by using Bernoulli collocation method. This attempt includes three chapters which are organized as follow. Chapter 1, is devoted to fundamental concepts of Bernoulli polynomials and expanded functions to its Bernoulli polynomials series. In chapter 2, Bernoulli matrix method is used for obtaining the numerical solution of systems of Volterra integral equations. Numerical solution of integral algebraic equations by using this method is presented in Chapter 3. To demonstrate the effectiveness of our work several numerical examples are presented at the end of each chapter.