en تخصصي Special پژوهشي Research In this paper, a two-dimensional multi-wavelet is constructed in terms of Chebyshev polynomials. The constructed multi-wavelet is an orthonormal basis for space. By discretizing two-dimensional Fredholm integral equation reduce to a algebraic system. The obtained system is solved by the Galerkin method in the subspace of by using two-dimensional multi-wavelet bases. Because the bases of subspaces are orthonormal, so the above mentioned system has a small dimension and also high accuracy in approximating solution of integral equations. For one-dimensional case, a similar works are done in [4, 5], which they have small dimension and high accuracy. In this article, we extend one-dimensional case to two-dimensional by extending and by choosing good functions on two axes. Numerical results show that the above mentioned method has a good accuracy. Two-Dimensional, Multi-Wavelet, Integral Equations, Galerkin, Chebyshev 45 54 http://jamlu.lahijan.iau.ir/browse.php?a_code=A-10-1-154&slc_lang=en&sid=1 M. Rabbani 1003194753284600516 1003194753284600516 Yes