|Boomer Manual and Download|
|PharmPK Listserv and other PK Resources|
|Previous Page||Course Index||Next Page|
Figure 10.3.1 Equations for the Fourth Order Runge-Kutta Method
Notice the time at which the first derivatives (f terms, i.e. differential equations) are determined. Calculations are made at the initial time, two at half the stepsize beyond the initial time and at the final time. These four calculations allow the use of larger overall stepsizes with good accuracy. The method can be represented graphically.
Figure 10.3.2 Plot of Cp versus time Illustrating the 4th order Runge-Kutta Method
Compare the accuracy using the fourth order Runge-Kutta with the accuracy achieved with Euler's method. As with the previous Euler's method example the initial value is 100 and the rate constant is 0.3 hr-1.
Table 10.3.1 Table Illustrating Accuracy Achieved with 4th order Runge-Kutta Method
This seems to be an efficient method of numerically integrating differential equations. However, what stepsize should we use? We could repeat the calculation with a smaller stepsize, say half the previous value. If the two estimates were within limits we could confirm that they were both satisfactory. The problem is that this would require 8 more calculations. That is the original four calculation plus the additional 8 for the re-calculation with the stepsize half the previous value. A solution to this inefficiency is the Runge-Kutta-Fehlberg method.
return to the Course index
Material on this website should be used for Educational or Self-Study Purposes Only
iBook and pdf versions of this material and other PK material is available
Copyright © 2001-2021 David W. A. Bourne (email@example.com)