Abstract
A method to compute explicit solutions of homogeneous triangular systems of first-order linear initial-value ordinary differential equations with constant coefficients is described. It is suitable for the limited case of well separated eigenvalues, or for multiple zero eigenvalues provided the entire column corresponding to a zero eigenvalue is zero. The solution for the case of constant inhomogeneity is described. The method requires only the computation of a constant matrix using a simple recurrence. Computing the solutions of the system from that matrix, for values of the independent variable, requires one to exponentiate only the diagonal of a matrix. It is not necessary to compute the exponential of a general triangular matrix. Although this work was motivated by a study of nuclear decay without fission or neutron absorption, which is used throughout as an example, it has wider applicability.
Supplemental Material
Available for Download
Software for Explicit Solutions of Triangular Systems of First-Order Linear Initial-Value Ordinary Differential Equations with Constant Coefficients
- Harry Bateman. 1910. Solution of a system of differential equations occurring in the theory of radioactive transformation. Proceedings of the Cambridge Philsophical Society 15 (1910), 423--427Google Scholar
- Jerzy Cetnar. 2006. General solution of Bateman equations for nuclear transmutations. Annals of Nuclear Energy 33 (2006), 640--645.Google ScholarCross Ref
- A. G. Groff. 1983. ORIGEN2: A versatile computer code for calculating the nuclide compositions and characteristics of nuclear materials. Nuclear Technology 62 (September 1983), 335--352.Google Scholar
- Maple. 2013. Maple 17 User Manual. Maplesoft division of Waterloo Maple Inc., Waterloo, ON Canada.Google Scholar
- Cleve Moler and Charles Van Loan. 2003. Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later. SIAM Rev. 45, 1 (2003), 3--49.Google ScholarDigital Library
Index Terms
Algorithm 982: Explicit Solutions of Triangular Systems of First-Order Linear Initial-Value Ordinary Differential Equations with Constant Coefficients
Recommendations
Remark on Algorithm 982: Explicit Solutions of Triangular Systems of First-order Linear Initial-value Ordinary Differential Equations with Constant Coefficients
Algorithm 982: Explicit solutions of triangular systems of first-order linear initial-value ordinary differential equations with constant coefficients provides an explicit solution for an homogeneous system, and a brief description of how to compute a ...
A new integration algorithm for ordinary differential equations based on continued fraction approximations
A new integration algorithm is found, and an implementation is compared with other programmed algorithms. The new algorithm is a step-by-step procedure for solving the initial value problem in ordinary differential equations. It is designed to ...
A simple approach to the particular solution of constant coefficient ordinary differential equations
AbstractAn eigenfunction approach to compute the particular solution of constant coefficient ordinary differential equations is presented. It is shown that the exponentials are the eigenfunctions of such equations. Solutions corresponding to ...
Comments