Diffractive properties of two dimensional relief gratings made with electromagnetic anisotropic medium are theoretically analysed by using rigorous modal theory, and the degradation is discussed. Reflection transmission coefficient matrix (RTCM) algorithm and correct Fourier factorization rules are adopted to improve stability, convergence and calculation efficiency of numerical computation. Numerical compution on gratings shows that the method is satisfying in credibility, stability and convergence.