This work shows a self-consistent numerical model of ion cyclotron resonant heating (ICRH) in a plasma of arbitrary (2D or 3D) geometry. Effects of anisotropy and full shaping on both MHD equilibrium and the dielectric tensor are explored. In particular, it is shown that the dielectric tensor becomes dependent on poloidal angle. Also, the anisotropic effects dominate the dielectric tensor. Using single particle simulations it is demonstrated that anisotropy develops due to ICRH. The changing of the equilibrium was successfully fed back into the equilibrium.