Approximate solution of linear integro-differential equations by using modified Taylor expansion method
|
|
- Lynne McDowell
- 5 years ago
- Views:
Transcription
1 ISSN , Engl, UK World Journal of Modelling Simulation Vol 9 (213) No 4, pp Approximate solution of linear integro-differential equations by using modified Taylor expansion method Jalil Rashidinia 1, Ali Tahmasebi 2 1 School of Mathematics, Iran University of Science & Technology, Narmak 16844, Tehran, Iran 2 School of Mathematics Computer Science, Damghan University, Damghan, Iran (Received December 5 212, Revised July , Accepted September ) Abstract In this study we developed modified Taylor expansion method for approximating the solution of linear Fredholm Volterra integro-differential equations Via Taylor s expansion of the unknown function at an arbitrary point, the integro-differential equations to be solved is approximately transformed into a system of linear equations for the unknown its derivatives which can be dealt with in an easy way This method gives a simple closed form solution for a linear integro-differential equation This method also enable us to control truncation error by adjusting the step size used in the numerical scheme Some numerical examples are provided to illustrate the accuracy of our approach Keywords: linear integro-differential equations, approximate solution, taylor expansion method 1 Introduction Mathematical modeling of real-life problems usually results in functional equations, eg partial differential equations, integral integro-differential equations, stochastic equations others In particular, integro-differential equations arise in fluid dynamics, biological models chemical kinetics The analytical solutions of some integro-differential equations cannot be found, thus numerical methods are required In recent years, the numerical methods for linear integro-differential equations have been extensively studied by many authors The iterated Galerkin methods have been proposed in [1] Furthermore, mixed interpolation collocation methods for first- second-order Volterra linear integro-differential equations have been suggested in [2] Hu [3] developed the interpolation collocation method for solving Fredholm linear integro-differential equations Rashed treated a special type of integro-differential equations with derivatives appearing in integrals [4] Taylor series expansion method is a powerful technique for solving integro-differential equations The Taylor expansion approach for solving Volterra integral equations has been presented by Kanwal Liu [5] The Taylor expansion method for approximate solution of a class of linear integro-differential equations presented by Li [6] also, a Taylor series method [7 1, 22] was presented by Szer et al for differential, integral integro-differential equations Bessel or Chebyshev polynomial approach are used to solve integro-differential equations in [11 16] Sahin in [17] applied collection method for numerical solution of Volterra integral equations with variable coefficients An improved Legre method [18] is used for solution of integro-differential equations A numerical method for the solution of integro-differential-difference equation with variable coefficients has been proposed by Saadatmi [19] In this paper, the Taylors expansion method [6] is modified such that by this modification the corresponding accuracy is drastically improved Also, the approximate method described here is very easy to implement Corresponding author Tel: address: rashidinia@iustacir Published by World Academic Press, World Academic Union
2 29 J Rashidinia & A Tahmasebi: Approximate solution of linear integro-differential equations 2 Volterra integro-differential equations Consider the Volterra integro-differential equations of the form x β i (x)y (i) (x) = f(x) + λ v k v (x, t)y(t)dt, x, t I = [a, b], (1) under the initial conditions η c ji y (i) (η) = µ j, j =, 1,, m 1, (2) y(x) is an unknown function to be determined, the known continuous functions β i (x), f(x) k v (x, t) are defined on the interval I Here the real coefficients λ v c ji are known constants η is a given point in the spatial domain of the problem Remark 1 The following partition is used throughout of the paper = {a = x θl < < x 1 < x = η < x 1 < < x θr = b}, is an equidistance partition on I = [a, b], h is a step size of the partition x θ = η + θh, for θ { θ r,, 1,, 1,, θ r } Theorem 1 Consider linear integro-differential Eq (1)- Eq (2) Let the functions β i (x), f(x), k v (x, t) the exact solution y(x) are sufficiently differentiable on the interval I k v (x, t) be a separable kernel Then there exist linear independent functions ϕ i (χ) constants c i for i =, 1,, α, such that y(χ) = y(χ+x θ ) is an exact solution of β i (χ + x θ )y (i) (χ) = f(χ + x θ ) + α c i ϕ i (χ) + λ v k v (χ + x θ, + x θ )y()d, (3) with initial conditions x = χ + x θ y (i) () = y (i) (x θ ), i =, 1,, m 1, (4) Proof The kernel in Eq (1) is separable, then we have k(x, t) = p u i (x)v i (t), (5) by substituting Eq (5) into Eq (3) then change of variables x = χ + x θ, t = + x θ, we have β i (χ + x θ )y (i) (χ + x θ ) = f(χ + x θ ) + d i = η x θ v i ( + x θ )y( + x θ )d p d i u i (χ + x θ ) + λ v k v (χ + x θ, + x θ )y( + x θ )d, Now, by simplifying classifying on the middle terms of RHS Eq (6), we can write (6) p d i u i (χ + x θ ) = α c i ϕ i (χ), (7) WJMS for contribution: submit@wjmsorguk
3 World Journal of Modelling Simulation, Vol 9 (213) No 4, pp ϕ i (χ), i =, 1,, α are known linear independent functions c i, i =, 1,, α are unknown constants Using Eq (7), by denoting y(χ) = y(χ + x θ ), we can rewrite Eq (6) Eq (2), in the form of Eq (3) Eq (4), respectively To obtain c i, by derivation on Eq (3) for α 1 times at χ =, the following algebraic system is obtained α c i ϕ i () = m g i () f(x θ ), α c i ϕ i () = m g i () f (x θ ) ρ(), α c i ϕ (α 1) i () = m g (α 1) i () f (α 1) (x θ ) ρ (α 2) (), g i (χ) = β i (χ + x θ )y (i) (χ) ρ(χ) = k v (χ + x θ, χ + x θ )y(χ) By solving algebraic system Eq (8), the unknown constants can be determined as follows c i = τ i (x θ, y(), y (),, y (m+α 1) ()), i =, 1,, α, (9) τ i are known functions y (i) () = y (i) (x θ ), i =, 1,, α For solving linear integro-differential equations Eq (1) Eq (2), by modified Taylor expansion method, we assume that y() in Eq (3) can be represented in terms of the n th -order Taylor expansion as y() = y(χ) + y (χ)( χ) + + y(n) (χ) ( χ) n + y(n+1) (ζ) n! (n + 1)! ( χ)n+1, (1) ζ is between χ It is readily shown that the Lagrange remainder y (n+1) (ζ) ( χ)n (n + 1)! is sufficiently small for a large enough n provided that y (n+1) (χ) is a uniformly bounded function for any n on the interval I, so that we neglect the remainder the truncated Taylor expansions y(χ) as y() n k= (8) y (k) (χ) ( χ) k (11) k! It is worth noting that the Lagrange remainder vanishes for a polynomial of degree equal to or less than n, implying that the above n th -order Taylor expansion is exact Also, n is chosen any positive integer such that n m Now, we substitute Taylor s expansion Eq (11) for y() into Eq (3) find k (χ)y(χ) + k 1 (χ)y (χ) + + k n (χ)y (n) (χ) = f (χ), (12) { βi (χ + x k i = θ ) λv χ k v(χ + x θ, + x θ )( χ) i d, λv χ k v(χ + x θ, + x θ )( χ) i d, i =, 1,, m i = m + 1, m + 2,, n, (13) f (χ) = f(χ + x θ ) + α c i ϕ i (χ) (14) Thus Eq (12) becomes a m th -order, linear, ordinary differential equation with variable coefficients for y(χ) its derivations up to n Instead of solving analytically ordinary differential equation, we will determine y(χ), y (χ),, y (n) (χ) by solving linear equations To this end, other n independent linear equations WJMS for subscription: info@wjmsorguk
4 292 J Rashidinia & A Tahmasebi: Approximate solution of linear integro-differential equations for y(χ), y (χ),, y (n) (χ) are needed This can be achieved by integrating both sides of Eq (3) with respect to χ from to χ changing the order of the integrations Accordingly, one can get that β i ( + x θ )y (i) ()d = f ()d + λ v k v (s + x θ, + x θ )dsy()d, (15) we have replaced variable s to χ for convenience Applying integration by part, accordingly, one can get that Thus β i ( + x θ )y (i) ()d =β i ( + x θ )y (i 1) () χ β i( + x θ )y (i 2) () χ + + ( 1) i 1 β (i 1) i β i ( + x θ )y (i) ()d = ( + ( + x θ )y() χ + ( 1)i β (i) i ( + x θ )y()d ( 1) k i 1 β (k i 1) k ( + x θ ))y (i) ()) χ ( 1) k β (k) k ( + x θ)y()d, (16) Applying the Taylor expansion again substituting Eq (11) for y() into Eq (16) gives k 1i (χ) = k= k 1 (χ)y(χ) + k 11 (χ)y (χ) + + k 1n (χ)y (n) (χ) = f 1 (χ), (17) ( 1) k i 1 β (k i 1) k (χ + x θ ) + 1 λ v ( 1) k β (k) ( + x θ)( χ) i d k= k v (s + x θ, + x θ )ds( χ) i d, i =, 1,, m 1, (18) k k 1i (χ) = 1 λ v ( 1) k β (k) ( + x θ)( χ) i d k= k k v (s + x θ, + x θ )ds( χ) i d, i = m, m + 1,, n, f 1 (χ) = ( m ) ( 1) k i 1 β (k i 1) k (x θ ) y (i) () + f ()d (19) By repeating the integration on both side of Eq (15) for r times r = 2, 3,, n, one can arrive that (χ ) r 1 β i ( + x θ )y (i) ()d = By using integration by part, we can obtain (χ ) r 1 f ()d + λ v k v (s + x θ, + x θ )(χ s) r 1 dsy()d (2) WJMS for contribution: submit@wjmsorguk
5 World Journal of Modelling Simulation, Vol 9 (213) No 4, pp (χ ) r 1 β i ( + x θ )y (i) ()d = ( 1) k k β k ( + x θ ) (k) (χ ) r 1 k y (k) ()d k= ( m + ( 1) k i 1 k i 1 β k ( + x θ )(χ ) r 1 ) k i 1 y() χ (21) k ri (χ) = Substituting Eq (11) for y() in the Eq (21), we obtain λ v k r (χ)y(χ) + k r1 (χ)y (χ) + + k rn (χ)y (n) (χ) = f r (χ), (22) ( 1) k i 1 k i 1 β k ( + x θ )(χ ) r 1 k i 1 =χ + 1 ( 1) k k β k ( + x θ )( χ) r 1 k d k v (s + x θ, + x θ )(χ s) r 1 ds( χ) i d, r = 1,, m 1, (23) k= f r (χ) = k ri (χ) = 1 λ v ( m ( 1) k k β k ( + x θ )( χ) r 1 k d k= k v (s + x θ, + x θ )(χ s) r 1 ds( χ) i d, i = m, m + 1,, n ( 1) k i 1 k i 1 β k ( + x θ )(χ ) r 1 k i 1 ) y() =η + (χ ) r 1 f ()d (24) Therefore, Eqs (12), (17) (22) form a system of n + 1 linear equation for n + 1 unknown functions y(χ), y (χ),, y (n) (χ) For simplicity, this system can be rewritten in an alternative compact form as K nn (χ)y n (χ) = F(χ), (25) K nn (χ) is an (n + 1) (n + 1) square matrix function, Y n (χ) F n (χ) are two vectors of length n + 1, these are defined as k (χ) k 1 (χ) k n (χ) k 1 (χ) k 11 (χ) k 1n (χ) K nn (χ) =, (26) k n (χ) k n1 (χ) k nn (χ) Y n (χ) = y(χ) y (χ) y (n) (χ), F n(χ) = Therefore, by well-known Crammer s formula, we can obtain f (χ) f 1 (χ) f n (χ) (27) y (i) (χ) = Ψ i (χ, x θ, c, c 1,, c α, y(), y (),, y () ()), i =, 1,, m + α 1, (28) which is Taylor expansion of solution Eq (3) Eq (4) for y(χ) its derivations up to m + α 1 at χ =, by change of variable χ = x x θ, we have WJMS for subscription: info@wjmsorguk
6 294 J Rashidinia & A Tahmasebi: Approximate solution of linear integro-differential equations y (i) (x) = Ψ i (x x θ, x θ, c, c 1,, c α, y(x θ ), y (x θ ),, y () (x θ )), i =, 1,, m + α 1, (29) which is Taylor expansion of solution Eq (1) Eq (2) for y(x) it s derivations up to m + α 1 at x = x θ Note that, in the Eq (28) Eq (29), in the any step of proposed method c i, i =, 1,, α, y (i) (), i =, 1,, m 1 are known values by using (9) 21 Algorithm of the approach In this section, we try to propose an algorithm on the basis of the above discussions suppose that we face with the linear Volterra integro-differential Eq (1) Eq (2), its kernel satisfies the conditions of Theorem 1 This algorithm is presented in two stages such as initialization step main steps Initialization step: Choose step size h > for equidistance partition on I Set c j =, j =, 1,, α y (j) () = y (j) (η), j =, 1,, m 1 Set θ = go to main steps Main steps: Step 1 Compute by Eq (29) the following approximate solution y (i) (x) = Ψ i (x x θ, x θ, c, c 1,, c α, y(), y (),, y () ()), i =, 1,, m + α 1, (3) which is Taylor expansion approximate Eq (1)- Eq (2) for y(x) it s derivations up to m+α 1 at x = x θ Go to Step 2 Step 2 Set θ = θ + 1, if θ > θ r 1, stop; Otherwise, using Eq (3) compute the approximate values y (i) () = y (i) (x θ ), i =, 1,, m + α 1, (31) which are the initial conditions Eq (1) Eq (2) at x = x θ Go to next Step Step 3 By Eq (31) compute the following approximate values c i = τ i (x θ, y(), y (),, y (m+α 1) ()), i =, 1,, α, (32) go to Step 1 Note that in the Step 2, if we replaced condition θ > θ l 1, with the condition θ > θ l + 1 also θ 1 to θ + 1, then by applying the above algorithm we obtain Taylor expansion approximate solution of Eq (1)- Eq (2) its derivations at x = x θ, θ = 1, 2,, θ l 3 Fredholm integro-differential equations In this section, our study focuses on a class of linear Fredholm integro-differential equations of the following form b β i (x)y (i) (x) = f(x) + λ f k f (x, t)y(t)dt, x, t I = [a, b], (33) subject to the initial conditions a c ji y (i) () = µ j, a η b, j =, 1,, m 1 (34) We assume that the desired solution y(t) is n + 1 times continually differentiable on the interval I, i e y C n+1 [I] Consequently, for y C n+1, y(t) can be represented in terms of the truncated Taylor series expansion as n y (k) (h) y(t) (t h) k, (35) k! k= WJMS for contribution: submit@wjmsorguk
7 World Journal of Modelling Simulation, Vol 9 (213) No 4, pp h is step size in the partition Substituting the approximate solution Eq (35) for y(t) into Eq (33), one can get a y(h) + a 1 y (h) + + a n y (n) (h) = f(h), (36) a i = β i (h) λ f λ f b a k f (h, t)(t h) i dt, i =, 1,, m, b a k f (h, t)(t h) i dt, i = m + 1, m + 2,, n (37) Using a procedure analogous to the previous section, we can obtain other n linear equations for y(h), y (h),, y (n) (h) as follows a ri = h h λ f a r y(h) + a r1 y (h) + + a rn y (n) (h) = f r (h), (38) ( 1) k i 1 k i 1 β k ()(h ) r 1 k i 1 =h + 1 η h η ( 1) k k (β k ()( h) r 1 ) k d k= k f (s, )(h s) r 1 ds( h) i d, r = 1,, m 1, (39) a ri = 1 h η λ f ( 1) k k (β k ()( h) r 1 ) k d k= h h η k f (s, )(h s) r 1 ds( h) i d, i = m, m + 1,, n, f r (h) = ( m ) ( 1) k i 1 k i 1 β k ()(h ) r 1 y() =η + k i 1 h η (h ) r 1 f()d (4) Therefore Eq (36) Eq (38) form a system of n + 1 linear equation for n + 1 unknown values y(h), y (h),, y (n) (h) By using Crammer s formula we can obtain is initial conditions Eq (1) Eq (2) at x = h Now let q k f = u i (x)v i (t), y i = y (i) (h), i =, 1,, q 1, q < n, (41) by using theorem 1, linear integro-differential equations Eq (1) converted to the following form with initial conditions β i (χ + x θ )y (i) (χ) = f(χ + x θ ) + y(χ) = y(χ + x θ ), x = χ + x θ q I i u i (χ + x θ ), (42) y (i) () = y (i) (x θ ), i =, 1,, m 1, (43) WJMS for subscription: info@wjmsorguk
8 296 J Rashidinia & A Tahmasebi: Approximate solution of linear integro-differential equations b I i = λ f v i (t)y(t)dt a is an unknown constant From (42) the following algebraic equation cab be drived b (χ)y(χ) + b 1 (χ)y (χ) + + b n (χ)y (n) (χ) = g (χ), (44) b i (χ) = { βi (χ + x θ ), i =, 1,, m,, i = m + 1, m + 2,, n, (45) g (χ) = f(χ + x θ ) + q I i u i (χ + x θ ) (46) Also by using a procedure analogous to the section 2, we can obtain other n linear equations for y(χ), y (χ),, y (n) (χ) as follows b r (χ)y(χ) + b r1 (χ)y (χ) + + b rn (χ)y (n) (χ) = g r (χ), (47) b ri (χ) = ( 1) k i 1 k i 1 β k ( + x θ )(χ ) r k i 1 =χ ( 1) k k β k ( + x θ )( χ) r 1 k d, r = 1,, m 1, (48) k= g r (χ) = b ri (χ) = 1 ( m ( 1) k k β k ( + x θ )( χ) r 1 k d, i = m, m + 1,, n, k= ( 1) k i 1 k i 1 β k ( + x θ )(χ ) r 1 k i 1 ) y() =η + (χ ) r 1 g ()d (49) Therefore Eq (44) Eq (47) form a system of n + 1 linear equation for n + 1 unknown functions y(χ), y (χ),, y (n) (χ) For simplicity, this system can be rewritten in an alternative compact form as B nn (χ)y n (χ) = F(χ), (5) B nn (χ) is an (n + 1) (n + 1) square matrix function, Y n (χ) F n (χ) are two vectors of length n + 1, these are defined as b (χ) b 1 (χ) b n (χ) b 1 (χ) b 11 (χ) b 1n (χ) B nn (χ) =, (51) b n (χ) b n1 (χ) b nn (χ) Y n (χ) = y(χ) y (χ) y (n) (χ), F n(χ) = g (χ) g 1 (χ) g n (χ) (52) WJMS for contribution: submit@wjmsorguk
9 World Journal of Modelling Simulation, Vol 9 (213) No 4, pp By well-known Crammer s formula, we have y (i) (χ) = Ψ i (χ, x θ, I, I 1,, I q, y(), y (),, y () ()), i =, 1,, m 1, (53) y (i) () = y (i) (x θ ), i =, 1,, m 1, {I k } q k= are unknown parameters will be determined Using Eq (41) Eq (53) the following algebraic system is yield Ψ i (, η + h, I, I 1,, I q, y(), y (),, y () ()) = y i, i =, 1,, q 1, (54) y (i) () = y (i) (x θ ), i =, 1,, m 1 By solving algebraic system Eq (54), the unknown constants can be determined as follows By substituting, the obtained values I k into Eq (53), we can obtain I k = γ k, k =, 1,, q (55) y (i) (χ) = Ψ i (χ, x θ, γ, γ 1,, γ q, y(), y (),, y () ()), i =, 1,, m 1, (56) which is Taylor expansion approximate solution Eq (42) Eq (43) for y(χ) its derivations up to m + α 1 at χ =, by change of variable χ = x x θ, we have y (i) (x) = Ψ i (x x θ, x θ, γ, γ 1,, γ q, y(), y (),, y () ()), i =, 1,, m 1, (57) which is Taylor expansion approximate solution Eq (33) Eq (34) for y(x) it s derivations up to m + α 1 at x = x θ 31 Algorithm of the approach In this section, we propose an algorithm on the basis of the above discussions This algorithm is presented in two stages such as initialization step main steps Initialization step: Choose step size h > for equidistance partition on I Set y (j) () = y (j) (), j =, 1,, m 1 Set θ = go to main steps Main steps: Step 1 By Eq (53) Eq (55), compute the following approximate solution y (i) (x) = Ψ i (x x θ, x θ, γ, γ 1,, γ q, y(), y (),, y () ()), i =, 1,, m 1, (58) which is Taylor expansion approximate Eq (33) Eq (34) for y(x) it s derivations up to m 1 at x = x θ Go to next Step Step 2 Set θ = θ + 1, if θ > θ r 1, stop; Otherwise, using Eq (58) compute the approximate values y (i) () = y (i) (x θ ), i =, 1,, m 1 (59) which are the initial conditions Eq (33) Eq (34) at x = x θ go to Step 1 Note that in the Step 2, if we replaced condition θ > θ l 1, with the condition θ > θ l + 1 also θ 1 to θ + 1, then by applying the above algorithm we obtain Taylor expansion approximate solution of Eq (33) Eq (34) its derivations at x = x θ, θ = 1, 2,, θ l WJMS for subscription: info@wjmsorguk
10 298 J Rashidinia & A Tahmasebi: Approximate solution of linear integro-differential equations 4 Numerical examples To illustrate the effectiveness of the proposed method we consider several test examples to carried out the applicability of our approach Example 1 Consider the first order linear integro-differential equation [22] y = 1 + λ x y(t)dt, 1 x, t 1 (6) We consider this Example in the two different cases: one is that λ = 1, y() = corresponding exact solutions y(x) = sin x other is λ = 1, y() = 1 corresponding exact solutions y(x) = e x In the case of λ = 1, y() =, we explain the steps of our method for n = 3 h = 1 as follows By using theorem 1 linear integro-differential equations convert to the following form with initial condition y = 1 + c y()d, (61) y () = y(x θ ), (62) c = 1 y () Let n = 4, in Eq (11), by solving corresponding algebraic system Eq (25) by using Eq (28), we can obtain ( χ 9 y(χ) = c χ9 144 χ c χ χ c χ χ c χ χ8 y() χ1 y() χ12 y() χ14 y() ) ( χ χ χ χ ( χ y 8 (χ) = c χ8 97 χ1 97c χ1 59 χ ) 1 (63) 59c χ χ c χ χ9 y() χ11 y() 112 ) ( χ13 y() χ χ χ χ 16 ) 1 (64) In the initialization step of algorithms (21), by setting c =, θ =, y() = into Eq (63), Eq (64), by using Eq (29), we have ( x 9 y(x) = 144 x x x 15 ) ( x x x x 16 ) 1 (65) ( x y 8 97 x1 (x) = ( x 8 ) 59 x x x x x ) 1, (66) WJMS for contribution: submit@wjmsorguk
11 World Journal of Modelling Simulation, Vol 9 (213) No 4, pp which is Taylor expansion approximation solution Eq (1) Eq (2) its derivative at x = In the next step, setting x = h, (θ = 1), into Eq (65) Eq (66), the following approximate values can be obtain y() = , y () = 9954, c = , y() = y(1) y () = y (1) Now again, by setting the obtained values c, y(), y () θ = 1, into Eq (63), Eq (64), by using Eq (29), Taylor approximation solution Eq (1) Eq (2) its derivative at x = h can be obtained Repeating the above procedure for θ = 2, ±3,, ±9, we can obtain Taylor approximate solution Eq (6) at x = θh in the case of λ = 1, y() = Consequently, by using similarly procedure, we can obtain the approximate solution Eq (6) for λ = 1, y() = 1 in the interval [ 1, 1] Also we employ the Taylor series method [6] to determine the three- five-order approximate solution for λ = 1, y() = 1 λ = 1, y() = 1 From numerical results in Tab 1 2, we can find that the accuracy of approximate solution by Modified Taylor series is more accurate than the solution by Taylor eries method, particularly at end points Also, Tab 1 Tab 2 show that proposed method can be used for broad intervals (67) Table 1 Absolute error of Example 1 (λ = 1, y() = ) x i Taylor series [6] n = 3 Present method n = 3 Taylor series [6] n = 5 Present method n = Table 2 Absolute error of Example 1 (λ = 1, y() = 1) x i Taylor series [6] n = 3 Present method n = 3 Taylor series [6] n = 5 Present method n = Example 2 Consider the second order linear integro-differential equation [6, 22] y + xy xy = e x 2 sin x + with initial conditions y() = 1, y (), which is the exact solution y(x) = e x 1 1 e t y(t) sin xdt, 1 x, t 1, (68) WJMS for subscription: info@wjmsorguk
12 3 J Rashidinia & A Tahmasebi: Approximate solution of linear integro-differential equations Our method is applied in this Example with h = 1 with order of approximate n = 4 The absolute error in the grid points are given also compared the absolute error in the solution with [6] are tabulated in Tab 3 From Tab 3, we fined that numerical results based on the present method are closer to the exact solution than those in [6] Table 3 Absolute error of Example 2 x i Taylor series [6] n = 4 Present method n = Example 3 Consider the following high order linear Fredholm integro-differential equation 1 y (5) x 2 y (3) y xy = x 2 cos x x sin x + y(t)dt, 1 x, t 1, (69) 1 with initial condition y() = y () = y (4) =, y () = 1, y (3) = 1 exact solution y(x) = sin x We approximate the solution by n = 5, h = 2 of our method Tab 4 shows that the numerical results of our method are better than results obtained in [6] Table 4 Absolute error of Example 3 x i Taylor method [6] n = 5 Present method n = Example 4 Finally, consider the following third order linear Voltera integro-differential equation y (3) xy (2) = 4 7 x9 8 5 x7 x 6 + 6x x x 2 t 3 y(t)dt, x, t 1, (7) with the initial conditions y() = 1, y () = 2, y (2) () = We find the approximate solution by the proposed method for n = 3 is the same as y(x) = x 3 + 2x + 1, the exact solution WJMS for contribution: submit@wjmsorguk
13 World Journal of Modelling Simulation, Vol 9 (213) No 4, pp Conclusion The approximation by Taylor series method in [7 1, 22] has the appropriate accuracy at the closed neighborhood of the point of expansion, but for the points far from points the accuracy of the approximation reduced considerably, to overcome such difficultly, we modified develop Taylor expansion approximation to obtain Taylor expansion approximate at grade points x = x i on the interval [a, b] The comparison of the results obtained by the presented method conventional Taylor series method reveals that ours method is very effective convenient The efficiency of the algorithms are verified by some test problems References [1] W Volk The iterated Galerkin methods for linear integro-differential equations Journal of Computational Applied Mathematics, 1988, 21: [2] H Brunner, A Makroglou, K Miller Mixed interpolation collocation methods for first second order Volterra integro-differential equations with periodic solution Applied numerical mathematics, 1997, 23: [3] Q Hu Interpolation correction for collocation solutions of Fredholm integro-differential equation Mathematics of Computation of the American Mathematical Society, 1998, 223: [4] M Rashed Numerical solution of a special type of integro-differential equations Applied mathematics computation, 23, 143: [5] R Kanwall, K Liu A Taylor expansion approach for solving integral equations International Journal of Mathematical Education in Science Technology, 1989, 2: [6] Y Huang, X Li Approximate solution of a class of linear integro-differential equations by Taylor expansion method International Journal of Computer Mathematics, 21, 87(6): [7] M Sezer A method for the approximate solution of the second-order linear differential equations in terms of Taylor polynomials International Journal of Mathematical Education in Science Technology, 1996, 27(6) [8] M Sezer Taylor polynomial solution of Volterra integral equations International Journal of Mathematical Education in Science Technology, 1994, 25: [9] S Yalcinbas, M Sezer The approximate solution of high-order linear Volterra-Fredholm integro-differential equations in terms of Taylor polynomials Applied Mathematics Computation, 2, 112: [1] S Yalcinbas Taylor polynomial solutions of nonlinear Volterra-Fredholm integral equations Applied Mathematics Computation, 22, 127: [11] S Yuzbasi, N Sahin, M Sezer Bessel polynomial solutions of high-order linear Volterra integro-differential equations Computers & Mathematics with Applications, 211, 62: [12] S Yuzbas, N Sahin, M Sezer Numerical solutions of systems of linear Fredholm integro-differential equations with Bessel polynomial bases Computers & Mathematics with Applications, 211, 61: [13] A Akyuz-Dascioglu, M Sezer Chebyshev polynomial solutions of systems of higher-order linear Fredholm- Volterra integro-differential equations Journal of the Franklin Institute, 25, 342: [14] A Akyuz-Dascioglu A chebyshev polynomial approach for linear Fredholm-Volterra integro-differential equations in the most general form Applied mathematics computation, 26, 181: [15] N Sahin, S Yuzbasi, M Sezer A Bessel Polynomial approach for solving general linear Fredholm integrodifferential-difference equations International Journal of Computer Mathematics, 211, 88: [16] Y Gamze, Y Suayip, S Mehmet A Chebyshev Method for a class of high-order linear Fredholm integrodifferential equations Advanced Research in Applied Math, 212, 4: [17] N Çahın, Ş Yüzbaşı, M Gülsu A collocation approach for solving systems of linear Volterra integral equations with variable coefficients Computers & Mathematics with Applications, 211, 62: [18] S Yuzbasi, M Sezer, B Kemanci Numerical solutions of integro-differential equations application of a population model with an improved Legendre method Applied Mathematical Modelling, 213, 37(4): [19] A Saadatmi, M Dehghan Numerical solution of the higher-order linear Fredholm integro-differentialdifference equation with variable coefficients Computers & Mathematics with Applications, 21, 59: [2] L Huang, X Li, Y Huang Approximate solution of Abel integral equation Computers & Mathematics with Applications, 28, 56: [21] X Li, L Huang, Y Huang A new Abel inversion by means of the integrals of an input function with noise Journal of Physics A: Mathematical Theoretical, 27, 4: [22] S Yalcinbas, M Sezer The approximate solution of high-order linear Volterra-Fredholm integro-differential equations in terms of Taylor polynomials Applied Mathematics Computation, 2, 112: [23] M Berenguer, M Fernández Muˆnoz, et al A sequential approach for solving the Fredholm integro-differential equation Applied Numerical Mathematics, 212, 62(4): WJMS for subscription: info@wjmsorguk
RESEARCH ARTICLE. Spectral Method for Solving The General Form Linear Fredholm Volterra Integro Differential Equations Based on Chebyshev Polynomials
Journal of Modern Methods in Numerical Mathematics ISSN: 29-477 online Vol, No, 2, c 2 Modern Science Publishers wwwm-sciencescom RESEARCH ARTICLE Spectral Method for Solving The General Form Linear Fredholm
More informationApproximation of Systems of Volterra Integro-Differential Equations Using the New Iterative Method
Approximation of Systems of Volterra Integro-Differential Equations Using the New Iterative Method Hassan IBRAHIM 1, Peter Vanenchii AYOO 2 1, 2 Department of Mathematics, Federal University Lafia, PMB
More informationA Taylor polynomial approach for solving differential-difference equations
Journal of Computational and Applied Mathematics 86 (006) 349 364 wwwelseviercom/locate/cam A Taylor polynomial approach for solving differential-difference equations Mustafa Gülsu, Mehmet Sezer Department
More informationNumerical Solution of Linear Fredholm Integro-Differential Equations by Non-standard Finite Difference Method
Available at http://pvamu.edu/aam Appl. Appl. Math. ISS: 1932-9466 Vol. 10, Issue 2 (December 2015), pp. 1019-1026 Applications and Applied Mathematics: An International Journal (AAM) umerical Solution
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 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 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 informationResearch Article A Matrix Method Based on the Fibonacci Polynomials to the Generalized Pantograph Equations with Functional Arguments
Advances in Mathematical Physics, Article ID 694580, 5 pages http://dx.doi.org/10.1155/2014/694580 Research Article A Matrix Method Based on the Fibonacci Polynomials to the Generalized Pantograph Equations
More informationApplication of Taylor-Padé technique for obtaining approximate solution for system of linear Fredholm integro-differential equations
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 1 (June 2017), pp. 392 404 Applications and Applied Mathematics: An International Journal (AAM) Application of Taylor-Padé
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 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 informationTaylor polynomial approach for systems of linear differential equations in normal form and residual error estimation
NTMSCI 3, No 3, 116-128 (2015) 116 New Trends in Mathematical Sciences http://wwwntmscicom Taylor polynomial approach for systems of linear differential equations in normal form and residual error estimation
More informationThe Müntz-Legendre Tau method for Weakly singular Volterra integral equations
ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 14 (218) No. 3, pp. 199-28 The Müntz-Legendre Tau method for Weakly singular Volterra integral equations Ali Tahmasbi 1 *, Behruz
More informationNumerical solution of system of linear integral equations via improvement of block-pulse functions
Journal of Mathematical Modeling Vol. 4, No., 16, pp. 133-159 JMM Numerical solution of system of linear integral equations via improvement of block-pulse functions Farshid Mirzaee Faculty of Mathematical
More informationComparison of He s variational iteration and Taylor Expansion Methods for the Solutions of Fractional Integro-Differential Equations
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 11, Issue 4 Ver. III (Jul - Aug. 2015), PP 20-26 www.iosrjournals.org Comparison of He s variational iteration and Taylor
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 informationModified Variational Iteration Method for the Multi-pantograph Equation with Convergence Analysis
Australian Journal of Basic and Applied Sciences, 5(5): 886-893, 0 ISSN 99-878 Modified Variational Iteration Method for the Multi-pantograph Equation with Convergence Analysis Mohsen Alipour, Kobra Karimi,
More informationCOLLOCATION METHOD FOR VOLTERRA-FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS
ISSN: 0976-2876 (Print) ISSN: 2250-0138(Online) COLLOCATION METHOD FOR VOLTERRA-FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS NEHZAT EBRAHIMI a1 AND JALIL RASHIDINIA b ab Department of Mathematics, Islamic Azad
More informationHybrid Functions Approach for the Fractional Riccati Differential Equation
Filomat 30:9 (2016), 2453 2463 DOI 10.2298/FIL1609453M Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Hybrid Functions Approach
More informationChebyshev polynomials for solving two dimensional linear and nonlinear integral equations of the second kind
Volume 31, N. 1, pp. 127 142, 2012 Copyright 2012 SBMAC ISSN 0101-8205 / ISSN 1807-0302 (Online) www.scielo.br/cam Chebyshev polynomials for solving two dimensional linear and nonlinear integral equations
More informationAnalytical Solution of BVPs for Fourth-order Integro-differential Equations by Using Homotopy Analysis Method
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.9(21) No.4,pp.414-421 Analytical Solution of BVPs for Fourth-order Integro-differential Equations by Using Homotopy
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 informationThe Approximate Solution of Non Linear Fredholm Weakly Singular Integro-Differential equations by Using Chebyshev polynomials of the First kind
AUSTRALIA JOURAL OF BASIC AD APPLIED SCIECES ISS:1991-8178 EISS: 2309-8414 Journal home page: wwwajbaswebcom The Approximate Solution of on Linear Fredholm Weakly Singular Integro-Differential equations
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 informationLaguerre polynomial solution of high- order linear Fredholm integro-differential equations
NTMSCI 4, No 2, 273-284 (216) 273 New Trends in Mathematical Sciences http://dxdoiorg/12852/ntmsci216218534 Laguerre polynomial solution of high- order linear Fredholm integro-differential equations Nurcan
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 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 informationChebyshev finite difference method for a Two Point Boundary Value Problems with Applications to Chemical Reactor Theory
Iranian Journal of Mathematical Chemistry, Vol. 3, o., February 22, pp. 7 IJMC Chebyshev finite difference method for a Two Point Boundary Value Problems with Applications to Chemical Reactor Theory ABBAS
More informationSOLUTION OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD
SOLUTION OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD R. C. Mittal 1 and Ruchi Nigam 2 1 Department of Mathematics, I.I.T. Roorkee, Roorkee, India-247667. Email: rcmmmfma@iitr.ernet.in
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 informationAn Approximate Solution for Volterra Integral Equations of the Second Kind in Space with Weight Function
International Journal of Mathematical Analysis Vol. 11 17 no. 18 849-861 HIKARI Ltd www.m-hikari.com https://doi.org/1.1988/ijma.17.771 An Approximate Solution for Volterra Integral Equations of the Second
More informationApplication of the Bernstein Polynomials. for Solving the Nonlinear Fredholm. Integro-Differential Equations
Journal of Applied Mathematics & Bioinformatics, vol., no.,, 3-3 ISSN: 79-66 (print), 79-6939 (online) International Scientific Press, Application of the Bernstein Polynomials for Solving the Nonlinear
More informationNumerical solution of delay differential equations via operational matrices of hybrid of block-pulse functions and Bernstein polynomials
Computational Methods for Differential Equations http://cmdetabrizuacir Vol, No, 3, pp 78-95 Numerical solution of delay differential equations via operational matrices of hybrid of bloc-pulse functions
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 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 informationNon-Polynomial Spline Solution of Fourth-Order Obstacle Boundary-Value Problems
Non-Polynomial Spline Solution of Fourth-Order Obstacle Boundary-Value Problems Jalil Rashidinia Reza Jalilian Abstract In this paper we use quintic non-polynomial spline functions to develop numerical
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 informationSMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract
SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS Kai Diethelm Abstract Dedicated to Prof. Michele Caputo on the occasion of his 8th birthday We consider ordinary fractional
More informationLinear Fredholm Integro-Differential Equation. of the Second Kind. By Khulood Nedal Iseed Thaher. Supervised Prof. Naji Qatanani
An-Najah National University Faculty of Graduate Studies Linear Fredholm Integro-Differential Equation of the Second Kind By Khulood Nedal Iseed Thaher Supervised Prof. Naji Qatanani This Thesis is Submitted
More informationMethod of Successive Approximations for Solving the Multi-Pantograph Delay Equations
Gen. Math. Notes, Vol. 8, No., January, pp.3-8 ISSN 9-784; Copyright c ICSRS Publication, www.i-csrs.org Available free online at http://www.geman.in Method of Successive Approximations for Solving the
More informationNumerical analysis of the one-dimensional Wave equation subject to a boundary integral specification
Numerical analysis of the one-dimensional Wave equation subject to a boundary integral specification B. Soltanalizadeh a, Reza Abazari b, S. Abbasbandy c, A. Alsaedi d a Department of Mathematics, University
More 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 informationApplication of homotopy perturbation method to non-homogeneous parabolic partial and non linear differential equations
ISSN 1 746-733, England, UK World Journal of Modelling and Simulation Vol. 5 (009) No. 3, pp. 5-31 Application of homotopy perturbation method to non-homogeneous parabolic partial and non linear differential
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More informationNumerical Solution of Fredholm and Volterra Integral Equations of the First Kind Using Wavelets bases
The Journal of Mathematics and Computer Science Available online at http://www.tjmcs.com The Journal of Mathematics and Computer Science Vol.5 No.4 (22) 337-345 Numerical Solution of Fredholm and Volterra
More informationNumerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method
Numerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method Y. Azari Keywords: Local RBF-based finite difference (LRBF-FD), Global RBF collocation, sine-gordon
More informationApplication of the Bernstein Polynomials for Solving Volterra Integral Equations with Convolution Kernels
Filomat 3:4 (216), 145 152 DOI 1.2298/FIL16445A Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Application of the Bernstein Polynomials
More informationACTA UNIVERSITATIS APULENSIS No 20/2009 AN EFFECTIVE METHOD FOR SOLVING FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS. Wen-Hua Wang
ACTA UNIVERSITATIS APULENSIS No 2/29 AN EFFECTIVE METHOD FOR SOLVING FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS Wen-Hua Wang Abstract. In this paper, a modification of variational iteration method is applied
More informationSolving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions.
Journal of Mathematical Modeling Vol 1, No 1, 213, pp 28-4 JMM Solving a class of nonlinear two-dimensional Volterra integral equations by using two-dimensional triangular orthogonal functions Farshid
More informationON PROJECTIVE METHODS OF APPROXIMATE SOLUTION OF SINGULAR INTEGRAL EQUATIONS. Introduction Let us consider an operator equation of second kind [1]
GEORGIAN MATHEMATICAL JOURNAL: Vol. 3, No. 5, 1996, 457-474 ON PROJECTIVE METHODS OF APPROXIMATE SOLUTION OF SINGULAR INTEGRAL EQUATIONS A. JISHKARIANI AND G. KHVEDELIDZE Abstract. The estimate for the
More informationSolving Fredholm integral equations with application of the four Chebyshev polynomials
Journal of Approximation Theory and Applied Mathematics, 204 Vol. 4 Solving Fredholm integral equations with application of the four Chebyshev polynomials Mostefa NADIR Department of Mathematics University
More informationResearch Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method
Applied Mathematics Volume 2012, Article ID 605741, 10 pages doi:10.1155/2012/605741 Research Article The Numerical Solution of Problems in Calculus of Variation Using B-Spline Collocation Method M. Zarebnia
More informationYURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL
Journal of Comput. & Applied Mathematics 139(2001), 197 213 DIRECT APPROACH TO CALCULUS OF VARIATIONS VIA NEWTON-RAPHSON METHOD YURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL Abstract. Consider m functions
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 informationODE Homework Series Solutions Near an Ordinary Point, Part I 1. Seek power series solution of the equation. n(n 1)a n x n 2 = n=0
ODE Homework 6 5.2. Series Solutions Near an Ordinary Point, Part I 1. Seek power series solution of the equation y + k 2 x 2 y = 0, k a constant about the the point x 0 = 0. Find the recurrence relation;
More informationCh 5.4: Regular Singular Points
Ch 5.4: Regular Singular Points! Recall that the point x 0 is an ordinary point of the equation if p(x) = Q(x)/P(x) and q(x)= R(x)/P(x) are analytic at at x 0. Otherwise x 0 is a singular point.! Thus,
More informationAdaptation of Taylor s Formula for Solving System of Differential Equations
Nonlinear Analysis and Differential Equations, Vol. 4, 2016, no. 2, 95-107 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/nade.2016.51144 Adaptation of Taylor s Formula for Solving System of Differential
More informationDesign of digital differentiators and integrators of order 1 2
ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 4 (2008) No. 3, pp. 182-187 Design of digital differentiators and integrators of order 1 2 B. T. Krishna, K. V. V. S. Reddy Department
More informationA New Numerical Scheme for Solving Systems of Integro-Differential Equations
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 1, No. 2, 213, pp. 18-119 A New Numerical Scheme for Solving Systems of Integro-Differential Equations Esmail Hesameddini
More informationResearch Article Quasilinearization Technique for Φ-Laplacian Type Equations
International Mathematics and Mathematical Sciences Volume 0, Article ID 975760, pages doi:0.55/0/975760 Research Article Quasilinearization Technique for Φ-Laplacian Type Equations Inara Yermachenko and
More informationUse of Bernstein Polynomials in Numerical Solutions of Volterra Integral Equations
Applied Mathematical Sciences, Vol. 2, 28, no. 36, 1773-1787 Use of Bernstein Polynomials in Numerical Solutions of Volterra Integral Equations Subhra Bhattacharya Department of Mathematics Jadavpur University,
More informationFractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials
Applied Mathematical Sciences, Vol. 5, 211, no. 45, 227-2216 Fractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials Z. Avazzadeh, B. Shafiee and G. B. Loghmani Department
More informationNew Iterative Method for Time-Fractional Schrödinger Equations
ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 9 2013) No. 2, pp. 89-95 New Iterative Method for Time-Fractional Schrödinger Equations Ambreen Bibi 1, Abid Kamran 2, Umer Hayat
More informationSolving Integral Equations by Petrov-Galerkin Method and Using Hermite-type 3 1 Elements
5 1 July 24, Antalya, Turkey Dynamical Systems and Applications, Proceedings, pp. 436 442 Solving Integral Equations by Petrov-Galerkin Method and Using Hermite-type 3 1 Elements M. Karami Department of
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 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 informationSection 5.2 Series Solution Near Ordinary Point
DE Section 5.2 Series Solution Near Ordinary Point Page 1 of 5 Section 5.2 Series Solution Near Ordinary Point We are interested in second order homogeneous linear differential equations with variable
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationFree convection modeling over a vertical flat plate embedded in saturated porous medium with a variable heat source and radiation flux
ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 9 (2013) No. 3, pp. 163-172 Free convection modeling over a vertical flat plate embedded in saturated porous medium with a variable
More informationPolynomial Solutions of the Laguerre Equation and Other Differential Equations Near a Singular
Polynomial Solutions of the Laguerre Equation and Other Differential Equations Near a Singular Point Abstract Lawrence E. Levine Ray Maleh Department of Mathematical Sciences Stevens Institute of Technology
More informationLCPI method to find optimal solutions of nonlinear programming problems
ISSN 1 746-7233, Engl, UK World Journal of Modelling Simulation Vol. 14 2018) No. 1, pp. 50-57 LCPI method to find optimal solutions of nonlinear programming problems M.H. Noori Skari, M. Ghaznavi * Faculty
More informationNUMERICAL METHOD FOR THE MIXED VOLTERRA-FREDHOLM INTEGRAL EQUATIONS USING HYBRID LEGENDRE FUNCTIONS
Conference Applications of Mathematics 215, in honor of the birthday anniversaries of Ivo Babuška (9), Milan Práger (85), and Emil Vitásek (85) J. Brandts, S. Korotov, M. Křížek, K. Segeth, J. Šístek,
More informationDifferential transformation method for solving one-space-dimensional telegraph equation
Volume 3, N 3, pp 639 653, 2 Copyright 2 SBMAC ISSN -825 wwwscielobr/cam Differential transformation method for solving one-space-dimensional telegraph equation B SOLTANALIZADEH Young Researchers Club,
More 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 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 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 Second Order Singularly Perturbed Differential Difference Equation with Negative Shift. 1 Introduction
ISSN 749-3889 print), 749-3897 online) International Journal of Nonlinear Science Vol.8204) No.3,pp.200-209 Numerical Solution of Second Order Singularly Perturbed Differential Difference Equation with
More informationMA22S3 Summary Sheet: Ordinary Differential Equations
MA22S3 Summary Sheet: Ordinary Differential Equations December 14, 2017 Kreyszig s textbook is a suitable guide for this part of the module. Contents 1 Terminology 1 2 First order separable 2 2.1 Separable
More informationNumerical Solutions of Volterra Integral Equations of Second kind with the help of Chebyshev Polynomials
Annals of Pure and Applied Mathematics Vol. 1, No. 2, 2012, 158-167 ISSN: 2279-087X (P), 2279-0888(online) Published on 16 November 2012 www.researchmathsci.org Annals of Numerical Solutions of Volterra
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 informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More 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 informationReview of Power Series
Review of Power Series MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Introduction In addition to the techniques we have studied so far, we may use power
More informationThe Trigonometric Cubic B-spline Algorithm for Burgers Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.4(017) No., pp.10-18 The Trigonometric Cubic B-spline Algorithm for Burgers Equation Idris Dag 1, Ozlem Ersoy Hepson,
More informationSolution of Quadratic Integral Equations by the Adomian Decomposition Method
Copyright 213 Tech Science Press CMES, vol.92, no.4, pp.369-385, 213 Solution of Quadratic Integral Equations by the Adomian Decomposition Method Shou-Zhong Fu 1, Zhong Wang 1 and Jun-Sheng Duan 1,2,3
More information(Received 10 December 2011, accepted 15 February 2012) x x 2 B(x) u, (1.2) A(x) +
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol13(212) No4,pp387-395 Numerical Solution of Fokker-Planck Equation Using the Flatlet Oblique Multiwavelets Mir Vahid
More informationNumerical Solution of One-dimensional Telegraph Equation using Cubic B-spline Collocation Method
2014 (2014) 1-8 Available online at www.ispacs.com/iasc Volume 2014, Year 2014 Article ID iasc-00042, 8 Pages doi:10.5899/2014/iasc-00042 Research Article Numerical Solution of One-dimensional Telegraph
More informationNumerical Solution of Volterra-Fredholm Integro-Differential Equation by Block Pulse Functions and Operational Matrices
Gen. Math. Notes, Vol. 4, No. 2, June 211, pp. 37-48 ISSN 2219-7184; Copyright c ICSRS Publication, 211 www.i-csrs.org Available free online at http://www.gean.in Nuerical Solution of Volterra-Fredhol
More informationFibonacci tan-sec method for construction solitary wave solution to differential-difference equations
ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 7 (2011) No. 1, pp. 52-57 Fibonacci tan-sec method for construction solitary wave solution to differential-difference equations
More informationCubic B-spline Collocation Method for Fourth Order Boundary Value Problems. 1 Introduction
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.142012 No.3,pp.336-344 Cubic B-spline Collocation Method for Fourth Order Boundary Value Problems K.N.S. Kasi Viswanadham,
More informationAcceleration of Levenberg-Marquardt method training of chaotic systems fuzzy modeling
ISSN 746-7233, England, UK World Journal of Modelling and Simulation Vol. 3 (2007) No. 4, pp. 289-298 Acceleration of Levenberg-Marquardt method training of chaotic systems fuzzy modeling Yuhui Wang, Qingxian
More informationAn Exponential Matrix Method for Numerical Solutions of Hantavirus Infection Model
Available at http://pvamu.edu/aam Appl. Appl. Math. ISS: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 99-115 Applications and Applied Mathematics: An International Journal (AAM) An Exponential Matrix Method
More informationJim Lambers MAT 772 Fall Semester Lecture 21 Notes
Jim Lambers MAT 772 Fall Semester 21-11 Lecture 21 Notes These notes correspond to Sections 12.6, 12.7 and 12.8 in the text. Multistep Methods All of the numerical methods that we have developed for solving
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 informationAdvanced Computational Fluid Dynamics AA215A Lecture 2 Approximation Theory. Antony Jameson
Advanced Computational Fluid Dynamics AA5A Lecture Approximation Theory Antony Jameson Winter Quarter, 6, Stanford, CA Last revised on January 7, 6 Contents Approximation Theory. Least Squares Approximation
More informationDirect method for variational problems by using hybrid of block-pulse and Bernoulli polynomials
Direct method for variational problems by using hybrid of block-pulse and Bernoulli polynomials M Razzaghi, Y Ordokhani, N Haddadi Abstract In this paper, a numerical method for solving variational problems
More informationOptimal Control of Nonlinear Systems Using the Shifted Legendre Polynomials
Optial Control of Nonlinear Systes Using Shifted Legendre Polynoi Rahan Hajohaadi 1, Mohaad Ali Vali 2, Mahoud Saavat 3 1- Departent of Electrical Engineering, Shahid Bahonar University of Keran, Keran,
More informationSemiorthogonal Quadratic B-Spline Wavelet Approximation for Integral Equations
Mathematical Sciences Vol. 3, No. (29) 99- Semiorthogonal Quadratic B-Spline Wavelet Approximation for Integral Equations Mohsen Rabbani a,, Nasser Aghazadeh b a Department of Mathematics, Islamic Azad
More informationEngg. Math. II (Unit-IV) Numerical Analysis
Dr. Satish Shukla of 33 Engg. Math. II (Unit-IV) Numerical Analysis Syllabus. Interpolation and Curve Fitting: Introduction to Interpolation; Calculus of Finite Differences; Finite Difference and Divided
More informationDifferentiable Functions
Differentiable Functions Let S R n be open and let f : R n R. We recall that, for x o = (x o 1, x o,, x o n S the partial derivative of f at the point x o with respect to the component x j is defined as
More information