Print anisotropy as Stevens for given N

[kbarros]

This adds a named parameter N to the function print_anisotropy_as_stevens(). Also includes various cleanups, especially improvements to printing of polynomials.