Abstract
A code based on the two-stage Gauss formula (order four) for second-order initial value problems of a special type is developed. This code can be used to obtain a low- to medium-precision integration for a wide range of problems in the class of oscillatory type, Hamiltonian problems, and time-dependent partial differential equations discretized in space by finite differences or finite elements. The iteration process used in solving for the stage values of the Gauss formula, the selection of the initial step size, and the choice of an appropriate local error estimator for determining the step size change according to a particular tolerance specified by the user are studied. Moreover, a global error estimate and a dense output at equidistant points in the integration interval are supplied with the code. Numerical experiments and some comparisons with certain standard codes on relevant test problems are also given.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
