Using the Smit - Suhl method, we have calculated the magnetic susceptibility tensor of magnetized ferromagnetic media with the uniaxial magnetocrystalline anisotropy. An analysis of the obtained formulas in the particular cases of a composite material with randomly oriented ferromagnetic inclusions of ellipsoidal shape has been carried out.