In paramagnetic materials with multiple unpaired electrons per magnetic center, even when the ground state has no first-order orbital angular momentum, structural anisotropy around the magnetic center gives rise to zero-field splitting and g-factor anisotropy, resulting in magnetic anisotropy. This paper presents an exact algebraic expression for the angle-dependent magnetization in the axially symmetric spin-quartet (S = 3/2) ground state without first-order orbital angular momentum, taking into account zero-field splitting and g-factor anisotropy. The use of this angle-dependent expression enables accurate and rapid simulations of magnetic anisotropy, particularly for powder averaging of magnetization. The study also includes simulations of anisotropic saturation magnetization curves, cross-sectional views of angle-dependent magnetization, and three-dimensional representations of magnetization. Furthermore, the powder average and arithmetic average of anisotropic magnetization have been compared, and the circumstances under which powder average calculations are necessary have been discussed.