In this paper we extend a method for iteratively improving slow manifolds so that it also can be used to approximate the fiber directions. The extended method is applied to general finite-dimensional real analytic systems where we obtain exponential estimates of the tangent spaces to the fibers. The method is demonstrated on the Michaelis--Menten--Henri model and the Lindemann mechanism. The latter example also serves to demonstrate the method on a slow-fast system in nonstandard slow-fast form. Finally, we extend the method further so that it also approximates the curvature of the fibers.
S I a M Journal on Applied Dynamical Systems, 2014, Vol 13, Issue 2