Solution of the Cauchy Problem for Ordinary Differential Equations Using the Collocation and Least Squares Method with the Pade Approximation

A new method for solving the Cauchy problem for an ordinary differential equation is proposed and implemented using the collocation and least squares method of increased accuracy. It is based on the derivation of an approximate nonlinear equation by a multipoint approximation of the problem under consideration. An approximate solution of the problem in the form of the Pade approximation is reduced to an iterative solution of the linear least squares problem with respect to the coefficients of the desired rational function. In the case of nonlinear differential equations, their preliminary linearization is applied. A significant superiority in accuracy of the method proposed in the paper for solving the problem over the accuracy of the NDSolve procedure in the Mathematica system is shown. The solution of a specific example shows the superiority in accuracy of the proposed method over the fourth-order Runge-Kutta method. Examples of solving the Cauchy problem for linear and non-linear equations with an accuracy close to the value of rounding errors during operations on a computer with numbers in the double format are given. It is shown that the accuracy of solving the problem essentially depends on the complexity of the behavior of the values of the right-hand side of the equation on a given interval. An example of constructing a spline from pieces of Pade approximants on partial segments into which a given segment is divided is given in the case when it is necessary to improve the accuracy of the solution.

Keywords

Cauchy problem; ordinary differential equation; Pade аpproximation; collocation and least squares method; high accuracy; Mathematica System.

References

1. Baker G.A.Jr. Essentials of Pade Approximants. New York, Academic Press, 1975.
2. Gonchar A.A., Rakhmanov E.A. On the Convergence of Joint Pade Approximations for Systems of Markov Type Functions. Proceedings of the Steklov Institute of Mathematics, 1981, vol. 157, pp. 31-50.
3. Gonchar A.A., Rakhmanov E.A. Equilibrium Distributions and the Rate of Rational Approximation of an Analytical Function. Mathematics of the USSR-Sbornik, 1989, vol. 62, no. 2, pp. 305-348. DOI: 10.1070/SM1989v062n02ABEH003242
4. Gonchar A.A., Novikova N.N., Khenkin G.M. Multipoint of Pade Approximants. Matematichesky Sbornik [Mathematical Collection], 1996, vol. 187, no. 12, pp. 57-86. (in Russian)
5. Suetin S.P. Pade Approximations and Effective Analytical Continuation of a Power Series. Russian Mathematical Surveys, 2002, vol. 57, no. 1, pp. 43-141. DOI:10.1070/RM2002v057n01ABEH000475
6. Gonchar A.A. Rational Approximations of Analytic Functions. Proceedings of the Steklov Institute of Mathematics, 2011, vol. 257, no. 2, pp. 44-57. DOI: 10.1134/S0081543811030047
7. Aptekarev A.I., Buslaev V.I., Martines-Finkelstein A., Suetin S.P. Pade Approximations, Continuous Fractions and Orthogonal Polynomials. Russian Mathematical Survey, 2011, vol. 66, no. 6, pp. 1049-1131. DOI: 10.1070/RM2011v066n06ABEH004770
8. Khovansky A.N. The Application of Continued Fractions and their Generalizations to Problems in Approximation Theory. P. Noordhoff, Groivingen, the Netherlands, 1963.
9. Velichko I.G., Tkachenko I.G., Balabanova V.V. Application of the Method of Continued Fractions to Obtain Approximations of the Case Solutions of Cauchy Problems for Differential Equations of the First Order. Bulletin of Science and Education in North-West Russia, 2015, vol. 1, no. 3, pp. 1-9. (in Russian)
10. Vishnevsky V.E., Zubov A.V., Ivanova O.A. [Pade Approximation of the Solution of the Cauchy Problem]. Vestnik Sankt-Peterburgskogo Universiteta, 2012, no. 4, pp. 3-17. (in Russian)
11. Isaev V.I., Shapeev V.P. High-Accuracy Versions of the Collocations and Least Squares Method for the Numerical Solution of the Navier-Stokes Equations. Computational Mathematics and Mathematical Physics, 2010, vol. 50, no. 10, pp. 1670-1681. DOI: 10.1134/S0965542510100040
12. Shapeev V.P., Vorozhtsov E.V. CAS Application to the Construction of the Collocations and Least Residuals Method for the Solution of the Burgers and Korteweg-de Vries-Burgers Equations. International Workshop on Computer Algebra in Scientific Computing, Warsaw, 2014, vol. 8660, pp. 432-446. DOI: 10.1007/978-3-319-10515-4-31.
13. Shapeev V.P., Belyaev V.A. [Variants of the Collocation Method and the Least Residuals of Increased Accuracy in a Region with a Curvilinear Boundary]. Vychislitel'nye Tekhnologii [Computational Technologies], 2016, vol. 21, no 5. pp. 95-110. (in Russian)
14. Belyaev V.A., Shapeev V.P. Versions of the Collocation and Least Squares Method for Solving Biharmonic Equations in Non-Canonical Domains. AIP Conference Proceedings, 2017, vol. 1893, no. 1, article ID: 030102. DOI: 10.1063/1.5007560.
15. Shapeev V.P., Belyaev V.A., Golushko S.K., Idimeshev S.V. New Possibilities and Applications of the Least Squares Collocation Method. EPJ Web of Conferences, 2018, vol. 173, article ID: 01012. DOI: 10.1051/epjconf/201817301012
16. Vorozhtsov E.V., Shapeev V.P. On the Efficiency of Combining Different Methods for Acceleration of Iterations at the Solution of PDEs by the Method of Collocations and Least Residuals. Applied Mathematics and Computation, 2019, vol. 363, pp. 1-19. DOI: 10.1016/j.amc.2019.124644
17. Shapeev V.P., Vorozhtsov E.V. High-Accuracy Numerical Solution of the Second-Kind Integral Equations in the Mathematica Environment. Journal of Multidisciplinary. Engineering Science and Technology, 2018, vol. 5, no. 12, pp. 9308-9319.
18. Dyakonov V.P. Mathematica 5/6/7. Polnoe rukovodstvo [Mathematica 5/6/7. Complete Guide]. Мoscow, DMK Press, 2010. (in Russian)