A general discussion of the simulation procedure of the full susceptibility tensor and isothermal magnetization pseudovector for compounds comprising weakly-interacting magnetic centers is presented. A single-crystal-sample as well as a powder-sample case are considered. The procedure is used to obtain explicit expressions for the full susceptibility tensor for spins S=1, 3/2, 2, and 5/2 for non-vanishing rhombic local anisotropy and any form of spectroscopic tensor.