The variational approach based on the Bogoliubov inequality for the free energy is used to study the three-dimensional anisotropic Heisenberg XXZ model with a crystal field. The magnetization and the phase diagrams are obtained as a function of the parameters of the Hamiltonian. Limiting cases,...