Haar wavelet method for nonlinear integro-differential equations
|
|
- Delphia Cox
- 5 years ago
- Views:
Transcription
1 Applied Mathematics and Computation 176 (6) Haar wavelet method for nonlinear integro-differential equations Ülo Lepik Institute of Applied Mathematics, University of Tartu, Liivi Str., 549 Tartu, Estonia Abstract A numerical method for solving nonlinear integral equations based on the Haar wavelets is presented. The method is applicable for Volterra integral equations and integro-differential equations; it can be used also for solving boundary value problems of ordinary differential equations. The efficiency of the proposed method is tested with the aid of four examples. High accuracy even for a small number of collocation points is stated. Ó 5 Elsevier Inc. All rights reserved. Keywords: Integro-differential equations; Numerical solution; Nonlinear; Haar wavelets; Collocation method 1. Introduction Many problems of theoretical physics and other disciplines lead to nonlinear Volterra integral equations or integro-differential equations. For solving these equations several numerical approaches have been proposed, an overview can be found in the monograph [1]. Beginning from 1991 the wavelet method has been applied for solving integral equations, a short survey on these papers can be found in []. The solutions are often quite complicated and the advantages of the wavelet method get lost, therefore any kind of simplifications are welcome. One possibility for it is to make use of the Haar wavelets, which are mathematically the most simple wavelets. The Haar wavelet method for solving linear integral equations of different type was proposed in [] and for nonlinear Fredholm integral equations in [3]. In the present article the Haar wavelets are applied for solving of nonlinear Volterra integral equations and integro-differential equations. The proposed method is based on the collocation technique; the wavelet coefficients are calculated iteratively making use of the Newton method. The main idea of the method is to double the number of collocation points at each iteration. With minor changes the proposed method is applicable also for solving two-point boundary value problems of ordinal differential equations. The method is tested by the help of four numerical examples, for which the exact solution is known. To make the article more readable in Section a short description on the Haar wavelets is added. address: ulo.lepik@ut.ee 96-33/$ - see front matter Ó 5 Elsevier Inc. All rights reserved. doi:1.116/j.amc.5.9.1
2 Ü. Lepik / Applied Mathematics and Computation 176 (6) Haar wavelet basis The Haar wavelet family is 8 >< 1 for t ½t ð1þ ; t ðþ Þ; h i ðtþ ¼ 1 for t ½t ðþ ; t ð3þ Þ; >: elsewhere. ð:1þ Here the notations t ð1þ ¼ k m ; tðþ ¼ k þ :5 m ; tð3þ ¼ k þ 1 ð:þ m are introduced. The integer m = j, j =,1,...,J indicates the level of the wavelet; k =,1,...,m 1 is the translation parameter. The integer J determines the maximal level of resolution. The index i is calculated from the formula i = m + k + 1. Here the minimal value is i = (then m =1,k = ); maximal value is i =Mwhere M = J. The index i = 1 corresponds to the scaling function 1 for 6 t 6 1; h 1 ðtþ ¼ ð:3þ elsewhere. We shall divide the interval t [,1] into M parts of equal length Dt = 1/(M); the grid points are s l ¼ðl 1ÞDt; l ¼ 1; ;...; M þ 1. ð:4þ Since in the subsequent sections the collocation method is used then we define also the collocation points t l ¼ðl :5ÞDt; l ¼ 1; ;...; M. ð:5þ Following Chen and Hsiao [4] the Haar coefficient matrix H is introduced; it is a M M matrix with the elements H(i,l)=h i (t l ). Let us integrate (.1) q i ¼ Z t h i ðtþdt. In the collocation points (.6) gets the form Q(i,l) =q i (t l ), where Q is a M M matrix. Chen and Hsiao [4] presented this matrix in the form Q = PH, where PH is interpreted as the product of the matrices P and H. Chen and Hsiao called P the operational matrix of integration; they shoved that for calculating this matrix of order l the following matrix equation holds: " P ;5l ð 1 P l ¼ lþh # ;5l ð 1 1. ð:7þ lþh ;5l It should be noted that calculations for H l, P l must be carried out only once; the obtained results are applicable for solving whatever problems. 3. Solving nonlinear integral equations by the Haar wavelet method ð:6þ Consider integro-differential equation of the type au ðxþþbuðxþ ¼ Kðx; t; uðtþ; u ðtþþdt þ f ðxþ for 6 x 6 1 and u () = u. Here a, b are constants, the functions K, f are prescribed. Satisfying (3.1) in the collocation points (.5) we obtain au ðx l Þþbuðx l Þ¼ where l =1,,...,M. l Kðx l ; t; uðtþ; u ðtþþdt þ f ðx l Þ; ð3:1þ ð3:þ
3 36 Ü. Lepik / Applied Mathematics and Computation 176 (6) The function u (t) is developed into the Haar series u ðtþ ¼ XM a i h i ðtþ. i¼1 By integrating (3.3) we obtain uðtþ ¼ XM a i q i ðtþþuðþ. i¼1 ð3:3þ ð3:4þ If (3.3) and (3.4) are substituted into (3.) and the integrations are carried out a system of M equations a i is obtained for calculating the wavelet coefficients. Since this system is nonlinear we must use some numerical approach. In the present article for this the Newton method is applied, which brings us to the equation where X M p¼1 ah p ðx l Þþbq p ðx l Þ ¼ au ðx l Þ buðx l Þþ ¼ oa p ou l l oa p dt Da p ou þ ou ¼ oa p ou oa p ou q pðtþþ ou h pðtþ. K dt þ f ðx l Þ; l ¼ 1; ;...; M; ð3:5þ The main problem is to evaluate the integrals in (3.5). This can be done in the following way. Let us denote l l uðlþ ¼ K dt; wðp; lþ ¼ dt ð3:7þ oa p and consider the subinterval t [s s,s s+1 ], s =1,,...,M, where s s is the sth grid point defined by (.4). In each subinterval ð3:6þ h p ðtþ ¼h p ðt s Þ¼const., q p ðtþ ¼q p ðs s Þþðt s s Þh p ðt s Þ. Here t s denotes the sth collocation point defined by (.5) and t s ¼ s s þ :5Dt. Since u (t)=u (t s ) = const., we obtain uðtþ ¼uðs s Þþðt s s Þu ðt s Þ; uðt s Þ¼uðs s Þþ:5u ðt s ÞDt. Next the following notations are introduced: DGðx l ; t s ; uðt s Þ; u ðt s ÞÞ ¼ D egðx l ; t l ; uðt l Þ; u ðt l ÞÞ ¼ DG u ¼ D eg u ¼ Z ssþ1 s s Z tl s l Z ssþ1 s s Z tl s l Kðx l ; t; uðtþ; u ðt s ÞÞdt; Kðx l ; t; uðtþ; u ðt l ÞÞdt; Z ssþ1 dt ¼ oa p s s ou q pðtþþ ou h pðtþ dt; oa p dt ¼ Z tl s l ou q pðtþþ ou h pðtþ dt. ð3:8þ ð3:9þ ð3:1þ ð3:11þ
4 Evaluating the integrals (3.7) for each subinterval t [s s,s s+1 ] and summing up the results we obtain uðlþ ¼ Xl 1 DG þ D eg; wðp; lþ ¼ Xl 1 DG u þ DG e u. ð3:1þ It is convenient to put our results into the matrix form. For this purpose we introduce the row vectors x =[x(l)], a =[a i ], Da =[Da i ], u =[u(t l )], u =[u (t l )], f =[f(x l )], u =[u(l)] and M M matrices H =[h p (x l )], Q =[g p (x l )]. Besides we denote S ¼ ah þ bq w; ð3:13þ F ¼ au bu þ u þ f. The matrix form of the system (3.5) is now DaS ¼ F ; ð3:14þ which has the solution Da ¼ FS 1. ð3:15þ The Volterra integral equation uðxþ ¼ Ü. Lepik / Applied Mathematics and Computation 176 (6) Kðx; t; uðtþþdt þ f ðxþ ð3:16þ is a special case of (3.5). Our results remain applicable also for the Volterra equation if out the following changes are carried out: (i) a =,b =1, ou ¼ ; (ii) Eqs. (3.3) and (3.4) are replaced with u = ah. 4. Calculation of the wavelet coefficients The wavelet coefficients a i can be calculated in the following way. Let us assume that the problem is solved for the level m = J 1, to which correspond M = J collocation points. The wavelet coefficients for this approximation are labeled as a ðmþ i ; i ¼ 1; ;...; M. Next the value of J is increased by one, thus the number of collocation points is doubled. The wavelet coefficients at the new level are estimated as ( ^a ðmþ1þ i ¼ aðmþ i for i ¼ 1; ;...; M; ð4:1þ for i ¼ M þ 1;...; M. The estimates ^u ðmþ1þ ; ð^u Þ ðmþ1þ are calculated according to (3.3) and (3.4). The values ^a ðmþ1þ i are corrected by (3.15) and we obtain a ðmþ1þ i ¼ ^a ðmþ1þ i þ Da ðmþ1þ i ; i ¼ 1; ;...; M. ð4:þ The corrected values for u (m+1),(u ) (m+1) are calculated again from (3.3) and (3.4). This cycle is repeated until the necessary exactness of the results is obtained. Assumption (4.1) is motivated by the fact that higher coefficients of the sequence a i are usually quite small. The question how to start this procedure arises. We recommend to take the starting solution in the form (u ) () = a, u() = a x + u() and satisfy (3.1) in the collocation point x =.5. This leads to the equation aa þ bð:5a þ uðþþ ¼ Z :5 from which the coefficient a can be evaluated. Kð:5; t; a t þ uðþ; a Þdt þ f ð:5þ ð4:3þ
5 38 Ü. Lepik / Applied Mathematics and Computation 176 (6) Estimates for the next step are ^a ð1þ 1 ¼ a ; ^a ð1þ ¼ and ^u ð1þ ðtþ ¼^a ð1þ 1 h 1ðtÞþ^a ð1þ h ðtþ ¼a ; ^u ð1þ ðtþ ¼^a ð1þ 1 q 1ðtÞþ^a ð1þ q ðtþþuðþ ¼a t þ uðþ. ð4:4þ These estimates are corrected by solving (3.15) for M = 1. Our calculations show that this simple method gives good results in many cases, but there may be problems for which the sequence of iterations does not converge (the convergence of the Newton method depends upon the successful choice of the initial approximation). In such cases for getting a convergent solution we could solve (3.1) directly for two or four collocation points. The other possibility is to apply the Newton method for a fixed number of collocation points more than once and only after that double the number of collocation points. 5. Error estimates With the purpose to demonstrate the applicability and efficiency of the proposed method in the following sections some tutorial examples, for which the exact solution is known, are solved. The error of the mth iteration can be estimated as eðmþ ¼ max 16i6M juðmþ ðx i Þ u ex ðx i Þj; ð5:1þ where u ex (x) is the exact solution. For determining the convergence rate of the solution the quantity qðmþ ¼ eðm 1Þ ; m ¼ ; 3; 4;... ð5:þ eðmþ is introduced. In most cases we do not know the exact solution. Here the error of the results can be estimated in the following way. We introduce the quantity SðmÞ ¼Dx ðmþ XM u ðmþ ðx i Þ i¼1 ð5:3þ where Dx (m) denotes the stepsize of the mth approximation. Since u (m) (x) is in each subinterval x [x i,x i+1 ]a linear function then S(m) has a distinct geometrical meaning: it is the area which lies in the interval x [, 1] underneath the curve ju (m) (x)j. For estimating the exactness of the solution the quantity Sðm þ 1Þ DðmÞ ¼ 1 SðmÞ ; ð5:4þ is introduced. Geometrically D(m) denotes the relative increase of the area underneath the curve ju (m) (x)j, x [, 1] for one iteration step. The convergence rate of the process can be estimated with the aid of the function Dðm 1Þ rðmþ ¼ ; m ¼ ; 3;... DðmÞ According to the criterion (5.1) the error estimate is based on a single collocation point, therefore it is a local criterion of convergence. The estimate (5.4) makes use of all collocation points therefore it is a integral criterion. 6. Examples Let us consider some test problems.
6 Example 6.1. Consider the Volterra equation 1 þ u ðtþ uðxþ ¼ ; ð6:1þ 1 þ t which has the exact solution u ex = x. The wavelet solution is sought in the form uðtþ ¼ XM i¼1 a i h i ðtþ. Since K ¼ 1þu, it follows from (3.11) and (3.1) that: 1þt uðlþ ¼ Xl 1 ½1 þ u ðt s ÞAðsÞŠ þ 1 þ u ðt l ÞaðlÞ; wðp; lþ ¼ Xl 1 uðt s Þh p ðt s ÞAðsÞþ½1þuðt l Þh p ðt l ÞŠAðlÞ; p; l ¼ 1; ;...; M; where AðsÞ ¼arctan s sþ1 arctan s s ; AðlÞ ¼arctan t l arctan s l. We start our solution with u = a = const. Satisfying (4.3) we get the equation a a arctan :5 þ 1 ¼ ; ð6:4þ which has two roots a ð1þ ¼ :6747 and a ðþ ¼ 1:481. It follows from the calculations that the value a ðþ brings to a nonconvergent iteration process and therefore we shall take a = Next two collocation points x 1 =.5 and x =.75 are taken and the estimates for a and u are ^a ¼ða ; Þ; ^u ¼ða ; a Þ. Correcting these values with the aid of (3.15) we obtain a (1) = (.53,.33) and u (1) = (.199,.846); the error function (5.1) is e (1) =.96. Results of the following iterations are shown in Table 1. Example 6.. Solve the Volterra equation of second kind Ü. Lepik / Applied Mathematics and Computation 176 (6) uðxþ ¼ 1 uðtþuðx tþdt þ 1 sin x; < x < 1. ð6:5þ The exact solution of (6.5) is u(x) =J 1 (x), where symbol J 1 denotes the Bessel function of order 1. In the present case ð6:þ ð6:3þ K ¼ 1 uðtþuðx tþ; ¼ 1 oa p ½uðx tþh pðx tþ uðtþh p ðtþš. ð6:6þ Table 1 Error estimates e, D and convergence rates q, r for Eq. (6.1) J M e q D r 1 4.7E 7.7E 8 1.6E E E E E E E E E E 5 3.9
7 33 Ü. Lepik / Applied Mathematics and Computation 176 (6) Table Error estimates e, D and convergence rates q, r for Eq. (6.5) J M e q D r E.E 8 1.E E E E E E E E E E Making use of (3.11) and (3.1) we obtain uðlþ ¼ Dt X l 1 uðt s Þuðx l t s Þ; wðp; lþ ¼ Dt X l 1 ½uðt s Þh p ðt s Þ uðx l t s Þh p ðx l t s ÞŠ. Eq. (3.14) get the form S = H, F = u + u + f, where f(l) =.5sinx l. The starting solution is taken again in the form u = a. Satisfying (4.3) we get for a the quadratic equation a 4a þ sin :5 ¼ ; which has the solutions a ð1þ ¼ :56; a ðþ ¼ 3:744. Since the value a ðþ leads to a nonconverging iteration process we take a =.56. The estimated solution for two collocation points is ^a 1 ¼ða ; Þ; ^u ¼ða ; a Þ. Correcting it with the aid of (3.15) we obtain a = (.4,.117), u = (.137,.357) with the error e =.1. Error estimates of the subsequent approximations are shown in Table. It follows from this table that the convergence rates are smaller as in the case of Example 6.1. Example 6.3. Consider the integro-differential equation u ðxþ ¼1 þ ð6:7þ uðtþu ðtþdt; 6 x 6 1; uðþ ¼. ð6:8þ The exact solution of (6.8) is p u ex ðxþ ¼ ffiffiffi x tan p ffiffi. We shall seek the wavelet solution in the form (3.3) and (3.4); from (3.11) and (3.1) we find uðlþ ¼Dt Xl 1 wðp; lþ ¼Dt Xl 1 uðt s Þu ðt s Þþ Dt uðt lþu ðt l Þ; ½uðt s Þh p ðt s Þþu ðt s Þq p ðt s ÞŠ þ Dt ½uðt lþh p ðt l Þþu ðt l Þq k ðt l ÞŠ. ð6:9þ ð6:1þ Eqs. (3.13) and (3.14) obtain the form DaS = F, where S = H w, F = u + u + E (symbol E denotes the M dimensional unit vector). For starting we take u = a, u = a x. Satisfying (4.3) for x =.5 the equation a 8a þ 8 ¼ is obtained. It has two roots from which we shall take a = 1.17 (the other root leads to a nonconverging process). Correcting the estimates ^a ¼ða ; Þ, ^u ¼ða ; a Þ, ^u ¼ :5ða ; 3a Þ with the aid of (3.15) we find a = (1.37,.19), u = (1.47,1.46), u = (.6,.88) the error of this approximation is e =.51. Results for the next approximations are presented in Table 3. This example is taken from the article [5] by Avudainayagam and Vani, who solved it by the wavelet- Galerkin method. For computing the integrals the connection coefficient method was used. This makes the solution quite complicated, besides it is applicable only in the case of quadratic nonlinearities of the type u,
8 Ü. Lepik / Applied Mathematics and Computation 176 (6) Table 3 Error estimates e, D and convergence rates q, r for Eq. (6.8) J M e q D r E 3 4.E 8 3.3E E E E E E E E E E uu or u. The same example was solved also in [6] with the aid of the AdomianÕs decomposition method. The authors of [6] denote that their solution is more simple and easy to use. We would like to note that its accuracy of the results may be insufficient. So it follows from Table 3 of the paper [6] that u = for x = 1 while the exact value is u = 1.85; consequently, the error is 9.%. 7. Boundary value problems of ODE The method of solution proposed in this article can be applied also for solving boundary value problems of ordinary differential equations. To illustrate this let us solve the two-point boundary value problem for the differential equation u ¼ Kðx; u; u Þ ð7:1þ with the boundary conditions u() = A, u(1) = B. By integrating (7.1) we get the integro-differential equation u ðxþ ¼ Kðt; uðtþ; u ðtþþdt þ u ðþ; 6 x 6 1. ð7:þ As before we seek its solution in the form (3.3) and (3.4). Since 1 for i ¼ 1; q i ð1þ ¼ ði ¼ 1;...; MÞ for i > 1; it follows from (3.4) that a 1 = B A, consequently, the first coefficient in (3.4) is specified. According to (3.3) u ðþ ¼ XM i¼1 a i h i ðþ; where 8 >< 1 if i ¼ 1; h i ðþ ¼ 1 if i ¼ J þ 1; J ¼ ; 1; ;... >: elsewhere. Due to the fact that a 1 is fixed and Da 1 = some changes in the system of Eqs. (3.14) must be executed. Let the symbol ~a denote the vector a for which the first component is deleted. In similar way the bs is a reduced matrix in which the first row and column are deleted. Instead of (3.13) (3.15) we get now D~a e S ¼ ef ; ð7:5þ where es ¼ eh w ~ ee ~ hðþ T ; ef ¼ ~u þ ~u ðþþ~u þ f ~. Here ee is a M 1 dimensional unit vector, and - denotes the Kronecker tensor product. ð7:3þ ð7:4þ
9 33 Ü. Lepik / Applied Mathematics and Computation 176 (6) The following solution proceeds according to the algorithm presented in Sections 3 and 4: we generate the vector a =(B A,a,...,a M ) and evaluate u(x), u (x) for the next approximation from (3.3) and (3.4). Example 7.1. Consider the boundary value problem u uu 1 ¼ qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 1 p ffiffiffi ; uðþ ¼1; uð1þ ¼ ; 4 ð1 þ xþ 3 which has the exact solution p u ¼ ffiffiffiffiffiffiffiffiffiffiffi 1 þ x. By integrating (3.4) we obtain u ðxþ ¼ uu dt þ u ðþ 1 ð1 þ xþþ 1 p ffiffiffiffiffiffiffiffiffiffi. ð7:7þ 1 þ x p It follows from boundary conditions that a 1 ¼ uð1þ uðþ ¼ ffiffiffi 1. We start our solution with evaluating the quantities u(l) and w(p,l). Making use of (3.11) and (3.1) we obtain uðlþ ¼Dt Xl 1 ð7:6þ uðt s Þu ðt s Þþ Dt u ðt l Þ uðt l Þ Dt 4 u ðt l Þ ; ð7:8þ wðp; lþ ¼Dt Xl 1 ½h p ðt s Þuðt s Þþq p ðt s Þu ðt s ÞŠ þ Dt ðh pðt l Þuðt l Þþq p ðt l Þu ðt l ÞÞ Dt h p ðt l Þu ðt l Þ. ð7:9þ Since now a 1 is fixed we take the starting solution in the form u ðtþ ¼a 1 þ a h ðtþ; ð7:1þ uðtþ ¼a 1 t þ a q ðtþþuðþ. Replacing (7.1) into (7.7) and satisfying this equation in the point x =.5 we get the quadratic equation a þ ð5 þ p ffiffi p Þa þ ffiffiffi 8 þ p ffiffiffi 7 ¼. 6 This equation has two roots from which fits to us the root a =.7. The estimates for the first approximation are ^a ¼ða 1 ; a Þ¼ð:414; :7Þ, ^u ¼ð:4943; :3539Þ, ^u ¼ð1:14; 1:336Þ. Correcting the values with the aid of (7.5), (7.8) and (7.9) we find a (1) = (.44,.56), u (1) = (1.1,1.33), u (1) = (.48,.368) with the error e () =.9. Results of the subsequent approximations are presented in Table 4. Table 4 Error estimates e, D and convergence rates q, r for Eq. (7.7) J M e q D r E 3 6.5E 8 1.1E E E E E E E E E E 6 3.8
10 Ü. Lepik / Applied Mathematics and Computation 176 (6) Conclusions A new method for numerical solution of Volterra integral-equations and integro-differential equations, which is based on the Haar wavelets, is proposed. Its applicability and efficiency is checked on four test problems. The benefits of the Haar wavelet approach are sparse matrices of representation, fast transformation and possibility of implementation of fast algorithms. Simplicity of our solution is due to a great extent to the assumption (4.1), based on the fact that usually higher wavelet coefficients a i are small. It should be mentioned that such an approach is not applicable in the case of the conventional piecewise constant approximations method. It follows from Tables 1 4 that the accuracy of the obtained solutions is quite high even if the number of collocation points is small (as a rule a four or eight point solution guarantees satisfactory exactness). By increasing the number of collocation points the error of the solution rapidly decreases. Approximation with the Haar wavelets is equivalent with the approximation for piecewise constant functions. Therefore the convergence rate for piecewise constant functions, which is O(M ) can be transferred to Haar wavelet approach. As to the Newton method then in the case of sufficiently good initial values it has also quadratic convergence. So it could be expected that in the case of our solution by doubling the number of collocation points the error function roughly decreases four times. This theoretical estimation is in general consistent with the data in Tables 1, 3 and 4. Deviations appear only for small values of J 6 3, where due to small number of calculation points some adaption takes place. For bigger values of J the convergence rates q and r come nearer to the theoretical value 4. Exceptional is Example 6., where convergence rates are considerably smaller. Acknowledgement Financial support from the Estonian Science Foundation under Grant ETF-54 is gratefully acknowledged. References [1] H. Brunner, Collocation method for Volterra integral and related functional equations, Cambridge Monograph on Applied and Computational Mathematics, Cambridge University Press, Cambridge, MA, 4. [] Ü. Lepik, E. Tamme, Application of the Haar wavelets for solution of linear integral equations, in: Dynamical Systems and Applications, Antalaya, 4, Proceedings, 5, pp [3] Ü. Lepik, E. Tamme, Solution of nonlinear integral equations via Haar wavelet method, Int. J. Wavelets, Multiresol. Inform. Process., submitted for publication. [4] C.F. Chen, C.H. Hsiao, Haar wavelet method for solving lumped and distributed-parameter systems, IEE Proc. Control Theory Appl. 144 (1997) [5] A. Avudainayagam, C. Vani, Wavelet-Galerkin method for integro-differential equations, Appl. Numer. Math. 3 () [6] S.M. El-Sayed, M.R. Abdel-Azis, A comparison of AdomianÕs decomposition method and wavelet-galerkin method for solving integro-differential equations, Appl. Math. Comput. 136 (3)
Stochastic approximations of perturbed Fredholm Volterra integro-differential equation arising in mathematical neurosciences
Applied athematics and Computation 6 (27) 73 2 www.elsevier.com/locate/amc Stochastic approximations of perturbed Fredholm Volterra integro-differential equation arising in mathematical neurosciences.
More informationPAijpam.eu NUMERICAL SOLUTION OF WAVE EQUATION USING HAAR WAVELET Inderdeep Singh 1, Sangeeta Arora 2, Sheo Kumar 3
International Journal of Pure and Applied Mathematics Volume 98 No. 4 25, 457-469 ISSN: 3-88 (printed version); ISSN: 34-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/.2732/ijpam.v98i4.4
More informationTwo numerical methods for solving a backward heat conduction problem
pplied Mathematics and Computation 79 (6) 37 377 wwwelseviercom/locate/amc Two numerical methods for solving a backward heat conduction problem Xiang-Tuan Xiong, Chu-Li Fu *, hi Qian epartment of Mathematics,
More informationNumerical solutions of fourth-order Volterra integro-differential equations by the Green s function and decomposition method
Math Sci (16) 115 166 DOI 117/s46-16-1- ORIGINAL RESEARCH Numerical solutions of fourth-order Volterra integro-differential equations by the Green s function and decomposition method Randhir Singh 1 Abdul-Majid
More informationNumerical comparison of methods for solving linear differential equations of fractional order
Chaos, Solitons and Fractals 31 (27) 1248 1255 www.elsevier.com/locate/chaos Numerical comparison of methods for solving linear differential equations of fractional order Shaher Momani a, *, Zaid Odibat
More informationNumerical solution of a Fredholm integro-differential equation modelling _ h-neural networks
Available online at www.sciencedirect.com Applied Mathematics and Computation 195 (2008) 523 536 www.elsevier.com/locate/amc Numerical solution of a Fredholm integro-differential equation modelling _ h-neural
More informationNewton-homotopy analysis method for nonlinear equations
Applied Mathematics and Computation 188 (2007) 1794 1800 www.elsevier.com/locate/amc Newton-homotopy analysis method for nonlinear equations S. Abbasbandy a, *, Y. Tan b, S.J. Liao b a Department of Mathematics,
More informationExistence of positive periodic solutions for a periodic logistic equation
Applied Mathematics and Computation 139 (23) 311 321 www.elsevier.com/locate/amc Existence of positive periodic solutions for a periodic logistic equation Guihong Fan, Yongkun Li * Department of Mathematics,
More informationInitial inverse problem in heat equation with Bessel operator
International Journal of Heat and Mass Transfer 45 (22) 2959 2965 www.elsevier.com/locate/ijhmt Initial inverse problem in heat equation with Bessel operator Khalid Masood *, Salim Messaoudi, F.D. Zaman
More informationA new modified Halley method without second derivatives for nonlinear equation
Applied Mathematics and Computation 189 (2007) 1268 1273 www.elsevier.com/locate/amc A new modified Halley method without second derivatives for nonlinear equation Muhammad Aslam Noor *, Waseem Asghar
More informationShooting Methods for Numerical Solution of Stochastic Boundary-Value Problems
STOCHASTIC ANALYSIS AND APPLICATIONS Vol. 22, No. 5, pp. 1295 1314, 24 Shooting Methods for Numerical Solution of Stochastic Boundary-Value Problems Armando Arciniega and Edward Allen* Department of Mathematics
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat xxx (2009) xxx xxx Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns A one-step optimal homotopy
More informationSynchronization of an uncertain unified chaotic system via adaptive control
Chaos, Solitons and Fractals 14 (22) 643 647 www.elsevier.com/locate/chaos Synchronization of an uncertain unified chaotic system via adaptive control Shihua Chen a, Jinhu L u b, * a School of Mathematical
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationNumerical solution of hyperbolic heat conduction in thin surface layers
International Journal of Heat and Mass Transfer 50 (007) 9 www.elsevier.com/locate/ijhmt Numerical solution of hyperbolic heat conduction in thin surface layers Tzer-Ming Chen * Department of Vehicle Engineering,
More informationMathematics (JUN11MPC301) General Certificate of Education Advanced Level Examination June Unit Pure Core TOTAL
Centre Number Candidate Number For Examiner s Use Surname Other Names Candidate Signature Examiner s Initials Mathematics Unit Pure Core 3 Monday 13 June 2011 General Certificate of Education Advanced
More informationMemory Effects Due to Fractional Time Derivative and Integral Space in Diffusion Like Equation Via Haar Wavelets
Applied and Computational Mathematics 206; (4): 77-8 http://www.sciencepublishinggroup.com/j/acm doi: 0.648/j.acm.206004.2 SSN: 2328-60 (Print); SSN: 2328-63 (Online) Memory Effects Due to Fractional Time
More information"APPENDIX. Properties and Construction of the Root Loci " E-1 K ¼ 0ANDK ¼1POINTS
Appendix-E_1 5/14/29 1 "APPENDIX E Properties and Construction of the Root Loci The following properties of the root loci are useful for constructing the root loci manually and for understanding the root
More informationAn inverse vibration analysis of a tower subjected to wind drags on a shaking ground
Applied Mathematical Modelling 26 (2002) 517 528 www.elsevier.com/locate/apm An inverse vibration analysis of a tower subjected to wind drags on a shaking ground C.C. Kang a, C.Y. Lo b, * a Department
More informationNUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING HAAR WAVELET OPERATIONAL MATRIX
Palestine Journal of Mathematics Vol. 6(2) (217), 515 523 Palestine Polytechnic University-PPU 217 NUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING HAAR WAVELET OPERATIONAL MATRIX Raghvendra
More informationNumerical Solution of Fourth Order Boundary-Value Problems Using Haar Wavelets
Applied Mathematical Sciences, Vol. 5, 20, no. 63, 33-346 Numerical Solution of Fourth Order Boundary-Value Problems Using Haar Wavelets Fazal-i-Haq Department of Maths/Stats/CS Khyber Pukhtoon Khwa Agricultral
More informationJournal of Computational Physics
Journal of Computational Physics 9 (010) 3019 3045 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp BDF-like methods for nonlinear
More informationA numerical analysis of a model of growth tumor
Applied Mathematics and Computation 67 (2005) 345 354 www.elsevier.com/locate/amc A numerical analysis of a model of growth tumor Andrés Barrea a, *, Cristina Turner b a CIEM, Universidad Nacional de Córdoba,
More informationThis article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author s institution, sharing
More informationA dynamic multiscale viscosity method for the spectral approximation of conservation laws
Comput. Methods Appl. Mech. Engrg. 195 (2006) 1778 1792 www.elsevier.com/locate/cma A dynamic multiscale viscosity method for the spectral approximation of conservation laws Assad A. Oberai *, John Wanderer
More informationSOLUTION OF DIFFERENTIAL EQUATIONS BASED ON HAAR OPERATIONAL MATRIX
Palestine Journal of Mathematics Vol. 3(2) (214), 281 288 Palestine Polytechnic University-PPU 214 SOLUTION OF DIFFERENTIAL EQUATIONS BASED ON HAAR OPERATIONAL MATRIX Naresh Berwal, Dinesh Panchal and
More informationAttitude measurement system based on micro-silicon accelerometer array
Chaos, Solitons and Fractals 29 (2006) 141 147 www.elsevier.com/locate/chaos Attitude measurement system based on micro-silicon accelerometer array Li Qin *, Wendong Zhang, Huixin Zhang, Weixing Xu Key
More informationGoal-oriented, model-constrained optimization for reduction of large-scale systems
Journal of Computational Physics 224 (27) 88 896 www.elsevier.com/locate/jcp Goal-oriented, model-constrained optimization for reduction of large-scale systems T. Bui-Thanh a, K. Willcox a, *, O. Ghattas
More informationCHAPTER III HAAR WAVELET METHOD FOR SOLVING FISHER S EQUATION
CHAPTER III HAAR WAVELET METHOD FOR SOLVING FISHER S EQUATION A version of this chapter has been published as Haar Wavelet Method for solving Fisher s equation, Appl. Math.Comput.,(ELSEVIER) 211 (2009)
More informationApplied Mathematical Modelling
Applied Mathematical Modelling 33 (29) 2293 2297 Contents lists available at ScienceDirect Applied Mathematical Modelling journal homepage: www.elsevier.com/locate/apm Stability criteria for a nonlinear
More informationProperties of BPFs for Approximating the Solution of Nonlinear Fredholm Integro Differential Equation
Applied Mathematical Sciences, Vol. 6, 212, no. 32, 1563-1569 Properties of BPFs for Approximating the Solution of Nonlinear Fredholm Integro Differential Equation Ahmad Shahsavaran 1 and Abar Shahsavaran
More informationJournal of Computational Physics
Journal of Computational Physics 229 (21) 8333 8363 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp Time-dependent generalized polynomial
More informationPDE: The Method of Characteristics Page 1
PDE: The Method of Characteristics Page y u x (x, y) + x u y (x, y) =, () u(, y) = cos y 2. Solution The partial differential equation given can be rewritten as follows: u(x, y) y, x =, (2) where = / x,
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 15 (21) 3965 3973 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns A similarity solution in
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat xxx ) xxx xxx Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: wwwelseviercom/locate/cnsns A note on the use of Adomian
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationDipartimento di Matematica, Politecnico di Torino, Torino 10129, Italy PLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [Camporesi, Roberto] On: 19 May 211 Access details: Access Details: [subscription number 936142537] Publisher Taylor & Francis Informa Ltd Registered in England and Wales
More informationA modified nodal scheme for the time-dependent, incompressible Navier Stokes equations
Journal of Computational Physics 87 (2003) 68 96 www.elsevier.com/locate/jcp A modified nodal scheme for the time-dependent, incompressible Navier Stokes equations Fei Wang, Rizwan-uddin * Department of
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 informationHigher-Order Differential Equations
CHAPTER Higher-Order Differential Equations. HIGHER-ORDER HOMOGENEOUS EQUATIONS WITH CONSTANT COEFFICIENTS. ðd 7D þ 0Þy 0 has auxiliary equation m 7m þ 00 which factors into ðm Þðm 5Þ 0 and so has roots
More informationVolterra integral equations solved in Fredholm form using Walsh functions
ANZIAM J. 45 (E) ppc269 C282, 24 C269 Volterra integral equations solved in Fredholm form using Walsh functions W. F. Blyth R. L. May P. Widyaningsih (Received 8 August 23; revised 6 Jan 24) Abstract Recently
More informationFast principal component analysis using fixed-point algorithm
Pattern Recognition Letters 28 (27) 1151 1155 www.elsevier.com/locate/patrec Fast principal component analysis using fixed-point algorithm Alok Sharma *, Kuldip K. Paliwal Signal Processing Lab, Griffith
More informationSequential dynamical systems over words
Applied Mathematics and Computation 174 (2006) 500 510 www.elsevier.com/locate/amc Sequential dynamical systems over words Luis David Garcia a, Abdul Salam Jarrah b, *, Reinhard Laubenbacher b a Department
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... find the approximate solutions of derivative (first- and second-order) and antiderivative (definite integral only). Numerical Differentiation
More informationApplied Mathematics Letters. Nonlinear stability of discontinuous Galerkin methods for delay differential equations
Applied Mathematics Letters 23 21 457 461 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Nonlinear stability of discontinuous Galerkin
More informationNon-linear dynamic analysis of beams with variable stiffness
Journal of Sound and Vibration 27 (24) 847 863 JOURNAL OF SOUND AND VIBRATION www.elsevier.com/locate/jsvi Non-linear dynamic analysis of beams with variable stiffness J.T. Katsikadelis*, G.C. Tsiatas
More informationA two level finite difference scheme for one dimensional PennesÕ bioheat equation
Applied Mathematics and Computation 171 (5) 3 331 wwwelseviercom/locate/amc A two level finite difference scheme for one dimensional PennesÕ bioheat equation Jennifer J Zhao a,1, Jun Zhang b, *,, Ning
More informationFinite difference heterogeneous multi-scale method for homogenization problems
Journal of Computational Physics 191 (2003) 18 39 www.elsevier.com/locate/jcp Finite difference heterogeneous multi-scale method for homogenization problems Assyr Abdulle a, *, Weinan E b a Computational
More informationNUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS
TWMS J Pure Appl Math V5, N2, 24, pp22-228 NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS S ASADI, AH BORZABADI Abstract In this paper, Haar wavelet benefits are applied to the delay
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 informationSeries solutions of non-linear Riccati differential equations with fractional order
Available online at www.sciencedirect.com Chaos, Solitons and Fractals 40 (2009) 1 9 www.elsevier.com/locate/chaos Series solutions of non-linear Riccati differential equations with fractional order Jie
More informationApproximate step response of a nonlinear hydraulic mount using a simplified linear model
Journal of Sound and Vibration 99 (007) 656 663 Short Communication JOURNAL OF SOUND AND VIBRATION Approximate step response of a nonlinear hydraulic mount using a simplified linear model Song He, Rajendra
More informationExperimental designs for precise parameter estimation for non-linear models
Minerals Engineering 17 (2004) 431 436 This article is also available online at: www.elsevier.com/locate/mineng Experimental designs for precise parameter estimation for non-linear models Z. Xiao a, *,
More informationRobust synchronization of Sprott circuits using sliding mode control
Chaos, Solitons and Fractals 3 (6) 11 18 www.elsevier.com/locate/chaos Robust synchronization of Sprott circuits using sliding mode control David I. Rosas Almeida *, Joaquín Alvarez, Juan Gonzalo Barajas
More informationA Hypothesis-Free Multiple Scan Statistic with Variable Window
Biometrical Journal 50 (2008) 2, 299 310 DOI: 10.1002/bimj.200710412 299 A Hypothesis-Free Multiple Scan Statistic with Variable Window L. Cucala * Institut de Math matiques de Toulouse, Universit Paul
More informationA supply chain under limited-time promotion: The effect of customer sensitivity
Available online at www.sciencedirect.com European Journal of Operational Research 188 (008) 73 9 O.R. Applications A supply chain under limited-time promotion: The effect of customer sensitivity Konstantin
More informationChebyshev pseudospectral method for computing numerical solution of convection diffusion equation
Available online at wwwsciencedirectcom Applied Mathematics and Computation (8) 537 546 wwwelseviercom/locate/amc Chebyshev pseudospectral method for computing numerical solution of convection diffusion
More informationEstimation of carbon transfer coefficients using Duke Forest free-air CO 2 enrichment data
Applied Mathematics and Computation 1 () 11 1 wwwelseviercom/locate/amc Estimation of carbon transfer coefficients using Duke Forest free-air CO enrichment data Luther White a, *, Yiqi Luo b a Department
More informationFinite difference/spectral approximations for the time-fractional diffusion equation q
Journal of Computational Physics 225 (2007) 533 552 www.elsevier.com/locate/jcp Finite difference/spectral approximations for the time-fractional diffusion equation q Yumin Lin, Chuanju Xu * School of
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationFlexibility-based upper bounds on the response variability of simple beams
Comput. Methods Appl. Mech. Engrg. 94 (25) 385 44 www.elsevier.com/locate/cma Flexibility-based upper bounds on the response variability of simple beams Vissarion Papadopoulos a, *, George Deodatis b,
More informationAcommon practice for computing elementary transcendental
678 IEEE TRANSACTIONS ON COMPUTERS, VOL. 53, NO. 6, JUNE 004 Near Optimality of Chebyshev Interpolation for Elementary Function Computations Ren-Cang Li Abstract A common practice for computing an elementary
More informationTail probability of random variable and Laplace transform
Applicable Analysis Vol. 84, No. 5, May 5, 499 5 Tail probability of random variable and Laplace transform KENJI NAKAGAWA* Department of Electrical Engineering, Nagaoka University of Technology, Nagaoka,
More informationJournal of Computational Physics
Journal of Computational Physics 230 (2011) 8134 8154 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp Insights from von Neumann
More informationQuarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations
MATEMATIKA, 2011, Volume 27, Number 2, 199 208 c Department of Mathematical Sciences, UTM Quarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations 1 E. Aruchunan
More informationNew design of estimators using covariance information with uncertain observations in linear discrete-time systems
Applied Mathematics and Computation 135 (2003) 429 441 www.elsevier.com/locate/amc New design of estimators using covariance information with uncertain observations in linear discrete-time systems Seiichi
More informationEigenvalue inverse formulation for optimising vibratory behaviour of truss and continuous structures
Computers and Structures 8 () 397 43 www.elsevier.com/locate/compstruc Eigenvalue inverse formulation for optimising vibratory behaviour of truss and continuous structures H. Bahai *, K. Farahani, M.S.
More informationu x + u y = x u . u(x, 0) = e x2 The characteristics satisfy dx dt = 1, dy dt = 1
Õ 83-25 Þ ÛÐ Þ Ð ÚÔÜØ Þ ÝÒ Þ Ô ÜÞØ ¹ 3 Ñ Ð ÜÞ u x + u y = x u u(x, 0) = e x2 ÝÒ Þ Ü ÞØ º½ dt =, dt = x = t + c, y = t + c 2 We can choose c to be zero without loss of generality Note that each characteristic
More informationHAAR AND LEGENDRE WAVELETS COLLOCATION METHODS FOR THE NUMERICAL SOLUTION OF SCHRODINGER AND WAVE EQUATIONS. H. Kheiri and H.
Acta Universitatis Apulensis ISSN: 1582-5329 No. 37/214 pp. 1-14 HAAR AND LEGENDRE WAVELETS COLLOCATION METHODS FOR THE NUMERICAL SOLUTION OF SCHRODINGER AND WAVE EQUATIONS H. Kheiri and H. Ghafouri Abstract.
More informationArbitrary order Krylov deferred correction methods for differential algebraic equations
Journal of Computational Physics 221 (2007) 739 760 www.elsevier.com/locate/jcp Arbitrary order Krylov deferred correction methods for differential algebraic equations Jingfang Huang *,1, Jun Jia 1, Michael
More informationRKC time-stepping for advection diffusion reaction problems
Journal of Computational Physics 21 (24) 61 79 www.elsevier.com/locate/jcp RKC time-stepping for advection diffusion reaction problems J.G. Verwer *, B.P. Sommeijer, W. Hundsdorfer CWI, P.O. Box 9479,
More informationA Haar Wavelet Study on Convective-Radiative Fin under Continuous Motion with Temperature-Dependent Thermal Conductivity
Engineering and Physical Sciences A Haar Wavelet Study on Convective-Radiative Fin under Continuous Motion with Temperature-Dependent Thermal Conductivity Adivi Sri Venkata RAVI KANTH 1,* and Niyan UDAY
More informationApplied Mathematics and Computation
Alied Mathematics and Comutation 18 (011) 944 948 Contents lists available at ScienceDirect Alied Mathematics and Comutation journal homeage: www.elsevier.com/locate/amc A generalized Birkhoff Young Chebyshev
More informationAPPLICATION OF HAAR WAVELETS IN SOLVING NONLINEAR FRACTIONAL FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS
APPLICATION OF HAAR WAVELETS IN SOLVING NONLINEAR FRACTIONAL FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS H. SAEEDI DEPARTMENT OF MATHEMATICS, SAHID BAHONAR UNIVERSITY OF KERMAN, IRAN, 76169-14111. E-MAILS:
More informationChaos synchronization of nonlinear Bloch equations
Chaos, Solitons and Fractal7 (26) 357 361 www.elsevier.com/locate/chaos Chaos synchronization of nonlinear Bloch equations Ju H. Park * Robust Control and Nonlinear Dynamics Laboratory, Department of Electrical
More informationA note on Two-warehouse inventory model with deterioration under FIFO dispatch policy
Available online at www.sciencedirect.com European Journal of Operational Research 190 (2008) 571 577 Short Communication A note on Two-warehouse inventory model with deterioration under FIFO dispatch
More informationA Numerical Method for Linear ODEs using Non-recursive Haar Connection Coefficients
International Journal of Computational Science and Matematics ISSN 0974-389 Volume, Number 3 (00), pp 49--440 International Researc Publication House ttp://wwwirpousecom A Numerical Metod for Linear ODEs
More informationSolving Nonlinear Two-Dimensional Volterra Integral Equations of the First-kind Using the Bivariate Shifted Legendre Functions
International Journal of Mathematical Modelling & Computations Vol. 5, No. 3, Summer 215, 219-23 Solving Nonlinear Two-Dimensional Volterra Integral Equations of the First-kind Using the Bivariate Shifted
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationAn analytic approach to solve multiple solutions of a strongly nonlinear problem
Applied Mathematics and Computation 169 (2005) 854 865 www.elsevier.com/locate/amc An analytic approach to solve multiple solutions of a strongly nonlinear problem Shuicai Li, Shi-Jun Liao * School of
More informationON THE NUMERICAL SOLUTION FOR THE FRACTIONAL WAVE EQUATION USING LEGENDRE PSEUDOSPECTRAL METHOD
International Journal of Pure and Applied Mathematics Volume 84 No. 4 2013, 307-319 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v84i4.1
More informationResearch Article Assessment of Haar Wavelet-Quasilinearization Technique in Heat Convection-Radiation Equations
Applied Computational Intelligence and So Computing, Article ID 45423, 5 pages http://.doi.org/.55/24/45423 Research Article Assessment of Haar Wavelet-Quasilinearization Technique in Heat Convection-Radiation
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 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 informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 6 (2) 63 75 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns The scaled boundary FEM for nonlinear
More informationGalerkin method for the numerical solution of the RLW equation using quintic B-splines
Journal of Computational and Applied Mathematics 19 (26) 532 547 www.elsevier.com/locate/cam Galerkin method for the numerical solution of the RLW equation using quintic B-splines İdris Dağ a,, Bülent
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 information(Received 13 December 2011, accepted 27 December 2012) y(x) Y (k) = 1 [ d k ] dx k. x=0. y(x) = x k Y (k), (2) k=0. [ d k ] y(x) x k k!
ISSN 749-3889 (print), 749-3897 (online) International Journal of Nonlinear Science Vol.6(23) No.,pp.87-9 Solving a Class of Volterra Integral Equation Systems by the Differential Transform Method Ercan
More informationA multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system
Chaos, Solitons and Fractals 30 (006) 197 03 www.elsevier.com/locate/chaos A multiple Riccati equations rational expansion method and novel solutions of the Broer Kaup Kupershmidt system Qi Wang a,c, *,
More informationJournal of Computational Physics
Journal of Computational Physics 8 (009) 5989 605 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp High-order compact schemes for
More informationNONSTATIONARY DISCRETE CHOICE: A CORRIGENDUM AND ADDENDUM. PETER C. B. PHILIPS, SAINAN JIN and LING HU COWLES FOUNDATION PAPER NO.
NONSTATIONARY DISCRETE CHOICE: A CORRIGENDUM AND ADDENDUM BY PETER C. B. PHILIPS, SAINAN JIN and LING HU COWLES FOUNDATION PAPER NO. 24 COWLES FOUNDATION FOR RESEARCH IN ECONOMICS YALE UNIVERSITY Box 2828
More informationSynchronization of self-sustained oscillators by common white noise
Physica A 351 (5) 16 13 www.elsevier.com/locate/physa Synchronization of self-sustained oscillators by common white noise D.S. Goldobin a,b, A.S. Pikovsky a, a Department of Physics, University of Potsdam,
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 informationApplied Mathematics and Computation
Applied Mathematics and Computation 245 (2014) 86 107 Contents lists available at ScienceDirect Applied Mathematics and Computation journal homepage: www.elsevier.com/locate/amc An analysis of a new family
More informationImproved Lebesgue constants on the triangle
Journal of Computational Physics 207 (2005) 625 638 www.elsevier.com/locate/jcp Improved Lebesgue constants on the triangle Wilhelm Heinrichs Universität Duisburg-Essen, Ingenieurmathematik (FB 10), Universitätsstrasse
More informationA null space method for solving system of equations q
Applied Mathematics and Computation 149 (004) 15 6 www.elsevier.com/locate/amc A null space method for solving system of equations q Pu-yan Nie 1 Department of Mathematics, Jinan University, Guangzhou
More informationLie group variational integrators for the full body problem
Comput. Methods Appl. Mech. Engrg. 96 (2007) 2907 2924 www.elsevier.com/locate/cma Lie group variational integrators for the full body problem Taeyoung Lee a, *,,2, Melvin Leok b,, N. Harris McClamroch
More informationDirect method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
Journal of Computational and Applied Mathematics 22 (28) 51 57 wwwelseviercom/locate/cam Direct method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions
More informationTHE ALGORITHM TO CALCULATE THE PERIOD MATRIX OF THE CURVE x m þ y n ¼ 1
TSUKUA J MATH Vol 26 No (22), 5 37 THE ALGORITHM TO ALULATE THE PERIOD MATRIX OF THE URVE x m þ y n ¼ y Abstract We show how to take a canonical homology basis and a basis of the space of holomorphic -forms
More informationHow do our representations change if we select another basis?
CHAPTER 6 Linear Mappings and Matrices 99 THEOREM 6.: For any linear operators F; G AðV Þ, 6. Change of Basis mðg FÞ ¼ mðgþmðfþ or ½G FŠ ¼ ½GŠ½FŠ (Here G F denotes the composition of the maps G and F.)
More information