The diffraction properties of electromagnetic anisotropic gratings with arbitrary profile, illuminated at an arbitrary incidence angle are studied, by generalizing Fourier series expansion method of inverse rule and recursive algorithm of reflection transmission coefficient matrix. The physical concepts are clear in this method, and its numerical calculation, compared with those of other methods, are stable, convergent and fast.