  This paper continues the studies of symbolic integration by focusing on the
stability problems on D-finite functions. We introduce the notion of stability
index in order to investigate the order growth of the differential operators
satisfied by iterated integrals of D-finite functions and determine bounds and
exact formula for stability indices of several special classes of differential
operators. With the basic properties of stability index, we completely solve
the stability problem on general hyperexponential functions.
