"FixedStep" Method for NDSolve
Introduction
It is often useful to carry out a numerical integration using fixed step sizes.
For example, certain methods such as
and
carry out a sequence of fixed-step integrations before combining the solutions to obtain a more accurate method with an error estimate that allows adaptive step sizes to be taken.
The method
allows any one-step integration method to be invoked using fixed step sizes.
This loads a package with some example problems and a package with some utility functions.
Examples
Define an example problem.
| Out[5]= |  |
This integrates a differential system using the method

with a fixed step size of

.
| Out[6]= |  |
Actually the

method has no adaptive step size control. Therefore, the integration is already carried out using fixed step sizes so the specification of

is unnecessary.
| Out[8]= |  |
Here are the step sizes taken by the method

for this problem.
| Out[10]= |  |
This specifies that fixed step sizes should be used for the method

.
| Out[12]= |  |
The option MaxStepFraction provides an absolute bound on the step size that depends on the integration interval.
Since the default value of
MaxStepFraction is

the step size in this example is bounded by one-tenth of the integration interval, which leads to using a constant step size of

.
| Out[15]= |  |
By setting the value of
MaxStepFraction to a different value, the dependence of the step size on the integration interval can be relaxed or removed entirely.
| Out[17]= |  |
Option Summary
| | |
| Method | None | specify the method to use with fixed step sizes |
Option of the method
.