In this paper, we determine analytically the demagnetization field within a uniformly magnetized ellipsoid embedded in an infinite anisotropic media. The anisotropic boundary condition consists of decomposing the external magnetization into an isotropic part parallel to the external field and an anisotropic part. The resulting demagnetization field within the ellipsoid is then calculated by superposing the demagnetization fields produced by the ellipsoid magnetization, the isotropic part of the external magnetization and field, and the anisotropic part of the external magnetization. Such computation is needed in developing accurate anisotropic multi-scale models of electrical steel sheets, for example.