Mathematics – Numerical Analysis
Scientific paper
2009-03-27
Numerical Functional Analysis and Optimization, Vol.32, No.5, 541-582, (2011)
Mathematics
Numerical Analysis
Errors concerning numbering of figures are fixed
Scientific paper
10.1080/01630563.2011.557917
Based on functional analysis, we propose an algorithm for finite-norm solutions of higher-order linear Fuchsian-type ordinary differential equations (ODEs) P(x,d/dx)f(x)=0 with P(x,d/dx):=[\sum_m p_m (x) (d/dx)^m] by using only the four arithmetical operations on integers. This algorithm is based on a band-diagonal matrix representation of the differential operator P(x,d/dx), though it is quite different from the usual Galerkin methods. This representation is made for the respective CONSs of the input Hilbert space H and the output Hilbert space H' of P(x,d/dx). This band-diagonal matrix enables the construction of a recursive algorithm for solving the ODE. However, a solution of the simultaneous linear equations represented by this matrix does not necessarily correspond to the true solution of ODE. We show that when this solution is an l^2 sequence, it corresponds to the true solution of ODE. We invent a method based on an integer-type algorithm for extracting only l^2 components. Further, the concrete choice of Hilbert spaces H and H' is also given for our algorithm when p_m is a polynomial or a rational function with rational coefficients. We check how our algorithm works based on several numerical demonstrations related to special functions, where the results show that the accuracy of our method is extremely high.
Hayashi Masahito
Sakaguchi Fuminori
No associations
LandOfFree
General theory for integer-type algorithm for higher order differential equations does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with General theory for integer-type algorithm for higher order differential equations, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and General theory for integer-type algorithm for higher order differential equations will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-21792