Optimally accurate Petrov-Galerkin method of finite elements

We perform analysis for a finite elements method applied to the singular self-adjoint problem. This method uses continuous piecewise polynomial spaces for the trial and the test spaces. We fit the trial polynomial space by piecewise exponentials and we apply so exponentially fitted Galerkin method to singular self-adjoint problem by approximating driving terms by Lagrange piecewise polynomials, linear, quadratic and cubic. We measure the erroe in max norm. We show that method is optimal of the first order in the error estimate. We also give numerical results for the Galerkin approximation.