In this article, we study the multivariate quantitative smooth approximation under differentiation of functions. The approximators here are multivariate neural network operators activated by the symmetrized and perturbed hyperbolic tangent activation function. All domains used here are infinite. The multivariate neural network operators are of quasi-interpolation type: the basic type, the Kantorovich type, and the quadrature type. We give pointwise and uniform multivariate approximations with rates. We finish with illustrations.