scholarly journals Solution of Higher-Order ODEs Using Backward Difference Method

2011 ◽  
Vol 2011 ◽  
pp. 1-18 ◽  
Author(s):  
Mohamed Bin Suleiman ◽  
Zarina Bibi Binti Ibrahim ◽  
Ahmad Fadly Nurullah Bin Rasedee

The current numerical technique for solving a system of higher-order ordinary differential equations (ODEs) is to reduce it to a system of first-order equations then solving it using first-order ODE methods. Here, we propose a method to solve higher-order ODEs directly. The formulae will be derived in terms of backward difference in a constant stepsize formulation. The method developed will be validated by solving some higher-order ODEs directly with constant stepsize. To simplify the evaluations of the integration coefficients, we find the relationship between various orders. The result presented confirmed our hypothesis.

2014 ◽  
Vol 2014 ◽  
pp. 1-13 ◽  
Author(s):  
Hazizah Mohd Ijam ◽  
Mohamed Suleiman ◽  
Ahmad Fadly Nurullah Rasedee ◽  
Norazak Senu ◽  
Ali Ahmadian ◽  
...  

We describe the development of a 2-point block backward difference method (2PBBD) for solving system of nonstiff higher-order ordinary differential equations (ODEs) directly. The method computes the approximate solutions at two points simultaneously within an equidistant block. The integration coefficients that are used in the method are obtained only once at the start of the integration. Numerical results are presented to compare the performances of the method developed with 1-point backward difference method (1PBD) and 2-point block divided difference method (2PBDD). The result indicated that, for finer step sizes, this method performs better than the other two methods, that is, 1PBD and 2PBDD.


Open Physics ◽  
2020 ◽  
Vol 18 (1) ◽  
pp. 738-750
Author(s):  
Saima Akram ◽  
Allah Nawaz ◽  
Thabet Abdeljawad ◽  
Abdul Ghaffar ◽  
Kottakkaran Sooppy Nisar

AbstractThis article concerns with the development of the number of focal values. We analyzed periodic solutions for first-order cubic non-autonomous ordinary differential equations. Bifurcation analysis for periodic solutions from a fine focus {\mathfrak{z}}=0 is also examined. In particular, we are interested to detect the maximum number of periodic solutions for various classes of higher order in which a given solution can bifurcate under perturbation of the coefficients. We calculate the maximum number of periodic solutions for different classes, namely, {C}_{10,5} and {C}_{12,6} with trigonometric coefficients, and they are found with nine and eight multiplicities at most. The classes {C}_{8,3} and {C}_{8,4} with algebraic coefficients have at most eight limit cycles. The new formula {\varkappa }_{10} is developed by which we succeeded to find highest known multiplicity ten for class {C}_{\mathrm{9,3}} with polynomial coefficient. Periodicity is calculated for both trigonometric and algebraic coefficients. Few examples are also considered to explain the applicability and stability of the methods presented.


2014 ◽  
Vol 2014 ◽  
pp. 1-10 ◽  
Author(s):  
Ahmad Fadly Nurullah Rasedee ◽  
Mohamed bin Suleiman ◽  
Zarina Bibi Ibrahim

The current numerical techniques for solving a system of higher order ordinary differential equations (ODEs) directly calculate the integration coefficients at every step. Here, we propose a method to solve higher order ODEs directly by calculating the integration coefficients only once at the beginning of the integration and if required once more at the end. The formulae will be derived in terms of backward difference in a constant step size formulation. The method developed will be validated by solving some higher order ODEs directly using variable order step size. To simplify the evaluations of the integration coefficients, we find the relationship between various orders. The results presented confirmed our hypothesis.


2006 ◽  
Vol 49 (2) ◽  
pp. 170-184
Author(s):  
Richard Atkins

AbstractThis paper investigates the relationship between a system of differential equations and the underlying geometry associated with it. The geometry of a surface determines shortest paths, or geodesics connecting nearby points, which are defined as the solutions to a pair of second-order differential equations: the Euler–Lagrange equations of the metric. We ask when the converse holds, that is, when solutions to a system of differential equations reveals an underlying geometry. Specifically, when may the solutions to a given pair of second order ordinary differential equations d2y1/dt2 = f (y, ẏ, t) and d2y2/dt2 = g(y, ẏ, t) be reparameterized by t → T(y, t) so as to give locally the geodesics of a Euclidean space? Our approach is based upon Cartan's method of equivalence. In the second part of the paper, the equivalence problem is solved for a generic pair of second order ordinary differential equations of the above form revealing the existence of 24 invariant functions.


1972 ◽  
Vol 39 (3) ◽  
pp. 689-695 ◽  
Author(s):  
W. W. Recker

The two-dimensional equations of magnetoelastodynamics are considered as a symmetric hyperbolic system of linear first-order partial-differential equations in three independent variables. The characteristic properties of the system are determined and a numerical method for obtaining the solution to mixed initial and boundary-value problems in plane magnetoelastodynamics is presented. Results on the von Neumann necessary condition are presented. Application of the method to a problem which has a known solution provides further numerical evidence of the convergence and stability of the method.


Sign in / Sign up

Export Citation Format

Share Document