The resolvent Rλ(A) of a complex r×r matrix A is an analytic function in any domain with empty intersection with the spectrum ΣA of A. The well known Taylor expansion of Rλ(A) in a neighborhood of any given λ0∉ΣA is modified taking into account that only the first powers of Rλ0(A) are linearly independent. The main tool in this framework is given by the multivariable polynomials depending on the invariants v1,v2,…,vr of Rλ(A) (m denotes the degree of the minimal polynomial). These functions are used in order to represent the coefficients of the subsequent powers of Rλ0(A) as a linear combination of the first m of them.