In this work we utilize the boundary integral equation and the Dual Reciprocity
Boundary Element Method (DRBEM) for the solution of the steady state
convection-diffusion-reaction equations with variable convective coefficients in twodimension.
The DRBEM is a numerical method to transform the domain integrals into
the boundary only integrals by using the fundamental solution of Helmholtz equation.
Some examples are calculated to confirm the accuracy of the approach. The results
obtained by the analytic solutions are in good agreement with ones provided by the