We present a simulation of the magnetization curves, energy, probability, and torque landscapes of uniaxial systems with up to five anisotropy constants. The total energy used in the simulation is the sum of the anisotropy and Zeeman energies. The exchange interaction is not considered in the present work in which we treat single-domain-particle systems within a classical mechanics-based model. Diverse features of the calculated magnetization curves are highlighted for the studied systems. These diverse features are strongly dependent on the sign and magnitude of the simulation parameters. The model is versatile enough to handle both hypothetical and real material systems, e.g. HoFe 11 Ti and Y 2 Co 17 .