Numerical modeling for anisotropic magnetic media including saturation effects