Theory and performance of a subroutine for solving Volterra Integral Equations

AbstractThe author considers Volterra Integral Equations of either of the two forms $$u(x) = f(x) + \int\limits_a^x {k(x - t)g(u(t))dt, a \leqslant } x \leqslant b,$$ wheref, k, andg are continuous andg satisfies a local Lipschitz condition, or $$u(x) = f(x) + \int\limits_a^x {\sum\limits_{j = 1}^m {c_j (x)g_j (t,u(t))dt} ,} $$ wheref,cj, andgj,j=1,2,...,m, are continuous and eachgj satisfies a local Lipschitz condition in its second variable. It is shown that in each case the respective integral equation can be solved by conversion to a system of ordinary differential equations which can be solved by referring to a described FORTRAN subroutine. Subroutine VE1.In the case of the convolution equation, it is shown how VE1 converts the equation via a Chebyshev expansion, and a theorem is proved, and implemented in VE1, wherein the solution error due to truncation of the expansion can be simultaneously computed at the discretion of the user. Some performance data are supplied and a comparison with other standard schemes is made. Detailed performance data and a program listing are available from the author.ZusammenfassungDer Verfasser betrachtet zwei Typen Volterrascher Integralgleichung; die erste besitzt die Form $$u(x) = f(x) + \int\limits_a^x {k(x - t)g(u(t))dt, a \leqslant } x \leqslant b,$$ wobeif, k, undg stetig sind, undg eine lokale Lipschitz-Bedingung erfüllt; die zweite Gleichung ist $$u(x) = f(x) + \int\limits_a^x {\sum\limits_{j = 1}^m {c_j (x)g_j (t,u(t))dt} ,} $$ wobeif,cj, undgj (j=1,...m) stetig sind und jede der Funktionengj eine lokale Lipschitz-Bedingung in den zweiten Argumenten erfüllt. Für beide Fälle wird gezeigt, daß die Integralgleichung durch Zurückführung auf ein System gewöhnlicher Differentialgleichungen gelöst werden kann, wobei letzteres mit einer beschriebenen FORTRAN Subroutine VE1 behandelt werden kann.