Summary: In a recent theoretical study [see the authors, Int. J. Plast. 22, No. 12, 2346–2365 (2006; Zbl 1229.74015)], a constitutive model for anisotropic elastoplasticity at finite strains has been developed. The model is based on the multiplicative decomposition of deformation gradient. The stored energy function as well as the flow rule have been considered as quadratic functions of their arguments. In both cases, the list of arguments is extended to include structural tensors which describe the anisotropy of the material response at hand. Nonlinear isotropic hardening is considered as well. In this paper, we present the integration of the constitutive law. The associative flow rule is integrated using the exponential map which preserves the plastic incompressibility condition. The numerical treatment of the problem is fully developed, and expressions related to the local iteration and the consistent tangent operator are considered in detail. It is shown that while the consistent linearisation of the model is quite complicated, it still can be achieved if various intriguing implicit relations are identified and correctly dealt with. Various numerical examples of three-dimensional deformations of whole structural components are presented. The examples clearly illustrate the influence of anisotropy on finite elastoplastic deformations.