In this dissertation, the solution of boundary value problems for materials with nonlocal behavior has been studied. Eringen's model of nonlocal integral elasticity is considered that it di?ers from the justify; LINE-HEIGHT: normal; TEXT-INDENT: 0in; MARGIN: 0in 0in 8pt; unicode-bidi: embed; DIRECTION: ltr" In the presented method, the solution is split into two parts, i.e. homogeneous and particular parts. Introduction of the EBFs into the homogeneous governing differential equations leads to a characteristic equation through which the exponents of the EBFs are defined. This combination of the bases satisfies the differential equation of equilibrium in a central zone. For many cases, the characteristic equation possesses some multiple roots. In such situations polynomial functions are added to EBFs. After selection of the EBFs as the bases for the approximation, the unknown coefficients of the series are determined by the satisfaction of the boundary conditions and the equilibrium equations on boundary zone, in an element based manner, through a weight integral technique. The particular part of the solution is constructed by a rather similar approach. In this thesis, EBFs are evaluated for 1D and 2D problems. The capabilities of the method have been shown by defining some norms of residuals and comparing the results with semi-analytical solutions. Key Words Exponential basis functions, Eringen's nonlocal elasticity, Fourier series, Chebyshev polynomials, weight integral