Multigrid method is an iterative method with high rate of convergence. Different size of grids will be used in this method. Indeed, in this method we start to solve the problem on coarse grid and smooth the error. This is called relaxation on the solution. Than we solve the problem on fine grid. We will continue this procedure on coarse and fine grid up to find an acceptable solution. This method can be used to solve differential and integral equations. In this article we will use this method to solve Fredholm integral equation of second kind. To do this we facilitated a FORTRAN computer code on which we can solve this kind of integral equation for different kernels.