A generalized Taylor method of order three for the solution of initial value problems in standard and infinity floating-point arithmetic

A well-known drawback of algorithms based on Taylor series formulae is that the explicit calculation of higher order derivatives formally is an over-elaborate task. To avoid the analytical computation of the successive derivatives, numeric and automatic differentiation are usually used. A recent alternative to these techniques is based on the calculation of higher derivatives by using the Infinity Computer-a new computational device allowing one to work numerically with infinities and infinitesimals. Two variants of a one-step multi-point method closely related to the classical Taylor formula of order three are considered. It is shown that the new formula is order three accurate, though requiring only the first two derivatives of y(t) (rather than three if compared with the corresponding Taylor formula of order three). To get numerical evidence of the theoretical results, a few test problems are solved by means of the new methods and the obtained results are compared with the performance of Taylor methods of order up to four.