Abstract. We present a new algorithm for computing covariant power expansions of tensor fields in generalized Riemannian normal coordinates, introduced in some neighborhood of a parallelized k-dimensional submanifold (k = 0, 1, . . . < n; the case k = 0 corresponds to a point), by transforming the expansions to the corresponding Taylor series. For an arbitrary real analytic tensor field, the coefficients of such series are expressed in terms of its covariant derivatives and covariant derivatives of the curvature and the torsion. The algorithm computes the corresponding Taylor polynomials of arbitrary orders for the field components and is applicable to connections that are, in general, nonmetric and not torsion-free. We show that this computational problem belongs to the complexity class LEXP.