RATIONAL CHEBYSHEV COLLOCATION METHOD FOR SOLVING NONLINEAR ORDINARY DIFFERENTIAL EQUATIONS OF LANE-EMDEN TYPE
|
|
- Marjorie Paula Shaw
- 5 years ago
- Views:
Transcription
1 INTERNATIONAL JOURNAL OF INFORMATON AND SYSTEMS SCIENCES Volume 6, Number 1, Pages c 2010 Institute for Scientific Computing and Information RATIONAL CHEBYSHEV COLLOCATION METHOD FOR SOLVING NONLINEAR ORDINARY DIFFERENTIAL EQUATIONS OF LANE-EMDEN TYPE KOUROSH PARAND AND MEHDI SHAHINI (Communicated by H. Azari) Abstract. Lane-Emden equation is a nonlinear singular equation that plays an important role in the astrophysics. In this paper, we have applied the collocation method based on rational Chebyshev functions to solve Lane-Emden type equations. The method reduces solving the nonlinear ordinary differential equation to solving a system of nonlinear algebraic equations. The comparison of the results with the other numerical methods shows that our solutions are highly accurate and the method is convergent. Key Words. Lane-Emden equation, rational Chebyshev, collocation method, nonlinear ODE, semi-infinite interval. 1. Introduction Many science and engineering problems arise in unbounded domains. Different spectral methods have been proposed for solving problems on unbounded domains. A direct approach is to use spectral method associated with orthogonal systems in unbounded domains, such as the Hermite spectral method and the Laguerre method [1, 2, 3, 4, 5, 6, 7]. The second approach is to reformulate original problems in unbounded domains to singular problems in bounded domains by variable transformations, and then using suitable Jacobi polynomials to approximate the resulting singular problems [8, 9, 10]. Another approach is replacing infinite domain with [ L, L] and semi-infinite interval with [0, L] by choosing L, sufficiently large. This method is named as domain truncation [11]. Another effective direct approach for solving such problems is based on rational approximations. Christov [12] and Boyd [13, 14] developed some spectral methods on unbounded intervals by using mutually orthogonal systems of rational functions. Boyd [14] defined a new spectral basis, named rational Chebyshev functions on the semi-infinite interval, by mapping to the Chebyshev polynomials. Guo et al. [15] introduced a new set of rational Legendre functions which are mutually orthogonal in L 2 (0, + ). They applied a spectral scheme using the rational Legendre functions for solving the Korteweg-de Vries equation on the half line. Boyd et al. [16] applied pseudospectral methods on a semi-infinite interval and compared rational Chebyshev, Laguerre and mapped Fourier sine. Received by the editors June 17, 2008 and, in revised form, April 7, Mathematics Subject Classification. 34B16, 34B40, 65M70. 72
2 RATIONAL CHEBYSHEV COLLOCATION TO SOLVE LANE-EMDEN EQUATIONS 73 Parand and Razzaghi [17, 18, 19] applied spectral method to solve nonlinear ordinary differential equations on semi-infinite intervals. Their approach was based on rational Tau method. They obtained the operational matrices of derivative and product of rational Chebyshev and Legendre functions and then they applied these matrices together with the Tau method to reduce the solution of these problems to the solution of a system of algebraic equations. In this paper a collocation method by rational Chebyshev functions is applied to solve Lane-Emden type equations on semi-infinite domains. Many problems in mathematical physics and astrophysics which occur on semiinfinite interval, are related to the diffusion of heat perpendicular to the parallel planes and can be modeled by the heat equation (1) x k d dx or equivalently ( x k dy dx ) + f(x)g(y) = h(x), x > 0, k > 0, (2) y + k x y + f(x)g(y) = h(x), x > 0, k > 0, where y is the temperature. For the steady-state case, and for k = 2, h(x) = 0, this equation is the generalized Emden-Fowler equation [20, 21, 22] given by (3) y + 2 x y + f(x)g(y) = 0, x > 0, subject to the conditions (4) y(0) = a, y (0) = b, where f(x) and g(y) are given functions of x and y, respectively. When f(x) = 1, (3) reduces to the Lane-Emden equation which, with specified g(y), was used to model several phenomena in mathematical physics and astrophysics such as the theory of stellar structure, the thermal behavior of a spherical cloud of gas, isothermal gas sphere and theory of thermionic currents. The case f(x) = 1 and g(y) = y m, is the standard Lane-Emden equation that corresponds to the polytropic models. This equation is one of the basic equations in the theory of stellar structure and has been the focus of many studies. This equation describes the temperature variation of a spherical gas cloud under the mutual attraction of its molecules and subject to the laws of classical thermodynamics. It also describes the variation of density as a function of the radial distance for a polytrope. It was first studied by the astrophysicists Jonathan Homer Lane and Robert Emden, which considered the thermal behavior of a spherical cloud of gas acting under the mutual attraction of its molecules and subject to the classical laws of thermodynamics. Let us consider a spherical cloud of gas and denote its hydrostatic pressure at a distance r from the center by P. Let M(r) be the mass of the sphere at radius r, Φ is the gravitational potential of the gas and g is the acceleration of the gravity. From these we have the usual equations (5) g = GM(r) r 2 = dφ dr, where G is the gravitational constant. Three conditions are assumed for the determination of Φ and P, namely, (6) dp = gρdr = ρdφ,
3 74 K. PARAND AND M. SHAHINI where ρ is the density of the gas; (7) 2 Φ = 4πGρ, in which 2 = d 2 /dr 2 + 2/r d/dr; (8) P = Kρ γ, where γ and K are empirical constants. From (6) and (8), together with the assumption that Φ = 0, when ρ = 0, it is readily found that ρ = LΦ m, where m = (γ 1) 1 and L = (m + 1)K m. By substituting this value of ρ in (7), we obtain (9) 2 Φ = a 2 Φ m, where a 2 = 4LG. If we now let Φ = Φ 0 y, where Φ 0 is the value of Φ at the center of the sphere and let x (10) r =, aφ 1 2 (m 1) 0 then we will have the Lane-Emden equation with g(y) = y m. The equation has the form (11) y + 2 x y + y m = 0, x > 0, where y and x are dimensionless variables, and m is an index related to the ratio of specific heats of the gas comprising the star. Physically interesting values of m lie in the interval [0, 5]. The conditions associated with (11) are (12) y(0) = 1, y (0) = 0, Exact solutions for Lane-Emden equation, i.e. (11) are known only for the values m = 0, 1, and 5, for other values of m the Lane-Emden equation is to be integrated numerically. The first zero of y gives the radius of the star, so y must be computed up to this zero. The gross properties of the star such as, mass, central pressure, binding energy, etc. can be computed through their relations to y. Inserting f(x) = 1 and g(y) = (y 2 C) 3 2 into (3) gives the white-dwarf equation (13) y + 2 x y + (y 2 C) 3 2 = 0, x > 0, which was introduced by [20] in his study of the gravitational potential of the degenerate white dwarf stars. The boundary conditions of this equation are the same as the Lane-Emden boundary conditions in (12). It is interesting to point out that setting C = 0 reduces the white-dwarf equation to Lane-Emden equation of index m = 3. Different approaches were used to solve Lane-Emden type equations. Bender et al. [23] proposed a perturbative technique for solving nonlinear differential equation such as Lane-Emden, consist of replacing nonlinear terms in the Lagrangian such as y n by y 1+δ and then treating δ as a small parameter. Shawagfeh [24] applied a nonperturbative approximate analytical solution for the Lane-Emden equation using the Adomian decomposition method. His solution was in the form of a power series. He used Padé approximation method to accelerate the convergence of the power series. Mandelzweig and Tabakin [25] used quasilinearization approach to solve Lane- Emden equation. This method approximates the solution of a nonlinear differential equation by treating the nonlinear terms as a perturbation about the linear ones,
4 RATIONAL CHEBYSHEV COLLOCATION TO SOLVE LANE-EMDEN EQUATIONS 75 and unlike the perturbation theories is not based on the existence of some kind of small parameters. Wazwaz [26] employed the Adomian decomposition method with an alternate framework designed to overcome the difficulty of the singular point. It was applied to the differential equations of Lane-Emden type. Further he used [27] the modified decomposition method for solving analytic treatment of nonlinear differential equations such as Lane-Emden equation. The modified method accelerates the rapid convergence of the series solution, dramatically reduces the size of work, and provides the solution by using few iterations only without any need to the so-called Adomian polynomials. Liao [28] provided an analytic algorithm for Lane-Emden type equations. This algorithm logically contains the well-known Adomian decomposition method. Different from all other analytic techniques, this algorithm itself provides us with a convenient way to adjust convergence regions even without Padé technique. By the semi-inverse method, He [29] obtained a variational principle for the Lane- Emden equation, which gives much numerical convenience when applying finite element methods or Ritz method. Parand and Razzaghi [19] presented a numerical technique to solve higher ordinary differential equations such as Lane-Emden. Their approach was based on a rational Legendre Tau method. In their work the operational matrices of the derivative and product of rational Legendre functions together with the Tau method were utilized to reduce the solution of these physical problems to the solution of systems of algebraic equations. Ramos [30, 31, 32] solved Lane-Emden equation by different methods. He presented linearization methods for singular initial-value problems in second-order ordinary differential equations such as Lane-Emden. These methods result in linear constant-coefficients ordinary differential equations which can be integrated analytically, thus they yield piecewise analytical solutions and globally smooth solutions [30]. Later, he developed piecewise-adaptive decomposition methods for the solution of nonlinear ordinary differential equations. Piecewise-decomposition methods provide series solutions in intervals which are subject to continuity conditions at the end points of each interval, and their adaption is based on the use of either a fixed number of approximants and a variable step size, a variable number of approximants and a fixed step size or a variable number of approximants and a variable step size [31]. In [32], series solutions of the Lane-Emden equation based on either a Volterra integral equation formulation or the expansion of the dependent variable in the original ordinary differential equation are presented and compared with series solutions obtained by means of integral or differential equations based on a transformation of the dependent variables. Yousefi [33] used integral operator and converted Lane-Emden equations to integral equations and then applied Legendre wavelet approximations. In his work properties of Legendre wavelet together with the Gaussian integration method are utilized to reduce the integral equations to the solution of algebraic equations. By his method the equation was formulated on [0, 1]. Chowdhury and Hashim [34] obtained analytical solutions of the generalized Emden-Fowler type equations in the second order ordinary differential equations by homotopy-perturbation method (HPM). This method is a coupling of the perturbation method and the homotopy method. The main feature of the HPM is that it deforms a difficult problem into a set of problems which are easier to solve. In
5 76 K. PARAND AND M. SHAHINI this work, HPM yields solutions in convergent series forms with easily computable terms. Aslanov [35] constructed a recurrence relation for the components of the approximate solution and investigated the convergence conditions for the Emden-Fowler type of equations. He improved the previous results on the convergence radius of the series solution. Dehghan and Shakeri [36] investigated Lane-Emden equation using the variational iteration method, and showed the efficiency and applicability of this procedure for solving this equation. Their technique did not require any discretization, linearization or small perturbations and therefore reduces the numerical computations a lot. Bataineh et al. [37] obtained analytical solutions of singular initial value problems (IVPs) of the Emden-Fowler type by the homotopy analysis method (HAM). Their solutions contained an auxiliary parameter which provided a convenient way of controlling the convergence region of the series solutions. It was shown that the solutions obtained by the Adomian decomposition method (ADM) and the homotopy-perturbation method (HPM) are only special cases of the HAM solutions. As one more step in this direction, we use rational Chebyshev collocation method to solve singular nonlinear differential equations, i.e. Lane-Emden and white-dwarf equations. The main point of our analysis lies in the fact that there is no reconstruction of the problem on the finite domain. We show that our results have good agreement with exact results, which demonstrate the viability of this approximate method. In this sense, this method has the potential to provide a wider applicability. On the other hand, the comparison of the results obtained by this method and others shows that this method provides more accurate and numerically stable solutions than those obtained by other methods. The organization of the paper is as follows: In Section 2 we explain the formulation of rational Chebyshev functions required for our subsequent development. In Section 3, after a short introduction to essentials of Lane-Emden equation, we summarize the application of rational Chebyshev functions method for solving Lane-Emden and white-dwarf equations. Then, a comparison is made with the existing methods in the literature. Section 4 is devoted to conclusions. 2. Rational Chebyshev interpolation In this section, at first, we introduced rational Chebyshev functions and expressed some basic properties of them. More, we presented approximating a function using Gauss-Radau integration and rational Chebyshev-Gauss-Radau points Rational Chebyshev functions. The well-known Chebyshev polynomials are orthogonal in the interval [-1,1] with respect to the weight function ρ(y) = 1/ 1 y 2 and can be determined with the help of the following recurrence formula: (14) T 0 (y) = 1, T 1 (y) = y, T n+1 (y) = 2yT n (y) T n 1 (y), n 1. The new basis functions, denoted by R n (x), are defined by [14, 11, 41] (15) R n (x) = T n (y) = cos(nt),
6 RATIONAL CHEBYSHEV COLLOCATION TO SOLVE LANE-EMDEN EQUATIONS 77 where L is a constant parameter and (16) (17) y = x L, y [ 1, 1], x + L ( ) x t = 2 cot 1, t [0, π]. L The constant parameter L sets the length scale of the mapping. boyd [38] offered guidelines for optimizing the map parameter L. R n (x) is the nth eigenfunction of the singular Sturm-Liouville problem x ( (18) (x + L) (x + L) xr L n(x)) + n 2 R n (x) = 0, x I = [0, ), n = 0, 1, 2,..., and by (14) satisfies in the following recurrence relation: (19) R n+1 (x) = 2 R 0 (x) = 1, R 1 (x) = x L x + L, ( ) x L R n (x) R n 1 (x), n 1. x + L 2.2. Function approximation. Let w(x) = L/( x(x + L)) denotes a nonnegative, integrable, real-valued function over the interval I = [0, ). We define (20) L 2 w(i) = {v : I R v is measurable and v w < }, where ( (21) v w = 0 ) 1 v(x) 2 2 w(x)dx, is the norm induced by the inner product of the space L 2 w(i), (22) < u, v > w = 0 u(x)v(x)w(x)dx. Thus {R n (x)} n 0 denotes a system which are mutually orthogonal under (22), i.e., (23) < R n, R m > w = c mπ 2 δ nm, with (24) c m = { 2 m = 0, 1 m 1. where δ nm is the Kronecker delta function. This system is complete in L 2 w(i). For any function u L 2 w(i) the following expansion holds (25) u(x) = with + k=0 a k R k (x), (26) a k = < u, R k > w R k 2. w The a k s are the discrete expansion coefficients associated with the family {R k (x)}.
7 78 K. PARAND AND M. SHAHINI 2.3. Rational Chebyshev interpolation approximation. Canuto et al. [39] and Gottlieb et al. [40] introduced Gauss-Radau integration. Further Guo et al. [41] introduced rational Chebyshev-Gauss-Radau points. Let (27) R N = span { R 0, R 1,..., R N }, and (28) y j = cos ( ) 2j + 1 π j = 0, 1,..., N, 2N + 1 are the N + 1 Chebyshev-Gauss-Radau points; thus we define (29) x j = L 1 + y j 1 y j j = 0, 1,..., N, which is named as rational Chebyshev-Gauss-Radau nodes. In fact, these points are zeros of the function R N+1 (x) + R N (x). Using Gauss-Radau integration we have: (30) (31) 0 u(x)w(x)dx = = 1 u(l 1 + y 1 1 y )ρ(y)dy N u(x j )w j u R 2N, j=0 where w j = 2π/(2N + 1), 0 j N 1 and w N = π/(2n + 1) are the corresponding weights with the N + 1 rational Chebyshev-Gauss-Radau nodes. In a collocation method the fundamental representation of a smooth function u on a semi-infinite interval is in terms of its values at discrete Gauss-Radau type points. Derivatives of the function are approximated by analytic derivatives of the interpolating function. The interpolating function is denoted by I N u. It is an element of R N and is defined as N (32) I N u(x) = a k R k (x), k=0 such that I N u(x j ) = u(x j ), 0 j N. I N u is the orthogonal projection of u upon R N with respect to the discrete inner product and discrete norm as: (33) < u, v > w,n = N u(x j )v(x j )w j, j=0 (34) u w,n =< u, u > 1/2 w,n. The Gauss-Radau integration formula implies that (35) < I N u, v > w,n =< u, v > w,n, u.v R 2N. 3. Illustrative examples In this section we apply the collocation method to find solutions for Lane-Emden and white-dwarf equations which satisfy the boundary conditions in (12). For a collocation approximation, by (32)-(35), we approximate y as the truncated series I N y. The solution I N y is represented by its values at the grid points x j introduced in (29). The grid values of I N y are related to the discrete expansion coefficients by (26). Thus, our goal is to find the coefficients a k, 0 k N.
8 RATIONAL CHEBYSHEV COLLOCATION TO SOLVE LANE-EMDEN EQUATIONS 79 Table 1. Comparison the first zero of y for m = 2, present method, perturbation method [23] and exact numerical values [42]. N L Present method Padé Exact value Table 2. Comparison the first zero of y for m = 3, present method, perturbation method [23] and exact numerical values [42]. N L Present method Padé Exact value Table 3. Comparison the first zero of y for m = 4, present method, perturbation method [23] and exact numerical values [42]. N L Present method Padé Exact value Solving Lane-Emden equation. The Lane-Emden equation is introduced in (11) with boundary conditions in (12). As mentioned before, it has exact solutions for m = 0, 1 and 5. Here, we solved it for m = 2, 3 and 4 by collocation method. By applying the above discussion for Lane-Emden equation and using rational Chebyshev-Gauss-Radau points in (29), we have (36) such that (37) (38) ( d 2 I N y dx di N y x dx + I Ny m) = 0, j = 0,..., N 1, x=x j I N y(x 0 ) = 1, di N y(x) dx = 0. x=x0 By omitting the first equation in (36) and with two boundary conditions in (37) and (38), we have N + 1 equations that generate a set of N + 1 nonlinear equations that can be solved by Newton method for unknown coefficients a k s. The first zero of y gives the radius of the star, so y must be computed up to this zero. The approximations of the first zero of y obtained by this method and perturbative technique [23] with exact numerical value [42] for m = 2 are listed in Table 1. Values of the first zero of y, for m = 3 and m = 4 are shown in Table 2 and Table 3, respectively. The comparison of the results with exact values, our solution is more accurate. Horedt [42] has given exact numerical values of y for some optional x. In Table 4 these values are compared with our results for m = 4 and N = 8. It shows that our results are highly accurate. The resulting graph of Lane-Emden equation for N = 7 is shown in Figure 1.
9 80 K. PARAND AND M. SHAHINI Table 4. Comparison of y(x), between present method and exact values given by Horedt [42], for m = 4. x Present method Exact value Figure 1. Lane-Emden graph obtained by present method (N = 7) for m = 2 (dotted line), m = 3 (solid line) and m = 4 (dashed line) y(x) x 3.2. Solving white-dwarf equation. Solving white-dwarf equation (13) is the same as solving the Lane-Emden equation. So we should find function I N y that satisfies ( d 2 I N y (39) dx ) di N y x dx + (I Ny 2 C) 3 2 = 0, j = 0,..., N 1, x=x j and boundary conditions (40) (41) I N y(x 0 ) = 1, di N y(x) dx = 0. x=x0 This function can be determined by its coefficients a k s. These coefficients can be found by solving the set of N + 1 nonlinear equations in (39), for j = 1,..., N 1, (40) and (41). In this paper, this equation is solved for C = 0.2, 0.4, 0.6, 0.8 and 1.0. The graph of white-dwarf equation for N = 6 and L = 0.4 is shown in Figure 2.
10 RATIONAL CHEBYSHEV COLLOCATION TO SOLVE LANE-EMDEN EQUATIONS 81 Figure 2. White-dwarf graph obtained by present method (N = 6, L = 0.4) for C = 0.2 (dashed line), C = 0.4 (dashed-dotted line), C = 0.6 (dotted line), C = 0.8 (spaced-dotted line) and C = 1.0 (solid line) y(x) x 4. Conclusion In the above discussion, we applied the collocation method to solve initial value problems, i.e. Lane-emden and white-dwarf. Lane-Emden equation occurs in the theory of stellar structure and describes the temperature variation of a spherical gas cloud. The white-dwarf equation appears in the gravitational potential of the degenerate white dwarf stars. The difficulty in this type of equations, due to the existence of singular point at x = 0, is overcomed here. In the Lane-Emden equation, the first zero of y is an important point of the function, so we have computed y to this zero. In this paper, this equation is solved for m = 2, 3 and 4, which does not have exact solutions. White-dwarf equation is solved for C = 0.2, 0.4, 0.6, 0.8 and 1.0. The validity of the method is based on the assumption that it converges by increasing the number of collocation points. Here, we tried to apply a method that give more accurate answers without reformulating the equation to bounded domains. Our goal has been achieved by approximations with a high degree of accuracy. References [1] Funaro, D., Computational aspects of pseudospectral Laguerre approximations, Appl. Numer. Math., 6 (1990) [2] Funaro, D. and Kavian, O., Computational aspects of pseudospectral Laguerre approximations, Math. Comp. 57 (1990) [3] Guo, B.Y., Error estimation of Hermite spectral method for nonlinear partial differential equations, Math. Comp., 68 (1999)
11 82 K. PARAND AND M. SHAHINI [4] Guo, B.Y. and Shen J., Laguerre-Galerkin method for nonlinear partial differential equations on a semi-infinite interval, Numer. Math., 86 (2000) [5] Maday, Y., Pernaud-Thomas, B. and Vandeven, H., Reappraisal of Laguerre type spectral methods, La Recherche Aerospatiale, 6 (1985) [6] Shen, J., Stable and efficient spectral methods in unbounded domains using Laguerre functions, SIAM J. Numer. Anal., 38 (2000) [7] Siyyam, H.I., Laguerre Tau methods for solving higher order ordinary differential equations, J. Comput. Anal. Appl., 3 (2001) [8] Guo, B.Y., Gegenbauer approximation and its applications to differential equations on the whole line, J. Math. Anal. Appl., 226 (1998) [9] Guo, B.Y., Jacobi spectral approximation and its applications to differential equations on the half line, J. Comput. Math., 18 (2000) [10] Guo, B.Y., Jacobi approximations in certain Hilbert spaces and their applications to singular differential equations, J. Math. Anal. Appl., 243 (2000) [11] Boyd, J.P., Chebyshev and Fourier Spectral Methods, second ed., Dover, New York, [12] Christov, C.I., A complete orthogonal system of functions in L 2 (, ) space. SIAM J. Appl. Math., 42 (1982) [13] Boyd, J.P., Spectral methods using rational basis functions on an infinite interval, J. Comput. Phys., 69 (1987) [14] Boyd, J.P., Orthogonal rational functions on a semi-infinite interval, J. Comput. Phys., 70 (1987) [15] Guo, B.Y., Shen, J. and Wang, Z.Q., A rational approximation and its applications to differential equations on the half line, J. Sci. Comput., 15 (2000) [16] Boyd, J.P., Rangan, C. and Bucksbaum, P.H., Pseudospectral methods on a semi-infinite interval with application to the hydrogen atom: a comparison of the mapped Fourier-sine method with Laguerre series and rational Chebyshev expansions, J. Comput. Phys., 188 (2003) [17] Parand, K. and Razzaghi, M., Rational Chebyshev Tau method for solving Volterra s population model, Appl. Math. Comput., 149 (2004) [18] Parand, K. and Razzaghi, M., Rational Chebyshev Tau method for solving higher-order ordinary differential equations, Int. J. Comput. Math., 81 (2004) [19] Parand, K. and Razzaghi, M., Rational Legendre approximation for solving some physical problems on semi-infinite intervals, Phys. Scr., 69 (2004) [20] Chandrasekhar, S., Introduction to the Study of Stellar Structure, Dover, New York, [21] Davis, H.T., 1962 Introduction to Nonlinear Differential and Integral Equations, Dover, New York, [22] Richardson, O.U., The Emission of Electricity from Hot Bodies, Longman s Green and Company, London, [23] Bender, C.M., Milton, K.A., Pinsky, S.S. and Simmons Jr. L.M., A new perturbative approach to nonlinear problems, J. Math. Phys., 30 (1989) [24] Shawagfeh, N.T., Nonperturbative approximate solution for Lane-Emden equation, J. Math. Phys., 34 (1993) [25] Mandelzweig V.B. and Tabakin F., Quasilinearization approach to nonlinear problems in physics with application to nonlinear ODEs, Comput. Phys. Commun., 141 (2001) [26] Wazwaz, A.M., A new algorithm for solving differential equations of Lane-Emden type, Appl. Math. Comput., 118 (2001) [27] Wazwaz, A.M., The modified decomposition method for analytic treatment of differential equationse, Appl. Math. Comput., 173 (2006) [28] Liao, S., A new analytic algorithm of Lane-Emden type equations, Appl. Math. Comput., 142 (2003) [29] He, J.H., Variational approach to the Lane-Emden equation, Appl. Math. Comput., 143 (2003) [30] Ramos, J.I., Linearization techniques for singular initial-value problems of ordinary differential equations, Appl. Math. Comput., 161 (2005) [31] Ramos, J.I., Piecewise-adaptive decomposition methods, Chaos Soliton. Frac., doi: /j.chaos [32] Ramos, J.I., Series approach to the Lane-Emden equation and comparison with the homotopy perturbation method, Chaos Soliton. Frac., 38 (2008) [33] Yousefi, S.A., Legendre wavelets method for solving differential equations of Lane-Emden type, Appl. Math. Comput., 181 (2006)
12 RATIONAL CHEBYSHEV COLLOCATION TO SOLVE LANE-EMDEN EQUATIONS 83 [34] Chowdhury, M.S.H. and Hashim I., Solutions of Emden-Fowler equations by homotopyperturbation method, Nonlinear Anal.-Real, doi: /j.nonrwa [35] Aslanov, A., Determination of convergence intervals of the series solutions of Emden-Fowler equations using polytropes and isothermal spheres, Phys. Lett. A, 372 (2008) [36] Dehghan, M. and Shakeri, F., Approximate solution of a differential equation arising in astrophysics using the variational iteration method, New Ast., 13 (2008) [37] Bataineh, A.S., Noorani, M.S.M. and Hashim I., Homotopy analysis method for singular IVPs of Emden-Fowler type, Nonlinear Sci. Numer. Simul., doi: /j.cnsns [38] Boyd, J.P., The optimization of convergence for Chebyshev polynomial methods in an unbounded domain, J. Comput. Phys., 45 (1982) [39] Canuto, C., Hussaini, M.Y., Quarteroni, A. and Zang, T.A., Spectral Methods in Fluid Dynamic, Springer-Verlag, New York, [40] Gottlieb, D., Hussaini, M.Y. and Orszag S., Theory and Applications of Spectral Methods in Spectral Methods for Partial Differential Equations Edited by R. Voigt, D. Gottlieb, M.Y. Hussaini, SIAM, Philadelphia, [41] Guo, B.Y., Shen, J. and Wang, Z.Q., Chebyshev rational spectral and pseudospectral methods on a semi-infinite interval, Int. J. Numer. Meth. Engng., 53 (2002) [42] Horedt, G.P., Polytropes Applications in Astrophysics and Related Fields, Kluwer Academic Publishers, Dordrecht, Department of Computer Sciences, Shahid Beheshti University, Tehran, Iran k parand@sbu.ac.ir and m.shahini@mail.sbu.ac.ir
An Elegant Perturbation Iteration Algorithm for the Lane-Emden Equation
Volume 32 - No.6, December 205 An Elegant Perturbation Iteration Algorithm for the Lane-Emden Equation M. Khalid Department of Mathematical Sciences Federal Urdu University of Arts, Sciences & Techonology
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 informationResearch Article A Nonclassical Radau Collocation Method for Nonlinear Initial-Value Problems with Applications to Lane-Emden Type Equations
Applied Mathematics Volume 2012, Article ID 103205, 13 pages doi:10.1155/2012/103205 Research Article A Nonclassical Radau Collocation Method for Nonlinear Initial-Value Problems with Applications to Lane-Emden
More informationarxiv: v2 [math-ph] 16 Aug 2010
Numerical approximations for population growth model by Rational Chebyshev and Hermite Functions collocation approach: A comparison K. Parand 1,, A. R. Rezaei, A. Taghavi Department of Computer Sciences,
More informationSolving Two Emden Fowler Type Equations of Third Order by the Variational Iteration Method
Appl. Math. Inf. Sci. 9, No. 5, 2429-2436 215 2429 Applied Mathematics & Information Sciences An International Journal http://d.doi.org/1.12785/amis/9526 Solving Two Emden Fowler Type Equations of Third
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 informationResearch Article Analysis of IVPs and BVPs on Semi-Infinite Domains via Collocation Methods
Journal of Applied Mathematics Volume 202, Article ID 696574, 2 pages doi:0.55/202/696574 Research Article Analysis of IVPs and BVPs on Semi-Infinite Domains via Collocation Methods Mohammad Maleki, Ishak
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 informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationA Computationally Hybrid Method for Solving a Famous Physical Problem on an Unbounded Domain
Commun. Theor. Phys. 71 (2019) 9 15 Vol. 71, No. 1, January 1, 2019 A Computationally Hybrid Method for Solving a Famous Physical Problem on an Unbounded Domain F. A. Parand, 1, Z. Kalantari, 2 M. Delkhosh,
More informationResearch Article A Coupled Method of Laplace Transform and Legendre Wavelets for Lane-Emden-Type Differential Equations
Journal of Applied Mathematics Volume 22, Article ID 6382, 6 pages doi:.55/22/6382 Research Article A Coupled Method of Laplace Transform and Legendre Wavelets for Lane-Emden-Type Differential Equations
More informationMethod for solving Lane-Emden type differential equations by Coupling of wavelets and Laplace transform. Jai Prakesh Jaiswal, Kailash Yadav 1
International Journal of Advances in Mathematics Volume 219, Number 1, Pages 15-26, 219 eissn 2456-698 c adv-math.com Method for solving Lane-Emden type differential equations by Coupling of wavelets and
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 New Analytic Solution to the Lane-Emden Equation of Index 2
Mathematical Problems in Engineering Volume 212, Article ID 614796, 2 pages doi:1.1155/212/614796 Research Article New Analytic Solution to the Lane-Emden Equation of Index 2 S. S. Motsa 1 and S. Shateyi
More informationEFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL EQUATIONS BASED ON THE GENERALIZED LAGUERRE POLYNOMIALS
Journal of Fractional Calculus and Applications, Vol. 3. July 212, No.13, pp. 1-14. ISSN: 29-5858. http://www.fcaj.webs.com/ EFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL
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 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 informationOn the Homotopy Perturbation Method and the Adomian Decomposition Method for Solving Abel Integral Equations of the Second Kind
Applied Mathematical Sciences, Vol. 5, 211, no. 16, 799-84 On the Homotopy Perturbation Method and the Adomian Decomposition Method for Solving Abel Integral Equations of the Second Kind A. R. Vahidi Department
More informationA Numerical Solution of Volterra s Population Growth Model Based on Hybrid Function
IT. J. BIOAUTOMATIO, 27, 2(), 9-2 A umerical Solution of Volterra s Population Growth Model Based on Hybrid Function Saeid Jahangiri, Khosrow Maleknejad *, Majid Tavassoli Kajani 2 Department of Mathematics
More informationSeries Solution of Weakly-Singular Kernel Volterra Integro-Differential Equations by the Combined Laplace-Adomian Method
Series Solution of Weakly-Singular Kernel Volterra Integro-Differential Equations by the Combined Laplace-Adomian Method By: Mohsen Soori University: Amirkabir University of Technology (Tehran Polytechnic),
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 informationComparison of homotopy analysis method and homotopy perturbation method through an evolution equation
Comparison of homotopy analysis method and homotopy perturbation method through an evolution equation Songxin Liang, David J. Jeffrey Department of Applied Mathematics, University of Western Ontario, London,
More informationarxiv: v1 [math.na] 24 Feb 2016
Newton type method for nonlinear Fredholm integral equations arxiv:16.7446v1 [math.na] 4 Feb 16 Mona Nabiei, 1 Sohrab Ali Yousefi. 1, Department of Mathematics, Shahid Beheshti University, G. C. P.O. Box
More informationhttp://www.springer.com/3-540-30725-7 Erratum Spectral Methods Fundamentals in Single Domains C. Canuto M.Y. Hussaini A. Quarteroni T.A. Zang Springer-Verlag Berlin Heidelberg 2006 Due to a technical error
More informationImproving the Accuracy of the Adomian Decomposition Method for Solving Nonlinear Equations
Applied Mathematical Sciences, Vol. 6, 2012, no. 10, 487-497 Improving the Accuracy of the Adomian Decomposition Method for Solving Nonlinear Equations A. R. Vahidi a and B. Jalalvand b (a) Department
More informationSolving Singular BVPs Ordinary Differential Equations by Modified Homotopy Perturbation Method
Journal of mathematics and computer Science 7 (23) 38-43 Solving Singular BVPs Ordinary Differential Equations by Modified Homotopy Perturbation Method Article history: Received March 23 Accepted Apri
More informationULTRASPHERICAL WAVELETS METHOD FOR SOLVING LANE-EMDEN TYPE EQUATIONS
ULTRASPHERICAL WAVELETS METHOD FOR SOLVING LANE-EMDEN TYPE EQUATIONS Y. H. YOUSSRI, W. M. ABD-ELHAMEED,, E. H. DOHA Department of Mathematics, Faculty of Science, Cairo University, Giza 63, Egypt E-mail:
More informationAn Efficient Numerical Solution of Nonlinear Hunter Saxton Equation
Commun. Theor. Phys. 67 (2017) 483 492 Vol. 67, No. 5, May 1, 2017 An Efficient Numerical Solution of Nonlinear Hunter Saxton Equation Kourosh Parand 1,2, and Mehdi Delkhosh 1 1 Department of Computer
More informationAn Implicit Method for Numerical Solution of Second Order Singular Initial Value Problems
Send Orders for Reprints to reprints@benthamscience.net The Open Mathematics Journal, 2014, 7, 1-5 1 Open Access An Implicit Method for Numerical Solution of Second Order Singular Initial Value Problems
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 informationThe Series Solution of Problems in the Calculus of Variations via the Homotopy Analysis Method
The Series Solution of Problems in the Calculus of Variations via the Homotopy Analysis Method Saeid Abbasbandy and Ahmand Shirzadi Department of Mathematics, Imam Khomeini International University, Ghazvin,
More informationNumerical Study of Generalized Three-Dimensional MHD Flow over a Porous Stretching Sheet
Journal of Applied Fluid Mechanics, Vol. 7, No. 3, pp. 473-483, 2014. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. Numerical Study of Generalized Three-Dimensional MHD Flow
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 informationarxiv: v1 [math.na] 8 Jan 2019
arxiv:190102503v1 [mathna] 8 Jan 2019 A Numerical Approach for Solving of Fractional Emden-Fowler Type Equations Josef Rebenda Zdeněk Šmarda c 2018 AIP Publishing This article may be downloaded for personal
More informationON THE EXPONENTIAL CHEBYSHEV APPROXIMATION IN UNBOUNDED DOMAINS: A COMPARISON STUDY FOR SOLVING HIGH-ORDER ORDINARY DIFFERENTIAL EQUATIONS
International Journal of Pure and Applied Mathematics Volume 105 No. 3 2015, 399-413 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v105i3.8
More informationApproximate Solution of BVPs for 4th-Order IDEs by Using RKHS Method
Applied Mathematical Sciences, Vol. 6, 01, no. 50, 453-464 Approximate Solution of BVPs for 4th-Order IDEs by Using RKHS Method Mohammed Al-Smadi Mathematics Department, Al-Qassim University, Saudi Arabia
More informationA Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions
Applied Mathematical Sciences, Vol. 5, 2011, no. 23, 1145-1152 A Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions Z. Avazzadeh
More informationAccurate spectral solutions of first- and second-order initial value problems by the ultraspherical wavelets-gauss collocation method
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 1, Issue (December 15), pp. 835 851 Applications and Applied Mathematics: An International Journal (AAM) Accurate spectral solutions
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationExponential function method for solving nonlinear ordinary differential equations with constant coefficients on a semi-infinite domain
Proc Indian Acad Sci (Math Sci) Vol 26, No, February 206, pp 79 97 c Indian Academy of Sciences Exponential function method for solving nonlinear ordinary differential equations with constant coefficients
More informationDIFFERENTIAL TRANSFORMATION METHOD FOR SOLVING DIFFERENTIAL EQUATIONS OF LANE-EMDEN TYPE
Mathematical and Computational Applications, Vol, No 3, pp 35-39, 7 Association for Scientific Research DIFFERENTIAL TRANSFORMATION METHOD FOR SOLVING DIFFERENTIAL EQUATIONS OF LANE-EMDEN TYPE Vedat Suat
More informationA NUMERICAL APPROACH TO SOLVE LANE-EMDEN-TYPE EQUATIONS BY THE FRACTIONAL ORDER OF RATIONAL BERNOULLI FUNCTIONS
A NUMERICAL APPROACH TO SOLVE LANE-EMDEN-TYPE EQUATIONS BY THE FRACTIONAL ORDER OF RATIONAL BERNOULLI FUNCTIONS K. PARAND 1,2, H. YOUSEFI 2, M. DELKHOSH 2 1 Department of Cognitive Modelling, Institute
More informationHomotopy perturbation method with Laplace Transform (LT HPM) for solving Lane Emden type differential equations (LETDEs)
DOI 10.1186/s40064-016-3487-4 RESEARCH Homotopy perturbation method with Laplace Transform LT HPM for solving Lane Emden type differential equations LETDEs Rajnee Tripathi and Hradyesh Kumar Mishra * Open
More informationSolution of an anti-symmetric quadratic nonlinear oscillator by a modified He s homotopy perturbation method
Nonlinear Analysis: Real World Applications, Vol. 10, Nº 1, 416-427 (2009) Solution of an anti-symmetric quadratic nonlinear oscillator by a modified He s homotopy perturbation method A. Beléndez, C. Pascual,
More informationApproximate Solution of an Integro-Differential Equation Arising in Oscillating Magnetic Fields Using the Differential Transformation Method
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 13, Issue 5 Ver. I1 (Sep. - Oct. 2017), PP 90-97 www.iosrjournals.org Approximate Solution of an Integro-Differential
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Mathematics 234 (21) 1364 1373 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More information1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel
1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel Xiaoyong Zhang 1, Junlin Li 2 1 Shanghai Maritime University, Shanghai,
More informationTHE DIFFERENTIAL TRANSFORMATION METHOD AND PADE APPROXIMANT FOR A FORM OF BLASIUS EQUATION. Haldun Alpaslan Peker, Onur Karaoğlu and Galip Oturanç
Mathematical and Computational Applications, Vol. 16, No., pp. 507-513, 011. Association for Scientific Research THE DIFFERENTIAL TRANSFORMATION METHOD AND PADE APPROXIMANT FOR A FORM OF BLASIUS EQUATION
More informationChapter 2 Analytical Approximation Methods
Chapter 2 Analytical Approximation Methods 2.1 Introduction As we mentioned in the previous chapter, most of the nonlinear ODEs have no explicit solutions, i.e., solutions, which are expressible in finite
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 informationExact Solutions of Fractional-Order Biological Population Model
Commun. Theor. Phys. (Beijing China) 5 (009) pp. 99 996 c Chinese Physical Society and IOP Publishing Ltd Vol. 5 No. 6 December 15 009 Exact Solutions of Fractional-Order Biological Population Model A.M.A.
More informationThe Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation
The Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation M. M. KHADER Faculty of Science, Benha University Department of Mathematics Benha EGYPT mohamedmbd@yahoo.com N. H. SWETLAM
More informationA Modification in Successive Approximation Method for Solving Nonlinear Volterra Hammerstein Integral Equations of the Second Kind
Journal of Mathematical Extension Vol. 8, No. 1, (214), 69-86 A Modification in Successive Approximation Method for Solving Nonlinear Volterra Hammerstein Integral Equations of the Second Kind Sh. Javadi
More informationThe variational iteration method for solving linear and nonlinear ODEs and scientific models with variable coefficients
Cent. Eur. J. Eng. 4 24 64-7 DOI:.2478/s353-3-4-6 Central European Journal of Engineering The variational iteration method for solving linear and nonlinear ODEs and scientific models with variable coefficients
More informationResearch Article New Algorithm for the Numerical Solutions of Nonlinear Third-Order Differential Equations Using Jacobi-Gauss Collocation Method
Mathematical Problems in Engineering Volume 20, Article ID 83728, 4 pages doi:0.55/20/83728 Research Article New Algorithm for the Numerical Solutions of Nonlinear hird-order Differential Equations Using
More informationHomotopy Perturbation Method for the Fisher s Equation and Its Generalized
ISSN 749-889 (print), 749-897 (online) International Journal of Nonlinear Science Vol.8(2009) No.4,pp.448-455 Homotopy Perturbation Method for the Fisher s Equation and Its Generalized M. Matinfar,M. Ghanbari
More informationNumerical solution of Troesch s problem using Christov rational functions
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 3, No. 4, 05, pp. 47-57 Numerical solution of Troesch s problem using Christov rational functions Abbas Saadatmandi Department
More informationThe Modified Variational Iteration Method for Solving Linear and Nonlinear Ordinary Differential Equations
Australian Journal of Basic and Applied Sciences, 5(10): 406-416, 2011 ISSN 1991-8178 The Modified Variational Iteration Method for Solving Linear and Nonlinear Ordinary Differential Equations 1 M.A. Fariborzi
More informationAn Analytic Study of the (2 + 1)-Dimensional Potential Kadomtsev-Petviashvili Equation
Adv. Theor. Appl. Mech., Vol. 3, 21, no. 11, 513-52 An Analytic Study of the (2 + 1)-Dimensional Potential Kadomtsev-Petviashvili Equation B. Batiha and K. Batiha Department of Mathematics, Faculty of
More informationNumerical Solution of Fredholm Integro-differential Equations By Using Hybrid Function Operational Matrix of Differentiation
Available online at http://ijim.srbiau.ac.ir/ Int. J. Industrial Mathematics (ISS 8-56) Vol. 9, o. 4, 7 Article ID IJIM-8, pages Research Article umerical Solution of Fredholm Integro-differential Equations
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationSolving systems of ordinary differential equations in unbounded domains by exponential Chebyshev collocation method
NTMSCI 1, No 1, 33-46 2016) 33 Journal of Abstract and Computational Mathematics http://wwwntmscicom/jacm Solving systems of ordinary differential equations in unbounded domains by exponential Chebyshev
More informationVariation of Parameters Method for Solving Fifth-Order. Boundary Value Problems
Applied Mathematics & Information Sciences 2(2) (28), 135 141 An International Journal c 28 Dixie W Publishing Corporation, U. S. A. Variation of Parameters Method for Solving Fifth-Order Boundary Value
More informationA numerical approximation to some specific nonlinear differential equations using magnus series expansion method
NTMSCI 4, No. 1, 125-129 (216) 125 New Trends in Mathematical Sciences http://dx.doi.org/1.2852/ntmsci.21611566 A numerical approximation to some specific nonlinear differential equations using magnus
More informationA Multistage Adomian Decomposition Method for Solving The Autonomous Van Der Pol System
Australian Journal of Basic Applied Sciences, 3(4): 4397-4407, 2009 ISSN 1991-8178 A Multistage Adomian Decomposition Method for Solving The Autonomous Van Der Pol System Dr. Abbas Y. AL_ Bayati Dr. Ann
More informationInstructions for Matlab Routines
Instructions for Matlab Routines August, 2011 2 A. Introduction This note is devoted to some instructions to the Matlab routines for the fundamental spectral algorithms presented in the book: Jie Shen,
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationHomotopy Perturbation Method for Computing Eigenelements of Sturm-Liouville Two Point Boundary Value Problems
Applied Mathematical Sciences, Vol 3, 2009, no 31, 1519-1524 Homotopy Perturbation Method for Computing Eigenelements of Sturm-Liouville Two Point Boundary Value Problems M A Jafari and A Aminataei Department
More informationMathematical Modeling using Partial Differential Equations (PDE s)
Mathematical Modeling using Partial Differential Equations (PDE s) 145. Physical Models: heat conduction, vibration. 146. Mathematical Models: why build them. The solution to the mathematical model will
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 informationNumerical solution of linear time delay systems using Chebyshev-tau spectral method
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 93-9466 Applications and Applied Mathematics: An International Journal (AAM) Vol., Issue (June 07), pp. 445 469 Numerical solution of linear time
More informationSpectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions
Applied Mathematical Sciences, Vol. 1, 2007, no. 5, 211-218 Spectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions M. Javidi a and A. Golbabai b a Department
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 informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationSolving Linear Time Varying Systems by Orthonormal Bernstein Polynomials
Science Journal of Applied Mathematics and Statistics 2015; 3(4): 194-198 Published online July 27, 2015 (http://www.sciencepublishinggroup.com/j/sjams) doi: 10.11648/j.sjams.20150304.15 ISSN: 2376-9491
More informationc 2003 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 41, No. 6, pp. 333 349 c 3 Society for Industrial and Applied Mathematics CONVERGENCE ANALYSIS OF SPECTRAL COLLOCATION METHODS FOR A SINGULAR DIFFERENTIAL EQUATION WEIZHANG HUANG,
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 informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
More informationOn The Exact Solution of Newell-Whitehead-Segel Equation Using the Homotopy Perturbation Method
On The Exact Solution of Newell-Whitehead-Segel Equation Using the Homotopy Perturbation Method S. Salman Nourazar, Mohsen Soori, Akbar Nazari-Golshan To cite this version: S. Salman Nourazar, Mohsen Soori,
More informationNUMERICAL SOLUTION FOR CLASS OF ONE DIMENSIONAL PARABOLIC PARTIAL INTEGRO DIFFERENTIAL EQUATIONS VIA LEGENDRE SPECTRAL COLLOCATION METHOD
Journal of Fractional Calculus and Applications, Vol. 5(3S) No., pp. 1-11. (6th. Symp. Frac. Calc. Appl. August, 01). ISSN: 090-5858. http://fcag-egypt.com/journals/jfca/ NUMERICAL SOLUTION FOR CLASS OF
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationMODIFIED LAGUERRE WAVELET BASED GALERKIN METHOD FOR FRACTIONAL AND FRACTIONAL-ORDER DELAY DIFFERENTIAL EQUATIONS
MODIFIED LAGUERRE WAVELET BASED GALERKIN METHOD FOR FRACTIONAL AND FRACTIONAL-ORDER DELAY DIFFERENTIAL EQUATIONS Aydin SECER *,Neslihan OZDEMIR Yildiz Technical University, Department of Mathematical Engineering,
More informationVARIATION OF PARAMETERS METHOD FOR SOLVING SIXTH-ORDER BOUNDARY VALUE PROBLEMS
Commun. Korean Math. Soc. 24 (29), No. 4, pp. 65 615 DOI 1.4134/CKMS.29.24.4.65 VARIATION OF PARAMETERS METHOD FOR SOLVING SIXTH-ORDER BOUNDARY VALUE PROBLEMS Syed Tauseef Mohyud-Din, Muhammad Aslam Noor,
More informationThe Homotopy Perturbation Method for Solving the Modified Korteweg-de Vries Equation
The Homotopy Perturbation Method for Solving the Modified Korteweg-de Vries Equation Ahmet Yildirim Department of Mathematics, Science Faculty, Ege University, 351 Bornova-İzmir, Turkey Reprint requests
More informationOrdinary Differential Equations
Ordinary Differential Equations Part 4 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Ordinary Differential Equations p. 1/23 Two-point Boundary Value Problems NRiC 17. BCs specified at two
More informationDiscrete Projection Methods for Integral Equations
SUB Gttttingen 7 208 427 244 98 A 5141 Discrete Projection Methods for Integral Equations M.A. Golberg & C.S. Chen TM Computational Mechanics Publications Southampton UK and Boston USA Contents Sources
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 the Coupled Klein-Gordon Schrödinger Equation Using the Modified Decomposition Method
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.4(2007) No.3,pp.227-234 Solution of the Coupled Klein-Gordon Schrödinger Equation Using the Modified Decomposition
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 informationSOLUTION OF NONLINEAR VOLTERRA-HAMMERSTEIN INTEGRAL EQUATIONS VIA SINGLE-TERM WALSH SERIES METHOD
SOLUTION OF NONLINEAR VOLTERRA-HAMMERSTEIN INTEGRAL EQUATIONS VIA SINGLE-TERM WALSH SERIES METHOD B. SEPEHRIAN AND M. RAZZAGHI Received 5 November 24 Single-term Walsh series are developed to approximate
More informationQuadrature approaches to the solution of two point boundary value problems
Quadrature approaches to the solution of two point boundary value problems Seth F. Oppenheimer Mohsen Razzaghi Department of Mathematics and Statistics Mississippi State University Drawer MA MSU, MS 39763
More informationLegendre Wavelets Based Approximation Method for Cauchy Problems
Applied Mathematical Sciences, Vol. 6, 212, no. 126, 6281-6286 Legendre Wavelets Based Approximation Method for Cauchy Problems S.G. Venkatesh a corresponding author venkamaths@gmail.com S.K. Ayyaswamy
More informationNUMERICAL SOLUTIONS OF SYSTEM OF SECOND ORDER BOUNDARY VALUE PROBLEMS USING GALERKIN METHOD
GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 37 (2017) 161-174 NUMERICAL SOLUTIONS OF SYSTEM OF SECOND ORDER BOUNDARY VALUE PROBLEMS USING GALERKIN METHOD Mahua Jahan Rupa 1 and Md. Shafiqul Islam 2
More informationTHE USE OF ADOMIAN DECOMPOSITION METHOD FOR SOLVING PROBLEMS IN CALCULUS OF VARIATIONS
THE USE OF ADOMIAN DECOMPOSITION METHOD FOR SOLVING PROBLEMS IN CALCULUS OF VARIATIONS MEHDI DEHGHAN AND MEHDI TATARI Received 16 March 25; Revised 9 August 25; Accepted 12 September 25 Dedication to Professor
More informationAND NONLINEAR SCIENCE SERIES. Partial Differential. Equations with MATLAB. Matthew P. Coleman. CRC Press J Taylor & Francis Croup
CHAPMAN & HALL/CRC APPLIED MATHEMATICS AND NONLINEAR SCIENCE SERIES An Introduction to Partial Differential Equations with MATLAB Second Edition Matthew P Coleman Fairfield University Connecticut, USA»C)
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationCURRICULUM VITAE. Ahmad Sami Bataineh. October 16, 2014
CURRICULUM VITAE Ahmad Sami Bataineh October 16, 2014 1 PERSONAL DATA Name : AHMAD SAMI BATAINEH Date & Place of Birth : 10 October 1981, Irbid, Jordan Sex : Male Marital Status : Married Nationality :
More informationA Variational Iterative Method for Solving the Linear and Nonlinear Klein-Gordon Equations
Applied Mathematical Sciences, Vol. 4, 21, no. 39, 1931-194 A Variational Iterative Method for Solving the Linear and Nonlinear Klein-Gordon Equations M. Hussain and Majid Khan Department of Sciences and
More informationWojciech Czernous PSEUDOSPECTRAL METHOD FOR SEMILINEAR PARTIAL FUNCTIONAL DIFFERENTIAL EQUATIONS
Opuscula Mathematica Vol. 30 No. 2 2010 http://dx.doi.org/10.7494/opmath.2010.30.2.133 Wojciech Czernous PSEUDOSPECTRAL METHOD FOR SEMILINEAR PARTIAL FUNCTIONAL DIFFERENTIAL EQUATIONS Abstract. We present
More information