On the numerical resolution of anisotropic equations with high order differential operators arising in plasma physics