A new numerical method for the fractional integral that only stores part history data is presented, and its discretization error is estimated. The method can be
Different from previous temporal evolution assumption, the spatially growing mode was employed to analyze the linear stability for the channel flow of fiber sus