Optimal Solution of Linear Ordinary Differential Equations by Conjugate Gradient Method

arXiv1805.01085MaRDI QIDQ6301100FDOQ6301100

Authors: Wenqiang Yang, Wenyuan Wu, Robert M. Corless 

Publication date: 2 May 2018

Abstract: Solving initial value problems and boundary value problems of Linear Ordinary Differential Equations (ODEs) plays an important role in many applications. There are various numerical methods and solvers to obtain approximate solutions represented by points. However, few work about optimal solution to minimize the residual can be found in the literatures. In this paper, we first use Hermit cubic spline interpolation at mesh points to represent the solution, then we define the residual error as the square of the L2 norm of the residual obtained by substituting the interpolation solution back to ODEs. Thus, solving ODEs is reduced to an optimization problem in curtain solution space which can be solved by conjugate gradient method with taking advantages of sparsity of the corresponding matrix. The examples of IVP and BVP in the paper show that this method can find a solution with smaller global error without additional mesh points.

This page was built for publication: Optimal Solution of Linear Ordinary Differential Equations by Conjugate Gradient Method

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6301100)