In this paper, we propose two arrays, containing the coefficients of fractional Adams-Bashforth and Adams-Moulton methods, and also recursive relations to produce the elements of these arrays. Then, we illustrate the application of these arrays in a suitable way to construct higher-order fractional linear multi-step methods in general form, with extended stability regions. The effectiveness of the new method is shown in comparison with some available previous results in an illustrative test problem.