A numerical algorithm to analyze the plane-wave three-dimensional oblique incidence on a strip grating is presented. Electromagnetic field is decomposed into vector Floquet harmonics of the E-type and H-type modes. To impose boundary conditions on the incident, reflected and transmitted waves, two integral equations of Fredholm of first kind are obtained. These equations are solved numerically with the standard Galerkin procedure, and the convergence of the algorithm is examined numerically. Since the superficial current near the edges of a conducting strip have been taken into account, the computational algorithm shows a fast convergence. Results are compared with other numerical results available in the literature to demonstrate the accuracy of the proposed method. Some numerical results are shown at normal and oblique incidence, to make clear the behaviour of the grids at different azimuth incidence angles.