A New Class of Highly Accurate Solvers for Ordinary Differential Equations

  • Authors:
  • Andreas Glaser;Vladimir Rokhlin

  • Affiliations:
  • Program in Applied Mathematics, Yale University, New Haven, USA 06511;Departments of Computer Science and Mathematics, Yale University, New Haven, USA 06511

  • Venue:
  • Journal of Scientific Computing
  • Year:
  • 2009

Quantified Score

Hi-index 0.00

Visualization

Abstract

We introduce a new class of predictor-corrector schemes for the numerical solution of the Cauchy problem for non-stiff ordinary differential equations (ODEs), obtained via the decomposition of the solutions into combinations of appropriately chosen exponentials; historically, such techniques have been known as exponentially fitted methods. The proposed algorithms differ from the classical ones both in the selection of exponentials and in the design of the quadrature formulae used by the predictor-corrector process. The resulting schemes have the advantage of significantly faster convergence, given fixed lengths of predictor and corrector vectors. The performance of the approach is illustrated via a number of numerical examples.