Abstract
In order to optimize the stability regions we construct families of k multi-step methods with free parameters and of order k+1. We exhibit parameters for which the real interval of stability is arbitrarly large and using the Flexible Tolerance Method we show how these parameters optimize - under a predefined criterion - the stability region. Computational experience, while limited, has indicated that the procedure is efficient