Numerical solution of system of linear integral equations via improvement of block-pulse functions
|
|
- Garry Price
- 6 years ago
- Views:
Transcription
1 Journal of Mathematical Modeling Vol. 4, No., 16, pp JMM Numerical solution of system of linear integral equations via improvement of block-pulse functions Farshid Mirzaee Faculty of Mathematical Sciences and Statistics, Malayer University, P.O. Box , Malayer, Iran Abstract. In this article, a numerical method based on improvement of block-pulse functions (IBPFs) is discussed for solving the system of linear Volterra and Fredholm integral equations. By using IBPFs and their operational matrix of integration, such systems can be reduced to a linear system of algebraic equations. An efficient error estimation and associated theorems for the proposed method are also presented. Some examples are given to clarify the efficiency and accuracy of the method. Keywords: system of linear integral equations, improvement of block-pulse functions, operational matrix, vector forms, error analysis. AMS Subject Classification: 45A5, 45D5, 45F5, 65R, 65D3. 1 Introduction Systems of linear integral equations and their solutions have great importance in science and engineering. Most physical problems, such as biological applications in population dynamics and genetics where impulses arise naturally or are caused by control, can be modeled by an integral equation or a system of these equations. The systems of integral equations are usually difficult to solve analytically, thus some numerical methods are applied to approximately solve them. Systems of linear Volterra and Fredholm inte- Corresponding author. Received: 18 April 16 / Revised: 1 August 16 / Accepted: September 16. c 16 University of Guilan
2 134 F. Mirzaee gral equations can be appeared in the following form, respectively and f(x) = g(x) + x f(x) = g(x) + k(x, y)f(y)dy, x D = [, 1), (1) k(x, y)f(y)dy, x D, () where f(x) is an unknown function and g(x) and k(x, y) are analytical function on D and D D, respectively, as follows f(x) = [f 1 (x), f (x),..., f m (x)] T, (3) g(x) = [g 1 (x), g (x),..., g m (x)] T, (4) k(x, y) = [k ij (x, y)], i, j = 1,,..., m. (5) Recently, several numerical methods such as the homotopy perturbation method [4, 5], the Lagrange method [3], the modified homotopy perturbation method [8], the rationalized Haar functions method [1, 11], the differential transformation method [1], the Tau method [16], the variational iteration method [19], the Legendre matrix method [], the Adomian method [5], the Galerkin method [1], the Bessel matrix method [1], and other methods [14 19] have been used for solving integral and integrodifferential equations systems. In this paper, IBPFs are introduced. Also, some theorems are proved for IBPFs method that show the results of this numerical expansions are more precise than the results of block pulse expansions. This functions are disjoint, orthogonal and complete. According to the disjointness of IBPFs, the joint terms will disappear in each subinterval when multiplication, division and some other operations are applied. Also, the orthogonality property of IBPFs will cause that the operational matrix to be a sparse matrix. The completeness of IBPFs guarantees that an arbitrary small mean square error can be obtained for a real bounded function, which has only a finite number of discontinuous point in the interval x [, 1), by increasing the number of terms in the improved block pulse series. The rest of paper is organized as follows: In Section, we describe IBPFs and their properties. In Section 3, we apply IBPFs for approximating the solution of system of linear Volterra and Fredholm integral equations. Convergence analysis is discussed in Section 4. Numerical results are given in Section 5 to illustrate the efficiency and accuracy of proposed method. Finally, Section 6 concludes the paper.
3 Numerical solution of system of linear integral equations 135 IBPFs and their properties.1 Definition of IBPFs An (n + 1)-set of IBPFs consists of (n + 1) functions which are defined over district D as follows { 1, x [, h φ (x) = ),, otherwise, φ i (x) = { 1, x [(i 1)h + h, ih + h ),, otherwise, i = 1,,..., n 1, φ n (x) = { 1, x [1 h, 1),, otherwise, where n is an arbitrary positive integer and h = 1 n. The IBPFs are disjoint φ i (x)φ j (x) = { φi (x), i = j,, otherwise, where i, j =, 1,..., n and are orthogonal to each other h 1, i = j {, n}, φ i (x)φ j (x)dx = h, i = j {1,,..., n 1},, otherwise, (6) (7) where x D.. Vector forms Consider the first (n + 1) terms of IBPFs and write them concisely as (n + 1)-vector Φ n (x) = [φ (x), φ 1 (x),..., φ n (x)] T, x D. (8) Eq. (6) implies that Φ n (x)φ T n (x) = φ (x)... φ 1 (x) = diag(φ n(x)). (9)... φ n (x)
4 136 F. Mirzaee Now suppose that X be an (n + 1)-vector. Hence by using Eq. (9) we obtain Φ n (x)φ T n (x)x = XΦ n (x), where X = diag(x) is an (n + 1) (n + 1) diagonal matrix..3 Operational matrix We have x { x, x [, h φ (y)dy = ), h, otherwise. Note that x = h 4, at mid-point of [, h ). x [, h ), by h 4. x φ i (x) = So we can approximate x, for, x [, (i 1)h + h ), x (i 1)h h, x [(i 1)h + h, ih + h ), h, otherwise, for i = 1,,..., n 1. Also, we can approximate x (i 1)h h, for x [(i 1)h + h, ih + h ), by h. x { x 1 + h φ n (x) =, x [1 h, 1), otherwise. So, we can approximate x 1 + h, for x [1 h, 1), by h 4. Therefore, the integration of the vector Φ n (x) defined in Eq. (8) can be approximated by where Using Eq. (9), we get x Φ n (y)dy P 1Φ n (x), (1) P 1 = h Φ n (y)φ T n (y)dy P, (11)
5 Numerical solution of system of linear integral equations 137 where P = h IBPFs expansions A continues function f(x) L (D) may be expanded by the IBPFs as n f(x) f n (x) = f i φ i (x) = Fn T Φ n (x) = Φ T n (x)f n, (1) i= where F n is an (n + 1) 1 vector given by F n = [f, f 1,..., f n ] T, and Φ n (x) is defined in Eq. (8), and f i is obtained as n h f(x)dx, i =, f i = n ih+ h f(x)dx, i = 1,,..., n 1, (i 1)h+ h n 1 h f(x)dx, i = n. (13) Similarly a function of two variables, k(x, y) L (D D) can be approximated by IBPFs as follows k(x, y) k n (x, y) = Φ T n (x)k n Φ n (y), (14) where Φ n (x) and Φ n (y) are IBPFs vector of dimension (n + 1), and K n = [k ij ] is the (n + 1) (n + 1) IBPFs coefficients matrix of k(x, y). 3 Method of solution 3.1 System of linear Volterra integral equations In this section, we solve system of linear Volterra integral equations of the form Eq. (1) by using IBPFs. We can rewrite Eq. (1) as follows f i (x) = g i (x) + m j=1 x k ij (x, y)f j (y)dy, i = 1,,..., m. (15)
6 138 F. Mirzaee We now approximate functions f i (x), g i (x) and k ij (x, y), i, j = 1,,..., m, by IBPFs as follows f i (x) Φ T n (x)f i, g i (x) Φ T n (x)g i, (16) k ij (x, y) Φ T n (x)k ij Φ n (y), where i, j = 1,,..., m, Φ n (x) is defined by (8), F i, G i and K ij are IBPFs coefficients of f i (x), g i (x) and k ij (x, y), respectively. By substituting Eq. (16) in Eq. (15), we get Φ T n (x)f i = Φ T n (x)g i + m j=1 x therefore by Eq. (9), we have Φ T n (x)f i = Φ T n (x)g i + From Eq.(1), we have where Φ T n (x)k ij Φ n (y)φ T n (y)f j dy, i = 1,,..., m, m x Φ T n (x)k ij diag(φ n (y))dyf j, i = 1,,..., m. j=1 Φ T n (x)f i = Φ T n (x)g i + m Φ T n (x)k ij AF j, i = 1,,..., m, (17) j=1 P 1 Φ n (x)... P 1 1 Φ n (x)... A =......,... P 1 n Φ n (x) and P 1 i is ith row of P 1. On the other hand, using (6), we have Kij P 1 Φ n (x) Kij 1P 1 1Φ n (x)... Kij np 1 nφ n (x) Φ T n (x)k ij A = Φ T Kij 1 n (x) P 1 Φ n (x) Kij 11P 1 1Φ n (x)... Kij 1nP 1 nφ n (x) Kij np 1 Φ n (x) Kij n1p 1 1Φ n (x)... Kij nnp 1 nφ n (x) [ n ] n n = Kij l P 1 l φ l (x) Kij l1 P 1 1l φ l (x)... Kij ln P 1 nl φ l (x) l= = Φ T n (x)b ij, l= l=
7 Numerical solution of system of linear integral equations 139 where B ij = K ij P 1 K 1 ij P K n ij P 1 n K 1 ij P 1 1 K 11 ij P K 1n ij P 1 n1. K n ij P 1 n..... K n1 ij P 1 1n... K nn ij P 1 nn. So, we can rewrite Eq. (17) as follows or Let m Φ T n (x)f i = Φ T n (x)g i + Φ T n (x)b ij F j, i = 1,,..., m, j=1 m F i = G i + B ij F j, i = 1,,..., m. (18) j=1 F = [F 1, F,..., F m ] T, G = [G 1, G,..., G m ] T, (19) B = [B ij ], i, j = 1,,..., m. Therefore, Eq. (18) can be written as F = G + BF. After solving above linear system, we can find F and accordingly find F i, i = 1,,..., m, so f i (x) Φ T n (x)f i. 3. System of linear Fredholm integral equations In this section, we solve system of linear Fredholm integral equations of the form Eq. () by using IBPFs. We can rewrite Eq. () as follows f i (x) = g i (x) + m j=1 Substituting Eq. (16) in Eq. (), we get Φ T n (x)f i = Φ T n (x)g i + m j=1 k ij (x, y)f j (y)dy, i = 1,,..., m. () Φ T n (x)k ij Φ n (y)φ T n (y)f j dy, i = 1,,..., m.
8 14 F. Mirzaee Using Eq. (11), we get Φ T n (x)f i = Φ T n (x)g i + m Φ T n (x)k ij P F j, i = 1,,..., m. j=1 Letting C = [K ij P ], from Eq.(19) we have F = G + CF. After solving above linear system, we can find F and accordingly find F i ; i = 1,,..., m, so f i (x) Φ T n (x)f i. 4 Convergence analysis In this section, we show that the method discussed in the previous section is convergent and its order of convergence is O( 1 n ). We define and ( f(x) = f(x) dx) 1, f(x) = ( m i=1 where f(x) L (D) and f(x) is defined in Eq. (3) and f i (x) ) 1, (1) and ( k(x, y) = i=1 j=1 k(x, y) dxdy) 1, m m k(x, y) = k ij (x, y) 1, () where k(x, y) L (D D) and k(x, y) is defined in Eq. (5). Also, for our purposes we will need the following theorems. Theorem 1. Let f(x) L (D) and f n (x) be the IBPFs expansion of f(x) that is defined as n f n (x) = f i φ i (x), i=
9 Numerical solution of system of linear integral equations 141 where f i ; i =, 1,..., n, are defined in Eq. (13). Then the criterion of this approximation is that the mean square error between f(x) and f n (x) in the interval x D (f(x) f n (x)) dx, achieves its minimum value and also we have f (x)dx = fi φ i (x). i= Proof. It is an immediate consequence of theorem which is proved in [9]. Theorem. Suppose f(x) is continuous on D, differentiable on (, 1), and there exists a positive scalar M such that f (x) M, for every x D. Then f(b) f(a) M b a, a, b D. Proof. See [18]. Theorem 3. Suppose f n (x) is the IBPFs expansions of f(x) that is defined as Eq. (1) and f(x) is differentiable on D such that f (x) M. Also, assume that e n (x) = f(x) f n (x), then e n (x) = O(h). Proof. Suppose x =, x i = ih h, i = 1,..., n and x n+1 = 1. We define the error between f(x) and its IBPFs expansion over every subinterval I i = [x i, x i+1 ) as follows e n,i (x) = f(x) f i (x), x I i, where i =, 1,..., n. By using mean value theorem for integral, we have e n, (x) = h e n,(x)dx = h where ξ I. Also, for i = 1,,..., n 1, we have e n,i (x) = ih+ h ih h e n,i(x)dx = where ξ i I i. Furthermore, we have e n,n (x) = 1 h e n,n(x)dx = ih+ h ih h 1 h (f(x) f ) dx = h (f(ξ ) f ), (f(x) f i ) dx = h (f(ξ i ) f i ), (f(x) f n ) dx = h (f(ξ n) f n ),
10 14 F. Mirzaee where ξ n I n. Using Eq. (13) and the mean value theorem, we have f i = n h f(x)dx = n h f(η ) = f(η ), i =, n ih+ h f(x)dx = nhf(η (i 1)h+ h i ) = f(η i ), i = 1,,..., n 1, n 1 h f(x)dx = n h f(η n) = f(η n ), i = n, where η i I i, i =, 1,..., n. From the above equations and Theorem, we get e n,i (x) = h (f(ξ ) f(η )) M h ξ η M h 3 8, i =, h (f(ξ i ) f(η i )) M h ξ i η i M h 3, i = 1,,..., n 1, h (f(ξ n) f(η n )) M h ξ n η n M h 3 8, i = n. (3) We have e n (x) = = e n(x)dx = i= ( n e n,i (x)) dx i= ( n ) e n,i(x) dx + i j e n,i (x)e n,j (x)dx. Since for i j, I i I j =, then e n (x) = ( n ) e n,i(x) dx = i= n e n,i (x). (4) i= Substituting Eq. (3) into Eq. (4), we get which completes the proof. e n (x) M h 3M h 3, 4 Suppose that e n(x) is the error between f(x) and its BPFs expansion. From [9], it is clear that e n (x) e n(x).
11 Numerical solution of system of linear integral equations 143 Lemma 1. Let g(x) be as defined in Eq. (4), g n (x) be the IBPFs of g(x) and e g (x) = g(x) g n (x). Then Proof. From Eq. (1), we have e g (x) = e g (x) = O(h). ( m i=1 g i (x) g n,i (x) ) 1, and from Theorem, g i (x) g n,i (x) C i h. Then e g (x) ( m which completes the proof. i=1 C i h ) 1 = ( m i=1 C i ) 1 h = Ch, (5) Theorem 4. Let k n (x, y) be the IBPFs expansions of k(x, y) defined as Eq. (14) and k(x, y) be differentiable on D D such that k (x, y) M. Also, assume that e n (x, y) = k(x, y) k n (x, y), then e n (x, y) = O(h). Proof. Suppose x = y =, x i = y i = ih h, i = 1,..., n and x n+1 = y n+1 = 1. We define the error between k(x, y) and its IBPFs expansion over every subinterval I i,j = [x i, x i+1 ) [y j, y j+1 ) as follows e n,ij (x, y) = k(x, y) k ij (x, y), x I i,j, i, j =, 1,..., n. By using the mean value theorem for integral and similar to the proof of Theorem 3, we get e n,ij (x, y) { M h 4 8, j =, n, 5M h 4 8, j = 1,,..., n 1, (6) for i =, n and e n,ij (x, y) { 5M h 4 8, j =, n, M h 4, j = 1,,..., n 1, (7)
12 144 F. Mirzaee for i = 1,,..., n 1. We have e n (x, y) = = e n(x, y)dxdy = n i= j= + i k j l n e n,ij(x, y) dxdy n i= j= e n,ij (x, y)e n,kl (x, y)dxdy. n e n,ij (x, y) Since for i k and j l, we have I i I k = and I j I l =, then e n (x, y) = = n i= j= i= j= n e n,ij(x, y) dxdy dxdy n n e n,ij (x, y). (8) Substituting Eqs. (6) and (7) into Eq. (8), we get e n (x, y) M h 3M h 3. Suppose e n(x, y) be the error between k(x, y) and its BPFs expansion. From [14], it is clear that e n (x, y) e n(x, y). Lemma. Let k(x, y) be as defined in Eq. (5), k n (x, y) be the IBPFs of k(x, y) and e k (x, y) = k(x, y) k n (x, y). Then Proof. From Eq. (), we have e k (x, y) = O(h). 1 m m e k (x, y) = k ij (x, y) k n,ij (x, y). i=1 j=1
13 Numerical solution of system of linear integral equations 145 From Theorem 3, we conclude that k ij (x, y) k n,ij (x, y) C ij h. Therefore m m e k (x, y) i=1 j=1 C ijh 1 m m = i=1 j=1 C ij 1 h = Ch. (9) Let the error of IBPFs be denoted by E n = f(x) f n (x), x D, where f(x) was defined in Eq. (3). Furthermore, assume the following hypotheses: (M1) Let f(x) N for x D; (M) Let k(x, y) N for (x, y) D D; (M3) According to lemma 1 and, let and E g = e g (x) Ch, E k = e k (x, y) C h, where C and C are coefficients defined in Eqs. (5) and (9) and g(x) and k(x, y) were defined in Eqs. (4) and (5), respectively; (M4) Let N + C h < 1. Theorem 5. Let f(x) and f n (x) be the exact and approximate solutions of Eq. (1) or Eq. (), respectively. Also assumptions (M1)-(M4) are satisfied. Then we have Proof. For the first case, from Eq. (1), we have f(x) f n (x) = g(x) g n (x) + and therefore E n (C + C N)h 1 N C h. (3) x (k(x, y)f(y) k n (x, y)f n (y)) dy, E n E g + x k(x, y)f(y) k n (x, y)f n (y).
14 146 F. Mirzaee It is clear that x 1, So E n E g + k(x, y)f(y) k n (x, y)f n (y). (31) Also for the second case, from Eq. (), we have f(x) f n (x) = g(x) g n (x) + and therefore (k(x, y)f(y) k n (x, y)f n (y)) dy, E n E g + k(x, y)f(y) k n (x, y)f n (y). So Eq. (31) is true in the both cases. Now, according to assumptions (M1)-(M3), we have k(x, y)f(y) k n (x, y)f n (y) k(x, y) E n + E k (E n + f(x) ) N E n + C h (E n + N). (3) Also from assumptions (M3), Eqs. (31) and (3), we have E n (C + C N)h + (N + C h)e n. Therefore according to (M4), Eq. (3) is satisfied and this completes the proof. Also we have E n = O(h). Lemma 3. Suppose f(x) and f n (x) are the exact and approximate solution of Eq. (1) or Eq. (), respectively, where f(x) was defined in Eq. (3) and f n (x) = [f 1,n (x), f,n (x),..., f m,n (x)] T. Then e i,n = f i (x) f i,n (x) = O(h). (33) Proof. From Theorem 4, we have E n Ch and according to Eq. (1) we have e i,n E n Ch.
15 Numerical solution of system of linear integral equations 147 Table 1: Numerical results for Example 1 (n = 8, 16). Nodes x Exact solution Present method n = 8 n = 16. (1.,.,.) (1.68,.948,.3384) (1.17,.4696,.167).1 (1.,.3,.1) (1.33,.37638,.15868) (1.8,.37535,.15686). (1.,.6,.8) (1.391,.75169,.37768) (1.87,.5687,.5843).3 (1.,.9,.48) (1.391,.75169,.37768) (1.15,.93799,.5857).4 (1.,3.,.7) (1.477,3.1733,.65979) (1.117,3.1558,.65713).5 (1.,3.5,1.) (1.61,3.5347,1.56) (1.15,3.586,1.131).6 (1.,3.8,1.3) (1.831,3.8835, ) (1.5,3.8763,1.488).7 (1.,4.1,1.68) (1.1195,4.5841, ) (1.44,4.6416, ).8 (1.,4.4,.8) (1.1195,4.5841, ) (1.359,4.4411,.13678).9 (1.,4.7,.5) (1.1799,4.6384,.463) (1.443,4.683,.411) 5 Numerical experiments In this section, several examples are presented to demonstrate the effectiveness of our approach. In this regard, we have reported the values of the exact solution f(x), the approximate solution f n (x); n = 8; 16; 3; 64 and the L -norms of errors that have been calculated by using ( E,n = f(x) f n (x) dx) 1, x D. at the selected points of the given interval in tables and figures. All computations is preformed in Matlab. In order to show the convergence rate of this numerical method, we introduce the following notation E n,n = log f(x) f n (x) f(x) f n (x). The following problems have been tested. Example 1. Consider the following system of linear Volterra integral equations [] x x f 1 (x) = 1 x x x5 1 3 x6 + yf (y)dy + y 3 f 3 (y)dy, x x f (x) = + 3x 1 x 1 4 x4 5 x5 + yf 1 (y)dy + y f 3 (y)dy, f 3 (x) = x + x x 3 1 x4 3 5 x5 + x x f 1 (y)dy + x y 3 f (y)dy, where x [, 1). The exact solution is f 1 (x) = 1, f (x) = + 3x and f 3 (x) = x+x. Table 1 and illustrates the comparison between the exact
16 148 F. Mirzaee solution and numerical solution given by the proposed method (IBPFs) for different values of n. Also, we compare the infinity-norm of absolute error computed by the present method, block pulse functions method (BPFs) [13] and rationalized Haar functions method (RHFs) [15] in Table 3. Furthermore, for this example, numerical convergence rate for proposed method is tabulated in Table 4. Graph of the absolute errors for f 1 (x), f (x) and f (x) are shown in Figures 1, and 3, respectively. Table : Numerical results for Example 1 (n = 3, 64). Nodes x Present method n = 3 n = 64. (1.4,.346,.797) (1.1,.117,.395).1 (1.19,.8133,.11148) (1.5,.817,.11137). (1.,.5659,.5797) (1.5,.694,.8569).3 (1.6,.9376,.58) (1.6,.8965,.47319).4 (1.31,3.1891,.73657) (1.8,3.1879,.73639).5 (1.37,3.51,1.33) (1.9,3.55,1.8).6 (1.47, ,1.998) (1.1,3.7813,1.9894).7 (1.61,4.691, ) (1.16,4.1948,1.696).8 (1.89, ,.1338) (1.1,4.3978,.6713).9 (1.13, ,.551) (1.31, ,.54917) Table 3: Approximate L -norm of absolute error for Example 1. Methods E,n Method of [15] n = 8 (1.5e-,1.7e-1,1.31e-1) n = 16 (3.e-3,5.68e-,6.11e-) n = 3 (6.93e-4,3.1e-,3.16e-) n = 64 (1.87e-4,.1e-,1.96e-) Method of [13] n = 8 (1.6e-,1.9e-1,1.16e-1) n = 16 (.67e-3,5.4e-,5.8e-) n = 3 (6.67e-4,.71e-,.9e-) n = 64 (1.67e-4,1.4e-,1.48e-) Present method n = 8 (8.89e-3,1.3e-1,1.9e-1) n = 16 (.44e-3,5.9e-,5.63e-) n = 3 (6.37e-4,.67e-,.86e-) n = 64 (1.63e-4,1.35e-,1.44e-)
17 Numerical solution of system of linear integral equations 149 Table 4: Numerical convergence rate for Example 1. m E n,n for f 1 (x) for f (x) for f 3 (x) Error for m=8 Error for m=16 Error for m=3 Error for m=64.1 Error x Figure 1: Absolute value of error, Example 1 for f 1 (x). Example. Consider the following system of linear Volterra integral equations [4]: x f 1 (x) = 1 x + (f 1 (y) + ye y f (y)) dy, f (x) = 1 + x + x ( ye y f 1 (y) f (y) ) dy, where x [, 1). The exact solution is f 1 (x) = e x and f (x) = e x. Table 5 illustrates the comparison between the exact solution and numerical solution given by the proposed method (IBPFs) for different values of n. Also, we compare the infinity-norm of absolute error computed by the present method, BPFs [4] and RHFs [15] method in Table 6. Furthermore, for this example, numerical convergence rate for proposed method is tabulated in Table 7. Graph of the absolute errors for f 1 (x) and f (x) are shown in Figures 4 and 5, respectively.
18 15 F. Mirzaee Error for m=8 Error for m=16 Error for m=3 Error for m=64 Error x Figure : Absolute value of error, Example 1 for f (x) Error for m=8 Error for m=16 Error for m=3 Error for m=64 Error x Figure 3: Absolute value of error, Example 1 for f 3 (x). Table 5: Numerical results for Example. Nodes x, Exact Present method n = 8 n = 16 n = 3 n = 64., (1.,1.) (1.36,.96939) (1.1596,.98454) (1.789,.993) (1.393,.9961).1, (1.1517,.9484) ( ,.877) (1.141,.8778) (1.118,.9839) (1.9984,.9931)., (1.14,.81873) (1.31,.763) (1.1968,.8151) (1.138,.8471) (1.9,.81371).3, ( ,.748) (1.31,.763) (1.396,.71933) ( ,.7491) ( ,.73977).4, (1.4918,.673) (1.516,.6693) ( ,.6797) ( ,.65795) (1.5114,.6619).5, (1.6487,.6653) (1.747,.57499) ( ,.58868) (1.674,.5975) (1.6649,.6165).6, (1.81,.54881) (.8,.49749) ( ,.5144) (1.8465,.54164) (1.856,.54681).7, (.1375,.49659) (.9115,.494) (.631,.486) (.531,.4911) (.3877,.4891).8, (.554,.44933) (.9115,.494) (.3475,.4191) (.9969,.4391) (.486,.4446).9, (.4596,.4657) (.61538,.36847) (.5313,.3997) (.5876,.3937) (.5165,.3979)
19 Numerical solution of system of linear integral equations 151 Table 6: Approximate L -norm of absolute error for Example. Methods E,n Method of [15] n = 8 (6.47e-,7.89e-) n = 16 (8.13e-,4.7e-) n = 3 (4.1e-,.7e-) n = 64 (.34e-,6.8e-) Method of [4] n = 8 (5.61e-,7.17e-) n = 16 (7.94e-,.89e-) n = 3 (3.71e-,1.35e-) n = 64 (.3e-,5.97e-3) Present method n = 8 (1.44e-1,4.e-) n = 16 (7.13e-,.17e-) n = 3 (3.54e-,1.13e-) n = 64 (1.77e-,5.74e-3) Table 7: Numerical convergence rate for Example. m for f 1 (x) for f (x) E n,n Example 3. Consider the following system of linear Fredholm integral equations [13]: f 1 (x) = 11 6 x (x + y)f 1 (y)dy f (x) = 5 4 x x xy f 1 (y)dy (x y )f (y)dy, x yf (y)dy, where x [, 1). The exact solution is f 1 (x) = x and f (x) = x. Table 8 illustrates the comparison between the exact solution and numerical solution given by the proposed method (IBPFs) for different values of n. Also,
20 15 F. Mirzaee Error for m=8 Error for m=16 Error for m=3 Error for m=64 Error x Figure 4: Absolute value of error, Example for f 1 (x)..1.8 Error for m=8 Error for m=16 Error for m=3 Error for m=64.6 Error x Figure 5: Absolute value of error, Example for f (x). we compare the infinity-norm of absolute error computed by the present method, BPFs [13] and RHFs [15] method in Table 9. Furthermore, for this example, numerical convergence rate for proposed method is tabulated in Table 1. Graph of the absolute errors for f 1 (x) and f (x) are shown in Figures 6 and 7 respectively. Example 4. Consider the following system of linear Fredholm integral equations [13]: f 1 (x) = e x + ex+1 1 x+1 e x y f 1 (y)dy e (x+)y f (y)dy, f (x) = e x + e x + ex+1 1 x+1 e xy f 1 (y)dy e x+y f (y)dy,
21 Numerical solution of system of linear integral equations Error for m=8 Error for m=16 Error for m=3 Error for m=64.6 Error x Figure 6: Absolute value of error, Example 3 for f 1 (x)..1.8 Error for m=8 Error for m=16 Error for m=3 Error for m=64.6 Error x Figure 7: Absolute value of error, Example 3 for f (x). Table 8: Numerical results for Example 3. Nodes x, Exact Present method n = 8 n = 16 n = 3 n = 64., (.,.) (.34,.133) (.1638,.33) (.81,.8) (.396,.).1, (.1,.1) (.1755,.174) (.1568,.1598) (.9393,.888) (.938,.881)., (.,.4) (.56,.645) (.18815,.3553) (.18767,.355) (.317,.418).3, (.3,.9) (.56,.645) (.3137,.987) (.3165,.9776) (.9691,.8816).4, (.4,.16) (.37697,.1433) (.37553,.1416) (.4638,.16515) (.468,.1657).5, (.5,.5) (.5168,.5189) (.545,.548) (.51,.51) (.53,.53).6, (.6,.36) (.6639,.397) (.6537,.39116) (.59385,.3567) (.59378,.3557).7, (.7,.49) (.7511,.5648) (.68784,.473) (.68759,.478) (.7315,.4944).8, (.8,.64) (.7511,.5648) (.8176,.6678) (.8157,.6631) (.79689,.6355).9, (.9,.81) (.87581,.768) (.875,.7668) (.963,.8146) (.966,.8133)
22 154 F. Mirzaee Table 9: Approximate L -norm of absolute error for Example 3. Methods E,n Method of [15] n = 8 (3.98e-,4.44e-) n = 16 (.1e-,.58e-) n = 3 (1.4e-,1.53e-) n = 64 (1.77e-,5.74e-3) Method of [13] n = 8 (3.61e-,4.16e-) n = 16 (1.8e-,.8e-) n = 3 (9.e-3,1.4e-) n = 64 (4.51e-3,5.1e-3) Present method n = 8 (3.44e-,3.88e-) n = 16 (1.76e-,.1e-) n = 3 (8.9e-3,1.e-) n = 64 (4.48e-3,5.16e-3) Table 1: Numerical convergence rate for Example 3. m for f 1 (x) for f (x) E n,n where x [, 1). The exact solution is f 1 (x) = e x and f (x) = e x. Table 11 illustrates the comparison between the exact solution and numerical solution given by the proposed method (IBPFs) for different values of n. Also, we compare the infinity-norm of absolute error computed by the present method, BPFs [13] and RHFs [15] method in Table 1. Furthermore, for this example, numerical convergence rate for proposed method is tabulated in Table 13. Graph of the absolute errors for f 1 (x) and f (x) are shown in Figures 8 and 9, respectively.
23 Numerical solution of system of linear integral equations 155 Table 11: Numerical results for Example 4. Nodes x, Exact Present method n = 8 n = 16 n = 3 n = 64., (1.,1.) (1.197,.975) (1.157,.9863) (1.73,.9961) (1.371,.996).1, (1.1517,.9484) (1.116,.88873) (1.199,.88415) (1.9747,.9193) (1.988,.916)., (1.14,.81873) (1.7114,.78496) (1.9,.8367) (1.537,.8945) (1.51,.8168).3, ( ,.748) (1.7114,.78496) (1.3633,.7333) (1.3659,.733) (1.3454,.7434).4, (1.4918,.673) ( ,.69334) ( ,.68889) (1.5,.66655) (1.593,.6665).5, (1.6487,.6653) ( ,.6143) ( ,.689) (1.6477,.6693) ( ,.6663).6, (1.81,.54881) (1.8579,.5495) (1.8641,.53677) (1.897,.5564) (1.815,.5535).7, (.1375,.49659) (.16,.47776) ( ,.5431) ( ,.531) (.1977,.49513).8, (.554,.44933) (.16,.47776) (.491,.44514) (.537,.44411) (.183,.4583).9, (.4596,.4657) (.38148,.4189) (.394,.418) (.4738,.4437) (.47471,.441) Table 1: Approximate L -norm of absolute error for Example 4. Methods E,n Method of [15] n = 8 (6.94e-,.85e-) n = 16 (3.46e-,1.33e-) n = 3 (1.88e-,6.17e-3) n = 64 (8.36e-3,3.1e-3) Method of [13] n = 8 (6.69e-,.45e-) n = 16 (3.6e-,1.e-) n = 3 (1.6e-,5.94e-3) n = 64 (8.7e-3,.98e-3) Present method n = 8 (6.6e-,.9e-) n = 16 (3.15e-,1.16e-) n = 3 (1.59e-,5.85e-3) n = 64 (8.e-3,.94e-3) Table 13: Numerical convergence rate for Example 4. m for f 1 (x) for f (x) E n,n
24 156 F. Mirzaee Error for m=8 Error for m=16 Error for m=3 Error for m=64.1 Error x Figure 8: Absolute value of error, Example 4 for f 1 (x) Error for m=8 Error for m=16 Error for m=3 Error for m=64 Error x Figure 9: Absolute value of error, Example 4 for f (x). 6 Conclusion The IBPFs, together with operational matrices of integration P 1 and P, are used to obtain the solution of system of linear Volterra and Fredholm integral equations. The present method reduces a system of linear Volterra and Fredholm integral equations into a system of algebraic equations. The matrices P 1 and P have many zeros, hence the method is computationally very attractive. Also, we have shown that our approach is convergence and its order of convergent is O( 1 n ). For the approximate solution of the given examples, we plot absolute value of error to discover applicability and accuracy of the proposed method. Comparisons of the results of applying IBPFs, BPFs and RHFs methods reveals that the new
25 Numerical solution of system of linear integral equations 157 technique is effective and convenient. This method can be easily extended and applied for solving systems of nonlinear Volterra and Fredholm integral equations. Acknowledgements The author appreciates the support provided by Malayer University and thank the associate editor and referees for their comments and suggestions which have improved the paper. References [1] A. Arikoglu and I. Ozkol, Solutions of integral and integro-differential equation systems by using differential transform method, Comput. Math. Appl. 56 (8) [] K.E. Atkinson, The numerical solution of integral equations of the second kind, Cambridge University Press, [3] E. Babolian, J. Biazar and A.R. Vahidi, The decomposition method applied to systems of Fredholm integral equations of the second kind, Appl. Math. Comput. 148 (4) [4] V. Balakumar and K. Murugesan, Numerical solution of systems of linear Volterra integral equations using block-pulse functions, Malaya J. Matematik (13) [5] J. Biazar, E. Babolian and R. Islam, Solution of a system of Volterra integral equations of the first kind by Adomian method, Appl. Math. Comput. 139 (3) [6] M.C. Bonis and C. Laurita, Numerical treatment of second kind Fredholm integral equations systems on bounded intervals, J. Comput. Appl. Math. 17 (8) [7] L.M. Delves and J.L. Mohamed, Computational methods for integral equations, Cambridge University Press, [8] M. Javidi, Modified homotopy perturbation method for solving system of linear Fredholm integral equations, Math. Comput. Model. 5 (9)
26 158 F. Mirzaee [9] Z. H. Jiang and W. Schaufelberger, Block Pulse functions and their applications in control systems, Spriger-Verlag, Berlin, 199. [1] K. Maleknejad and F. Mirzaee, Numerical solution of linear Fredholm integral equations system by rationalized Haar functions method, Int. J. Comput. Math. 8 (3) [11] K. Maleknejad, F. Mirzaee and S. Abbasbandy, Solving linear integrodifferential equations system by using rationalized Haar functions method, Appl. Math. Comput. 155 (4) [1] K. Maleknejad and M. Tavassoli Kajani, Solving linear integrodifferential equation system by Galerkin methods with hybrid functions, Appl. Math. Comput. 159 (4) [13] K. Maleknejad, M. Shahrezaee and H. Khatami, Numerical solution of integral equations system of the second kind by block-pulse functions, Appl. Math. Comput. 166 (5) [14] K. Maleknejad, S. Sohrabi and B. Baranji, Application of D-BPFs to nonlinear integral equations, Commun. Nonlinear Sci. Numer. Simulat. 15 (1) [15] F. Mirzaee, Numerical computational solution of the linear Volterra integral equations system via rationalized Haar functions, J. King Saud Uni. Sci. (1) [16] J. Pour-Mahmoud, M.Y. Rahimi-Ardabili and S. Shahmorad, Numerical solution of the system of Fredholm integro-differential equations by the Tau method, Appl. Math. Comput. 168 (5) [17] J. Rashidinia and M. Zarebnia, Convergence of approximate solution of system of Fredholm integral equations, J. Math. Anal. Appl. 333 (7) [18] W. Rudin, Principles of Mathematical Analysis, McGraw-Hill, Singapore, [19] J. Saberi-Nadjafi and M. Tamamgar, The variational iteration method, a highly promising method for solving the system of integro differential equations, Comput. Math. Appl. 56 (8) [] R.K. Saeed and C.S. Ahmed, Numerical Solution of the System of linear Volterra Integral Equations of the Second Kind using Monte- Carlo Method, J. Appl. Sci. Res. 4 (8)
27 Numerical solution of system of linear integral equations 159 [1] N. Sahin, S. Yuzba and M. Gulsu, A collocation approach for solving systems of linear Volterra integral equations with variable coefficients, Appl. Math. Comput. 6 (11) [] H.H. Sorkun and S. Yalinbas, Approximate solutions of linear Volterra integral equation systems with variable coefficients, Appl. Math. Model. 34 (1) [3] E. Yusufoglu, Numerical expansion methods for solving systems of linear integral equations using interpolation and quadrature rules, Int. J. Comput. Math. 84 (7) [4] E. Yusufoglu, An efficient algorithm for solving integro-differential equations system, Appl. Math. Comput. 19 (7) [5] E. Yusufoglu, A homotopy perturbation algorithm to solve a system of Fredholm-Volterra type integral equations, Math. Comput. Model. 47 (8)
Solving 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 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 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 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 informationBLOCK-PULSE FUNCTIONS AND THEIR APPLICATIONS TO SOLVING SYSTEMS OF HIGHER-ORDER NONLINEAR VOLTERRA INTEGRO-DIFFERENTIAL EQUATIONS
Electronic Journal of Differential Equations, Vol. 214 (214), No. 54, pp. 1 9. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu BLOCK-PULSE FUNCTIONS
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 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 informationOperational Matrices for Solving Burgers' Equation by Using Block-Pulse Functions with Error Analysis
Australian Journal of Basic and Applied Sciences, 5(12): 602-609, 2011 ISSN 1991-8178 Operational Matrices for Solving Burgers' Equation by Using Block-Pulse Functions with Error Analysis 1 K. Maleknejad,
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 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 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 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 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 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 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 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 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 informationApproximate solution of linear integro-differential equations by using modified Taylor expansion method
ISSN 1 746-7233, Engl, UK World Journal of Modelling Simulation Vol 9 (213) No 4, pp 289-31 Approximate solution of linear integro-differential equations by using modified Taylor expansion method Jalil
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 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 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 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 informationImproving homotopy analysis method for system of nonlinear algebraic equations
Journal of Advanced Research in Applied Mathematics Vol., Issue. 4, 010, pp. -30 Online ISSN: 194-9649 Improving homotopy analysis method for system of nonlinear algebraic equations M.M. Hosseini, S.M.
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 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 informationImproving the Accuracy of the Adomian Decomposition Method for Solving Nonlinear Equations
Applied Mathematical Sciences, Vol. 6, 2012, no. 10, 487-497 Improving the Accuracy of the Adomian Decomposition Method for Solving Nonlinear Equations A. R. Vahidi a and B. Jalalvand b (a) Department
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More informationNEW ITERATIVE METHODS BASED ON SPLINE FUNCTIONS FOR SOLVING NONLINEAR EQUATIONS
Bulletin of Mathematical Analysis and Applications ISSN: 181-191, URL: http://www.bmathaa.org Volume 3 Issue 4(011, Pages 31-37. NEW ITERATIVE METHODS BASED ON SPLINE FUNCTIONS FOR SOLVING NONLINEAR EQUATIONS
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 informationError formulas for divided difference expansions and numerical differentiation
Error formulas for divided difference expansions and numerical differentiation Michael S. Floater Abstract: We derive an expression for the remainder in divided difference expansions and use it to give
More informationHomotopy perturbation method for solving hyperbolic partial differential equations
Computers and Mathematics with Applications 56 2008) 453 458 wwwelseviercom/locate/camwa Homotopy perturbation method for solving hyperbolic partial differential equations J Biazar a,, H Ghazvini a,b a
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 informationNumerical Solution of Nonlocal Parabolic Partial Differential Equation via Bernstein Polynomial Method
Punjab University Journal of Mathematics (ISSN 116-2526) Vol.48(1)(216) pp. 47-53 Numerical Solution of Nonlocal Parabolic Partial Differential Equation via Bernstein Polynomial Method Kobra Karimi Department
More informationExact Solutions for Systems of Volterra Integral Equations of the First Kind by Homotopy Perturbation Method
Applied Mathematical Sciences, Vol. 2, 28, no. 54, 2691-2697 Eact Solutions for Systems of Volterra Integral Equations of the First Kind by Homotopy Perturbation Method J. Biazar 1, M. Eslami and H. Ghazvini
More informationDifferential Transform Method for Solving. Linear and Nonlinear Systems of. Ordinary Differential Equations
Applied Mathematical Sciences, Vol 5, 2011, no 70, 3465-3472 Differential Transform Method for Solving Linear and Nonlinear Systems of Ordinary Differential Equations Farshid Mirzaee Department of Mathematics
More informationNumerical Solution of Linear Fredholm Fuzzy Integral Equation of the Second Kind by Block-pulse Functions
Australian Journal of Basic and Applied Sciences, 3(3): 2637-2642, 2009 ISSN 1991-8178 Numerical Solution of Linear Fredholm Fuzzy Integral Equation of the Second Kind by Block-pulse Functions 1 2 3 M.
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 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 informationMax-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig
Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig A CAS Wavelet Method for Solving Nonlinear Fredholm Integro- Differential Equations of Fractional Order by Habib Allah Saeedi, Nasibeh
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 informationSolving Integral Equations of the Second Kind by Using Wavelet Basis in the PG Method
5 1 July 24, Antalya, Turkey Dynamical Systems and Applications, Proceedings, pp. 515 52 Solving Integral Equations of the Second Kind by Using Wavelet Basis in the PG Method K. Maleknejad Department of
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
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 informationBernstein operational matrices for solving multiterm variable order fractional differential equations
International Journal of Current Engineering and Technology E-ISSN 2277 4106 P-ISSN 2347 5161 2017 INPRESSCO All Rights Reserved Available at http://inpressco.com/category/ijcet Research Article Bernstein
More informationACTA UNIVERSITATIS APULENSIS No 18/2009 NEW ITERATIVE METHODS FOR SOLVING NONLINEAR EQUATIONS BY USING MODIFIED HOMOTOPY PERTURBATION METHOD
ACTA UNIVERSITATIS APULENSIS No 18/2009 NEW ITERATIVE METHODS FOR SOLVING NONLINEAR EQUATIONS BY USING MODIFIED HOMOTOPY PERTURBATION METHOD Arif Rafiq and Amna Javeria Abstract In this paper, we establish
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 informationAn Efficient Numerical Method for Solving. the Fractional Diffusion Equation
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 1-12 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 An Efficient Numerical Method for Solving the Fractional
More informationarxiv: v1 [math.na] 19 May 2015
arxiv:505.04855v [math.na] 9 May 205 Approximation solution of two-dimensional linear stochastic Volterra integral equation by applying the Haar wavelet M. Fallahpour, M. Khodabin 2, K. Maleknejad 3 Department
More informationImprovements in Newton-Rapshon Method for Nonlinear Equations Using Modified Adomian Decomposition Method
International Journal of Mathematical Analysis Vol. 9, 2015, no. 39, 1919-1928 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2015.54124 Improvements in Newton-Rapshon Method for Nonlinear
More 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 informationDepartment of Applied Mathematics and Theoretical Physics. AMA 204 Numerical analysis. Exam Winter 2004
Department of Applied Mathematics and Theoretical Physics AMA 204 Numerical analysis Exam Winter 2004 The best six answers will be credited All questions carry equal marks Answer all parts of each question
More information6.1 Matrices. Definition: A Matrix A is a rectangular array of the form. A 11 A 12 A 1n A 21. A 2n. A m1 A m2 A mn A 22.
61 Matrices Definition: A Matrix A is a rectangular array of the form A 11 A 12 A 1n A 21 A 22 A 2n A m1 A m2 A mn The size of A is m n, where m is the number of rows and n is the number of columns The
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 informationA Maple program for computing Adomian polynomials
International Mathematical Forum, 1, 2006, no. 39, 1919-1924 A Maple program for computing Adomian polynomials Jafar Biazar 1 and Masumeh Pourabd Department of Mathematics, Faculty of Science Guilan University
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 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 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 informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationA Numerical Solution of Volterra s Population Growth Model Based on Hybrid Function
IT. J. BIOAUTOMATIO, 27, 2(), 9-2 A umerical Solution of Volterra s Population Growth Model Based on Hybrid Function Saeid Jahangiri, Khosrow Maleknejad *, Majid Tavassoli Kajani 2 Department of Mathematics
More 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 informationAn approximation to the solution of parabolic equation by Adomian decomposition method and comparing the result with Crank-Nicolson method
International Mathematical Forum, 1, 26, no. 39, 1925-1933 An approximation to the solution of parabolic equation by Adomian decomposition method and comparing the result with Crank-Nicolson method J.
More informationA Maple Program for Solving Systems of Linear and Nonlinear Integral Equations by Adomian Decomposition Method
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 29, 1425-1432 A Maple Program for Solving Systems of Linear and Nonlinear Integral Equations by Adomian Decomposition Method Jafar Biazar and Masumeh
More informationLECTURE 16 GAUSS QUADRATURE In general for Newton-Cotes (equispaced interpolation points/ data points/ integration points/ nodes).
CE 025 - Lecture 6 LECTURE 6 GAUSS QUADRATURE In general for ewton-cotes (equispaced interpolation points/ data points/ integration points/ nodes). x E x S fx dx hw' o f o + w' f + + w' f + E 84 f 0 f
More information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
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 informationWavelets in Scattering Calculations
Wavelets in Scattering Calculations W. P., Brian M. Kessler, Gerald L. Payne polyzou@uiowa.edu The University of Iowa Wavelets in Scattering Calculations p.1/43 What are Wavelets? Orthonormal basis functions.
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 informationA new approach to solve fuzzy system of linear equations by Homotopy perturbation method
Journal of Linear and Topological Algebra Vol. 02, No. 02, 2013, 105-115 A new approach to solve fuzzy system of linear equations by Homotopy perturbation method M. Paripour a,, J. Saeidian b and A. Sadeghi
More informationNumerical analysis of the one-dimensional Wave equation subject to a boundary integral specification
Numerical analysis of the one-dimensional Wave equation subject to a boundary integral specification B. Soltanalizadeh a, Reza Abazari b, S. Abbasbandy c, A. Alsaedi d a Department of Mathematics, University
More informationNumerical Solutions of Volterra Integral Equations Using Galerkin method with Hermite Polynomials
Proceedings of the 3 International Conference on Applied Mathematics and Computational Methods in Engineering Numerical of Volterra Integral Equations Using Galerkin method with Hermite Polynomials M.
More informationA Product Integration Approach Based on New Orthogonal Polynomials for Nonlinear Weakly Singular Integral Equations
Acta Appl Math (21) 19: 861 873 DOI 1.17/s144-8-9351-y A Product Integration Approach Based on New Orthogonal Polynomials for Nonlinear Weakly Singular Integral Equations M. Rasty M. Hadizadeh Received:
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 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 informationSolution of Linear System of Partial Differential Equations by Legendre Multiwavelet Andchebyshev Multiwavelet
International Journal of Scientific and Innovative Mathematical Research (IJSIMR) Volume 2, Issue 12, December 2014, PP 966-976 ISSN 2347-307X (Print) & ISSN 2347-3142 (Online) www.arcjournals.org Solution
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 informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
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 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 informationdy = f( x) dx = F ( x)+c = f ( x) dy = f( x) dx
Antiderivatives and The Integral Antiderivatives Objective: Use indefinite integral notation for antiderivatives. Use basic integration rules to find antiderivatives. Another important question in calculus
More informationarxiv: v1 [math.na] 20 May 2013
A numerical method based on the reproducing kernel Hilbert space method for the solution of fifth-order boundary-value problems arxiv:1305.4445v1 [math.na] 0 May 013 Mustafa Inc, Ali Akgül and Mehdi Dehghan
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 informationElectromagnetic wave scattering by many small bodies and creating materials with a desired refraction coefficient
168 Int. J. Computing Science and Mathematics, Vol. Electromagnetic wave scattering by many small bodies and creating materials with a desired refraction coefficient A.G. Ramm Mathematics epartment, Kansas
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
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 informationSolutions: Problem Set 4 Math 201B, Winter 2007
Solutions: Problem Set 4 Math 2B, Winter 27 Problem. (a Define f : by { x /2 if < x
More informationNumerical methods. Examples with solution
Numerical methods Examples with solution CONTENTS Contents. Nonlinear Equations 3 The bisection method............................ 4 Newton s method.............................. 8. Linear Systems LU-factorization..............................
More informationA New Method for Solving General Dual Fuzzy Linear Systems
Journal of Mathematical Extension Vol. 7, No. 3, (2013), 63-75 A New Method for Solving General Dual Fuzzy Linear Systems M. Otadi Firoozkooh Branch, Islamic Azad University Abstract. According to fuzzy
More informationImproved Newton s method with exact line searches to solve quadratic matrix equation
Journal of Computational and Applied Mathematics 222 (2008) 645 654 wwwelseviercom/locate/cam Improved Newton s method with exact line searches to solve quadratic matrix equation Jian-hui Long, Xi-yan
More informationExistence of Positive Periodic Solutions of Mutualism Systems with Several Delays 1
Advances in Dynamical Systems and Applications. ISSN 973-5321 Volume 1 Number 2 (26), pp. 29 217 c Research India Publications http://www.ripublication.com/adsa.htm Existence of Positive Periodic Solutions
More informationReal Analysis Qualifying Exam May 14th 2016
Real Analysis Qualifying Exam May 4th 26 Solve 8 out of 2 problems. () Prove the Banach contraction principle: Let T be a mapping from a complete metric space X into itself such that d(tx,ty) apple qd(x,
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 informationWORKSHEET FOR THE PUTNAM COMPETITION -REAL ANALYSIS- lim
WORKSHEET FOR THE PUTNAM COMPETITION -REAL ANALYSIS- INSTRUCTOR: CEZAR LUPU Problem Let < x < and x n+ = x n ( x n ), n =,, 3, Show that nx n = Putnam B3, 966 Question? Problem E 334 from the American
More informationExistence and Convergence Results for Caputo Fractional Volterra Integro-Differential Equations
J o u r n a l of Mathematics and Applications JMA No 41, pp 19-122 (218) Existence and Convergence Results for Caputo Fractional Volterra Integro-Differential Equations Ahmed A. Hamoud*, M.Sh. Bani Issa,
More informationAbsolute value equations
Linear Algebra and its Applications 419 (2006) 359 367 www.elsevier.com/locate/laa Absolute value equations O.L. Mangasarian, R.R. Meyer Computer Sciences Department, University of Wisconsin, 1210 West
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
More informationWORKSHEET FOR THE PUTNAM COMPETITION -REAL ANALYSIS- lim
WORKSHEET FOR THE PUTNAM COMPETITION -REAL ANALYSIS- INSTRUCTOR: CEZAR LUPU Problem Let < x < and x n+ = x n ( x n ), n =,, 3, Show that nx n = Putnam B3, 966 Question? Problem E 334 from the American
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 information