This paper is concerned with a Taylor series-based continuation algorithm, the so-called Asymptotic Numerical Method (ANM). It describes a generic continuation procedure to apply the ANM principle at best, in other words, that presents a high level of genericity without paying the price of this genericity by low computing performances. The way to quadratically recast a system of equation is now part of the method itself, and the way to handle elementary transcendental function is detailed with great attention. A sparse tensorial formalism is introduced for the internal representation of the system, which, when combined with a block condensation technique, provides a good computational efficiency of the ANM. Three examples are developed to show the performance and the versatility of the implementation of the continuation tool. Its robustness and its accuracy are explored. Finally, the potentiality of this method for complex nonlinear finite element analysis is enlightened by treating 2D elasticity problems with geometrical nonlinearities. KEYWORDS asymptotic numerical method, continuation, finite element method, nonlinear systems, quadratic recast, Taylor series *The term "numerical" stands in the name of the method because of this finite element discretization. Int J Numer Methods Eng. 2019;119:261-280.wileyonlinelibrary.com/journal/nme