An Embedded Fourth Order Method for Solving Structurally Partitioned Systems of Ordinary Differential Equations
|
|
- Alan Hood
- 6 years ago
- Views:
Transcription
1 Applied Mathematical Sciences, Vol. 9, 0, no. 97, HIKARI Ltd, An Embedded Fourth Order Method for Solving Structurally Partitioned Systems of Ordinary Differential Equations Igor V. Olemskoy and Alexey S. Eremin Faculty of Applied Mathematics and Control Processes Saint-Petersburg State University Universitetskii prospekt, Petergof, Saint-Petersburg, Russia, 9804 Copyright 0 Igor V. Olemskoy and Alexey S. Eremin. This article is distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract An embedded one-step numerical method for structurally partitioned systems of ordinary differential equations (ODEs) is considered. Twoparametric families of methods of order four with automatic step-control are constructed for systems of ODEs of first and second order. The methods have fewer stages than classic Runge Kutta methods. Mathematics Subject Classification: L0, L0 Keywords: initial value problem, embedded methods, partitioned systems, control term Introduction In the paper [] there was presented the algorithm of partitioning a system of ordinary differential equations (ODEs) z s = ϕ s (x, z,..., z M ), s =,..., M, () and the way of its reducing to the canonical form of structurally partitioned systems y 0 = f 0 (x, y 0,..., y n ), () y i = f i (x, y 0,..., y i, y l+,..., y n ), i =,..., l, () y j = f j (x, y 0,..., y j ), j = l +,..., n, (4)
2 4844 Igor V. Olemskoy and Alexey S. Eremin x [X 0, X k ] R, y s : [X 0, X k ] R rs, s = 0,..., n, n f 0 : [X 0, X k ] R M R r 0, r s = M, s=0 f i : [X 0, X k ] R M ri R r i, r i = f j : [X 0, X k ] R M ˆrj R r j, ˆr j = l r s, i =,..., l, s=i n r s, j = l +,..., n. The groups of equations () and (4) are structurally identical. Equations in them are ordered in such way that in every equation the right-hand side depends only on the previous functions of the same group. The group () contains all the equations that cannot be placed into () or (4), so it can be called general group. Sometimes it can be missing, as well as the group (). In works [, 9, 0] a generalisation of explicit Runge Kutta methods (RKs) for structurally partitioned systems (SPS) of ODEs () (4) was suggested. The effectiveness of the methods presented there is provided by the fact, that the number of stages m for the general group () is equal to that for classic RKs (m = Q, q 4), but for the groups () and (4), taking in account their special structure, it can be reduced while keeping the method s order. Moreover, when the general group is absent the number of necessary stages can be reduced even more significantly. Thus, a fifth order method for a system with only groups () and (4) was constructed with four stages [9] (while classic RKs have at least six). The considered generalisation of RKs is based on usage of the system s special structure in the algorithm. It is natural to apply the idea to embedded RKs. Dormand Prince type methods [] are widely used nowadays. Combining their approach with the First Same As Last (FSAL) idea the last stage of a step is the first of the next step [, 4] fast explicit one-step methods with stable step-size control were constructed. For systems with groups () and (4) an embedded Dormand Prince type method of order five with automatic step-size control was constructed with five stages []. s=j The integration method In the paper an explicit embedded fourth order method with automatic stepsize control for systems { y i = f i (x, y,..., y i, y l+,..., y n ), i =,..., l, y j () = f j (x, y,..., y j ), j = l +,..., n, n
3 An embedded fourth order method for SPS of ODEs 484 is constructed. Such systems appear, for example, in problems of optimal control and stabilisation [7, 8], medical modelling [], celestial mechanics or high-energy physics. For instance, large-scale oscillations of systems with rotational symmetry [] are described (after some reducing) with equations x x = K, z z = L, (x + z) (x + z) K = x, L = z. (x + z) (x + z) Here x is a dimensionless moment of inertia about the rotational axis; z is a dimensionless moment of inertia about the equatorial plane; K is a kinetic energy of motions parallel to the equatorial plane; and L is a kinetic energy of vertical motions. Change of variables ξ = x, ξ = z, ξ = x, ξ 4 = z, ξ = K, ξ = L reduces the system to the form (), which is an SPS without the general group. The approximated solution of () is found as m s y s (x + h) z s = y s (x) + b sv k sv (h), s =,..., n, (7) v= m i m = m =... = m l m l+ = m =... = m n m j, where k sv k sv (h) are calculated in the strict order as k iv = hf i (x + c iv h, y (x) + k, k,..., k n, k, k,..., v a ivw k w,..., y i (x) + w= w= v v a i,v,i,w k i,w, ) v y l+ (x) + a i,v,l+,w k l+,w,..., y n (x) + a ivnw k nw, (8) w= k jv = hf j (x + c jv h, y (x) + w= v a jvw k w,..., y j (x) + w= () ) v a j,v,j,w k j,w, w= c i = 0, i =,..., l, c j 0, j = l +,..., n. The method will be referred as structural. For fast and effective practical implementation of numerical methods for ODEs one needs a step-size control algorithm. Here the embedded Dormand Prince type method of order four, corresponding to the form (7) (8), is constructed on base of three-stage fourth order methods from []. In Dormand Prince methods two approximations of different orders are found simultaneously; the higher order result is used in further computations, and the lower
4 484 Igor V. Olemskoy and Alexey S. Eremin order result is used only for the step-size control (technically, it is not the error being estimated, but the last term equal in the Taylor series of the solution and the approximation). In [] a number of important relations between the method s parameters was obtained. They combine the parameters a spvw, b sp and c sp into quite few groups and can be referred as group equalities: c v =... = c lv =: c iv, c l+,v =... = c nv =: c jv, b v =... = b lv =: b iv, b l+,v =... = b nv =: b jv, a ivw =... = a i,v,i,w =: a ivîw, a i,v,l+,w =... = a ivnw =: a ivjw, a jvw =... = a jvlw =: a jviw, a j,v,l+,w =... = a j,v,j,w =: a jvĵw. (9) Here and further the indices i, j, î, ĵ of the parameters c iv, c jv, b iv, b jv, a ivîη, a ivjw, a jviw, a jvĵw correspond to the group attribute, not to a certain number. Namely, the index at first place connects parameters to the system () structural groups i to the first and j to the second; and the third index shows which argument of f i or f j in (8) the parameter corresponds i and î up to l-th function, j and ĵ from l + to n. The use of restrictions (9) significantly simplifies the design of order conditions and methods construction and realisation. We use the slightly modified Butcher tableau to present the embedded structural method (see table ). The control term estimation is done with the same k iv and k jv values as the result but with d iv and d jv parameters m s z s = y s (x) + d sv k sv. (0) The parameters of an m = (m i, m j )-stage embedded method (7) (8) of order p and estimator order q must provide that y s (x + h) z s = O(h p+ ) and y s (x + h) z s = O(h q+ ) for all s =,..., n. We consider methods where p > q. The difference z s z s is used for the step-size control. The method s (7) (8), (0) abbreviation RKSp(q)(m i, m j )F or RKSp(q)mF tells about its type RK (Runge Kutta); its generalisation to systems with special structure () S; its main order p and the estimator order (q); number of stages m if m i = m j = m or (m i, m j ); and about the FSAL idea used F. We construct fourth order methods with third order estimators. FSAL idea is also used. v=
5 An embedded fourth order method for SPS of ODEs 4847 Table : Butcher tableaux for an embedded structural method c is a isîv a isjv b is d is b i d i c i a iî a iî a ij b i d i c i a iî a iî a iî a ij a ij b i d i c imi a imi î a imi î... a imi îm i a imi j a imi j... a imi jm i b imi d imi c js a jsiv a jsîv b js d js c j a ji a jĵ b j d j c j a ji a ji a jĵ a jĵ b j d j c j a ji a ji a ji a jĵ a jĵ a jĵ b j d j c jmj a jmj i a jmj i... a jmj im j a jmj ĵ a jmj ĵ... a jmj ĵm j b jmj d jmj New methods will be compared to the classic four-stages Runge Kutta method ( The Runge Kutta method) with five-stages third-order estimator, named RK4()T in []; and to the same Runge Kutta method but with an embedded second-order estimator (named in Russian mathematical tradition Egorov control term after Prof. Vsevolod A. Egorov []) RK4(). Order conditions The order conditions for the fourth order method (7) (8), (0) with use of (9) and standard simplifying conditions w a swve = c sw () e= form the system of nonlinear algebraic equations: w= b sw c v sw =, v = 0,...,, () v +
6 4848 Igor V. Olemskoy and Alexey S. Eremin b iw c v iw w= w a iwîr c ir = r= w b iw c v iw a iwjr c jr = w= r=, v = 0,, () ( + v), v = 0,, (4) ( + v) b iw w w= r= a iwîr c ir =, b iw w w= r= a iwjr c jr = () w b iw a iwîr w= r= e= w r b iw a iwjr w= r= e= r a irîe c ie = 4, a jrie c ie = 4, b iw w r a iwîr w= r= e= w b iw a iwjr w= r= e= a irje c je = 4, () r a jrĵe c je = 4, (7) b jw c v jw w= b jw c v jw w= w a jwir c ir = r= w a jwĵr c jr = r=, v = 0,, (8) ( + v), v = 0,, (9) ( + v) b jw w w= r= w b jw a jwir w= r= e= w r b jw a jwĵr w= r= e= a jwir c ir =, r a irîe c ie = 4, a jrie c ie = 4, b jw w w= r= b jw a jwĵr c jr =, (0) w r a jwir w= r= e= w r b jw a jwĵr w= r= e= a irje c je = 4, () a jrĵe c je = 4, () 4 w= d sw c v sw =, v = 0,...,, () v + 4 w= 4 d iw d jw w r= w w= r= a iwîr c ir =, a jwir c ir =, 4 w= 4 a iwjr c jr =, (4) d iw w d jw r= w w= r= a jwĵr c ir =. ()
7 An embedded fourth order method for SPS of ODEs 4849 Table : RKS4()4F c is a isîv a isjv b is d is ξ + ξ 8 ξ 0 0 ξ c js a jsiv a jsĵv b js d js η η 0 η 0 0 η The system consists of three mutually connected blocks of equations. First are the conditions of the fourth order 8 equations () (). Second gives the third order of the estimator 0 equations () (). The last block () consists of 8 restrictions to provide the FSAL idea for the method (7) (8), (0): a i4iw = a j4iw = b iw, a i4jw = a j4jw = b jw, w =,...,, c i4 =, c j4 = + σ, a i4i4 = a i4j4 = 0, a j4i4 = a j4j4 = σ 0. () The number of unknown parameters is. Theorem. There exists an embedded RKS4()4F method (7) (8), (0) that satisfies conditions (9), (). The theorem is proved by solving the system () (). The fourth order method is determined uniquely; the third order estimator exist only when σ = / and has two free parameters d i4 0 and d j4 0. The table contains the parameters of RKS4()4F. The constructed method has better number of stages to order ratio than RK4()T and RK4(). Often the implementation of Runge Kutta methods is done so, that the step-size is changed every step. If so, k j4 values of RKS4()4F cannot be used as k j of the next step. It reduces the effectiveness of the method, but still it excels both RK4()T and RK4().
8 480 Igor V. Olemskoy and Alexey S. Eremin There is an option to control the step-size on base of error estimation of only first l functions (which can work fine when two groups correspond to two different physical values or are scaled differently []). The fourth stage for the second group of () can be omitted then. Such method RKS4(4)(4,)F (Table ) in fact takes only three stages per step, while RK4()T and RK4() both take four. Table : RKS4()(4,)F c is a isîv a isjv b is d is ξ + ξ 8 ξ 0 0 ξ c js a jsiv a jsĵv b js Testing The constructed methods RKS4()4F and RKS4()(4,)F with parameters ξ = η = were compared to RK4()T and RK4() solving the problem x = x + y, x(0) = x (0) =, y = x y, y(0) =, y (0) =, for t [0, π]. The general solution of the problem is x(t) = C cos(αt) + C sin(αt) + C cos(βt) + C 4 sin(βt), y(t) = C sin(αt) + C cos(αt) C sin(βt) + C 4 cos(βt), α = 4, β = +. 4 (7)
9 An embedded fourth order method for SPS of ODEs 48 Table 4: Methods comparison for the problem (7) lg( δ ) RKS4()4F RKS4()(4,)F RK4()T RK4() N h N g + N f N h N g + N f N h N g + N f N h N g + N f Since the practical effectiveness of Runge Kutta methods depends on the algorithm of step-size choosing, and all the methods we compare have the same order properties, all were implemented within a common algorithmic driver. The new step-size was determined as h new = c h(tol/err) /(q+) with the guarantee factor c = 0.9. A series of computations with global error δ varying from 0 to 0. was held. All four methods passed the convergence test, but it should be noted, that the step-size choosing algorithm was more sensitive to the solution behaviour in case of the new methods RKS4()4F, RKS4()(4,)F. In the table 4 number of necessary steps N h and the number of right-hand side evaluations N g + N f (N g in the first equation and N f in the second) over the whole integration interval are shown for the desired norm of the global error. The results confirm the theoretical expectations. For the same global error RKS4()4F and RKS4()(4,)F takes less steps and less right-hand side evaluations than RK4()T and RK4() do. And vice versa for the same amount of computations RK4()T and RK4() have larger global error. References [] O. B. Arushanyan, S. F. Zaletkin, Chislennoe resheniye differentsialnyh uravnenii na Fortrane [Numerical Solution of Ordinary Diffential Equations with Fortran], Moscow, Moscow State University, 990 (in Russian).
10 48 Igor V. Olemskoy and Alexey S. Eremin [] J. R. Dormand, P. J. Prince, New Runge Kutta algorithms for numerical simulation in dynamical astronomy, Celestial Mech., 8 (978),. [] J. R. Dormand, P. J. Prince, A family of embedded Runge Kutta formulae, J. Comput. Appl. Math., (980), 9. [4] J. R. Dormand, M. E. A. El-Mikkawy, P. J. Prince, Families of Runge Kutta Nyström formulae, J. Numer. Analys., 7 (987), 0. [] E. Hairer, S. P. Nørsett, G. Wanner. Solving Ordinary Differential Equations I: Nonstiff Problems, ed, Springer-Verlag Berlin Heidelberg, 008. [] S. A. Kutuzov, I. V. Olemskoy, L. P. Ossipkov, V. N. Starkov, Matematicheskiye metody issledovaniya kosmicheskih sistem [Mathematical Methods of Studying Space Systems], Saint-Petersburg, Saint-Petersburg State University, 00 (in Russian). [7] A. N. Kvitko, On a control problem, Differential Equations, (004), [8] A. N. Kvitko, A method for solving boundary value problems for nonlinear control systems in the class of discrete controls, Differential Equations, (008), [9] I. V. Olemskoy, Fifth-order four-stage method for numerical integration of special systems, Comp. Math. and Math. Phys., 8 (00), 4. [0] I. V. Olemskoy, Structural approach to the design of explicit one-stage methods, Comp. Math. and Math. Phys., 7 (00), [] I. V. Olemskoy, A fifth-order five-stage embedded method of the Dormand Prince type, Comp. Math. and Math. Phys., 7 (00), 8 9. [] I. V. Olemskoy, Modifikatsiya algoritma vydeleniya strukturnyh osobennostei [The modification of the algorithm for determining structural peculiarities], Vestnik Sankt-Peterburgskogo Universiteta, Series 0, (00), 4 4 (in Russian). [] I. V. Zhukova, E. P. Kolpak, Yu. E. Balykina, Mathematical model of growing tumor, Applied mathematical sciences, 0 (04), Received: June, 0; Published: July, 0
A New Embedded Phase-Fitted Modified Runge-Kutta Method for the Numerical Solution of Oscillatory Problems
Applied Mathematical Sciences, Vol. 1, 16, no. 44, 157-178 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.1988/ams.16.64146 A New Embedded Phase-Fitted Modified Runge-Kutta Method for the Numerical Solution
More informationRemark on the Sensitivity of Simulated Solutions of the Nonlinear Dynamical System to the Used Numerical Method
International Journal of Mathematical Analysis Vol. 9, 2015, no. 55, 2749-2754 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2015.59236 Remark on the Sensitivity of Simulated Solutions of
More informationarxiv: v1 [math.na] 31 Oct 2016
RKFD Methods - a short review Maciej Jaromin November, 206 arxiv:60.09739v [math.na] 3 Oct 206 Abstract In this paper, a recently published method [Hussain, Ismail, Senua, Solving directly special fourthorder
More informationInvestigation on the Most Efficient Ways to Solve the Implicit Equations for Gauss Methods in the Constant Stepsize Setting
Applied Mathematical Sciences, Vol. 12, 2018, no. 2, 93-103 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.711340 Investigation on the Most Efficient Ways to Solve the Implicit Equations
More informationA Study on Linear and Nonlinear Stiff Problems. Using Single-Term Haar Wavelet Series Technique
Int. Journal of Math. Analysis, Vol. 7, 3, no. 53, 65-636 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.988/ijma.3.3894 A Study on Linear and Nonlinear Stiff Problems Using Single-Term Haar Wavelet Series
More informationA SYMBOLIC-NUMERIC APPROACH TO THE SOLUTION OF THE BUTCHER EQUATIONS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 17, Number 3, Fall 2009 A SYMBOLIC-NUMERIC APPROACH TO THE SOLUTION OF THE BUTCHER EQUATIONS SERGEY KHASHIN ABSTRACT. A new approach based on the use of new
More informationRunge Kutta Collocation Method for the Solution of First Order Ordinary Differential Equations
Nonlinear Analysis and Differential Equations, Vol. 4, 2016, no. 1, 17-26 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/nade.2016.5823 Runge Kutta Collocation Method for the Solution of First
More informationExact Solutions for a Fifth-Order Two-Mode KdV Equation with Variable Coefficients
Contemporary Engineering Sciences, Vol. 11, 2018, no. 16, 779-784 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ces.2018.8262 Exact Solutions for a Fifth-Order Two-Mode KdV Equation with Variable
More informationIMPLICIT INTERVAL MULTISTEP METHODS FOR SOLVING THE INITIAL VALUE PROBLEM
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 8 (1), 17-30 (2002) IMPLICIT INTERVAL MULTISTEP METHODS FOR SOLVING THE INITIAL VALUE PROBLEM MAŁGORZATA JANKOWSKA 1, ANDRZEJ MARCINIAK 1,2 1 Poznań University
More informationOn the Equation of Fourth Order with. Quadratic Nonlinearity
International Journal of Mathematical Analysis Vol. 9, 015, no. 5, 659-666 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ijma.015.5109 On the Equation of Fourth Order with Quadratic Nonlinearity
More informationDynamic Model of Space Robot Manipulator
Applied Mathematical Sciences, Vol. 9, 215, no. 94, 465-4659 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.56429 Dynamic Model of Space Robot Manipulator Polina Efimova Saint-Petersburg
More informationA NOTE ON EXPLICIT THREE-DERIVATIVE RUNGE-KUTTA METHODS (ThDRK)
BULLETIN OF THE INTERNATIONAL MATHEMATICAL VIRTUAL INSTITUTE ISSN 303-4874 (p), ISSN (o) 303-4955 www.imvibl.org / JOURNALS / BULLETIN Vol. 5(015), 65-7 Former BULLETIN OF THE SOCIETY OF MATHEMATICIANS
More informationPeriodic and Soliton Solutions for a Generalized Two-Mode KdV-Burger s Type Equation
Contemporary Engineering Sciences Vol. 11 2018 no. 16 785-791 HIKARI Ltd www.m-hikari.com https://doi.org/10.12988/ces.2018.8267 Periodic and Soliton Solutions for a Generalized Two-Mode KdV-Burger s Type
More informationFourth-order symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: a review
Fourth-order symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: a review G. Vanden Berghe and M. Van Daele Vakgroep Toegepaste Wiskunde en Informatica, Universiteit Gent, Krijgslaan
More informationSolving Homogeneous Systems with Sub-matrices
Pure Mathematical Sciences, Vol 7, 218, no 1, 11-18 HIKARI Ltd, wwwm-hikaricom https://doiorg/112988/pms218843 Solving Homogeneous Systems with Sub-matrices Massoud Malek Mathematics, California State
More informationConvex Sets Strict Separation in Hilbert Spaces
Applied Mathematical Sciences, Vol. 8, 2014, no. 64, 3155-3160 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.44257 Convex Sets Strict Separation in Hilbert Spaces M. A. M. Ferreira 1
More informationTwo Optimized Runge-Kutta Methods for the Solution of the Schrödinger Equation
MATCH Communications in Mathematical and in Computer Chemistry MATCH Commun. Math. Comput. Chem. (8) 7-77 ISSN - Two Optimized Runge-Kutta Methods for the Solution of the Schrödinger Equation T.V. Triantafyllidis,
More informationLinearization of Two Dimensional Complex-Linearizable Systems of Second Order Ordinary Differential Equations
Applied Mathematical Sciences, Vol. 9, 2015, no. 58, 2889-2900 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.4121002 Linearization of Two Dimensional Complex-Linearizable Systems of
More informationEmbedded 5(4) Pair Trigonometrically-Fitted Two Derivative Runge- Kutta Method with FSAL Property for Numerical Solution of Oscillatory Problems
Embedded 5(4) Pair Trigonometrically-Fitted Two Derivative Runge- Kutta Method with FSAL Property for Numerical Solution of Oscillatory Problems N. SENU, N. A. AHMAD, F. ISMAIL & N. BACHOK Institute for
More informationRecurrence Relations between Symmetric Polynomials of n-th Order
Applied Mathematical Sciences, Vol. 8, 214, no. 15, 5195-522 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.214.47525 Recurrence Relations between Symmetric Polynomials of n-th Order Yuriy
More informationRiesz Representation Theorem on Generalized n-inner Product Spaces
Int. Journal of Math. Analysis, Vol. 7, 2013, no. 18, 873-882 HIKARI Ltd, www.m-hikari.com Riesz Representation Theorem on Generalized n-inner Product Spaces Pudji Astuti Faculty of Mathematics and Natural
More informationSymmetric Properties for the (h, q)-tangent Polynomials
Adv. Studies Theor. Phys., Vol. 8, 04, no. 6, 59-65 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/astp.04.43 Symmetric Properties for the h, q-tangent Polynomials C. S. Ryoo Department of Mathematics
More informationSymplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: revisited
Symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: revisited G. Vanden Berghe and M. Van Daele Vakgroep Toegepaste Wiskunde en Informatica, Universiteit Gent, Krijgslaan 8
More informationOn Positive Stable Realization for Continuous Linear Singular Systems
Int. Journal of Math. Analysis, Vol. 8, 2014, no. 8, 395-400 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.4246 On Positive Stable Realization for Continuous Linear Singular Systems
More informationFour Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential Equations
International journal of scientific and technical research in engineering (IJSTRE) www.ijstre.com Volume Issue ǁ July 206. Four Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential
More informationEfficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators
MATEMATIKA, 8, Volume 3, Number, c Penerbit UTM Press. All rights reserved Efficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators Annie Gorgey and Nor Azian Aini Mat Department of Mathematics,
More informationKKM-Type Theorems for Best Proximal Points in Normed Linear Space
International Journal of Mathematical Analysis Vol. 12, 2018, no. 12, 603-609 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2018.81069 KKM-Type Theorems for Best Proximal Points in Normed
More informationStability Analysis and Numerical Solution for. the Fractional Order Biochemical Reaction Model
Nonlinear Analysis and Differential Equations, Vol. 4, 16, no. 11, 51-53 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.1988/nade.16.6531 Stability Analysis and Numerical Solution for the Fractional
More informationQuadratic SDIRK pair for treating chemical reaction problems.
Quadratic SDIRK pair for treating chemical reaction problems. Ch. Tsitouras TEI of Chalkis, Dept. of Applied Sciences, GR 34400 Psahna, Greece. I. Th. Famelis TEI of Athens, Dept. of Mathematics, GR 12210
More informationof a Two-Operator Product 1
Applied Mathematical Sciences, Vol. 7, 2013, no. 130, 6465-6474 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2013.39501 Reverse Order Law for {1, 3}-Inverse of a Two-Operator Product 1 XUE
More informationZ. Omar. Department of Mathematics School of Quantitative Sciences College of Art and Sciences Univeristi Utara Malaysia, Malaysia. Ra ft.
International Journal of Mathematical Analysis Vol. 9, 015, no. 46, 57-7 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ijma.015.57181 Developing a Single Step Hybrid Block Method with Generalized
More informationNote About a Combinatorial Sum
Int. J. Contemp. Math. Sciences, Vol. 8, 203, no. 8, 349-353 HIKARI Ltd, www.m-hiari.com Note About a Combinatorial Sum Laurenţiu Modan Spiru Haret University, Academy of Economic Studies Department of
More informationExplicit Expressions for Free Components of. Sums of the Same Powers
Applied Mathematical Sciences, Vol., 27, no. 53, 2639-2645 HIKARI Ltd, www.m-hikari.com https://doi.org/.2988/ams.27.79276 Explicit Expressions for Free Components of Sums of the Same Powers Alexander
More informationPredictor-Corrector Finite-Difference Lattice Boltzmann Schemes
Applied Mathematical Sciences, Vol. 9, 2015, no. 84, 4191-4199 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.5138 Predictor-Corrector Finite-Difference Lattice Boltzmann Schemes G. V.
More informationAn Efficient Multiscale Runge-Kutta Galerkin Method for Generalized Burgers-Huxley Equation
Applied Mathematical Sciences, Vol. 11, 2017, no. 30, 1467-1479 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.7141 An Efficient Multiscale Runge-Kutta Galerkin Method for Generalized Burgers-Huxley
More informationarxiv: v1 [math.na] 15 Nov 2008
arxiv:0811.2481v1 [math.na] 15 Nov 2008 A Phase-Fitted Runge-Kutta-Nyström method for the Numerical Solution of Initial Value Problems with Oscillating Solutions Abstract D. F. Papadopoulos a, Z. A. Anastassi
More informationPoincaré`s Map in a Van der Pol Equation
International Journal of Mathematical Analysis Vol. 8, 014, no. 59, 939-943 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ijma.014.411338 Poincaré`s Map in a Van der Pol Equation Eduardo-Luis
More informationOn Symmetric Bi-Multipliers of Lattice Implication Algebras
International Mathematical Forum, Vol. 13, 2018, no. 7, 343-350 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/imf.2018.8423 On Symmetric Bi-Multipliers of Lattice Implication Algebras Kyung Ho
More informationVariational Theory of Solitons for a Higher Order Generalized Camassa-Holm Equation
International Journal of Mathematical Analysis Vol. 11, 2017, no. 21, 1007-1018 HIKAI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2017.710141 Variational Theory of Solitons for a Higher Order Generalized
More informationONE- AND TWO-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 5,53-65 (1999) ONE- AND TWO-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE ANDRZEJ MARCINIAK 1, BARBARA SZYSZKA 2 1 Institute of Computing Science,
More informationBoundary Value Problem for Second Order Ordinary Linear Differential Equations with Variable Coefficients
International Journal of Mathematical Analysis Vol. 9, 2015, no. 3, 111-116 HIKARI Ltd, www.m-hikari.com http://d.doi.org/10.12988/ijma.2015.411353 Boundary Value Problem for Second Order Ordinary Linear
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationOn Nonlinear Methods for Stiff and Singular First Order Initial Value Problems
Nonlinear Analysis and Differential Equations, Vol. 6, 08, no., 5-64 HIKARI Ltd, www.m-hikari.com https://doi.org/0.988/nade.08.8 On Nonlinear Methods for Stiff and Singular First Order Initial Value Problems
More informationSome Reviews on Ranks of Upper Triangular Block Matrices over a Skew Field
International Mathematical Forum, Vol 13, 2018, no 7, 323-335 HIKARI Ltd, wwwm-hikaricom https://doiorg/1012988/imf20188528 Some Reviews on Ranks of Upper Triangular lock Matrices over a Skew Field Netsai
More informationIdentities of Symmetry for Generalized Higher-Order q-euler Polynomials under S 3
Applied Mathematical Sciences, Vol. 8, 204, no. 3, 559-5597 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.2988/ams.204.4755 Identities of Symmetry for Generalized Higher-Order q-euler Polynomials under
More informationDynamical Analysis of a Harvested Predator-prey. Model with Ratio-dependent Response Function. and Prey Refuge
Applied Mathematical Sciences, Vol. 8, 214, no. 11, 527-537 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/12988/ams.214.4275 Dynamical Analysis of a Harvested Predator-prey Model with Ratio-dependent
More informationStarting Methods for Two-Step Runge Kutta Methods of Stage-Order 3 and Order 6
Cambridge International Science Publishing Cambridge CB1 6AZ Great Britain Journal of Computational Methods in Sciences and Engineering vol. 2, no. 3, 2, pp. 1 3 ISSN 1472 7978 Starting Methods for Two-Step
More informationGeneralized RK Integrators for Solving Ordinary Differential Equations: A Survey & Comparison Study
Global Journal of Pure and Applied Mathematics. ISSN 7-78 Volume, Number 7 (7), pp. Research India Publications http://www.ripublication.com/gjpam.htm Generalized RK Integrators for Solving Ordinary Differential
More informationSolution of Nonlinear Fractional Differential. Equations Using the Homotopy Perturbation. Sumudu Transform Method
Applied Mathematical Sciences, Vol. 8, 2014, no. 44, 2195-2210 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.4285 Solution of Nonlinear Fractional Differential Equations Using the Homotopy
More informationTHREE- AND FOUR-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 6, 41-59 (2000) THREE- AND FOUR-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE KAROL GAJDA 1, ANDRZEJ MARCINIAK 2, BARBARA SZYSZKA 1 1 Institute of
More informationRestrained Weakly Connected Independent Domination in the Corona and Composition of Graphs
Applied Mathematical Sciences, Vol. 9, 2015, no. 20, 973-978 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.4121046 Restrained Weakly Connected Independent Domination in the Corona and
More informationResearch Article P-Stable Higher Derivative Methods with Minimal Phase-Lag for Solving Second Order Differential Equations
Hindawi Publishing Corporation Journal of Applied Mathematics Volume 2011, Article ID 407151, 15 pages doi:10.1155/2011/407151 Research Article P-Stable Higher Derivative Methods with Minimal Phase-Lag
More informationOn a Certain Representation in the Pairs of Normed Spaces
Applied Mathematical Sciences, Vol. 12, 2018, no. 3, 115-119 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.712362 On a Certain Representation in the Pairs of ormed Spaces Ahiro Hoshida
More informationCaristi-type Fixed Point Theorem of Set-Valued Maps in Metric Spaces
International Journal of Mathematical Analysis Vol. 11, 2017, no. 6, 267-275 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2017.717 Caristi-type Fixed Point Theorem of Set-Valued Maps in Metric
More informationReducing round-off errors in symmetric multistep methods
Reducing round-off errors in symmetric multistep methods Paola Console a, Ernst Hairer a a Section de Mathématiques, Université de Genève, 2-4 rue du Lièvre, CH-1211 Genève 4, Switzerland. (Paola.Console@unige.ch,
More informationA Generalization of p-rings
International Journal of Algebra, Vol. 9, 2015, no. 8, 395-401 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ija.2015.5848 A Generalization of p-rings Adil Yaqub Department of Mathematics University
More informationInterval Images Recognition and Fuzzy Sets
International Mathematical Forum, Vol. 9, 2014, no. 19, 917-921 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2014.4490 Interval Images Recognition and Fuzzy Sets G. Sh. Tsitsiashvili, Yu.
More informationMetric Analysis Approach for Interpolation and Forecasting of Time Processes
Applied Mathematical Sciences, Vol. 8, 2014, no. 22, 1053-1060 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.312727 Metric Analysis Approach for Interpolation and Forecasting of Time
More informationSYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS
BIT 0006-3835/00/4004-0726 $15.00 2000, Vol. 40, No. 4, pp. 726 734 c Swets & Zeitlinger SYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS E. HAIRER Section de mathématiques, Université
More informationOn a 3-Uniform Path-Hypergraph on 5 Vertices
Applied Mathematical Sciences, Vol. 10, 2016, no. 30, 1489-1500 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.512742 On a 3-Uniform Path-Hypergraph on 5 Vertices Paola Bonacini Department
More informationGraceful Labeling for Complete Bipartite Graphs
Applied Mathematical Sciences, Vol. 8, 2014, no. 103, 5099-5104 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.46488 Graceful Labeling for Complete Bipartite Graphs V. J. Kaneria Department
More informationThe family of Runge Kutta methods with two intermediate evaluations is defined by
AM 205: lecture 13 Last time: Numerical solution of ordinary differential equations Today: Additional ODE methods, boundary value problems Thursday s lecture will be given by Thomas Fai Assignment 3 will
More informationSymmetric Identities of Generalized (h, q)-euler Polynomials under Third Dihedral Group
Applied Mathematical Sciences, vol. 8, 2014, no. 145, 7207-7212 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.49701 Symmetric Identities of Generalized (h, )-Euler Polynomials under
More informationStability Analysis of Plankton Ecosystem Model. Affected by Oxygen Deficit
Applied Mathematical Sciences Vol 9 2015 no 81 4043-4052 HIKARI Ltd wwwm-hikaricom http://dxdoiorg/1012988/ams201553255 Stability Analysis of Plankton Ecosystem Model Affected by Oxygen Deficit Yuriska
More informationModeling of Inflation of Liquid Spherical Layers. under Zero Gravity
Applied Mathematical Sciences Vol. 7 013 no. 138 6867-6876 HIKAI Ltd www.m-hikari.com http://dx.doi.org/10.1988/ams.013.310570 Modeling of Inflation of Liquid Spherical Layers under Zero Gravity G. I.
More informationA Note of the Strong Convergence of the Mann Iteration for Demicontractive Mappings
Applied Mathematical Sciences, Vol. 10, 2016, no. 6, 255-261 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.511700 A Note of the Strong Convergence of the Mann Iteration for Demicontractive
More informationDistribution Solutions of Some PDEs Related to the Wave Equation and the Diamond Operator
Applied Mathematical Sciences, Vol. 7, 013, no. 111, 5515-554 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ams.013.3844 Distribution Solutions of Some PDEs Related to the Wave Equation and the
More informationMorphisms Between the Groups of Semi Magic Squares and Real Numbers
International Journal of Algebra, Vol. 8, 2014, no. 19, 903-907 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ija.2014.212137 Morphisms Between the Groups of Semi Magic Squares and Real Numbers
More informationObservation Impact Assessment for Dynamic. Data-Driven Coupled Chaotic System
Applied Mathematical Sciences, Vol. 10, 2016, no. 45, 2239-2248 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.65170 Observation Impact Assessment for Dynamic Data-Driven Coupled Chaotic
More informationAn Application of Fibonacci Sequence on Continued Fractions
International Mathematical Forum, Vol. 0, 205, no. 2, 69-74 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.2988/imf.205.42207 An Application of Fibonacci Sequence on Continued Fractions Ali H. Hakami
More informationAn Approximate Solution for Volterra Integral Equations of the Second Kind in Space with Weight Function
International Journal of Mathematical Analysis Vol. 11 17 no. 18 849-861 HIKARI Ltd www.m-hikari.com https://doi.org/1.1988/ijma.17.771 An Approximate Solution for Volterra Integral Equations of the Second
More informationOptimal Homotopy Asymptotic Method for Solving Gardner Equation
Applied Mathematical Sciences, Vol. 9, 2015, no. 53, 2635-2644 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.52145 Optimal Homotopy Asymptotic Method for Solving Gardner Equation Jaharuddin
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationThe Law of Luminous Intensity Variation and Technical Vision
Adv. Studies Theor. Phys. Vol. 7 2013 no. 8 349-354 HIKARI Ltd www.m-hikari.com The Law of Luminous Intensity Variation and Technical Vision Anna Gorbenko Department of Intelligent Systems and Robotics
More informationExplicit Runge Kutta Pairs with Lower Stage-Order *
myjournal manuscript No. (will be inserted by the editor) Explicit Runge Kutta Pairs with Lower Stage-Order * J.H. Verner Dedicated to John C. Butcher in celebration of his eightieth birthday October,
More informationAdvanced Studies in Theoretical Physics Vol. 8, 2014, no. 22, HIKARI Ltd,
Advanced Studies in Theoretical Physics Vol. 8, 204, no. 22, 977-982 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.2988/astp.204.499 Some Identities of Symmetry for the Higher-order Carlitz Bernoulli
More informationReview Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester
HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE
More informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
More informationBackward error analysis
Backward error analysis Brynjulf Owren July 28, 2015 Introduction. The main source for these notes is the monograph by Hairer, Lubich and Wanner [2]. Consider ODEs in R d of the form ẏ = f(y), y(0) = y
More informationSecond Hankel Determinant Problem for a Certain Subclass of Univalent Functions
International Journal of Mathematical Analysis Vol. 9, 05, no. 0, 493-498 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/ijma.05.55 Second Hankel Determinant Problem for a Certain Subclass of Univalent
More informationOn Geometric Hyper-Structures 1
International Mathematical Forum, Vol. 9, 2014, no. 14, 651-659 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2014.312232 On Geometric Hyper-Structures 1 Mashhour I.M. Al Ali Bani-Ata, Fethi
More informationImproved Starting Methods for Two-Step Runge Kutta Methods of Stage-Order p 3
Improved Starting Methods for Two-Step Runge Kutta Methods of Stage-Order p 3 J.H. Verner August 3, 2004 Abstract. In [5], Jackiewicz and Verner derived formulas for, and tested the implementation of two-step
More informationNumerical Investigation of the Time Invariant Optimal Control of Singular Systems Using Adomian Decomposition Method
Applied Mathematical Sciences, Vol. 8, 24, no. 2, 6-68 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.2988/ams.24.4863 Numerical Investigation of the Time Invariant Optimal Control of Singular Systems
More informationDetermination of Young's Modulus by Using. Initial Data for Different Boundary Conditions
Applied Mathematical Sciences, Vol. 11, 017, no. 19, 913-93 HIKARI Ltd, www.m-hikari.com https://doi.org/10.1988/ams.017.7388 Determination of Young's Modulus by Using Initial Data for Different Boundary
More informationNonexistence of Limit Cycles in Rayleigh System
International Journal of Mathematical Analysis Vol. 8, 014, no. 49, 47-431 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ijma.014.4883 Nonexistence of Limit Cycles in Rayleigh System Sandro-Jose
More informationA New Block Method and Their Application to Numerical Solution of Ordinary Differential Equations
A New Block Method and Their Application to Numerical Solution of Ordinary Differential Equations Rei-Wei Song and Ming-Gong Lee* d09440@chu.edu.tw, mglee@chu.edu.tw * Department of Applied Mathematics/
More informationSecure Weakly Convex Domination in Graphs
Applied Mathematical Sciences, Vol 9, 2015, no 3, 143-147 HIKARI Ltd, wwwm-hikaricom http://dxdoiorg/1012988/ams2015411992 Secure Weakly Convex Domination in Graphs Rene E Leonida Mathematics Department
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationOn Numerical Solutions of Systems of. Ordinary Differential Equations. by Numerical-Analytical Method
Applied Mathematical Sciences, Vol. 8, 2014, no. 164, 8199-8207 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ams.2014.410807 On Numerical Solutions of Systems of Ordinary Differential Equations
More informationThe Modified Adomian Decomposition Method for. Solving Nonlinear Coupled Burger s Equations
Nonlinear Analysis and Differential Equations, Vol. 3, 015, no. 3, 111-1 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/nade.015.416 The Modified Adomian Decomposition Method for Solving Nonlinear
More informationOrder-theoretical Characterizations of Countably Approximating Posets 1
Int. J. Contemp. Math. Sciences, Vol. 9, 2014, no. 9, 447-454 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijcms.2014.4658 Order-theoretical Characterizations of Countably Approximating Posets
More informationDetermining the Rolle function in Lagrange interpolatory approximation
Determining the Rolle function in Lagrange interpolatory approimation arxiv:181.961v1 [math.na] Oct 18 J. S. C. Prentice Department of Pure and Applied Mathematics University of Johannesburg South Africa
More informationStatic Model of Decision-making over the Set of Coalitional Partitions
Applied Mathematical ciences, Vol. 8, 2014, no. 170, 8451-8457 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ams.2014.410889 tatic Model of Decision-maing over the et of Coalitional Partitions
More informationOn the satellite s electrodynamic attitude stabilization
On the satellite s electrodynamic attitude stabilization Alexander Yu. Aleksandrov Kirill A. Antipov Alexey A. Tikhonov aatikhonov@rambler.ru Abstract The paper deals with a satellite in a circular near-earth
More informationThe Improved Arithmetic-Geometric Mean Inequalities for Matrix Norms
Applied Mathematical Sciences, Vol 7, 03, no 9, 439-446 HIKARI Ltd, wwwm-hikaricom The Improved Arithmetic-Geometric Mean Inequalities for Matrix Norms I Halil Gumus Adıyaman University, Faculty of Arts
More informationImprovements in Newton-Rapshon Method for Nonlinear Equations Using Modified Adomian Decomposition Method
International Journal of Mathematical Analysis Vol. 9, 2015, no. 39, 1919-1928 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2015.54124 Improvements in Newton-Rapshon Method for Nonlinear
More informationThe Spin Angular Momentum. for the Celestial Objects
Adv. Studies Theor. Phys., Vol. 7, 2013, no. 11, 541-547 HIKARI Ltd, www.m-hikari.com The Spin Angular Momentum for the Celestial Objects Shubhen Biswas G.P.S.H.School, P.O.Alaipur, Pin.-741245(W.B), India
More informationStrong Stability of Singly-Diagonally-Implicit Runge-Kutta Methods
Strong Stability of Singly-Diagonally-Implicit Runge-Kutta Methods L. Ferracina and M. N. Spijker 2007, June 4 Abstract. This paper deals with the numerical solution of initial value problems, for systems
More informationLeft R-prime (R, S)-submodules
International Mathematical Forum, Vol. 8, 2013, no. 13, 619-626 HIKARI Ltd, www.m-hikari.com Left R-prime (R, S)-submodules T. Khumprapussorn Department of Mathematics, Faculty of Science King Mongkut
More informationExponential Integrators
Exponential Integrators John C. Bowman (University of Alberta) May 22, 2007 www.math.ualberta.ca/ bowman/talks 1 Exponential Integrators Outline Exponential Euler History Generalizations Stationary Green
More information