Invoking the dsolve function with the numeric option and output=procedurelist (default), output=listprocedure or output=operator causes dsolve to return procedure(s) to numerically solve the input ODE system.

Prior to returning this procedure, the input system is converted to a first order system (by introduction of new variables if necessary), then solved with respect to those first order derivatives giving a system of the form:

During the course of the numerical solution of the ODE system, it is necessary to evaluate the functions $\mathrm{f1}$, $\mathrm{f2}$ above for different values of $t$, $x$, $y$, to obtain the values of the derivatives $x\'$, $y\'$.

The option maxfun=n in the call to dsolve,numeric gives an upper limit of $n$ on the total number of evaluations performed on $\mathrm{f1}$, $\mathrm{f2}$ for any single call to the procedure returned by dsolve. Note that the evaluation of both $\mathrm{f1}$, $\mathrm{f2}$ at one set of values of $t$, $x$, $y$ is considered to be a single evaluation.

This option is provided as a means of limiting the amount of work performed on any single call to the solution procedure, and as a means of preventing the computation from running forever when it encounters a singularity in the solution of the ODE.

When the maxfun function evaluation limit is exceeded, the solution procedure halts with an error of the form Error, (in <proc>) too many function evaluations in <method>.

This option is available for the rkf45, ck45, rosenbrock, dverk78, lsode, gear, and classical methods.

For the dverk78, lsode, and gear methods, this option is disabled by default. For the rkf45, ck45, and rosenbrock methods, the default setting is maxfun=30000, while for the classical methods, the default setting is maxfun=50000.

A setting of maxfun=0 disables the option.

This option is useful for some variable step methods, as without this limit it is possible for a calculation to run forever, and never get past a certain point in the integration (if the step size were, for example, halving at each step of the numerical integration). This can happen when singularities are present. It is also useful when the function to be evaluated is large and complex or is expensive to evaluate (say, for example, a procedure call that does a numerical sum).

$\mathrm{DE}≔\mathrm{diff}\left(y\left(t\right)\&comma;t\right)=\cos \left(t^2\right)y\left(t\right)$

$\mathrm{DE}≔\frac{\&DifferentialD;}{\&DifferentialD;t}y\left(t\right)=\cos \left(t^2\right)y\left(t\right)$

$\mathrm{ysol}≔\mathrm{dsolve}\left(\left\{\mathrm{DE}\&comma;y\left(0\right)=1\right\}\&comma;\mathrm{numeric}\&comma;\mathrm{method}=\mathrm{rkf45}\right)$

$\mathrm{ysol}≔\mathbf{proc}\left(\mathrm{x\_rkf45}\right)...\mathbf{end\; proc}$

$\mathrm{ysol}\left(5\right)$

$\left[t=5.\&comma;y\left(t\right)=1.84313059855996\right]$

$\mathrm{ysol}\left(200\right)$

Error, (in ysol) cannot evaluate the solution further right of 50.435782, maxfun limit exceeded (see ?dsolve,maxfun for details)

$\mathrm{ysol}≔\mathrm{dsolve}\left(\left\{\mathrm{DE}\&comma;y\left(0\right)=1\right\}\&comma;\mathrm{numeric}\&comma;\mathrm{method}=\mathrm{rkf45}\&comma;\mathrm{maxfun}=500000\right)$

$\mathrm{ysol}≔\mathbf{proc}\left(\mathrm{x\_rkf45}\right)...\mathbf{end\; proc}$

$\mathrm{ysol}\left(200\right)$

$\left[t=200.\&comma;y\left(t\right)=1.87563347676955\right]$

$\mathrm{dsys}≔\left\{\mathrm{diff}\left(y\left(t\right)\&comma;t\right)=\frac{{y\left(t\right)}^2}{1-100y\left(t\right)}\&comma;y\left(0\right)=\frac{1}{200}\right\}\&colon;$

$\mathrm{ysol}≔\mathrm{dsolve}\left(\mathrm{dsys}\&comma;\mathrm{numeric}\&comma;\mathrm{method}=\mathrm{rkf45}\&comma;\mathrm{abserr}=1.\times {10}^{\mathrm{-7}}\&comma;\mathrm{relerr}=1.\times {10}^{\mathrm{-7}}\right)$

$\mathrm{ysol}≔\mathbf{proc}\left(\mathrm{x\_rkf45}\right)...\mathbf{end\; proc}$

$\mathrm{ysol}\left(31\right)$

Error, (in ysol) cannot evaluate the solution further right of 30.685303, maxfun limit exceeded (see ?dsolve,maxfun for details)

$\mathrm{ysol}≔\mathrm{dsolve}\left(\mathrm{dsys}\&comma;\mathrm{numeric}\&comma;\mathrm{method}=\mathrm{rkf45}\&comma;\mathrm{abserr}=1.\times {10}^{\mathrm{-7}}\&comma;\mathrm{relerr}=1.\times {10}^{\mathrm{-7}}\&comma;\mathrm{range}=0..31\right)$

$\mathrm{ysol}≔\mathbf{proc}\left(\mathrm{x\_rkf45}\right)...\mathbf{end\; proc}$

$\mathrm{ysol}≔\mathrm{dsolve}\left(\mathrm{dsys}\&comma;\mathrm{numeric}\&comma;\mathrm{method}=\mathrm{rkf45}\&comma;\mathrm{abserr}=1.\times {10}^{\mathrm{-7}}\&comma;\mathrm{relerr}=1.\times {10}^{\mathrm{-7}}\&comma;\mathrm{maxfun}=100000\right)$

$\mathrm{ysol}≔\mathbf{proc}\left(\mathrm{x\_rkf45}\right)...\mathbf{end\; proc}$

$\mathrm{ysol}\left(31\right)$

Error, (in ysol) cannot evaluate the solution further right of 30.685337, maxfun limit exceeded (see ?dsolve,maxfun for details)