A fourth-order finite difference scheme for the numerical solution of 1D linear hyperbolic equation
|
|
- Willa McBride
- 5 years ago
- Views:
Transcription
1 203 (203) - Available online at Volume 203, Year 203 Article ID cna-0048, Pages doi:0.5899/203/cna-0048 Research Article A fourth-order finite difference scheme for the numerical solution of D linear hyperbolic equation Akbar Mohebbi Department of Applied Mathematics, Faculty of Mathematical Science, University of Kashan, Kashan, Iran Copyright 203 c Akbar Mohebbi. This is an open access article 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 paper, a high-order and unconditionally stable difference method is proposed for the numerical solution of onespace dimensional linear hyperbolic equation. We apply a compact finite difference approximation of fourth-order for discretizing spatial derivative of this equation and a Padé approximation of fifth-order for the resulting system of ordinary differential equations. It is shown through analysis that the proposed scheme is unconditionally stable. This new method is easy to implement, produces very accurate results and needs short CPU time. Some numerical examples are included to demonstrate the validity and applicability of the technique. We compare the numerical results of this paper with the numerical results of some methods in the literature. Keywords: Padé approximation, compact finite difference scheme, unconditionally stable, linear hyperbolic equation, telegraph equation, high accuracy. Introduction The hyperbolic partial differential equations model the vibrations of structures (e.g. buildings, beams and machines) and are the basis for fundamental equations of atomic physics. In this paper we consider the second order one-dimensional linear hyperbolic equation 2 u u (x,t) + 2α t2 t (x,t) + β 2 u(x,t) = 2 u (x,t) + f (x,t), x2 (x,t) (0,L) (0,T ], (.) with initial conditions u(x,0) = ϕ (x), u t (x,0) = ϕ 2(x), (.2) and boundary conditions u(0,t) = g 0 (t), Corresponding author. address: a mohebbi@kashanu.ac.ir, mohebbi.kashanu@gmail.com, Tel:
2 Page 2 of u(l,t) = g (t), t 0. (.3) For α > 0 and β = 0, Eq. (.) represents a damped wave equation and for α > β > 0, it is called Telegraph equation. We assume that α and β are positive. Eq. (.) referred to as second-order telegraph equation with constant coefficients, models mixture between diffusion and wave propagation by introducing a term that accounts for effects of finite velocity to standard heat or mass transport equation [5]. However, Eq. (.) is commonly used in signal analysis for transmission and propagation of electrical signals [24] and also has applications in other fields [6, 5, 26]. In recent years, much attention has been given in the literature to the development, analysis and implementation of stable methods for the numerical solution of second-order hyperbolic equations, see for example [4, 6, 22, 27]. These methods are conditionally stable. Mohanty [6, 7, 8] did nice investigations on the one space-dimensional hyperbolic equations. In [7], Mohanty carried over a new technique to solve the linear one-space-dimensional hyperbolic equation (.) which is unconditionally stable and is of second-order accurate in both time and space components. Also this author proposed in [8] a three level implicit unconditionally stable difference scheme of second-order accurate in both time and spatial components for the solution of Eq. (.) with variable coefficients that fictitious points are not needed at each time step along the boundary. Abbasbandy and Roohani in [] applied the well-known homotopy analysis method (HAM) and an interesting iterative algorithm for solving the telegraph equation with an integral condition. A numerical scheme is developed in [5] to solve the one-dimensional hyperbolic telegraph equation using collocation points and approximating the solution using a thin plate splines radial basis function. Another numerical method is presented in [6] to solve the one-dimensional hyperbolic telegraph equation using Chebyshev cardinal functions. Authors of [] applied the Rothe-Wavelet method to the solution of Eq. (.) and described in slightly different manner for the discretization of this equation. A numerical scheme based on the shifted Chebyshev tau method is proposed in [26] to solve Eq. (.). The combination of fourth order compact finite difference and collocation method for solving Eq. (.) is given in [20]. Authors of [7] used a numerical method based on the boundary integral equation and dual reciprocity method. These authors used time stepping scheme to deal with the time derivative and different types of radial basis functions to approximate functions in the dual reciprocity method. The variational iteration method has been applied successfully for the solution linear, variable coefficient, fractional derivative and multi space telegraph equations in [8]. In [5] a new spline difference scheme based on quadratic spline interpolations in space direction and finite difference discretization in time direction is presented for solving the Eq. (.). This method is second and fourth-order accurate in time and space components respectively and is unconditionally stable. Rashidinia et al. [25] developed two conditionally stable three level difference schemes of orders O(τ 2 +h 2 ) and O(τ 2 +h 4 ) based on non-polynomial cubic spline functions for the solution of one-dimensional second-order non-homogeneous hyperbolic equation. Gao and Chi in [2] proposed an explicit difference scheme for the numerical solution of Eq. (.). This method is unconditionally stable and has fifth and second order of accuracy in time and space directions respectively. A three level compact difference scheme which is unconditionally stable but is second and fourth order accurate in time and space directions is given in [9]. The aim of this paper is to introduce an unconditionally stable numerical method for the Eq. (.) which has high-order of accuracy in both time and space directions. We first discretize the spatial derivative of equation with a fourth-order compact finite difference and then apply Padé approximation method of fifth-order for the resulting system of ordinary differential equations. The outline of this paper is as follows. In Section 2, we apply a fourth-order compact finite difference scheme for the second order spatial derivative of partial differential equation (.) and obtain a system of ordinary differential equations (ODEs). Then we replace the matrix exponential in recurrence relation with the family of Padé approximations of fifth order to solve the resulting system of ODEs. In Section 3, we prove the unconditional stability property of proposed method. In section 4, the numerical results of applying the method of this article on three test problems for equation (.) are presented. Also in this section a comparison with the methods of [9, 2, 5, 9] is given. Finally a conclusion is drawn in Section 5. 2 Proposition of new method For a positive integer n let h = L n t. So we define denote the step size of spatial variable, x, and t for step size of time variable,
3 Page 3 of x i = ih, i = 0,,2,...,n, t k = k t, k = 0,,2,.... For derivation of the method, we first discretize Eq. (.) in space with a fourth order compact difference scheme to obtain a system of ODEs with unknown function at each spatial grid point. The fourth-order discretization of equation can be written as follows [20] or d 2 y = k(x), 0 < x < L, (2.4) dx2 k i = y xx,i = δ xx y i 2 δ xxk i h 2 + O(h 4 ), 0k i + k i+ + k i 2 = y i+ 2y i + y i h 2 + O(h 4 ), (2.5) where k i = k(x i ), y i = y(x i ) and δ xx y i = y i+ 2y i +y i. Now we rewrite (.) as follows h 2 2 u u (x,t) + 2α t2 t (x,t) + β 2 u(x,t) f (x,t) = 2 u (x,t). (2.6) x2 If we put u t (x,t) = v(x,t) then (2.6) can be written as { ut (x,t) = v(x,t), v t (x,t) + 2α v(x,t) + β 2 u(x,t) f (x,t) = u xx (x,t). If we discretize spatial derivative of (2.7) in each grid point by the fourth order scheme (2.5), we will obtain the following relation u i(t) = v i (t), 2 [0(v i(t) + 2αv i (t) + β 2 u i (t) f i (t)) + (v i+(t) + 2αv i+ (t) + β 2 u i+ (t) f i+ (t)) + (v i (t) + 2αv i (t) + β 2 u i (t) f i (t))] = h 2 (u i (t) 2u i (t) + u i+ (t)), (2.7) i =,2,...,n, (2.8) where u i (t) = u(x i,t), v i (t) = v(x i,t), u i (t) = u t(x i,t), v i (t) = v t(x i,t), and f i (t) = f (x i,t). The initial and boundary conditions (.2) and (.3) for equations (2.8) can be used as follows u i (0) = ϕ (x i ), v i (0) = ϕ 2 (x i ), u 0 (t) = g 0 (t), v 0 (t) = g 0(t), v 0(t) = g 0(t), u n (t) = g (t), v n (t) = g (t), v n(t) = g (t). (2.9) We assume that g 0 and g have second derivatives respect to time component. If we put U(t) = [u (t),u 2 (t),...,u n (t),v (t),v 2 (t),...,v n (t)] T and write equation (2.9) for each interior grid point, we obtain the system of 2n 2 ordinary differential equations A du(t) dt U(0) = U 0, = B U(t) +C(t), (2.0)
4 Page 4 of in which B = [ In O A = O A A = ] [ O In, B = B B ], C(t) = [0,0,...,0,C (t)] T, , 5β β 2 h h 2 β β 2 h β 2 h h 2 0 β β 2 h β 2 h h β h 2 5β h 2 β h β h 2 5β h 2 β h β h 2 5β h 2, C (t) = 2 f 0(t) f (t) + 2 f 2(t) 2 ( v 0 (t) + 2αv 0 (t) + β 2 u 0 (t) ) + u h 2 0 (t) 2 f (t) f 2(t) + 2 f 3(t). 2 f n 3(t) f n 2(t) + 2 f n (t) 2 f n 2(t) f n (t) + 2 f n(t) 2 ( v n (t) + 2αv n (t) + β 2 u n (t) ) + u h 2 n (t), and B 2 = 2αA and O is zero matrix of size n. Since matrix A is a diagonally dominant matrix we can conclude that A is a diagonally dominant matrix so it is invertible and we can write (2.0) as follows du(t) dt U(0) = U 0. = A B U(t) + A C(t), For solving system of ODEs (2.) we replace the matrix exponential in recurrence relation where DU(t) = du(t) dt (2.) U(t + t) = e td U(t), t = 0, t,2 t,..., (2.2) is given by (2.), with the family of Padé approximations. Using (2,2) Padé approximation, i.e. we can write (2.2) as follows e z = 2 + 6z + z2 2 6z + z 2 + O(z5 ), [2 6 td + ( td) 2 ]U(t + t) = [2 + 6 td + ( td) 2 ]U(t). (2.3)
5 Page 5 of Not that D 2 U(t) = d2 U(t) dt 2 is given by where G(t) = A C(t). Using (2.), (2.3) and (2.4) we can write d 2 U(t) dt 2 = (A B) 2 U(t) + A B G(t) + G (t), (2.4) U k+ = M N U k + M [(6 ti t 2 A B)G k+ + (6 ti + t 2 A B)G k + t 2 (G k G k+ )], (2.5) where U k = U(t k ), G k = G(t k ), G k = G (t k ), M = [2I 6 ta B + ( ta B) 2 ], N = [2I + 6 ta B + ( ta B) 2 ] and I is identity matrix of order 2n 2. Now we prove that the method is unconditionally stable. 3 Stability analysis For homogeneous boundary conditions, the proposed method (2.5) can be written as U k+ = φ U k, k = 0,,2,..., (3.6) where the amplification matrix is given by φ = [2I 6 ta B + ( ta B) 2 ] [2I + 6 ta B + ( ta B) 2 ]. (3.7) Let λ be an eigenvalue of A B and z = tλ. For unconditional stability of the new method it is necessary that the absolute values of the eigenvalues of amplification matrix (φ) be less than or equal to one, i. e z + z 2 2 6z + z 2. (3.8) Lemma 3.. ([2]) The relation (3.8) holds iff z be in the left-half complex plane. According to Lemma 3. if z or equivalently λ be in the left-half complex plane, the relation (3.8) holds. So we should prove that each eigenvalue of A B has negative real part. Lemma 3.2. ([2]) If the Hermitian matrix Q = [q i j ] 2n 2 2n 2 be strictly diagonally dominant with positive (negative) diagonal elements then Q is positive (negative) definite matrix. Lemma 3.3. Each eigenvalue of matrix B has negative real part. Proof. Let η = a + bi be an eigenvalue of B with corresponding nonzero eigenvector ξ = [ ] T ξ ξ 2, i.e. [ ][ ] [ ] O In ξ ξ = η. (3.9) B B 2 From (3.9) we have ξ 2 ξ 2 ξ 2 = ηξ, B ξ + B 2 ξ 2 = ηξ 2, (3.20) therefore we have B ξ + ηb 2 ξ = η 2 ξ or (B + ηb 2 η 2 I)ξ = 0. If ξ = 0 then from (3.20) we see that ξ = 0 which is impossible. So ξ 0 and we should have det(b + ηb 2 η 2 I) = 0 or det(η 2 I ηb 2 B ) = 0. Now we suppose that real part of η (i.e. a) is nonnegative. If we have b = 0 then (η 2 I ηb 2 B ) = (a 2 I ab 2 B ) which is positive definite and contradiction with det(η 2 I ηb 2 B ) = 0 (note that from Lemma 3.2, B and B 2 are negative definite matrices). So we should have b 0. In this case we put R = (η 2 I ηb 2 B )i = ((a 2 b 2 )I ab 2 B )i b(2ai B 2 ).
6 Page 6 of Since a is nonnegative, the symmetric part of R, i.e. 2 (R + R ) = b(2ai B 2 ) is positive or negative definite matrix depending to sign of b. This gives that R is nonsymmetric positive or negative definite [3] and so it is nonsingular or det(η 2 I ηb 2 B ) 0 which is contradiction. So the real part of η should be negative. Definition 3.. The matrix A is called generalized positive definite (SPD) if its symmetric part is positive definite, and is called generalized negative definite (GND) if A is generalized positive definite. Lemma 3.4. ([0]) All eigenvalues of matrix Q have negative real parts, if and only if there exist a generalized negative-definite matrix G and a symmetric positive-definite matrix L such that Q = GL. Theorem 3.. Each eigenvalue of matrix A B has negative real part. Proof. From Lemma 3, all eigenvalues of matrix B have negative real parts so from Lemma 4 there exists a SPD matrix P in which G = BP is GND matrix. From Lemma 5 and regarding to this point that A is SPD matrix and putting L = P A we can conclude that all eigenvalues of matrix Q = GL = BPP A = BA have negative real parts. Note that eigenvalues of A B and BA are similar. This theorem shows that each eigenvalue of matrix A B is in left-half complex plane and we can conclude that the proposed method is unconditionally stable. 4 Numerical results In this section we present the numerical results of the new method on three test problems. We performed our computations using Matlab 7 software on a PC with Intel Core 2 Duo, 2.8 GHz CPU and 2 GB RAM. We tested the accuracy and stability of the method presented in this paper by performing the mentioned method for different values of t and h. The L and Root-Mean-Square (RMS) errors obtained by new method are shown. Also we calculated the computational orders of the method presented in this article (denoted by C-Order) with the following formula log( E E 2 ) log( h h 2 ), in which E and E 2 are errors correspond to grids with mesh size h and h 2 respectively. Table : Absolute errors obtained for Test problem with h = π/30, t = 0. x Method of [2] Method of [5] Present method CPU time (s) T = π/ π/ π/ π/ π/ T = 2 π/ π/ π/ π/ π/
7 Communications in Numerical Analysis Page 7 of 4. Test Problem We consider the following hyperbolic equation 2u u 2u (x,t) + 4 (x,t) + 2u(x,t) = 2 (x,t), 2 t t x where the exact solution is given by u(x,t) = e t sin(x). 0 x π, t > 0, Figure : Surface plots of approximate solutions for Test problems (left panel) and 2 (right panel). The boundary and initial conditions can be obtained from exact solution. We compare the numerical results of new method presented in this paper with the results of [2, 5]. In Table we compare the absolute error obtained in solving Test problem with h = π /30, t = 0., T = and T = 2. As we see the new method has better results. Table 2 shows the L and Root-Mean-Square (RMS) errors, C-Order and CPU time of new method for solving Test problem with T =, t = 0.02 and several values of h. As we see the method achieve an accuracy of order 0 9 in under 0.6 seconds. Also the computational orders in Table 2 shows the fourth-order accuracy of method in spatial variable. Table 2: Numerical results obtained for Test problem with t = 0.02 h RMS L C-Order CPU time(s) π /5 π /0 π /20 π /40 π / Test Problem 2 We consider the following hyperbolic equation 2u u 2u u = 6x(x )(5x2 5x + )e t, t2 t x2 0 x, t > 0,
8 Page 8 of where the exact solution is given by u(x,t) = x 3 (x ) 3 e t. The boundary and initial conditions can be obtained from exact solution. We compare the numerical results of proposed method with the results of [9, 9]. Table 3 shows the L and Root-Mean-Square (RMS) errors, C-Order and CPU time of new method for solving Test problem 2 with T =, t = /20 and several values of h. As we see the method achieve an accuracy of order 0 0 in under second. Also the computational orders in Table 3 shows the fourth-order accuracy of method in spatial variable. Table 3: Numerical results obtained for Test problem 2 with t = /20 h RMS L C-Order CPU time(s) / / / / / / In Table 4 we compare the absolute error obtained in solving Test problem 2 with t = /64, T = and several values of h. As we see the new method has better results. Also Figure shows the surface plots of approximate solutions of Test problems and 2. Table 4: Absolute errors obtained for Test problem 2 with t = /64 h Method of [9] Method of [9] Present method CPU time(s) / / / / / Test Problem 3 We consider the following hyperbolic equation 2 u u + 2α t2 t + β 2 u = 2 u x 2 + β 2 sin(x)cos(t) 2α sin(x)sin(t), 0 x 2, t > 0, where the exact solution is given by u(x,t) = sin(x)cos(t). The boundary and initial conditions can be obtained from exact solution. This test problem is given in [20]. We show the RMS error obtained in solving this test problem at T = with h = 0.05, t = 0.02 and several values of α and β in Table 5. Table 6 presents the numerical results in solving Test problem 3 with α = 0, β = 5, h = /80 and several values of t. As we see with the above parameters the new method achieve an accuracy of order 0 9 in under 7 seconds. Figure 2 shows the surface plot of approximate solution of Test problem 3.
9 Page 9 of Figure 2: Surface plot of approximate solution for Test problems 3. Table 5: RMS error obtained for Test problem 3 with h = 0.05 and t = 0.02 T=0.5 T =.0 T = 2.0 T = 5.0 α = 0, β = α = 20, β = α = 50, β = α = 00, β = α = 00, β = Table 6: Numerical results obtained for Test problem 2 with α = 0, β = 5 and h = /80 t L RMS C-Order CPU time(s) / / / / / Conclusion In this paper we proposed a class of high order compact schemes for solving the D linear hyperbolic equation. We combined a high-order compact finite difference scheme of fourth-order to approximate the spatial derivative and a fifth order Padé approximation for the time integration of the resulted linear system of ordinary differential equations. The proposed method for solving the mentioned equation has high order of accuracy and is unconditionally stable. Computational experiments confirmed the unconditional stability and high accuracy of the proposed method and presented that the new method can achieve good accuracy in short CPU time.
10 Page 0 of Acknowledgements The author would like to thank the referee for his (or her) comments. The author is also thankful to the University of Kashan for the financial support (Grant No /) of this work. References [] S. Abbasbandy, H. Roohani Ghehsareh, A new semi-analytical solution of the telegraph equation with integral condition, Z. Naturforsch. 66a (20) [2] O. Axelsson, Iterative solution methods, New York, Cambridge University Press, (996). [3] A. Bhaya, E. Kaszkurewicz, R. Santos, Characterizations of classes of stable matrices, Lin. Alg. Appl, 374 (2003) [4] Ali J. Chamkha, Heat and Mass Transfer from MHD Flow over a Moving Permeable Cylinder with Heat Generation or Absorption and Chemical Reaction, Communications in Numerical Analysis, 20 (20) [5] M. Dehghan, A. Shokri, A numerical method for solving the hyperbolic telegraph equation, Numer. Methods Partial Differential Equations, 24 (2008) [6] M. Dehghan, M. Lakestani, The use of Chebyshev cardinal functions for solution of the second-order onedimensional telegraph equation, Numer. Methods Partial Differential Equations, 25 (2009) [7] M. Dehghan, A. Ghesmati, Solution of the second-order one-dimensional hyperbolic telegraph equation by using the dual reciprocity boundary integral equation (DRBIE) method, Eng. Anal. Bound. Elem, 34 (200) [8] M. Dehghan, S. A. Yousefi, A. Lotfi, The use of He s variational iteration method for solving the telegraph and fractional telegraph equations, Int. J. Numer. Meth. Biomed. Eng, 27 (20) [9] H. Ding, Y. Zhang, A new unconditionally stable compact difference scheme of O(τ 2 + h 4 ) for the D linear hyperbolic equation, Appl. Math. Comput, 207 (2009) [0] G. R. Duan, R. J. Patton, A note on Hurwitz stability of matrices, Automatica, 34 (998) [] M. S. El-Azab, M. El-Ghamel, A numerical algorithm for the solution of telegraph equations, Appl. Math. Comput, 90 (2007) [2] F. Gao, C. M. Chi, Unconditionally stable difference schemes for a one-space-dimensional linear hyperbolic equation, Appl. Math. Comput, 87 (2007) [3] D. Harville, Matrix Algebra from a Statistician s Perspective, Springer, New York, (997).
11 Page of [4] W. Liao, M. Dehghan, A. Mohebbi, Direct numerical method for an inverse problem of a parabolic partial differential equation, J. Comput. Appl. Math, 232 (2009) [5] H. W. Liu, L. B. Liu, An unconditionally stable spline difference scheme of O(k 2 + h 4 ) for solving the second order D linear hyperbolic equation, Math. Comput. Model, 49 (2009) [6] R. K. Mohanty, M. K. Jain, k. George, On the use of high order difference methods for the system of one space second order non-linear hyperbolic equations with variable coefficients, J. Comp. Appl. Math, 72 (996) [7] R. K. Mohanty, An unconditionally stable difference scheme for the one-space-dimensional linear hyperbolic equation, Appl. Math. Lett, 7 (2004) [8] R. K. Mohanty, An unconditionally stable finite difference formula for a linear second order one space dimensional hyperbolic equation with variable coefficients, Appl. Math. Comput, 65 (2005) [9] R. K. Mohanty, New unconditionally stable difference schemes for the solution of multi-dimensional telegraphic equations, Int. J. Compu. Math, 86 (2009) [20] A. Mohebbi, M. Dehghan, High order compact solution of the one-space-dimensional linear hyperbolic equation, Numer. Methods Partial Differential Equations, 24 (2008) [2] A. Mohebbi, M. Dehghan, High-order compact solution of the one-dimensional heat and advection-diffusion equations, Appl. Math. Model, 34 (200) [22] A. Mohebbi, M. Dehghan, High-order scheme for determination of a control parameter in an inverse problem from the over-specified data, Comput. Phys. Commun, 8 (200) [23] K. Parand, S. Kazem, A.R. Rezaei, A numerical study on reaction-diffusion problem using radial basis functions, Communications in Numerical Analysis, 20 (20) -. [24] D. M. Pozar, Microwave engineering, Addison-Wesley, (990). [25] J. Rashidinia, R. Mohammadi, R. Jalilian, Spline methods for the solution of hyperbolic equation with variable coefficients, Numer. Methods Partial Differential Equations, 32 (2006) -9. [26] A. Saadatmandi, M. Dehghan, Numerical solution of hyperbolic telegraph equation using the Chebyshev Tau method, Numer. Methods Partial Differential Equations, 26 (200) [27] E. H. Twizell, An explicit difference method for the wave equation with extended stability range, BIT, 9 (979)
A Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional Hyperbolic Telegraph Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(01) No.3,pp.59-66 A Differential Quadrature Algorithm for the Numerical Solution of the Second-Order One Dimensional
More informationNumerical Solution of One-dimensional Telegraph Equation using Cubic B-spline Collocation Method
2014 (2014) 1-8 Available online at www.ispacs.com/iasc Volume 2014, Year 2014 Article ID iasc-00042, 8 Pages doi:10.5899/2014/iasc-00042 Research Article Numerical Solution of One-dimensional Telegraph
More informationA Numerical Study of One-Dimensional Hyperbolic Telegraph Equation
Journal of Mathematics and System Science 7 (2017) 62-72 doi: 10.17265/2159-5291/2017.02.003 D DAVID PUBLISHING A Numerical Study of One-Dimensional Hyperbolic Telegraph Equation Shaheed N. Huseen Thi-Qar
More informationDifferential transformation method for solving one-space-dimensional telegraph equation
Volume 3, N 3, pp 639 653, 2 Copyright 2 SBMAC ISSN -825 wwwscielobr/cam Differential transformation method for solving one-space-dimensional telegraph equation B SOLTANALIZADEH Young Researchers Club,
More informationResearch Article The Spectral Method for Solving Sine-Gordon Equation Using a New Orthogonal Polynomial
International Scholarly Research Network ISRN Applied Mathematics Volume, Article ID 4673, pages doi:.54//4673 Research Article The Spectral Method for Solving Sine-Gordon Equation Using a New Orthogonal
More informationNumerical analysis of the one-dimensional Wave equation subject to a boundary integral specification
Numerical analysis of the one-dimensional Wave equation subject to a boundary integral specification B. Soltanalizadeh a, Reza Abazari b, S. Abbasbandy c, A. Alsaedi d a Department of Mathematics, University
More informationNumerical solution of General Rosenau-RLW Equation using Quintic B-splines Collocation Method
Available online at www.ispacs.com/cna Volume 2012, Year 2012 Article ID cna-00129, 16 pages doi:10.5899/2012/cna-00129 Research Article Numerical solution of General Rosenau-RLW Equation using Quintic
More informationFusion Higher -Order Parallel Splitting Methods for. Parabolic Partial Differential Equations
International Mathematical Forum, Vol. 7, 0, no. 3, 567 580 Fusion Higher -Order Parallel Splitting Methods for Parabolic Partial Differential Equations M. A. Rehman Department of Mathematics, University
More informationNumerical studies of non-local hyperbolic partial differential equations using collocation methods
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 6, No. 3, 2018, pp. 326-338 Numerical studies of non-local hyperbolic partial differential equations using collocation methods
More informationResearch Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method
Applied Mathematics Volume 2012, Article ID 605741, 10 pages doi:10.1155/2012/605741 Research Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method M. Zarebnia
More informationIn the present work, we consider the singularly perturbed Volterra integro-differential equations (SVIDE) x K(x, t, ε, y(t)) dt, x I = [0, X], (1)
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.12(211) No.4,pp.43-441 An Approximation Algorithm for the Solution of the Singularly Perturbed Volterra Integro-differential
More informationADOMIAN-TAU OPERATIONAL METHOD FOR SOLVING NON-LINEAR FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS WITH PADE APPROXIMANT. A. Khani
Acta Universitatis Apulensis ISSN: 1582-5329 No 38/214 pp 11-22 ADOMIAN-TAU OPERATIONAL METHOD FOR SOLVING NON-LINEAR FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS WITH PADE APPROXIMANT A Khani Abstract In this
More informationExistence of solution and solving the integro-differential equations system by the multi-wavelet Petrov-Galerkin method
Int. J. Nonlinear Anal. Appl. 7 (6) No., 7-8 ISSN: 8-68 (electronic) http://dx.doi.org/.75/ijnaa.5.37 Existence of solution and solving the integro-differential equations system by the multi-wavelet Petrov-Galerkin
More informationSolution of Non Linear Singular Perturbation Equation. Using Hermite Collocation Method
Applied Mathematical Sciences, Vol. 7, 03, no. 09, 5397-5408 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/ams.03.37409 Solution of Non Linear Singular Perturbation Equation Using Hermite Collocation
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 informationCubic B-spline collocation method for solving time fractional gas dynamics equation
Cubic B-spline collocation method for solving time fractional gas dynamics equation A. Esen 1 and O. Tasbozan 2 1 Department of Mathematics, Faculty of Science and Art, Inönü University, Malatya, 44280,
More informationApplication of fractional sub-equation method to the space-time fractional differential equations
Int. J. Adv. Appl. Math. and Mech. 4(3) (017) 1 6 (ISSN: 347-59) Journal homepage: www.ijaamm.com IJAAMM International Journal of Advances in Applied Mathematics and Mechanics Application of fractional
More informationFinite Difference Method of Fractional Parabolic Partial Differential Equations with Variable Coefficients
International Journal of Contemporary Mathematical Sciences Vol. 9, 014, no. 16, 767-776 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.1988/ijcms.014.411118 Finite Difference Method of Fractional Parabolic
More informationCubic spline Numerov type approach for solution of Helmholtz equation
Journal of Linear and Topological Algebra Vol. 03, No. 01, 2014, 47-54 Cubic spline Numerov type approach for solution of Helmholtz equation J. Rashidinia a, H. S. Shekarabi a and M. Aghamohamadi a a Department
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationA parallel algorithm for the heat equation with derivative boundary conditions
International Mathematical Forum, 2, 2007, no. 12, 565-574 A parallel algorithm for the heat equation with derivative boundary conditions M. Akram University College of Information Technology University
More informationAdomian Decomposition Method with Laguerre Polynomials for Solving Ordinary Differential Equation
J. Basic. Appl. Sci. Res., 2(12)12236-12241, 2012 2012, TextRoad Publication ISSN 2090-4304 Journal of Basic and Applied Scientific Research www.textroad.com Adomian Decomposition Method with Laguerre
More informationThe comparison of optimal homotopy asymptotic method and homotopy perturbation method to solve Fisher equation
Computational Methods for Differential Equations http://cmdetabrizuacir Vol 4, No, 206, pp 43-53 The comparison of optimal homotopy asymptotic method and homotopy perturbation method to solve Fisher equation
More informationA Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations
Mathematics A Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations Mohamed Meabed KHADER * and Ahmed Saied HENDY Department of Mathematics, Faculty of Science,
More informationResearch Article A New Method for Riccati Differential Equations Based on Reproducing Kernel and Quasilinearization Methods
Abstract and Applied Analysis Volume 0, Article ID 603748, 8 pages doi:0.55/0/603748 Research Article A New Method for Riccati Differential Equations Based on Reproducing Kernel and Quasilinearization
More informationA Parallel Algorithm for the Inhomogeneous Advection Equation
International Mathematical Forum, 3, 2008, no. 10, 463-472 A Parallel Algorithm for the Inhomogeneous Advection Equation M. Akram PUCIT, University of the Punjab, Old Campus Lahore-54000, Pakistan m.akram@pucit.edu.pk,
More informationAnalytical solution for determination the control parameter in the inverse parabolic equation using HAM
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 2 (December 2017, pp. 1072 1087 Applications and Applied Mathematics: An International Journal (AAM Analytical solution
More informationAn Exponential High-Order Compact ADI Method for 3D Unsteady Convection Diffusion Problems
An Exponential High-Order Compact ADI Method for 3D Unsteady Convection Diffusion Problems Yongbin Ge, 1 Zhen F. Tian, 2 Jun Zhang 3 1 Institute of Applied Mathematics and Mechanics, Ningxia University,
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
More informationNumerical Solution of Two-Dimensional Volterra Integral Equations by Spectral Galerkin Method
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 159-174 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 Numerical Solution of Two-Dimensional Volterra
More informationRunge-Kutta Method for Solving Uncertain Differential Equations
Yang and Shen Journal of Uncertainty Analysis and Applications 215) 3:17 DOI 1.1186/s4467-15-38-4 RESEARCH Runge-Kutta Method for Solving Uncertain Differential Equations Xiangfeng Yang * and Yuanyuan
More informationNumerical solution of Maxwell equations using local weak form meshless techniques
Journal of mathematics and computer science 13 2014), 168-185 Numerical solution of Maxwell equations using local weak form meshless techniques S. Sarabadan 1, M. Shahrezaee 1, J.A. Rad 2, K. Parand 2,*
More informationResearch Article Numerical Solution of the Inverse Problem of Determining an Unknown Source Term in a Heat Equation
Applied Mathematics Volume 22, Article ID 39876, 9 pages doi:.55/22/39876 Research Article Numerical Solution of the Inverse Problem of Determining an Unknown Source Term in a Heat Equation Xiuming Li
More informationNumerical solution of fourth order parabolic partial dierential equation using parametric septic splines
Hacettepe Journal of Mathematics and Statistics Volume 5 20, 07 082 Numerical solution of fourth order parabolic partial dierential equation using parametric septic splines Arshad Khan and Talat Sultana
More informationOn the Inverting of a General Heptadiagonal Matrix
British Journal of Applied Science & Technology 18(5): 1-, 2016; Article no.bjast.313 ISSN: 2231-0843, NLM ID: 101664541 SCIENCEDOMAIN international www.sciencedomain.org On the Inverting of a General
More informationResearch Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
Mathematical Problems in Engineering Volume 1, Article ID 693453, 1 pages doi:11155/1/693453 Research Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method
More informationSolution of Differential Equations of Lane-Emden Type by Combining Integral Transform and Variational Iteration Method
Nonlinear Analysis and Differential Equations, Vol. 4, 2016, no. 3, 143-150 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/nade.2016.613 Solution of Differential Equations of Lane-Emden Type by
More informationResearch Article Solution of the Porous Media Equation by a Compact Finite Difference Method
Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 2009, Article ID 9254, 3 pages doi:0.55/2009/9254 Research Article Solution of the Porous Media Equation by a Compact Finite Difference
More informationThe Sinc-Collocation Method for Solving the Telegraph Equation
Te Sinc-Collocation Metod for Solving te Telegrap Equation E. Hesameddini *1, E. Asadolaifard Department of Matematics, Faculty of Basic Sciences, Siraz University of Tecnology, Siraz, Iran *1 esameddini@sutec.ac.ir;
More informationChebyshev finite difference method for a Two Point Boundary Value Problems with Applications to Chemical Reactor Theory
Iranian Journal of Mathematical Chemistry, Vol. 3, o., February 22, pp. 7 IJMC Chebyshev finite difference method for a Two Point Boundary Value Problems with Applications to Chemical Reactor Theory ABBAS
More informationPotential Symmetries and Differential Forms. for Wave Dissipation Equation
Int. Journal of Math. Analysis, Vol. 7, 2013, no. 42, 2061-2066 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2013.36163 Potential Symmetries and Differential Forms for Wave Dissipation
More informationChebyshev finite difference method for solving a mathematical model arising in wastewater treatment plants
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 6, No. 4, 2018, pp. 448-455 Chebyshev finite difference method for solving a mathematical model arising in wastewater treatment
More informationA new approach to solve fuzzy system of linear equations by Homotopy perturbation method
Journal of Linear and Topological Algebra Vol. 02, No. 02, 2013, 105-115 A new approach to solve fuzzy system of linear equations by Homotopy perturbation method M. Paripour a,, J. Saeidian b and A. Sadeghi
More informationDiagonalizing Hermitian Matrices of Continuous Functions
Int. J. Contemp. Math. Sciences, Vol. 8, 2013, no. 5, 227-234 HIKARI Ltd, www.m-hikari.com Diagonalizing Hermitian Matrices of Continuous Functions Justin Cyr 1, Jason Ekstrand, Nathan Meyers 2, Crystal
More informationApplication of linear combination between cubic B-spline collocation methods with different basis for solving the KdV equation
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 4, No. 3, 016, pp. 191-04 Application of linear combination between cubic B-spline collocation methods with different basis
More informationEXACT SOLUTIONS OF NON-LINEAR FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS BY FRACTIONAL SUB-EQUATION METHOD
THERMAL SCIENCE, Year 15, Vol. 19, No. 4, pp. 139-144 139 EXACT SOLUTIONS OF NON-LINEAR FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS BY FRACTIONAL SUB-EQUATION METHOD by Hong-Cai MA a,b*, Dan-Dan YAO a, and
More informationAnalysis of Fractional Nonlinear Differential Equations Using the Homotopy Perturbation Method
Analysis of Fractional Nonlinear Differential Equations Using the Homotopy Perturbation Method Mehmet Ali Balcı and Ahmet Yıldırım Ege University, Department of Mathematics, 35100 Bornova-İzmir, Turkey
More informationA Smoothing Newton Method for Solving Absolute Value Equations
A Smoothing Newton Method for Solving Absolute Value Equations Xiaoqin Jiang Department of public basic, Wuhan Yangtze Business University, Wuhan 430065, P.R. China 392875220@qq.com Abstract: In this paper,
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationNumerical study of time-fractional hyperbolic partial differential equations
Available online at wwwisr-publicationscom/jmcs J Math Computer Sci, 7 7, 53 65 Research Article Journal Homepage: wwwtjmcscom - wwwisr-publicationscom/jmcs Numerical study of time-fractional hyperbolic
More informationSolving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions.
Journal of Mathematical Modeling Vol 1, No 1, 213, pp 28-4 JMM Solving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions Farshid
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationDifferential Transform Method for Solving the Linear and Nonlinear Westervelt Equation
Journal of Mathematical Extension Vol. 6, No. 3, (2012, 81-91 Differential Transform Method for Solving the Linear and Nonlinear Westervelt Equation M. Bagheri Islamic Azad University-Ahar Branch J. Manafianheris
More informationA Third-degree B-spline Collocation Scheme for Solving a Class of the Nonlinear Lane-Emden Type Equations
Iranian Journal of Mathematical Sciences and Informatics Vol. 12, No. 2 (2017), pp 15-34 DOI: 10.7508/ijmsi.2017.2.002 A Third-degree B-spline Collocation Scheme for Solving a Class of the Nonlinear Lane-Emden
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 informationLecture 1 INF-MAT3350/ : Some Tridiagonal Matrix Problems
Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems Tom Lyche University of Oslo Norway Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems p.1/33 Plan for the day 1. Notation
More informationAn efficient algorithm on timefractional. equations with variable coefficients. Research Article OPEN ACCESS. Jamshad Ahmad*, Syed Tauseef Mohyud-Din
OPEN ACCESS Research Article An efficient algorithm on timefractional partial differential equations with variable coefficients Jamshad Ahmad*, Syed Tauseef Mohyud-Din Department of Mathematics, Faculty
More informationarxiv: v1 [math.na] 21 Oct 2014
Computing Symmetric Positive Definite Solutions of Three Types of Nonlinear Matrix Equations arxiv:1410.5559v1 [math.na] 21 Oct 2014 Negin Bagherpour a, Nezam Mahdavi-Amiri a, a Department of Mathematical
More informationA Point Interpolation Meshless Method for the. Numerical Solution of the Singularly Perturbed. Integral and Integro-differential Equations
Int. Journal of Math. Analysis, Vol. 7, 2013, no. 13, 643-656 HIKARI Ltd, www.m-hikari.com A Point Interpolation Meshless Method for the Numerical Solution of the Singularly Perturbed Integral and Integro-differential
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 informationNumerical Solution of Nonlocal Parabolic Partial Differential Equation via Bernstein Polynomial Method
Punjab University Journal of Mathematics (ISSN 116-2526) Vol.48(1)(216) pp. 47-53 Numerical Solution of Nonlocal Parabolic Partial Differential Equation via Bernstein Polynomial Method Kobra Karimi Department
More informationBenha University Faculty of Science Department of Mathematics. (Curriculum Vitae)
Benha University Faculty of Science Department of Mathematics (Curriculum Vitae) (1) General *Name : Mohamed Meabed Bayuomi Khader *Date of Birth : 24 May 1973 *Marital Status: Married *Nationality : Egyptian
More informationSUMMATION-BY-PARTS IN TIME: THE SECOND DERIVATIVE
Department of Mathematics SUMMATION-BY-PARTS IN TIME: THE SECOND DERIVATIVE Jan Nordström and Tomas Lundquist LiTH-MAT-R--2014/11--SE Department of Mathematics Linköping University S-581 83 Linköping SUMMATION-BY-PARTS
More informationApplication of Semiorthogonal B-Spline Wavelets for the Solutions of Linear Second Kind Fredholm Integral Equations
Appl Math Inf Sci 8, No, 79-84 (4) 79 Applied Mathematics & Information Sciences An International Journal http://dxdoiorg/78/amis/8 Application of Semiorthogonal B-Spline Wavelets for the Solutions of
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationExact Solutions of the Generalized- Zakharov (GZ) Equation by the Infinite Series Method
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 05, Issue (December 010), pp. 61 68 (Previously, Vol. 05, Issue 10, pp. 1718 175) Applications and Applied Mathematics: An International
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationNUMERICAL SOLUTION OF FRACTIONAL RELAXATION OSCILLATION EQUATION USING CUBIC B-SPLINE WAVELET COLLOCATION METHOD
italian journal of pure and applied mathematics n. 36 2016 (399 414) 399 NUMERICAL SOLUTION OF FRACTIONAL RELAXATION OSCILLATION EQUATION USING CUBIC B-SPLINE WAVELET COLLOCATION METHOD Raghvendra S. Chandel
More informationBulletin of the. Iranian Mathematical Society
ISSN: 1017-060X (Print) ISSN: 1735-8515 (Online) Bulletin of the Iranian Mathematical Society Vol. 41 (2015), No. 5, pp. 1259 1269. Title: A uniform approximation method to solve absolute value equation
More informationNew Approach of ( Ǵ/G ) Expansion Method. Applications to KdV Equation
Journal of Mathematics Research; Vol. 6, No. ; ISSN 96-9795 E-ISSN 96-989 Published by Canadian Center of Science and Education New Approach of Ǵ/G Expansion Method. Applications to KdV Equation Mohammad
More informationA generalization of the Gauss-Seidel iteration method for solving absolute value equations
A generalization of the Gauss-Seidel iteration method for solving absolute value equations Vahid Edalatpour, Davod Hezari and Davod Khojasteh Salkuyeh Faculty of Mathematical Sciences, University of Guilan,
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
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 information(Received 10 December 2011, accepted 15 February 2012) x x 2 B(x) u, (1.2) A(x) +
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol13(212) No4,pp387-395 Numerical Solution of Fokker-Planck Equation Using the Flatlet Oblique Multiwavelets Mir Vahid
More informationThe Representation of Energy Equation by Laplace Transform
Int. Journal of Math. Analysis, Vol. 8, 24, no. 22, 93-97 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.2988/ijma.24.442 The Representation of Energy Equation by Laplace Transform Taehee Lee and Hwajoon
More informationResearch Article He s Variational Iteration Method for Solving Fractional Riccati Differential Equation
International Differential Equations Volume 2010, Article ID 764738, 8 pages doi:10.1155/2010/764738 Research Article He s Variational Iteration Method for Solving Fractional Riccati Differential Equation
More informationHOMOTOPY ANALYSIS METHOD FOR SOLVING KDV EQUATIONS
Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (21), 89 98 HOMOTOPY ANALYSIS METHOD FOR SOLVING KDV EQUATIONS Hossein Jafari and M. A. Firoozjaee Abstract.
More informationDomain decomposition schemes with high-order accuracy and unconditional stability
Domain decomposition schemes with high-order accuracy and unconditional stability Wenrui Hao Shaohong Zhu March 7, 0 Abstract Parallel finite difference schemes with high-order accuracy and unconditional
More informationNumerical methods for a fractional diffusion/anti-diffusion equation
Numerical methods for a fractional diffusion/anti-diffusion equation Afaf Bouharguane Institut de Mathématiques de Bordeaux (IMB), Université Bordeaux 1, France Berlin, November 2012 Afaf Bouharguane Numerical
More informationExact Solutions for a Class of Singular Two-Point Boundary Value Problems Using Adomian Decomposition Method
Applied Mathematical Sciences, Vol 6, 212, no 122, 697-618 Exact Solutions for a Class of Singular Two-Point Boundary Value Problems Using Adomian Decomposition Method Abdelhalim Ebaid 1 and Mona D Aljoufi
More informationA simple local variational iteration method for solving nonlinear Lane-Emden problems
A simple local variational iteration method for solving nonlinear Lane-Emden problems Asghar Ghorbani a,, Mojtaba Bakherad b a Department of Applied Mathematics, Faculty of Mathematical Sciences, Ferdowsi
More informationQuintic B-Spline Galerkin Method for Numerical Solutions of the Burgers Equation
5 1 July 24, Antalya, Turkey Dynamical Systems and Applications, Proceedings, pp. 295 39 Quintic B-Spline Galerkin Method for Numerical Solutions of the Burgers Equation İdris Dağ 1,BülentSaka 2 and Ahmet
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
More informationResearch Article An Operator-Difference Method for Telegraph Equations Arising in Transmission Lines
Discrete Dynamics in Nature and Society Volume 2, Article ID 565, 7 pages doi:.55/2/565 Research Article An Operator-Difference Method for Telegraph Equations Arising in Transmission Lines Mehmet Emir
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 informationSolving Non-Homogeneous Coupled Linear Matrix Differential Equations in Terms of Matrix Convolution Product and Hadamard Product
Journal of Informatics and Mathematical Sciences Vol. 10, Nos. 1 & 2, pp. 237 245, 2018 ISSN 0975-5748 (online); 0974-875X (print) Published by RGN Publications http://www.rgnpublications.com http://dx.doi.org/10.26713/jims.v10i1-2.647
More informationResearch Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation with Inconsistent Initial and Boundary Conditions
International Mathematics and Mathematical Sciences Volume 212, Article ID 82197, 13 pages doi:1.1155/212/82197 Research Article L-Stable Derivative-Free Error-Corrected Trapezoidal Rule for Burgers Equation
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationHigh-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
More informationConformable variational iteration method
NTMSCI 5, No. 1, 172-178 (217) 172 New Trends in Mathematical Sciences http://dx.doi.org/1.2852/ntmsci.217.135 Conformable variational iteration method Omer Acan 1,2 Omer Firat 3 Yildiray Keskin 1 Galip
More informationNumerical solution of optimal control problems by using a new second kind Chebyshev wavelet
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 4, No. 2, 2016, pp. 162-169 Numerical solution of optimal control problems by using a new second kind Chebyshev wavelet Mehdi
More informationProperties for the Perron complement of three known subclasses of H-matrices
Wang et al Journal of Inequalities and Applications 2015) 2015:9 DOI 101186/s13660-014-0531-1 R E S E A R C H Open Access Properties for the Perron complement of three known subclasses of H-matrices Leilei
More informationMath 331 Homework Assignment Chapter 7 Page 1 of 9
Math Homework Assignment Chapter 7 Page of 9 Instructions: Please make sure to demonstrate every step in your calculations. Return your answers including this homework sheet back to the instructor as a
More informationEE263: Introduction to Linear Dynamical Systems Review Session 5
EE263: Introduction to Linear Dynamical Systems Review Session 5 Outline eigenvalues and eigenvectors diagonalization matrix exponential EE263 RS5 1 Eigenvalues and eigenvectors we say that λ C is an eigenvalue
More informationSolution of Fourth Order Boundary Value Problems by Numerical Algorithms Based on Nonpolynomial Quintic Splines
Journal of Numerical Mathematics and Stochastics, 4(1) : 13-25, 2012 http://www.jnmas.org/jnmas4-2.pdf JNM@S Euclidean Press, LLC Online: ISSN 2151-2302 Solution of Fourth Order Boundary Value Problems
More informationLyapunov stability ORDINARY DIFFERENTIAL EQUATIONS
Lyapunov stability ORDINARY DIFFERENTIAL EQUATIONS An ordinary differential equation is a mathematical model of a continuous state continuous time system: X = < n state space f: < n! < n vector field (assigns
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 informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
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