A Parameter-uniform Method for Two Parameters Singularly Perturbed Boundary Value Problems via Asymptotic Expansion
|
|
- Clarissa Hoover
- 6 years ago
- Views:
Transcription
1 Appl. Math. Inf. Sci. 7, No. 4, (2013) 1525 Applied Mathematics & Information Sciences An International Journal A Parameter-uniform Method for Two Parameters Singularly Perturbed Boundary Value Problems via Asymptotic Expansion D. Kumar 1,, A. S. Yadaw 2 and M. K. Kadalbajoo 3 1 Department of Mathematics, B.I.T.S. Pilani, , India 2 Department of Pharmacology & Systems Biology Mount Sinai School of Medicine, New York, NY Department of Mathematics & Statistics, I.I.T. Kanpur, , India Received: 11 Dec. 2012, Revised: 13 Jan. 2013, Accepted: 4 Feb Published online: 1 Jul Abstract: An approximate method for two parameters singularly perturbed boundary value problems having boundary layers at both end points is given. The method is motivated by the asymptotic behavior of the solution. In the outer region the solution of the problem is approximated by the zeroth order asymptotic expansion while in the inner region the solution of the problem is obtained by using B- spline collocation method. The method is iterated on the transition point of the boundary layer region. To demonstrate the applicability of the method two test examples are considered. Keywords: Singular perturbation, two parameters, boundary layer, asymptotic expansion, B-spline collocation method. 1. Introduction The boundary value problems for ordinary differential equations in which one or more small positive parameters multiplying the derivatives arise in the field of physics and applied mathematics. The problems with one small positive parameter multiplying to highest derivative have considered by many authors [1]. Motivated by asymptotic expansion [2, 3] a second order boundary value problem with two small parameters multiplying to the highest and second highest derivative is considered. This type of problems arise in chemical reactor theory, engineering, biology, lubrication theory etc. Consider the problem Ly(x) ǫy (x) µa(x)y (x)+b(x)y(x) = f(x), (1) x [0,1], with y(0) = α, y(1) = β, (2) where ǫ and µ are small positive parameters satisfying ǫ/µ 2 0 as µ 0. The functions a(x), b(x) and f(x) are sufficiently smooth satisfying a(x) a > 0, b(x) b > 0. In 1967, Malley gave the asymptotic solution of two parameters singularly perturbed boundary value problems and demonstrate the roll of ǫ and µ on the solution. After three decades, various mathematician [4, 5, 6] gave the numerical solution of two parameters problems. Malley [2,3,7,8,9] examined the nature of asymptotic solution of the continuous problem where the ratio of µ 2 to ǫ was identified as significant. In [10, 11], the standard upwind finite difference operator on two different choices of Shishkin mesh was shown to be parameter-uniform of first order. In [12] parameter-uniform methods on a uniform mesh were constructed. Vulanović [13], used the higher order finite difference scheme on a piecewise uniform mesh both of Shishkin and Bakhavalov type for solving quasi-linear boundary value problems with small parameters. 2. Solution of the problem 2.1. Asymptotic solution Consider the asymptotic expansion of (1) with (2) of the form y(x,ǫ,µ) = (y 0 +(ǫ/µ)y 1 +(ǫ/µ) 2 y 2 +O((ǫ/µ) 3 )) (3) Corresponding author dkumar2845@gmail.com
2 1526 D. Kumar et al. : A Parameter-uniform Method for Two Parameters... +(z 0 +(ǫ/µ)z 1 +(ǫ/µ) 2 z 2 +O((ǫ/µ) 3 )) (4) = (y 0 +z 0 )+(ǫ/µ)(y 1 +z 1 ) (5) +(ǫ/µ) 2 (y 2 +z 2 )+O((ǫ/µ) 3 ) (6) = W 0 +(ǫ/µ)w 1 +(ǫ/µ) 2 W 2 +O((ǫ/µ) 3 ), (7) where the zeroth order asymptotic expansion W 0 is given by W 0 (x) = y 0 + z 0, where y 0 is the solution of the reduced problem µa(x)y 0(x)+b(x)y 0 (x) = f(x), y 0 (1) = β, (8) and y 1 (x) is the solution given by ǫy 1(x) µa(x)y 1(x)+b(x)y 1 (x) = µy 0(x), (9) y 1 (0) = 0, y 1 (1) = (ǫ/µ) 1 z 0 (1), (10) and z 0 is a layer correction given by d 2 z 0 /dτ 2 a(0)dz 0 /dτ = 0, z 0 (0) = α y 0 (0),z 0 ( ) = 0, (11) where τ = x/(ǫ/µ). The solution of Eq. (8) is given by ( 1 ) y 0 (x) = exp [ β x 1 x b(s) µa(s) ds f(s) µa(s) exp ( 1 x and the solution of Eq. (11) is given by z 0 (x) = (α y 0 (0))exp( a(0)x/(ǫ/µ)). ) b(s) µa(s) ds (12) ] ds (13), Thus the zeroth order asymptotic expansion W 0 (x) is obtained. It is easy to see that the zeroth order asymptotic expansion W 0 (x) of the problem (1) with (2) satisfies the following inequality: Theorem 1. The zeroth order asymptotic expansion W 0 (x) of the solutiony(x) of (1) with (2) satisfies y(x) W 0 (x) C(ǫ/µ), where ǫ/µ 2 0 asµ 0. Also the solution y(x) of (1) with (2) and its derivatives satisfy the following inequality: Theorem 2. Suppose a(x), b(x) and f(x) are sufficiently smooth and has derivatives at least of order k then we have y (i) (x) C[1+µ i exp( ν(1 x)/µ) (14) where b(x) a(x) ν. +ǫ i exp( αx/(ǫ/µ))], i = 0,1,...,k, (15) 2.2. Discrete solution We divide the interval [0, 1] into three non overlapping subintervals[0, k(ǫ/µ)], (k(ǫ/µ), 1 kµ) and[1 kµ, 1] where k be any positive integer such that kµ 1. The original problem is then divided into three equivalent problems (two inner regions and one outer region problem). To obtain the boundary conditions at the transition points k(ǫ/µ) and (1 kµ), we use the zeroth order asymptotic expansion W 0. Thus if W 0 (k(ǫ/µ)) = γ 1 and W 0 (1 kµ) = γ 2 then the three problems are P 1 :The left inner region problem (x [0,k(ǫ/µ)]) ǫy (x) µa(x)y (x)+b(x)y(x) = f(x), (16) y(0) = α, y(k(ǫ/µ)) = γ 1. (17) P 2 :The outer region problem (x (k(ǫ/µ),1 kµ)) ǫy (x) µa(x)y (x)+b(x)y(x) = f(x), (18) y(k(ǫ/µ)) = γ 1, y(1 kµ) = γ 2. (19) P 3 :The right inner region problem (x [1 kµ,1]) ǫy (x) µa(x)y (x)+b(x)y(x) = f(x), (20) y(1 kµ) = γ 2, y(1) = β. (21) After solving the inner and outer region problems, we combine their solutions to obtain an approximate solution of the problem (1) with (2) in the whole interval[0,1]. We change the value of k until the solution profiles do not differ much from iteration to iteration. For this we use the absolute error criteria y (m+1) (x) y (m) (x) σ, where σ is prescribed tolerance error bound. We use B-spline collocation method in the inner region and forp 1 we have x 0 = 0, x N = kǫ/µ, h 1 = kǫ/µn, x i = x 0 + ih 1, i = 1,2,...,N and for P 3 we have x 0 = 1 kµ, x N = 1, h 2 = kµ/n, x i = x 0 +ih 2, i = 1,2,...,N. For 0 < l 1 < l 2 < 1, we define L 2 [l 1,l 2 ] a vector space of all the square integrable function on [l 1,l 2 ] and let X be the linear subspace of L 2 [l 1,l 2 ]. Now define for i = 0,1,2,...,N (x x i 2 ) 3, x i 2 x x i 1, h 3 +3h 2 (x x i 1 )+3h(x x i 1 ) 2 B i (x) = 1 3(x x i 1 ) 3, x i 1 x x i, h 3 h 3 +3h 2 (x i+1 x)+3h(x i+1 x) 2 (22) 3(x i+1 x) 3, x i x x i+1, (x i+2 x) 3, x i+1 x x i+2, 0, otherwise. Let π = {x 0,x 1,...,x N } be the partition of [l 1,l 2 ]. We introduce four additional knots x 2 < x 1 < x 0 and x N+2 > x > x N. It is easy to check that each of the function B i (x) is twice continuously differentiable on the entire real line. Also 4, ifi = j, B i (x j ) = 1, ifi j = ±1, (23) 0, ifi j = ±2,
3 Appl. Math. Inf. Sci. 7, No. 4, (2013) / and that B i (x) = 0 for x x i+2 and x x i 2. We can also see that 0, if i = j, B i(x j ) = ± 3 h, if i j = ±1, (24) 0, if i j = ±2, and B i (x j ) = 12 h 2, if i = j, 6 h, if i j = ±1, 2 0, if i j = ±2. (25) Each B i (x) is also a piecewise cubic at π and B i (x) X. Let Π = {B 1,B 0,B 1,...,B } and let Φ 3 (π) is span ofπ. ThenΠ is linearly independent on[l 1,l 2 ], thus Φ 3 (π) is(n+3)-dimensional. In factφ 3 (π) is a subspace of X. Let L be a linear operator whose domain is X and whose range is also inx. Now we define S(x) = c i B i (x), (26) where c i are unknown real coefficients. Here we have introduced two extra cubic B-splines, B 1 and B to satisfy the boundary conditions. Thus LS(x i ) = f(x i ), 0 i N, (27) and S(x 0 ) = α, S(x N ) = β. (28) On solving the collocation equations (27) and putting the values of B-spline functions B i and of derivatives at mesh points, we obtain a system of (N +1) linear equations in (N +3) unknowns ( 6ǫ+3µa i h+b i h 2 )c i 1 +(12ǫ+4b i h 2 )c i +( 6ǫ 3µa i h+b i h 2 )c i+1 = f i h 2. (29) The given boundary conditions become c 1 +4c 0 +c 1 = α, c N 1 +4c N +c = β. (30) Thus the Eqs. (29) and (30) lead to an (N + 3) (N + 3) system with(n +3) unknownsc 1,c 0,...,c N,c. Now eliminating c 1 from first equation of (29) and from first equation of (30) (36ǫ 12µa 0 h)c 0 +( 6µa 0 h)c 1 = h 2 f 0 α( 6ǫ+3µa 0 h+b 0 h 2 ). (31) Similarly, eliminating c from the last equation of (29) and from second equation of (30) (6µa N h)c N 1 +(36ǫ+12µa N h)c N = h 2 f N β( 6ǫ 3µa N h+b N h 2 ). (32) Coupling equations (31) and (32) with the second through (N 1)st equations of (29). Thus by the elimination of c 1 and c, we lead to a system Tx N = d N of (N + 1) linear equations in (N + 1) unknowns, where x N = (c 0,c 1,...,c N ) t and d N = (h 2 f 0 α( 6ǫ + 3µa 0 h + b 0 h 2 ),h 2 f 1,...,h 2 f N 1,h 2 f N β( 6ǫ 3µa N h + b N h 2 )). The elements t i,j of the tridiagonal matrixt are given by 36ǫ 12µa 0 h, i = j = 0, 6µa 0 h, i = 0,j = 1, 6ǫ+3µa i h+b i h 2, i = j +1,j = 0(1)N 2, 12ǫ+4b i h 2, i = j = 1(1)N 1, t i,j = (33) 6ǫ 3µa i h+b i h 2, i = j 1,j = 2(1)N, 6µa N h, i = N,j = N 1, 36ǫ+12µa N h, i = j = N, 0, i j > 1. It is easy to see that the matrix T is strictly diagonally dominant and hence nonsingular. Since T is nonsingular, we can solve the systemtx N = d N forc 0,c 1,...,c N and substitute into the boundary conditions (30) to obtain c 1 and c. Lemma 1. The B-splinesB 1,B 0,...,B defined in equation (22), satisfy the inequality B i (x) 10, 0 x 1. Proof. We know that B i(x) B i(x). At any node x i, we have B i = B i 1 + B i + B i+1 = 6 < 10. Also we have B i (x) 4 and B i 1 (x) 4, x i 1 x x i. Similarly B i 2 (x) 1 and B i+1 1, x i 1 x x i. Now for any point x i 1 x x i we have B i(x) = B i 2 + B i 1 + B i + B i Theorem 3. Let S(x) be the collocation approximation from the space of cubic splines Φ 3 (π) to the solution y(x) of the boundary value problem (1) (2). If f C 2 [l 1,l 2 ], then the parameter-uniform error estimate is given by y(x) S(x) Ch 2, where C is a positive constant independent of ǫ and N. Proof. To estimate the error y(x) S(x), let Y n be the unique spline interpolate from Φ 3 (π) to the solution y(x) of our boundary value problem (1) (2). If f(x) C 2 [l 1,l 2 ] then y(x) C 4 [l 1,l 2 ] and it follows from De Boor-Hall error estimates [14] that D j (y(x) Y n ) γ j h 4 j, j = 0,1,2, (34) where h is uniform mesh spacing and γ j s are constants independent ofhand N. Let Y n (x) = p i B i (x). (35)
4 1528 D. Kumar et al. : A Parameter-uniform Method for Two Parameters... It follows immediately from the estimates (34) that LS(x i ) LY n (x i ) = f(x i ) LY(x i ) +Ly(x i ) Ly(x i ) λh 2, (36) where λ = ǫγ 2 + µγ 1 a(x) h + γ 0 b(x) h 2. Also LS(x i ) = Ly(x i ) = f(x i ). Let LY n (x i ) = ˆf n (x i ) for i = 1,2,...,N and ˆf n = (ˆf n (x 0 ), ˆf n (x 1 ),..., ˆf n (x N )) t. Now from the system Tx N = d N and (34), it is clear that the ith coordinate [T(x N y N )] i of T(x N y N ), where y N = (p 0,p 1,...,p N ) t, satisfies the inequality [T(x N y N )] i λh 4. (37) Since (Tx N ) i = h 2 f(x i ) and (Ty N ) i = h 2 ˆf n (x i ) for i = 0,1,2,...,N. The ith coordinate of [T(x N y N )] is theith equation ( 6ǫ+3µa i h+b i h 2 )δ i 1 +(12ǫ+4b i h 2 )δ i +( 6ǫ 3µa i h+b i h 2 )δ i+1 = ξ i, 1 i N 1,(38) where δ i = p i q i, 1 i N + 1 and ξ i = h 2 (f(x i ) ˆf n (x i )), 1 i N 1. Clearly ξ i λh 4. Let ξ = max ξ i. Also consider 1 i N 1 δ = (δ 1,δ 0,...,δ ) t, then we define e i = δ i and ẽ = max e i. Now equation (38) becomes 1 i N (12ǫ+4b i h 2 )δ i ξ i ξ i +(6ǫ b i h 2 )(δ i 1 +δ i+1 ) +3µa i h(δ i+1 δ i 1 ), 1 i N 1. (39) Taking absolute values with sufficiently small h, we obtain (12ǫ+4b i h 2 )e i ξ +2ẽ(6ǫ+3µa i h b i h 2 ). (40) Since 0 < a a(x) and 0 < b b(x), we get (12ǫ+4b h 2 )e i ξ +2ẽ(6ǫ+3µa h b h 2 ) (41) ξ +2ẽ(6ǫ+3µa h+b h 2 ). (42) In particular (12ǫ+4b h 2 )ẽ ξ +2ẽ(6ǫ+3µa h+b h 2 ). (43) Solving forẽ, we obtain(2b h 2 6µa h)ẽ ξ λh 4 or ẽ λh 3 2b h 6µa. (44) Now to estimate e 1,e 0,e N and e, we observe that the first and last equation of the the systemt(x N y N ) = h 2 ( f n ˆf n ) where f n = (f 0,f 1,...,f N ), gives e 0 and e N 2λb h 5 (36ǫ 12a hµ)(2b h 6µa ), (45) 2λb h 5 (36ǫ+12a hµ)(2b h 6µa ). (46) Now e 1 and e can be evaluated using boundary conditions δ 1 = (4δ 0 δ 1 ) and δ = (4δ N δ N 1 ) λh 3 ( 2b h 2 +9ǫ 3a ) hµ e 1 (2b h 6µa ) 9ǫ+3a, (47) hµ and e λh 3 (2b h 6µa ) ( 2b h 2 +9ǫ 3a hµ 9ǫ+3a hµ ). (48) Using valueλ = ǫγ 2 +µγ 1 b(x) h+γ 0 b(x) h 2 and since there exists a constant C such that e = max 1 i {e i} Ch 2. (49) The above inequality enables us to estimate S(x) Y n (x), and therefore y(x) S(x). In particular S(x) Y n (x) = (q i p i )B i (x). (50) Thus S(x) Y n (x) max q i p i B i (x). (51) Combining equations (49), (51) and using Lemma 1 we obtain S Y n Ch 2. Since y Y N γ 0 h 4 and y S y Y N + Y N S, we obtain y S Ch 2. Combining the results we get the required estimate. 3. Test examples and numerical results To demonstrate the efficiency of the method, two numerical examples are considered. Since the exact solution of the considered problems are given so the maximum absolute errors are estimated by using E N,ǫ = max 0 i N yn i S N i, where y N i is the exact solution and S N i is the computed solution. If the exact solution is not known then the maximum absolute errors can be obtained by using the double mesh principle. Example 1. Consider the boundary value problem ǫy (x) µy (x)+y(x) = x, y(0) = 1, y(1) = 0. The exact solution of the problem is given by y(x) = (1+µ)+(1 µ)em2 e m2 e m1 e m1x (52) + (1+µ)+(1 µ)em1 e m1 e m2 e m2x +x+µ, (53)
5 Appl. Math. Inf. Sci. 7, No. 4, (2013) / Table 1: Maximum absolute error at nodal points for Example 1 for ǫ = 10 6 and µ = 10 2 k Nodes E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E-11 O max E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E+00
6 1530 D. Kumar et al. : A Parameter-uniform Method for Two Parameters... Table 2: Maximum absolute error at nodal points for Example 2 for ǫ = 10 6 and µ = 10 2 k Nodes E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E+00 O max E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E+00
7 Appl. Math. Inf. Sci. 7, No. 4, (2013) / where m 1 = µ µ 2 +4ǫ 2ǫ, m 2 = µ+ µ 2 +4ǫ 2ǫ. For this problem the zeroth order asymptotic expansion W 0 is given by W 0 (x) = x exp((x 1)/µ) (54) + (1 exp( 1/µ))exp( x/(ǫ/µ)). (55) Example 2. Consider the boundary value problem ǫy (x) 2µy (x)+4y(x) = 1, y(0) = 0, y(1) = 1. The exact solution of the problem is given by 3+e m2 y(x) = 4(e m1 e m2 ) em1x 3+em1 4(e m1 e m2 ) em2x where m 1 = µ µ 2 +4ǫ ǫ, m 2 = µ+ µ 2 +4ǫ ǫ. For this example the zeroth order asymptotic expansion W 0 is given by W 0 (x) = (1/4)[1 exp( 2x/(ǫ/µ))] (56) + (3/4)[exp(2x/µ) exp( 2x/(ǫ/µ))] exp( 2/µ).(57) 4. Conclusion In this paper, we present an approximate method based on asymptotic expansion for two parameters singularly perturbed boundary value problems. The B-spline method is used in the inner region and for rest of the region we use the zeroth order asymptotic expansion approximation. To demonstrate the applicability of the method two numerical examples have been considered. For inner regions the absolute errors at nodal points are presented in the tables while for the outer region the maximum absolute error denoted by O max is presented at the middle of the table. It can be seen from the tables that the numerical solutions are very closed to exact solutions. The solution profiles for the considered examples for a fix µ and different values of ǫ are given in Figures 1 and ε = 10 3 ε = 10 4 ε = ε = 10 2 ε = 10 3 ε = Figure 2: Solution profile for Example 2 for µ = 10 4 using N = 32. References [1] S.M. Roberts, J. Math. Anal. Appl. 113, 411 (1988). [2] R.E. O Malley Jr., J. Math. Anal. Appl. 19, 291 (1967). [3] R.E. O Malley Jr., J. Math. Mech. 16, 1143 (1967). [4] small H.G. Roos and Z. Uzelac, Comput. Methods Appl. Math. 3, 443 (2003). [5] J.L. Gracia, E. O Riordan and M.L. Pickett, Appl. Numer. Math. 56, 962 (2006). [6] T. Linßand H.G. Roos, J. Math. Anal. Appl. 289, 355 (2004). [7] R.E. O Malley Jr., J. Math. Mech. 18, 835 (1969). [8] R.E. O Malley Jr., Introduction to Singular Perturbations (Academic Press, New York, 1974). [9] R.E. O Malley Jr., Singular Perturbation Methods for Ordinary Differential Equations (Springer, New York, 1990). [10] E. O Riordan, M.L. Pickett and G. I. Shishkin, Comput. Methods Appl. Math. 3, 424 (2003). [11] E. O Riordan, M.L. Pickett and G.I. Shishkin, Math. Comp. 75, 1135 (2006). [12] G.I. Shishkin and V.A. Titov, Chisl. Metody Mekh. Sploshn. Sredy 2, 145 (1976). [13] R. Vulanović, Computing 67, 287 (2001). [14] C. De Boor, A Practical Guide to Splines (Springer-Verlag, New York, 1978) Figure 1: Solution profile for Example 1 for µ = 10 4 using N = 32.
8 1532 D. Kumar et al. : A Parameter-uniform Method for Two Parameters... D. Kumar is presently employed as an assistant professor in the department of Mathematics at BITS Pilani, India. He obtained his PhD from IIT Kanpur (INDIA). He is an active researcher coupled with the teaching experience in various Institutes in India. He has published more than 10 research articles in reputed international journals of mathematical and engineering sciences. A. S. Yadaw is presently a postdoc fellow in the department of Pharmacology & Systems Biology Mount Sinai School of Medicine, New York. He obtained his PhD from IIT Kanpur (INDIA). He is an active researcher and has published more than 6 research articles in reputed international journals of mathematical and engineering sciences. M. K. Kadalbajoo is presently employed as a professor in the department of Mathematics at IIT Kanpur, India. He obtained his PhD from IIT Bombay (INDIA). He is an active researcher coupled with the vast (40 years) teaching experience at IIT Bombay and IIT Kanpur in India. He has been an invited speaker of number of conferences and has published more than 200 (Two hundred) research articles in reputed international journals of mathematical and engineering sciences.
Finite Difference, Finite element and B-Spline Collocation Methods Applied to Two Parameter Singularly Perturbed Boundary Value Problems 1
European Society of Computational Methods in Sciences and Engineering (ESCMSE) Journal of Numerical Analysis, Industrial and Applied Mathematics (JNAIAM) vol. 5, no. 3-4, 2011, pp. 163-180 ISSN 1790 8140
More informationNUMERICAL SOLUTION OF GENERAL SINGULAR PERTURBATION BOUNDARY VALUE PROBLEMS BASED ON ADAPTIVE CUBIC SPLINE
TWMS Jour. Pure Appl. Math., V.3, N.1, 2012, pp.11-21 NUMERICAL SOLUTION OF GENERAL SINGULAR PERTURBATION BOUNDARY VALUE PROBLEMS BASED ON ADAPTIVE CUBIC SPLINE R. MOHAMMADI 1 Abstract. We use adaptive
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 informationQuintic B-spline method for numerical solution of fourth order singular perturbation boundary value problems
Stud. Univ. Babeş-Bolyai Math. 63(2018, No. 1, 141 151 DOI: 10.24193/subbmath.2018.1.09 Quintic B-spline method for numerical solution of fourth order singular perturbation boundary value problems Ram
More informationA COMPUTATIONAL TECHNIQUE FOR SOLVING BOUNDARY VALUE PROBLEM WITH TWO SMALL PARAMETERS
Electronic Journal of Differential Equations, Vol. 2013 (2013), No. 30, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu A COMPUTATIONAL
More informationHIGH ORDER VARIABLE MESH EXPONENTIAL FINITE DIFFERENCE METHOD FOR THE NUMERICAL SOLUTION OF THE TWO POINTS BOUNDARY VALUE PROBLEMS
JOURNAL OF THE INTERNATIONAL MATHEMATICAL VIRTUAL INSTITUTE ISSN (p) 2303-4866, ISSN (o) 2303-4947 www.imvibl.org /JOURNALS / JOURNAL Vol. 8(2018), 19-33 DOI: 10.7251/JMVI1801019P Former BULLETIN OF THE
More informationCubic Splines MATH 375. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Cubic Splines
Cubic Splines MATH 375 J. Robert Buchanan Department of Mathematics Fall 2006 Introduction Given data {(x 0, f(x 0 )), (x 1, f(x 1 )),...,(x n, f(x n ))} which we wish to interpolate using a polynomial...
More informationAvailable online at ScienceDirect. Procedia Engineering 127 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 127 (2015 ) 424 431 International Conference on Computational Heat and Mass Transfer-2015 Exponentially Fitted Initial Value
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 informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More information1 Piecewise Cubic Interpolation
Piecewise Cubic Interpolation Typically the problem with piecewise linear interpolation is the interpolant is not differentiable as the interpolation points (it has a kinks at every interpolation point)
More informationAn Exponentially Fitted Non Symmetric Finite Difference Method for Singular Perturbation Problems
An Exponentially Fitted Non Symmetric Finite Difference Method for Singular Perturbation Problems GBSL SOUJANYA National Institute of Technology Department of Mathematics Warangal INDIA gbslsoujanya@gmail.com
More informationNumerical integration of singularly perturbed delay differential equations using exponential integrating factor
MATHEMATICAL COMMUNICATIONS 51 Math. Commun. 017), 51 64 Numerical integration of singularly perturbed delay differential equations using exponential integrating factor Lakshmi Sirisha Challa and Yanala
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationNON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE
Global and Stocastic Analysis Vol. 4 No. 1, January 2017, 1-10 NON STANDARD FITTED FINITE DIFFERENCE METHOD FOR SINGULAR PERTURBATION PROBLEMS USING CUBIC SPLINE K. PHANEENDRA AND E. SIVA PRASAD Abstract.
More informationParameter robust methods for second-order complex-valued reaction-diffusion equations
Postgraduate Modelling Research Group, 10 November 2017 Parameter robust methods for second-order complex-valued reaction-diffusion equations Faiza Alssaedi Supervisor: Niall Madden School of Mathematics,
More informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
4 Interpolation 4.1 Polynomial interpolation Problem: LetP n (I), n ln, I := [a,b] lr, be the linear space of polynomials of degree n on I, P n (I) := { p n : I lr p n (x) = n i=0 a i x i, a i lr, 0 i
More informationCubic Splines. Antony Jameson. Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305
Cubic Splines Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305 1 References on splines 1. J. H. Ahlberg, E. N. Nilson, J. H. Walsh. Theory of
More informationFitted operator finite difference scheme for a class of singularly perturbed differential difference turning point problems exhibiting interior layer
Fitted operator finite difference scheme for a class of singularly perturbed differential difference turning point problems exhibiting interior layer Dr. Pratima Rai Department of mathematics, University
More informationG-SPLINE INTERPOLATION FOR APPROXIMATING THE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS USING LINEAR MULTI- STEP METHODS
Journal of Al-Nahrain University Vol.0(), December, 00, pp.8- Science G-SPLINE INTERPOLATION FOR APPROXIMATING THE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS USING LINEAR MULTI- STEP METHODS Osama H.
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 informationNUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS
TWMS J Pure Appl Math V5, N2, 24, pp22-228 NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS S ASADI, AH BORZABADI Abstract In this paper, Haar wavelet benefits are applied to the delay
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More information3.1 Interpolation and the Lagrange Polynomial
MATH 4073 Chapter 3 Interpolation and Polynomial Approximation Fall 2003 1 Consider a sample x x 0 x 1 x n y y 0 y 1 y n. Can we get a function out of discrete data above that gives a reasonable estimate
More informationOn solving linear systems arising from Shishkin mesh discretizations
On solving linear systems arising from Shishkin mesh discretizations Petr Tichý Faculty of Mathematics and Physics, Charles University joint work with Carlos Echeverría, Jörg Liesen, and Daniel Szyld October
More informationA coupled system of singularly perturbed semilinear reaction-diffusion equations
A coupled system of singularly perturbed semilinear reaction-diffusion equations J.L. Gracia, F.J. Lisbona, M. Madaune-Tort E. O Riordan September 9, 008 Abstract In this paper singularly perturbed semilinear
More informationSPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM. Mirjana Stojanović University of Novi Sad, Yugoslavia
GLASNIK MATEMATIČKI Vol. 37(57(2002, 393 403 SPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM Mirjana Stojanović University of Novi Sad, Yugoslavia Abstract. We introduce piecewise interpolating
More informationNew Class of Boundary Value Problems
Inf. Sci. Lett. 1 No. 2, 67-76 (2012) Information Science Letters An International Journal 67 @ 2012 NSP Natural Sciences Publishing Cor. New Class of Boundary Value Problems Abdon Atangana Institute for
More informationSPLINE INTERPOLATION
Spline Background SPLINE INTERPOLATION Problem: high degree interpolating polynomials often have extra oscillations. Example: Runge function f(x = 1 1+4x 2, x [ 1, 1]. 1 1/(1+4x 2 and P 8 (x and P 16 (x
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationLecture 10 Polynomial interpolation
Lecture 10 Polynomial interpolation Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn
More informationLecture 1 INF-MAT3350/ : Some Tridiagonal Matrix Problems
Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems Tom Lyche University of Oslo Norway Lecture 1 INF-MAT3350/4350 2007: Some Tridiagonal Matrix Problems p.1/33 Plan for the day 1. Notation
More informationNumerische Mathematik
Numer. Math. (997) 76: 479 488 Numerische Mathematik c Springer-Verlag 997 Electronic Edition Exponential decay of C cubic splines vanishing at two symmetric points in each knot interval Sang Dong Kim,,
More informationCHAPTER 10 Shape Preserving Properties of B-splines
CHAPTER 10 Shape Preserving Properties of B-splines In earlier chapters we have seen a number of examples of the close relationship between a spline function and its B-spline coefficients This is especially
More informationMaria Cameron Theoretical foundations. Let. be a partition of the interval [a, b].
Maria Cameron 1 Interpolation by spline functions Spline functions yield smooth interpolation curves that are less likely to exhibit the large oscillations characteristic for high degree polynomials Splines
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 informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationNUMERICAL STUDY OF CONVECTION DIFFUSION PROBLEM IN TWO- DIMENSIONAL SPACE
IJRRAS 5 () November NUMERICAL STUDY OF CONVECTION DIFFUSION PROBLEM IN TWO- DIMENSIONAL SPACE. K. Sharath Babu,* & N. Srinivasacharyulu Research Scholar, Department of Mathematics, National Institute
More informationChordal cubic spline interpolation is fourth order accurate
Chordal cubic spline interpolation is fourth order accurate Michael S. Floater Abstract: It is well known that complete cubic spline interpolation of functions with four continuous derivatives is fourth
More informationMultiple-Scales Method and Numerical Simulation of Singularly Perturbed Boundary Layer Problems
Appl. Math. Inf. Sci., No. 3, 9-7 (6) 9 Applied Mathematics & Infmation Sciences An International Journal http://dx.doi.g/.8576/amis/33 Multiple-Scales Method and Numerical Simulation of Singularly Perturbed
More informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
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 informationNumerical Solution of One-dimensional Telegraph Equation using Cubic B-spline Collocation Method
2014 (2014) 1-8 Available online at www.ispacs.com/iasc Volume 2014, Year 2014 Article ID iasc-00042, 8 Pages doi:10.5899/2014/iasc-00042 Research Article Numerical Solution of One-dimensional Telegraph
More informationA uniformly convergent difference scheme on a modified Shishkin mesh for the singularly perturbed reaction-diffusion boundary value problem
Journal of Modern Methods in Numerical Mathematics 6:1 (2015, 28 43 A uniformly convergent difference scheme on a modified Shishkin mesh for the singularly perturbed reaction-diffusion boundary value problem
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 information290 J.M. Carnicer, J.M. Pe~na basis (u 1 ; : : : ; u n ) consisting of minimally supported elements, yet also has a basis (v 1 ; : : : ; v n ) which f
Numer. Math. 67: 289{301 (1994) Numerische Mathematik c Springer-Verlag 1994 Electronic Edition Least supported bases and local linear independence J.M. Carnicer, J.M. Pe~na? Departamento de Matematica
More informationIterative Domain Decomposition Methods for Singularly Perturbed Nonlinear Convection-Diffusion Equations
Iterative Domain Decomposition Methods for Singularly Perturbed Nonlinear Convection-Diffusion Equations P.A. Farrell 1, P.W. Hemker 2, G.I. Shishkin 3 and L.P. Shishkina 3 1 Department of Computer Science,
More informationA NOTE ON BLASIUS TYPE BOUNDARY VALUE PROBLEMS. Grzegorz Andrzejczak, Magdalena Nockowska-Rosiak, and Bogdan Przeradzki
Opuscula Math. 33, no. 1 213, 5 17 http://dx.doi.org/1.7494/opmath.213.33.1.5 Opuscula Mathematica A NOTE ON BLASIUS TYPE BOUNDARY VALUE PROBLEMS Grzegorz Andrzejczak, Magdalena Nockowska-Rosiak, and Bogdan
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationExact Solutions for a Class of Singular Two-Point Boundary Value Problems Using Adomian Decomposition Method
Applied Mathematical Sciences, Vol 6, 212, no 122, 697-618 Exact Solutions for a Class of Singular Two-Point Boundary Value Problems Using Adomian Decomposition Method Abdelhalim Ebaid 1 and Mona D Aljoufi
More informationA Singularly Perturbed Convection Diffusion Turning Point Problem with an Interior Layer
A Singularly Perturbed Convection Diffusion Turning Point Problem with an Interior Layer E. O Riordan, J. Quinn School of Mathematical Sciences, Dublin City University, Glasnevin, Dublin 9, Ireland. Abstract
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationUNIFORMLY CONVERGENT 3-TGFEM VS LSFEM FOR SINGULARLY PERTURBED CONVECTION-DIFFUSION PROBLEMS ON A SHISHKIN BASED LOGARITHMIC MESH
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 4, Number 4, Pages 35 334 c 23 Institute for Scientific Computing and Information UNIFORMLY CONVERGENT 3-TGFEM VS LSFEM FOR SINGULARLY
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationA PARAMETER ROBUST NUMERICAL METHOD FOR A TWO DIMENSIONAL REACTION-DIFFUSION PROBLEM
A PARAMETER ROBUST NUMERICAL METHOD FOR A TWO DIMENSIONAL REACTION-DIFFUSION PROBLEM C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN Abstract. In this paper a singularly perturbed reaction diffusion partial
More informationRobust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations
Robust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with
More informationOn a max norm bound for the least squares spline approximant. Carl de Boor University of Wisconsin-Madison, MRC, Madison, USA. 0.
in Approximation and Function Spaces Z Ciesielski (ed) North Holland (Amsterdam), 1981, pp 163 175 On a max norm bound for the least squares spline approximant Carl de Boor University of Wisconsin-Madison,
More informationx n+1 = x n f(x n) f (x n ), n 0.
1. Nonlinear Equations Given scalar equation, f(x) = 0, (a) Describe I) Newtons Method, II) Secant Method for approximating the solution. (b) State sufficient conditions for Newton and Secant to converge.
More informationA NUMERICAL APPROACH FOR SOLVING A CLASS OF SINGULAR BOUNDARY VALUE PROBLEMS ARISING IN PHYSIOLOGY
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 8, Number 2, Pages 353 363 c 2011 Institute for Scientific Computing and Information A NUMERICAL APPROACH FOR SOLVING A CLASS OF SINGULAR
More informationQUINTIC SPLINE SOLUTIONS OF FOURTH ORDER BOUNDARY-VALUE PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 5, Number, Pages 0 c 2008 Institute for Scientific Computing Information QUINTIC SPLINE SOLUTIONS OF FOURTH ORDER BOUNDARY-VALUE PROBLEMS
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 informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationNumerical Integration for Multivariable. October Abstract. We consider the numerical integration of functions with point singularities over
Numerical Integration for Multivariable Functions with Point Singularities Yaun Yang and Kendall E. Atkinson y October 199 Abstract We consider the numerical integration of functions with point singularities
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationA Nodal Spline Collocation Method for the Solution of Cauchy Singular Integral Equations 1
European Society of Computational Methods in Sciences and Engineering (ESCMSE) Journal of Numerical Analysis, Industrial and Applied Mathematics (JNAIAM) vol. 3, no. 3-4, 2008, pp. 211-220 ISSN 1790 8140
More informationAn interface problem with singular perturbation on a subinterval
Xie Boundary Value Problems 2014, 2014:201 R E S E A R C H Open Access An interface problem with singular perturbation on a subinterval Feng Xie * * Correspondence: fxie@dhu.edu.cn Department of Applied
More informationACCURATE NUMERICAL METHOD FOR BLASIUS PROBLEM FOR FLOW PAST A FLAT PLATE WITH MASS TRANSFER
FDS-2000 Conference, pp. 1 10 M. Sapagovas et al. (Eds) c 2000 MII ACCURATE UMERICAL METHOD FOR BLASIUS PROBLEM FOR FLOW PAST A FLAT PLATE WITH MASS TRASFER B. GAHA 1, J. J. H. MILLER 2, and G. I. SHISHKI
More informationCharacterizations of Regular Semigroups
Appl. Math. Inf. Sci. 8, No. 2, 715-719 (2014) 715 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.12785/amis/080230 Characterizations of Regular Semigroups Bingxue
More informationIntro Polynomial Piecewise Cubic Spline Software Summary. Interpolation. Sanzheng Qiao. Department of Computing and Software McMaster University
Interpolation Sanzheng Qiao Department of Computing and Software McMaster University July, 2012 Outline 1 Introduction 2 Polynomial Interpolation 3 Piecewise Polynomial Interpolation 4 Natural Cubic Spline
More informationCubic B-spline Collocation Method for One-Dimensional Heat Equation
Pure and Applied Mathematics Journal 27; 6(): 5-58 http://www.sciencepublishinggroup.com/j/pamj doi:.648/j.pamj.276.7 ISSN: 2326-979 (Print); ISSN: 2326-982 (Online) Cubic B-spline Collocation Method for
More informationTotal least squares. Gérard MEURANT. October, 2008
Total least squares Gérard MEURANT October, 2008 1 Introduction to total least squares 2 Approximation of the TLS secular equation 3 Numerical experiments Introduction to total least squares In least squares
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 informationAn Asymptotic Semi-Analytical Method for Third- Order Singularly Perturbed Boundary Value Problems
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 12, Number 1 (2016), pp. 253-263 Research India Publications http://www.ripublication.com An Asymptotic Semi-Analytical Method for
More informationA complex geometric proof of Tian-Yau-Zelditch expansion
A complex geometric proof of Tian-Yau-Zelditch expansion Zhiqin Lu Department of Mathematics, UC Irvine, Irvine CA 92697 October 21, 2010 Zhiqin Lu, Dept. Math, UCI A complex geometric proof of TYZ expansion
More informationSolution of a Fourth Order Singularly Perturbed Boundary Value Problem Using Quintic Spline
International Mathematical Forum, Vol. 7, 202, no. 44, 279-290 Solution of a Fourth Order Singularly Perturbed Boundary Value Problem Using Quintic Spline Ghazala Akram and Nadia Amin Department of Mathematics
More information1. Introduction We consider the following self-adjoint singularly perturbed boundary value problem: ɛu + p(x)u = q(x), p(x) > 0,
TWMS J. Pure Appl. Math. V., N., 00, pp. 6-5 NON-POLYNOMIAL SPLINE APPROXIMATIONS FOR THE SOLUTION OF SINGULARLY-PERTURBED BOUNDARY VALUE PROBLEMS J. RASHIDINIA, R. MOHAMMADI Abstract. We consider the
More informationParameter-Uniform Numerical Methods for a Class of Singularly Perturbed Problems with a Neumann Boundary Condition
Parameter-Uniform Numerical Methods for a Class of Singularly Perturbed Problems with a Neumann Boundary Condition P. A. Farrell 1,A.F.Hegarty 2, J. J. H. Miller 3, E. O Riordan 4,andG.I. Shishkin 5 1
More informationCubic Spline. s(x) = s j (x) = a j + b j (x x j ) + c j (x x j ) 2 + d j (x x j ) 3, j = 0,... n 1 (1)
Cubic Spline Suppose we are given a set of interpolating points (x i, y i ) for i = 0, 1, 2, n We seek to construct a piecewise cubic function s(x) that has the for x [[x j, x j+1 ] we have: s(x) = s j
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationNumerical Solutions of the Combined KdV-MKdV Equation by a Quintic B-spline Collocation Method
Appl. Math. Inf. Sci. Lett. 4 No. 1 19-24 (2016) 19 Applied Mathematics & Information Sciences Letters An International Journal http://d.doi.org/18576/amisl/040104 Numerical Solutions of the Combined KdV-MKdV
More informationThe Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods
The Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods by Hae-Soo Oh Department of Mathematics, University of North Carolina at Charlotte, Charlotte, NC 28223 June
More informationOn A Comparison of Numerical Solution Methods for General Transport Equation on Cylindrical Coordinates
Appl. Math. Inf. Sci. 11 No. 2 433-439 (2017) 433 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.18576/amis/110211 On A Comparison of Numerical Solution Methods
More informationLecture 1 INF-MAT : Chapter 2. Examples of Linear Systems
Lecture 1 INF-MAT 4350 2010: Chapter 2. Examples of Linear Systems Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo August 26, 2010 Notation The set of natural
More informationHigh Accuracy Finite Difference Approximation to Solutions of Elliptic Partial Differential Equations
Purdue University Purdue e-pubs Department of Computer Science Technical Reports Department of Computer Science 1977 High Accuracy Finite Difference Approximation to Solutions of Elliptic Partial Differential
More informationInterior Layers in Singularly Perturbed Problems
Interior Layers in Singularly Perturbed Problems Eugene O Riordan Abstract To construct layer adapted meshes for a class of singularly perturbed problems, whose solutions contain boundary layers, it is
More informationThere is a unique function s(x) that has the required properties. It turns out to also satisfy
Numerical Analysis Grinshpan Natural Cubic Spline Let,, n be given nodes (strictly increasing) and let y,, y n be given values (arbitrary) Our goal is to produce a function s() with the following properties:
More informationIntro Polynomial Piecewise Cubic Spline Software Summary. Interpolation. Sanzheng Qiao. Department of Computing and Software McMaster University
Interpolation Sanzheng Qiao Department of Computing and Software McMaster University January, 2014 Outline 1 Introduction 2 Polynomial Interpolation 3 Piecewise Polynomial Interpolation 4 Natural Cubic
More informationMath 246, Spring 2011, Professor David Levermore
Math 246, Spring 2011, Professor David Levermore 7. Exact Differential Forms and Integrating Factors Let us ask the following question. Given a first-order ordinary equation in the form 7.1 = fx, y, dx
More informationSome notes on Chapter 8: Polynomial and Piecewise-polynomial Interpolation
Some notes on Chapter 8: Polynomial and Piecewise-polynomial Interpolation See your notes. 1. Lagrange Interpolation (8.2) 1 2. Newton Interpolation (8.3) different form of the same polynomial as Lagrange
More informationSolving Linear Sixth-Order Boundary Value Problems by Using Hyperbolic Uniform Spline Method
International Journal of Mathematical Modelling & Computations Vol. 03, No. 03, 013, 169-180 Solving Linear Sixth-Order Boundary Value Problems by Using Hyperbolic Uniform Spline Method J. Dabounou a,
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 informationFinite Element Methods
Solving Operator Equations Via Minimization We start with several definitions. Definition. Let V be an inner product space. A linear operator L: D V V is said to be positive definite if v, Lv > for every
More informationNumerical Methods for Two Point Boundary Value Problems
Numerical Methods for Two Point Boundary Value Problems Graeme Fairweather and Ian Gladwell 1 Finite Difference Methods 1.1 Introduction Consider the second order linear two point boundary value problem
More informationConservation laws for the geodesic equations of the canonical connection on Lie groups in dimensions two and three
Appl Math Inf Sci 7 No 1 311-318 (013) 311 Applied Mathematics & Information Sciences An International Journal Conservation laws for the geodesic equations of the canonical connection on Lie groups in
More informationCAAM 454/554: Stationary Iterative Methods
CAAM 454/554: Stationary Iterative Methods Yin Zhang (draft) CAAM, Rice University, Houston, TX 77005 2007, Revised 2010 Abstract Stationary iterative methods for solving systems of linear equations are
More informationLinear Algebra Practice Problems
Linear Algebra Practice Problems Math 24 Calculus III Summer 25, Session II. Determine whether the given set is a vector space. If not, give at least one axiom that is not satisfied. Unless otherwise stated,
More informationInterpolation. Chapter Interpolation. 7.2 Existence, Uniqueness and conditioning
76 Chapter 7 Interpolation 7.1 Interpolation Definition 7.1.1. Interpolation of a given function f defined on an interval [a,b] by a polynomial p: Given a set of specified points {(t i,y i } n with {t
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationCHAPTER 4. Interpolation
CHAPTER 4 Interpolation 4.1. Introduction We will cover sections 4.1 through 4.12 in the book. Read section 4.1 in the book on your own. The basic problem of one-dimensional interpolation is this: Given
More informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More information