A New Embedded Phase-Fitted Modified Runge-Kutta Method for the Numerical Solution of Oscillatory Problems
|
|
- Joseph Lawrence
- 5 years ago
- Views:
Transcription
1 Applied Mathematical Sciences, Vol. 1, 16, no. 44, HIKARI Ltd, A New Embedded Phase-Fitted Modified Runge-Kutta Method for the Numerical Solution of Oscillatory Problems F. A. Fawzi Department of Mathematics, Universiti Putra Malaysia, 434 UPM Serdang Selangor, Malaysia & Department of Mathematics Faculty of Computer Science and Mathematics Tikrit University, Iraq N. Senu, F. Ismail and Z. A. Majid Department of Mathematics and Institute for Mathematical Research Universiti Putra Malaysia, 434 UPM Serdang Selangor, Malaysia Copyright c 16 F. A. Fawzi, N. Senu, F. Ismail and Z. A. Majid. 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 In this work, a new way for constructing an efficiently modified Runge-Kutta (RK) method to solve first-order ordinary differential equations with oscillatory solutions is provided. The proposed method solves the first-order ODEs by first converting the second order ODEs to an equivalent first-order ODEs. The method of the embedded has algebraic orders five and four. The numerical results of the new method have been compared with those of existing methods and showed that the new method is more efficient. Mathematics Subject Classification: 65L5, 65L
2 158 F. A. Fawzi, N. Senu, F. Ismail and Z. A. Majid Keywords: Embedded Runge Kutta Methods; Variable Step-size; Oscillatory Problems 1 Introduction The last few decades, a lot of research has been done in the area of numerical solution of first-order initial value problems of the form: where y (x) = f(x, y), y(x ) = y, y (x ) = y, x [a, b] (1) y(x) = [y 1 (x), y (x),..., y s (x)] T f(x, y) = [f 1 (x, y), f (x, y),..., f s (x, y)] T y is a given vector of initial condition and their solution is oscillating. This type of problem appears often in scientific areas such as astronomy, quantum mechanics, mechanics, and electronics. A lot of researches have been conducted for the construction of numerical methods that are both accurate and computationally efficient to solve the oscillatory problems. Moreover, they are a lot of numerical integration of adapting kind which is based on significant properties like exponentially-fitted, trigonometrically-fitted and phase-fitted [1-4]. Van de Vyver [5] presented an embedded phase-fitted modified RK method for the numerical solution of the Schrödinger equation. Fang et al. [6] constructed new embedded pair explicit RK methods with (FSAL) properties. Recently, Fang et al. [7] presented a new phase-fitted modified RK pair for the numerical solution of the radial Schrödinger equation. Simos and Vigo Aguiar [8] derived a modified RK method with phase-lag for the numerical solution of the Schrödinger equation. Also, they were concerned with phase-fitted in [9] based on the fifth - order method of Dormand and Prince. We used the phase-fitted technique to determine g. In [8, 9] the constant step size mode for phase-fitted modified Runge-Kutta methods have been used. Our purpose is to extend this technique for the development of variable step size algorithm. The aim of this paper is to develop a new efficient embedded phase-fitted modified Runge-Kutta method. In Section the preliminaries on the phase properties of the explicit modified RK method is discussed. In Section 3 we construct the new pair of explicit modified Runge-Kutta method and the stability and error analysis are also provided. The efficiency of the new method, when compared with the existing method, is given in Section 4.
3 A new embedded phase-fitted modified Runge-Kutta method 159 Phase Properties of Modified Runge-Kutta Method An explicit m-stage modified Runge-Kutta method formula is given by y n+1 = y n + h m b i f(x n + c i h, Y i ), i=1 i 1 Y i = g i y n + h a ij f(x n + c i h, Y j ), i = 1,..., m. () j=1 where a ij, c i, b i, i = 1,..., m are constant, h is the step size and the parameters g i, i = 1,..., m are even function of v = hw. The method is determined by means of Butcher tableau (see Table 1) [1] Table 1: The modified Runge-Kutta method c g a 1 c 3 g 3 a 31 a c s g s a s1 a s... a s,s 1 b 1 b... b s 1 b s The method is stated to be explicit when i j where a ij = and implicit otherwise. In RK method the embedded pair q(p) is based on the RK method (c, A, b) of order q and another RK method (c, A, b ) of order p < q. An embedded pair is characterized by Butcher tableau c g A An embedded pair of explicit modified Runge-Kutta method is used in a variable step size algorithm because they provide a cheap error estimate. b T b T
4 16 F. A. Fawzi, N. Senu, F. Ismail and Z. A. Majid From embedded method we obtain an estimated error. EST n+1 = y n+1 y n+1 (3) For the numerical integration of the equation (1) we used stepsize control procedure which is given by [11]: if EST n+1 < T ol 1, h n+1 = h n, if T ol 1 EST n+1 < T ol, h n+1 = h n, if EST n+1 T ol, h n+1 = hn and repeat the step. where T ol is the requested local error. It should be noted that the q th-order approximation y n is used as the initial value for the (n+1) th step, that mean the embedded pair is applied in local extrapolation mode or higher order mode. To develop the new method we utilize the test equation with analytical solution y = iwy, w ɛ R. (4) y(x n ) = y e iwx. (5) Based on reasons described by Houwen and Sommeijer [1], we shall confine our considerations to the homogeneous phase-lag. Applying the method () to the test equation (4) yields y n+1 = R(iH)y n, H = wh. (6) Definition.1 An explicit m-stage modified RK, (presented in Table 1) the quantities: φ(h) = H arg[r(ih)], α(h) = 1 R(iH) are called the dispersion or phase error or phase-lag and the amplification error respectively. If φ(h) = O(H q+1 ), and α(h) = O(H r+1 ) then the method is said to be phase-lag order q and dissipative order r. For explicit modified Runge-Kutta method R(iH) = A m (H ) + ihb m (H ), (7) where A m and B m are polynomials in H. From equation (7) it follows that [ φ(h) = H tan H B ] m(h ), (8) A m (H ) α(h) = 1 A m(h ) + H B m(h ). (9)
5 A new embedded phase-fitted modified Runge-Kutta method Analysis of Stability An m-stage modified Runge-Kutta method () is applied to equation (4), we obtain y n+1 = y n + HBY, (1) where and From (11), we have Y = y n G + HAY (11) Y = [Y 1, Y,..., Y s ], G = [g 1, g,..., g s ] B = [b 1 b... b s ] T, H = hw (I HA)Y = y n G, Y = (I HA)y n G, Y n+1 = y n + HB(I HA) 1 y n G (1) substituting equation (1) into equation (1), we obtain: where y n+1 = R(H)y n, is the stability function of the method. R(H) = 1 + HB(1 HA) 1 G (13) Definition. A modified Runge-kutta method is said to be absolutely stable if R(H) < 1 3 Construction of the New Method We consider an embedded pair of explicit Runge-Kutta which denoted by Butcher tableau which is given in [1]:
6 16 F. A. Fawzi, N. Senu, F. Ismail and Z. A. Majid 1 5 Table : Butcher Tableau RK method In this study, our focus is to have modified RK method with the phase-lag of order infinity. Therefore, the coefficients of this method is depending on the product H = wh. To get zero-phase error (φ(h) = ) for the fifth order method then the relation (9) become [ ] Bm (H ) tan(h) = H (14) A m (H ) = H [( g ) H H + 1 ] H H4 1 (15) H Solving equation (15) for g, we obtain g = 1 15H 4 tan(h) + 3 tan(h)h6 1 tan(h)h tan(h)h +5H 5 4H 3 + 4H (16) As H, g have the following Taylor expansions: g = H H H H H H (17)
7 A new embedded phase-fitted modified Runge-Kutta method 163 The same case for the fourth-order method tan(h) = ( H4 + ( 11g ) H + 19g ) H6 + ( 7 g ) H4 + ( 41 g ) 683 H. (18) g 3 = 9 tan(h) + 79 tan(h)h tan(h)h tan(h)h +... H(31H 3 tan(h) 41 tan(h)h + H 95) (19) As H, g 3 have the following Taylor expansions: g 3 = H H H H H H H () This new method is denoted as PHRK5(4). The stability function of the higher order of the new method is given by R(H) = 1 + H + 1 H H H H H H H H H H (1) It is not possible to find the stability region for infinity series of coefficients. Here, three different finite degrees of coefficients are considered, which are order four, six and eight. we obtained the stability interval to the new method of three different finite degrees of coefficients is ( 3.5, ) for order four, (.8, ) for order six and (.5, ) for order eight. If we continue to cut up to order ten the stability interval is (.4, ) that mean when we increase the order of coefficients the stability interval will be smaller. For the low order of the new method, we follow the same procedure then the stability function is given by R(H) = 1 + H + 1 H H H H H H H H H1
8 164 F. A. Fawzi, N. Senu, F. Ismail and Z. A. Majid H H H H () The stability region of the new method is obtained by equating the stability function R(H) with e Iθ and then solve for H using maple package. i.e R(H) = e iθ = cos(θ) + i sin(θ). The stability region for the high order method is shown in Figure 1. Figure 1: The stability region for the fifth order of the PHRK5(4) method for different order of the coefficients
9 A new embedded phase-fitted modified Runge-Kutta method 165 The stability region for the low order method is shown in Figure. Figure : The stability region for the fourth order of the PHRK5(4) method for different order of the coefficients The local truncation error analysis (LTE) of the new method based on the Taylor series expansion of the differences y n+1 and y(x n = h) will be considered The LTE of the high order is given by LT E = y n+1 y(x n + h) (3) LT E = h 6 [ 1 1 (f xf y f xyy y + f x f xy f yy y f y f xyy y ) 1 8 (f yf xy f yy y () f y f xyy y () ) (f xf y f yyy y () f x f xy f xy f xxy y f yy f xyy y () +f x f y f xxy f xy f xyy y () f yy f xyy y (3) f x f yyy () + f xf y f yy ) 1 7 (f yyf y f xx y +f xy f xxx f y f xxx f 3 y f xx f 4 y f xx f 4 y f x + f 5 y y + f xy f yyy y (3) + f y f xy f xx +f yy f xxx y + f yy f yyy y (4) ) 1 36 (f y f yy y f 3 y f xy y ) + f yy f xxx y + f yy f yyy y (4) )
10 166 F. A. Fawzi, N. Senu, F. Ismail and Z. A. Majid 1 36 (f y f yy y fy 3 f xy y ) f yfxyy 1 18 (f xfy f xy fy f yyy y (3) )] + O(h 7 ) (4) From equation (4), it is clear that the order of the new method is five because all the terms of h lower than h 6 are vanished. and the LTE for the low order: LT E = h 5 11 [ 7 (f xxxx+f yyyy y (4) ) f yf xxy y (f xf xyy y +f xy f yy y () ) 11 9 f xf y f xy (f y f xy y + f yy f xx y + f xxxy + y + f xyyy y (3) + f xx f xy +f yyy (3) ) f yf xyy y () (f xf yyy y () + f xxyy y () + f x f xxy ) f y f yy y () f xf y f yy y f xyy f xf yy 11 3 (f yf xxx + fy f xx ) (f y 3 f x + fy 4 y )] + O(h 6 ) (5) From equation (5), it is clear that the order of the new method is four because all the terms of h lower than h 5 are vanished. 4 Numerical Experiments In this section, we will apply the new method to solve different problems. The following explicit MRK method are selected for the numerical comparison. 4.1 Comparison with fixed step-size methods RK5: The phase-fitted fifth-order derived in this paper. PLRK4: The modified RK method derived by Simos and Vigo-Aguiar [8]. RK5B: Fifth-order six-stage RK method derived by Sakas and Simos [13]. RK4M: Fourth-order five-stage RK method given in Butcher [1]. RK4Z: Fourth-order five-stage RK method given in Hairer et. al [14].
11 A new embedded phase-fitted modified Runge-Kutta method 167 Problem 1: ( Homogeneous) [15] y 1 = y, y 1 (x) = 1, y = 64 y 1, y (x) = Theoretical solution : y 1 (x) = 1 sin(8x) + cos(8x) 4 y (x) = cos(8x) 8 sin(8x) Problem : ( Inhomogeneous) [1] y 1 = y, y 1 () = 1 y = v y 1 + (v 1) sin(x), y () = v + 1 Estimated frequency: v = 1 Theoretical solution : y 1 (x) = cos(vx) + sin(vx) + sin(x) y (x) = v sin(vx) + v cos(vx) + cos(x) Problem 3: ( Periodic orbit system) [3] y 1 = y, y 1 (x) = 1, y (x) = y = y cos(x) y 3 = y 4, y 3 (x) =, y 4 (x) =.9995 y 4 = y sin(x) Theoretical solution: y 1 (x) = cos(x) +.5x sin(x) y (x) = sin(x).5x cos(x) y 3 (x) = sin(x) +.5x cos(x) y 4 (x) = cos(x) +.5x sin(x)
12 168 F. A. Fawzi, N. Senu, F. Ismail and Z. A. Majid Problem 4: ( Nonlinear system) [17] y 1 = y 3, y 1 () = 1 y 3 = y 1 ( ) 3, y 3 () = y 1 + y y = y 4, y () = y 4 = y ( ) 3, y 4 () = 1 y 1 + y Theoretical solution : y 1 (x) = cos(x) y (x) = sin(x) y 3 (x) = sin(x) y 4 (x) = cos(x) Problem 5: ( Inhomogeneous) [18] y 1 = y, y 1 () = 1 y = y 1 + x, y () = Theoretical solution: y 1 = sin(x) + cos(x) + x y = cos(x) sin(x) + 1
13 A new embedded phase-fitted modified Runge-Kutta method log1(maxerr) 1 3 RK5 RK5B 4 PLRK4 RK4M RK4Z log 1 (STEP-SIZE) Figure 3: Comparison for RK5, RK5B, PLRK4, RK4M and RK4Z problem 1 with b=1 1 log1(maxerr) 1 RK5 3 RK5B PLRK4 RK4M RK4Z log 1 (STEP-SIZE) Figure 4: Comparison for RK5, RK5B, PLRK4, RK4M and RK4Z problem with b=1
14 17 F. A. Fawzi, N. Senu, F. Ismail and Z. A. Majid 1 RK5 RK5B PLRK4 RK4M RK4Z log1(maxerr) log 1 (STEP-SIZE) Figure 5: Comparison for RK5, RK5B, PLRK4, RK4M and RK4Z problem 3 with b=1 1 log1(maxerr) RK5 RK5B 6 PLRK4 RK4M RK4Z log 1 (STEP-SIZE) Figure 6: Comparison for RK5, RK5B, PLRK4, RK4M and RK4Z problem 4 with b=1
15 A new embedded phase-fitted modified Runge-Kutta method 171 log1(maxerr) RK5 RK5B PLRK4 RK4M RK4Z log 1 (STEP-SIZE) Figure 7: Comparison for RK5, RK5B, PLRK4, RK4M and RK4Z problem 5 with b=1
16 17 F. A. Fawzi, N. Senu, F. Ismail and Z. A. Majid 4. Comparison with variable step-size methods The methods we choose for comparison are as follows: PHMRKD5(4): The new pair of phase-fitted modified RK method presented in this paper. MODPHARK5: The phase-fitted modified RK pair given by Yonglei et al. [7]. MODRK54: The modified RK pair given by Van de Vyver [5]. RK5(4)F: The RK pair of Fehlbreg given in Butcher [1]. RK5(4)D: The RK pair of Dormand given by Hairer et al. [14]. Problem 6: ( Inhomogeneous system) y 1 = y 3, y 1 () = 1.1 y = y 4, y () = 1 Theoretical solution : y 3 = 16y e 1x + 1e 1x, y 3 () = 1 y 4 = 16y + 16e 1x + 1e 1x, y 4 () = 9.6 Source : Moo. et al. [17]. y 1 =.1 cos(4x) + e 1x y =.1 sin(4x) + e 1x y 3 =.4 sin(4x) 1e 1x y 4 =.4 cos(4x) 1e 1x Problem 7: ( Nonlinear system) Theoretical solution : y 1 = y 3, y 1 () = 1 y = y 4, y () = y 3 = 4x y y 1, y 3 () = y 1 + y y 4 = 4x y Source : Yonglei et al. [19]. y 1, y 4 () = y 1 + y3 y 1 (x) = cos(x ) y (x) = sin(x ) y 3 (x) = x sin(x ) y 4 (x) = x cos(x )
17 A new embedded phase-fitted modified Runge-Kutta method 173 log1(maxerr) PHMRKD5(4) RK5(4)D RK5(4)F MODRK54 MODPHARK5(4) log 1 (Function Evaluations) Figure 8: Accuracy curves for embedded 5(4) of MRK methods for Problem1 1 log1(maxerr) PHMRKD5(4) RK5(4)D RK5(4)F MODRK54 MODPHARK5(4) log 1 (Function Evaluations) Figure 9: Accuracy curves for embedded 5(4) of MRK methods for Problem
18 174 F. A. Fawzi, N. Senu, F. Ismail and Z. A. Majid log1(maxerr) PHMRKD5(4) RK5(4)D RK5(4)F MODRK54 MODPHARK5(4) log 1 (Function Evaluations) Figure 1: Accuracy curves for embedded 5(4) of MRK methods for Problem4 log1(maxerr) PHMRKD5(4) 1 RK5(4)D RK5(4)F MODRK54 MODPHARK5(4) log 1 (Function Evaluations) Figure 11: Accuracy curves for embedded 5(4) of MRK methods for Problem5
19 A new embedded phase-fitted modified Runge-Kutta method log1(maxerr) PHMRKD5(4) RK5(4)D RK5(4)F MODRK54 MODPHARK5(4) log 1 (Function Evaluations) 4 Figure 1: Accuracy curves for embedded 5(4) of MRK methods for Problem6 1 log1(maxerr) PHMRKD5(4) RK5(4)D RK5(4)F MODRK54 MODPHARK5(4) log 1 (Function Evaluations) Figure 13: Accuracy curves for embedded 5(4) of MRK methods for Problem7
20 176 F. A. Fawzi, N. Senu, F. Ismail and Z. A. Majid For comparison purpose, in analyzing the numerical results, methods with the same order will be compared. The results are plotted in Figures 3, 4, 5, 6 and 7 we present efficiency curves where the common logarithm of the maximum global error throughout the integration versus computational cost measured by h step-size. From the Figures 8-13, we present efficiency curves where the common logarithm of the maximum global error throughout the integration versus computational cost measured by number of function evaluation and we observed that the new PHMRKD5(4) method is more efficient for integration first-order differential equations possessing an oscillatory solution compared with other methods which are RK5(4)D, RK5(4)F, MODRK54 and MOD- PHARK5(4). 5 Conclusion In this paper, we derived a new embedded phase-fitted RK method for solving first-order ordinary differential equations with oscillatory problems. The methods of the embedded have algebraic orders five and four. The phase-fitting technique we used is an extension of the idea from [8]. At first, we introduced for fixed step-size then numerical findings indicate that the proposed method is more efficient than the existing method. Lastly, we used variable step-size become clear that the new embedded is more robustness than the existing method for solving oscillatory problems. References [1] T. E. Simos, Exponentially fitted Runge-Kutta methods for the numerical solution of the Schrödinger equation and related problems, Comput. Mater. Sci., 18 (), [] T. E. Simos, An exponentially-fitted Runge-Kutta method for the numerical integration of initial-value problems with periodic or oscillating solutions, Comput. Phys. Commun., 115 (1998), [3] G. Vanden Berghe, H. De Meyer, M. Van Daele, T. Van Hecke, Exponentially fitted explicit Runge-Kutta methods, Comput. Phys. Commun., 13 (1999), [4] Y. Fang, X. Wu, A trigonometrically fitted explicit Numerov-type method for second-order initial value problems with oscillating solutions, Appl.
21 A new embedded phase-fitted modified Runge-Kutta method 177 Numer. Math., 58 (8), [5] H. Van de Vyver, An embedded 5(4)pair of modified explicit Runge-Kutta methods for the numerical solution of the Schrödinger equation, International Journal of Modern Physics C, 16 (5), [6] Y. Fang, Y. Song and X. Wu, New embedded pairs of explicit Runge- Kutta methods with FSAL properties adapted to the numerical integration of oscillatory problems, Phys. Lett. A, 37 (8), [7] Y. Fang, X. You and Q. Ming, A new phase-fitted modified Runge-Kutta pair for the numerical solution of the radial Schrödinger equation, Appl. Math. Comp., 4 (13), [8] T. E. Simos and J. Vigo Aguiar, A modified Runge-Kutta method with phase-lag of order infinity for the numerical solution of the Schrödinger equation and related problems, Comput. Chem., 5 (1), [9] T. E. Simos and J. Vigo Aguiar, A modified phase-fitted Runge-Kutta method for the numerical solution of the Schrödinger equation, Journal of Mathematical Chemistry, 3 (1), [1] J. C. Butcher, Numerical Methods for Ordinary Differential Equations, John Wiley & Sons Ltd, Great Britain, 8. [11] A. D. Raptis and J. R. Cash, A variable step method for the numerical integration of the one-dimensional Schrödinger equation, J. Comput. Phys. Commum., 36 (1985), [1] P. J. van der Houwen and B. P. Sommeijer, Explicit Runge-Kutta (- Nyström) methods with reduced phase errors for computing oscillating solutions, SIAM Journal on Numerical Analysis, 4 (1987), [13] D. P. Sakas and T. E. Simos, A Fifth Algebraic Order Trigonometrically- Fitted Modified Runge-Kutta Zonneveld Method for the Numerical Solution of Orbital problems, Math. Comp. Mod., 4 (5),
22 178 F. A. Fawzi, N. Senu, F. Ismail and Z. A. Majid [14] E. Hairer, S. P. Nørsett, G. Wanner, Solving Ordinary Differential Equations I, Nonstiff Problems, Springer, Berlin, [15] M. M. Chawla, P. S. Rao, High-accuracy P-stable methods for y = f(t, y), IMA Journal of Numerical Analysis, 5 (1985), [16] D. G. Bettis, Runge-Kutta algorithms for oscillatory problems, Zeitschrift für Angewandte Mathematik und Physik, 3 (1979), [17] K. W. Moo, N. Senu, F. Ismail and M. Suleiman, New phase-fitted and amplification-fitted fourth-order and fifth-order Runge-Kutta-Nyström methods for oscillatory problems, Abstract and Applied Analysis, 13 (13), [18] R. C. Allen, Jr. and G. M. Wing, An invariant imbedding algorithm for the solution of inhomogeneous linear two-point boundary value problems, J. Computational Physics, 14 (1974), [19] Y. Fang, X. You and Q. Ming, Trigonometrically fitted two-derivative Runge-Kutta method for solving oscillatory differential equations, Numer. Algor., 65 (14), [] J. Vigo-Aguiar and T. E Simos, Review of multistep methods for the numerical solution of the radial Schrödinger equation, Int. J. Quant. Chem., 13 (5), [1] B. Paternoster, Runge-Kutta(-Nyström) methods for ODEs with periodic solutions based on trigonometric polynomials, Applied Numerical Mathematics, 8 (1998), [] Y. Zhang, H. Che, Y. Fang and X. You, A new trigonometrically fitted two-derivative Runge-Kutta method for the numerical solution of the Schrödinger equation and related problems, Journal of Applied Mathematics, 13 (13), [3] Firas Adel, Norazak Senu, Fudziah Ismail and Zanariah Abd. Majid, New Efficient Phase-Fitted and Amplification-Fitted Runge-Kutta Method for Oscillatory Problems, International Journal of Pure and Applied Mathematics, 17 (16), Received: April 7, 16; Published: July 6, 16
Embedded Explicit Two-Step Runge-Kutta-Nyström for Solving Special Second-Order IVPs
Applied Mathematical Sciences, Vol. 8, 14, no. 94, 4655 4673 HIKARI Ltd, www.mhikari.com http://dx.doi.org/1.1988/ams.14.4643 Embedded Explicit TwoStep RungeKuttaNyström for Solving Special SecondOrder
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 informationA NEW INTEGRATOR FOR SPECIAL THIRD ORDER DIFFERENTIAL EQUATIONS WITH APPLICATION TO THIN FILM FLOW PROBLEM
Indian J. Pure Appl. Math., 491): 151-167, March 218 c Indian National Science Academy DOI: 1.17/s13226-18-259-6 A NEW INTEGRATOR FOR SPECIAL THIRD ORDER DIFFERENTIAL EQUATIONS WITH APPLICATION TO THIN
More informationExponentially-Fitted Runge-Kutta Nystrom Method of Order Three for Solving Oscillatory Problems ABSTRACT 1. INTRODUCTION
Malaysian Journal of Mathematical Sciences 8(S): 7-4 (04) Special Issue: International Conference on Mathematical Sciences and Statistics 0 (ICMSS0) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage:
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 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 informationA Modified Runge-Kutta-Nyström Method by using Phase Lag Properties for the Numerical Solution of Orbital Problems
Appl. Math. Inf. Sci., No., - () Applied Mathematics & Information Sciences An International Journal A Modified Runge-Kutta-Nyström Method by using Phase Lag Properties for the Numerical Solution of Orbital
More informationAdebayo O. Adeniran1, Saheed O. Akindeinde2 and Babatunde S. Ogundare2 * Contents. 1. Introduction
Malaya Journal of Matematik Vol. No. 73-73 8 https://doi.org/.37/mjm/ An accurate five-step trigonometrically-fitted numerical scheme for approximating solutions of second order ordinary differential equations
More informationImproved Extended Runge-Kutta-like Method for Solving First Order IVPs
Indian Journal of Science and Technology, Vol 9(48), DOI: 0.7485/ijst/06/v9i48/0938, December 06 ISSN (Print) : 0974-6846 ISSN (Online) : 0974-5645 Improved Extended Runge-Kutta-like Method for Solving
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 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 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 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 informationResearch Article Diagonally Implicit Block Backward Differentiation Formulas for Solving Ordinary Differential Equations
International Mathematics and Mathematical Sciences Volume 212, Article ID 767328, 8 pages doi:1.1155/212/767328 Research Article Diagonally Implicit Block Backward Differentiation Formulas for Solving
More informationSymplectic Runge-Kutta-Nyström Methods with Phase-Lag Oder 8 and Infinity
Appl. Math. Inf. Sci. 9, No. 3, 1105-1112 (2015) 1105 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.12785/amis/090302 Symplectic Runge-Kutta-Nyström Methods with
More informationAn Embedded Fourth Order Method for Solving Structurally Partitioned Systems of Ordinary Differential Equations
Applied Mathematical Sciences, Vol. 9, 0, no. 97, 484-48 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/ams.0.98 An Embedded Fourth Order Method for Solving Structurally Partitioned Systems of Ordinary
More informationFifth-Order Improved Runge-Kutta Method With Reduced Number of Function Evaluations
Australian Journal of Basic and Applied Sciences, 6(3): 9-5, 22 ISSN 99-88 Fifth-Order Improved Runge-Kutta Method With Reduced Number of Function Evaluations Faranak Rabiei, Fudziah Ismail Department
More informationOn the reliability and stability of direct explicit Runge-Kutta integrators
Global Journal of Pure and Applied Mathematics. ISSN 973-78 Volume, Number 4 (), pp. 3959-3975 Research India Publications http://www.ripublication.com/gjpam.htm On the reliability and stability of direct
More informationExponentially Fitted Symplectic Runge-Kutta-Nyström methods
Appl. Math. Inf. Sci. 7, No. 1, 81-85 (2013) 81 Applied Mathematics & Information Sciences An International Journal Exponentially Fitted Symplectic Runge-Kutta-Nyström methods Th. Monovasilis 1, Z. Kalogiratou
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 informationDirect Explicit Integrators of RK Type for Solving Special Fourth-Order Ordinary Differential Equations With An Application
Global Journal of Pure and Applied Mathematics. ISSN 0973-768 Volume, Number 6 (06), pp. 4687 475 Research India Publications http://www.ripublication.com/gjpam.htm Direct Explicit Integrators of RK Type
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 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 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 informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
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 informationOne-Step Hybrid Block Method with One Generalized Off-Step Points for Direct Solution of Second Order Ordinary Differential Equations
Applied Mathematical Sciences, Vol. 10, 2016, no. 29, 142-142 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.6127 One-Step Hybrid Block Method with One Generalized Off-Step Points for
More informationSains Malaysiana 47(9)(2018):
Sains Malaysiana 47(9)(2018): 2223 2230 http://dx.doi.org/10.17576/jsm-2018-4709-33 Block Hybrid Method with Trigonometric-Fitting for Solving Oscillatory Problems (Kaedah Blok Hibrid dengan Penyuaian-Trigonometri
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 informationTrigonometrically tted predictor corrector methods for IVPs with oscillating solutions
Journal of Computational and Applied Mathematics 58 (23) 35 44 www.elsevier.com/locate/cam Trigonometrically tted predictor corrector methods for IVPs with oscillating solutions G.Psihoyios a, T.E. Simos
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 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 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 informationAn Optimized Symmetric 8-Step Semi-Embedded Predictor-Corrector Method for IVPs with Oscillating Solutions
Appl. Math. Inf. Sci., No., - () Applied Mathematics & Information Sciences An International Journal An Optimized Symmetric -Step Semi-Embedded Predictor-Corrector Method for IVPs with Oscillating Solutions
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 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 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 Two-step Iterative Method Free from Derivative for Solving Nonlinear Equations
Applied Mathematical Sciences, Vol. 8, 2014, no. 161, 8021-8027 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.49710 A Two-step Iterative Method Free from Derivative for Solving Nonlinear
More informationDeferred correction based on exponentially-fitted mono-implicit Runge-Kutta methods
Deferred correction based on exponentially-fitted mono-implicit Runge-Kutta methods M. Van Daele Department of Applied Mathematics, Computer Science and Statistics Ghent University 3th International Conference
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 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 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 informationEXAMPLE OF ONE-STEP METHOD
EXAMPLE OF ONE-STEP METHOD Consider solving y = y cos x, y(0) = 1 Imagine writing a Taylor series for the solution Y (x), say initially about x = 0. Then Y (h) = Y (0) + hy (0) + h2 2 Y (0) + h3 6 Y (0)
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 informationResearch Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
Applied Mathematics Volume 013, Article ID 591636, 5 pages http://dx.doi.org/10.1155/013/591636 Research Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
More informationOn the classification of certain curves up to projective tranformations
On the classification of certain curves up to projective tranformations Mehdi Nadjafikhah Abstract The purpose of this paper is to classify the curves in the form y 3 = c 3 x 3 +c 2 x 2 + c 1x + c 0, with
More informationAn Implementation of QSAOR Iterative Method for Non-Homogeneous Helmholtz Equations
International Journal of Contemporary Mathematical Sciences Vol. 11, 2016, no. 2, 85-96 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijcms.2016.51047 An Implementation of QSAOR Iterative Method
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 informationA Family of Optimal Multipoint Root-Finding Methods Based on the Interpolating Polynomials
Applied Mathematical Sciences, Vol. 8, 2014, no. 35, 1723-1730 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.4127 A Family of Optimal Multipoint Root-Finding Methods Based on the Interpolating
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 informationP-stable Exponentially fitted Obrechkoff methods
P-stable Exponentially fitted Obrechkoff methods Marnix Van Daele, G. Vanden Berghe Marnix.VanDaele@UGent.be Vakgroep Toegepaste Wiskunde en Informatica Universiteit Gent SciCADE 7 p. /4 Outline Introduction
More informationAnother Sixth-Order Iterative Method Free from Derivative for Solving Multiple Roots of a Nonlinear Equation
Applied Mathematical Sciences, Vol. 11, 2017, no. 43, 2121-2129 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.76208 Another Sixth-Order Iterative Method Free from Derivative for Solving
More informationConstruction and Implementation of Optimal 8 Step Linear Multistep
Page 138 ORIGINAL RESEARCH Construction and Implementation of Optimal 8 Step Linear Multistep method Bakre Omolara Fatimah, Awe Gbemisola Sikirat and Akanbi Moses Adebowale Department of Mathematics, Lagos
More informationExplicit One-Step Methods
Chapter 1 Explicit One-Step Methods Remark 11 Contents This class presents methods for the numerical solution of explicit systems of initial value problems for ordinary differential equations of first
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationSYMMETRY REDUCTION AND NUMERICAL SOLUTION OF A THIRD-ORDER ODE FROM THIN FILM FLOW
Mathematical and Computational Applications,Vol. 15, No. 4, pp. 709-719, 2010. c Association for Scientific Research SYMMETRY REDUCTION AND NUMERICAL SOLUTION OF A THIRD-ORDER ODE FROM THIN FILM FLOW E.
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 informationAdaptation of Taylor s Formula for Solving System of Differential Equations
Nonlinear Analysis and Differential Equations, Vol. 4, 2016, no. 2, 95-107 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/nade.2016.51144 Adaptation of Taylor s Formula for Solving System of Differential
More informationExtended Runge Kutta-like formulae
Applied Numerical Mathematics 56 006 584 605 www.elsevier.com/locate/apnum Extended Runge Kutta-like formulae Xinyuan Wu a Jianlin Xia b a State Key Laboratory for Novel Software Technology Department
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 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 informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
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 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 informationSymmetry Reduction of Chazy Equation
Applied Mathematical Sciences, Vol 8, 2014, no 70, 3449-3459 HIKARI Ltd, wwwm-hikaricom http://dxdoiorg/1012988/ams201443208 Symmetry Reduction of Chazy Equation Figen AÇIL KİRAZ Department of Mathematics,
More informationButcher tableau Can summarize an s + 1 stage Runge Kutta method using a triangular grid of coefficients
AM 205: lecture 13 Last time: ODE convergence and stability, Runge Kutta methods Today: the Butcher tableau, multi-step methods, boundary value problems Butcher tableau Can summarize an s + 1 stage Runge
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 informationDownloaded 12/15/15 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 40, No. 4, pp. 56 537 c 2002 Society for Industrial and Applied Mathematics PREDICTOR-CORRECTOR METHODS OF RUNGE KUTTA TYPE FOR STOCHASTIC DIFFERENTIAL EQUATIONS KEVIN BURRAGE
More informationApplication of the perturbation iteration method to boundary layer type problems
DOI 10.1186/s40064-016-1859-4 RESEARCH Open Access Application of the perturbation iteration method to boundary layer type problems Mehmet Pakdemirli * *Correspondence: mpak@cbu.edu.tr Applied Mathematics
More informationOn a Diophantine Equation 1
International Journal of Contemporary Mathematical Sciences Vol. 12, 2017, no. 2, 73-81 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ijcms.2017.728 On a Diophantine Equation 1 Xin Zhang Department
More informationSolving Zhou Chaotic System Using Fourth-Order Runge-Kutta Method
World Applied Sciences Journal 21 (6): 939-944, 2013 ISSN 11-4952 IDOSI Publications, 2013 DOI: 10.529/idosi.wasj.2013.21.6.2915 Solving Zhou Chaotic System Using Fourth-Order Runge-Kutta Method 1 1 3
More informationExact Solutions for the Nonlinear (2+1)-Dimensional Davey-Stewartson Equation Using the Generalized ( G. )-Expansion Method
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Published by Canadian Center of Science and Education Exact Solutions for the Nonlinear +-Dimensional Davey-Stewartson Equation
More informationThe Solution of Weakly Nonlinear Oscillatory Problems with No Damping Using MAPLE
World Applied Sciences Journal (): 64-69, 0 ISSN 88-495 IDOSI Publications, 0 DOI: 0.589/idosi.wasj.0..0.09 The Solution of Weakly Nonlinear Oscillatory Problems with No Damping Using MAPLE N. Hashim,
More informationMini project ODE, TANA22
Mini project ODE, TANA22 Filip Berglund (filbe882) Linh Nguyen (linng299) Amanda Åkesson (amaak531) October 2018 1 1 Introduction Carl David Tohmé Runge (1856 1927) was a German mathematician and a prominent
More informationODE Runge-Kutta methods
ODE Runge-Kutta methods The theory (very short excerpts from lectures) First-order initial value problem We want to approximate the solution Y(x) of a system of first-order ordinary differential equations
More informationOn parallel solution of ordinary differential equations
On parallel solution of ordinary differential equations arxiv:1601.02245v1 [math.na] 10 Jan 2016 Alejandra Gaitán Montejo 1, Octavio A. Michel-Manzo 2 and César A. Terrero-Escalante 2 1 Department of Mathematics,
More informationExponential Integrators
Exponential Integrators John C. Bowman and Malcolm Roberts (University of Alberta) June 11, 2009 www.math.ualberta.ca/ bowman/talks 1 Outline Exponential Integrators Exponential Euler History Generalizations
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationExponentially Fitted Error Correction Methods for Solving Initial Value Problems
KYUNGPOOK Math. J. 52(2012), 167-177 http://dx.doi.org/10.5666/kmj.2012.52.2.167 Exponentially Fitted Error Correction Methods for Solving Initial Value Problems Sangdong Kim and Philsu Kim Department
More information-Stable Second Derivative Block Multistep Formula for Stiff Initial Value Problems
IAENG International Journal of Applied Mathematics, :3, IJAM 3_7 -Stable Second Derivative Bloc Multistep Formula for Stiff Initial Value Problems (Advance online publication: 3 August ) IAENG International
More informationSolution for a non-homogeneous Klein-Gordon Equation with 5th Degree Polynomial Forcing Function
Advanced Studies in Theoretical Physics Vol., 207, no. 2, 679-685 HIKARI Ltd, www.m-hikari.com https://doi.org/0.2988/astp.207.7052 Solution for a non-homogeneous Klein-Gordon Equation with 5th Degree
More informationSolution of First Order Initial Value Problem by Sixth Order Predictor Corrector Method
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 6 (2017), pp. 2277 2290 Research India Publications http://www.ripublication.com/gjpam.htm Solution of First Order Initial
More informationA Compact Fourth-Order Finite Difference Scheme for Unsteady Viscous Incompressible Flows
Journal of Scientific Computing, Vol. 16, No. 1, March 2001 ( 2001) A Compact Fourth-Order Finite Difference Scheme for Unsteady Viscous Incompressible Flows Ming Li 1 and Tao Tang 2 Received January 23,
More informationGraded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo
2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical
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 informationYou may not use your books, notes; calculators are highly recommended.
Math 301 Winter 2013-14 Midterm 1 02/06/2014 Time Limit: 60 Minutes Name (Print): Instructor This exam contains 8 pages (including this cover page) and 6 problems. Check to see if any pages are missing.
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 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 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 informationResearch Article An Accurate Block Hybrid Collocation Method for Third Order Ordinary Differential Equations
Applied Mathematics, Article ID 549597, 9 pages http://dx.doi.org/1.1155/14/549597 Research Article An Accurate Bloc Hybrid Collocation Method for Third Order Ordinary Differential Equations Lee Ken Yap,
More informationKrylov Implicit Integration Factor Methods for Semilinear Fourth-Order Equations
mathematics Article Krylov Implicit Integration Factor Methods for Semilinear Fourth-Order Equations Michael Machen and Yong-Tao Zhang * Department of Applied and Computational Mathematics and Statistics,
More informationThe Generalized Viscosity Implicit Rules of Asymptotically Nonexpansive Mappings in Hilbert Spaces
Applied Mathematical Sciences, Vol. 11, 2017, no. 12, 549-560 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.718 The Generalized Viscosity Implicit Rules of Asymptotically Nonexpansive
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 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 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 informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationA Note on Product Range of 3-by-3 Normal Matrices
International Mathematical Forum, Vol. 11, 2016, no. 18, 885-891 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2016.6796 A Note on Product Range of 3-by-3 Normal Matrices Peng-Ruei Huang
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 informationParallel general linear methods for stiff ordinary differential and differential algebraic equations
~-~H-~-~ ~ APPLIED ~ NUMERICAL MATHEMATICS ELSEVIER Applied Numerical Mathematics 17 (1995) 213-222 Parallel general linear methods for stiff ordinary differential and differential algebraic equations
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More information