Multistep Methods

These methods use the numerical solution at tn and also at earlier times, tn—1, tn—2, etc., to obtain the solution at tn+1. (In contrast, methods like [19] use the solution only at the time tn to obtain the solution at t„+1 and hence are called single-step methods.) The idea behind multistep methods is to use the solution computed at those earlier steps to predict not only the slope, given by the RHS of [1], of the solution at tn, but also the curvature (the second derivative) and possibly higher-order derivatives of the solution. This allows one to approximate it at tn + 1 with an accuracy higher than that achieved by [2]. For example, the formula for a second-order-accurate, two-step method for [1] can be derived from the Taylor expansion of the same order (see [4]):

T2 T

U«+1 = U„ + T U'„ + — U"„ = U„ + Tf + ^ (fn-fn- 1) [20]

In deriving [20], one uses U'„ = fn and its corollary,

U"„ = (fn — fn_1)/t + O(t), and omits terms of order

0(t3) and higher. To start this method, one uses f0 from the initial condition and f1 found by some single-step method. Another well-known two-step second-order method is a so-called leap-frog method:

(It should be noted that this method is unstable for the test equation [10] for any Re(A) < 0; its stability region is the segment along the imaginary A-axis shown in Figure 1.) Formulas for higher-order multistep methods, known as Adams methods, can be found in most textbooks. The advantage of multistep methods over the single-step RK methods is that the latter require at least m function evaluations per step for an mth-order-accurate RK method (e.g., [19] requires four function evaluations), whereas a multistep method can achieve the same accuracy with only one function evaluation for any order m. Thus, multistep methods are faster than the RK ones. The main disadvantage of the multistep methods is that it is difficult to make them use adaptive step size, because their formulas are inherently based on the assumption that all steps have the same size t. Another disadvantage of multistep methods is that they have smaller stability regions, which shrink with increasing the method's order. Therefore, currently these methods are not widely used in commercial software, where the adaptive embedded RK methods are used instead.

Was this article helpful?

0 0
Trash To Cash

Trash To Cash

This book will surely change your life due to the fact that after reading this book and following through with the steps that are laid out for you in a clear and concise form you will be earning as much as several thousand extra dollars a month,  as you can see by the cover of the book we will be discussing how you can make cash for what is considered trash by many people, these are items that have value to many people that can be sold and help people who need these items most.

Get My Free Ebook

Post a comment