in this paper an exponentially fitted finite difference method is presented for solving singularly perturbed two-point boundary value problems with the boundary layer. a fitting factor is introduced and the model equation is discretized by a finite difference scheme on an uniform mesh. thomas algorithm is used to solve the tri-diagonal system. the stability of the algorithm is investigated. it ...