This paper presents a method to reduce time complexity of the computation of higher–order tensor lines. The method can be applied to higher–order tensors and the spherical harmonics representation, both widely used in medical imaging. It is based on a gradient descend technique and integrates well into fiber tracking algorithms. Furthermore, the method improves the angular resolution in contrast to discrete sampling methods which is especially important to tractography, since there, small errors accumulate fast and make the result unusable. Our implementation does not interpolate derived directions but works directly on the interpolated tensor information. The specific contribution of this paper is a fast algorithm for tracking lines tensor fields of arbitrary order that increases angular resolution compared to previous approaches.