1. Introduction We consider the following self-adjoint singularly perturbed boundary value problem: ɛu + p(x)u = q(x), p(x) > 0,
|
|
- Jeffry Ball
- 6 years ago
- Views:
Transcription
1 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 self-adjoint singularly perturbed two-point boundary value problems. We know that the numerical methods for solution of such problems based on nonpolynomial spline in grid points, can produce the fourth order method only. So that we look for an alternative to obtain higher order methods. We develop the non-polynomial spline in off-step points to rise the order of accuracy. Based on such spline, the purposed new methods are fourth, sixth and eighth-order accurate. These methods are applicable to problems both in singular and non-singular cases. The convergence analysis of the new eight-order method is proved. We applied the presented methods to test problems which have been solved by other existing methods in our references, for comparison of our methods with the existing methods. Numerical results are given to illustrate the efficiency of our methods. Keywords: Self-adjoint singularly-perturbed boundary value problem, Non-polynomial spline, Convergence analysis, Numerical conclusion. AMS Subject Classification: 4B5; F05; 65D0; 65L.. Introduction We consider the following self-adjoint singularly perturbed boundary value problem: ɛu + p(x)u = q(x), p(x) > 0, u(a) = η, u(b) = η, () where η, η are given constants and ɛ is a small positive parameter such that 0 < ɛ and p(x), q(x) are small bounded real functions. This problems occur naturally in various fields of science and engineering, for example, combustion, nuclear engineering, control theory, elasticity, fluid mechanics, fluid dynamics, quantum mechanics, optimal control, chemical-reactor theory, hydrodynamics, convection-diffusion process, geophysics, etc. A few notable examples are boundary-layer problems, WKB Theory, the modeling of steady and unsteady viscous flow problems with large Reynolds number and convective heat transport problems with large Peclet number. Secondly, the occurrence of sharp boundary-layers as ɛ, the coefficient of highest derivative, approaches zero creates difficulty for most standard numerical schemes. The application of splines for the numerical solution of singularly-perturbed boundary- value problems has been described in many papers []-[5], [7],[8],[]-[6]. Recently Aziz and khan [, ] and Khan et. al. [6] solve this problem by cubic spline in comparison, quintic spline and sextic spline with assuming p(x) = constant. A generalized scheme based on quartic non-polynomial spline functions was proposed by Tirmizi et. al. [6]. In this paper non-polynomial sextic spline relations have been derived using off-step points. In Section, we derive the formulation of our spline function approximation. We present the School of Mathematics, Iran University of Science and Technology, Tehran, Iran, rashidinia@iust.ac.ir Department of Mathematics, University of Neyshabour, Neyshabour, Iran, rez.mohammadi@gmail.com Manuscript received 8 November
2 J. RASHIDINIA, R. MOHAMMADI: NON-POLYNOMIAL SPLINE APPROXIMATIONS FOR... 7 spline relations to be used for discretization of the given problem (). In section, we present the formulation of our method. We drive boundary formulas in section 4. In section 5, convergence of the new eight-order method is stablished. Finally, in section 6, numerical evidence is included to demonstrate the efficiency of the method.. Non-polynomial spline functions We introduce the set of grid points in the interval [a, b] x 0 = a, x i = a + (i )h, h = b a N, i =,,..., N, x N = b. Definition. A non-polynomial spline function S l (x), interpolating to a function u(x) on [a, b] defined as: () In each subinterval [x l, x l+ ], S l (x) is a polynomial of degree at most six. () The first-fifth derivatives of S l (x) are continuous on [a, b]. () S l (x ) = u(x ), l = ()N. For each segment [x, x l+ ], l =,,..., N the non-polynomial spline S l (x), is define as S l (x) = a l + b l (x x l ) + c l (x x l ) + d l (x x l ) + e l (x x l ) 4 + +f l sin τ(x x l ) + g l cos τ(x x l ), l = 0,,,..., N, () where a l, b l, c l, d l, e l, f l and g l are constants and τ is free parameter. We further require that the values of the first-, second-, third-, fourth- and fifth-order derivatives be the same for the pair of segments that join at each point (x l, u l ). To derive expression for the coefficients of () in terms of u we first define: F and F l+ (i) S l (x ) = u,, (ii) S l (x l+ ) = u l+ (iii) S l (x ) = m, (iv) S l (x ) = M, (v) S l (x l+ ) = M l+, (vi) S (4) l (x ) = F, (vii) S (4) l (x l+ ) = F l+. From algebraic manipulation we get the following expression: a l =, u l+, m, M, M l+, 96τ 4 θ[ cos θ + θ( + θ ) sin θ] {[7θ(8 + θ ) + θ(9 + 4θ + 5θ 4 ) cos θ+ +( 5 48θ 7θ 4 + θ 6 ) sin θ]f [6θ(9 + 4θ + 5θ 4 ) + 44θ(8 + θ ) cos θ+ +( 04 96θ + 4θ 4 + θ 6 ) sin θ]f l+ + τ 4 {6h sin θ( θ + 5θ 4 )m + 4h [6θ( cos θ) + (9 + 6θ + θ 4 ) sin θ]m + h [44θ( cos θ) + (84 4θ θ 4 ) sin θ]m l+ + +[ 5( cos θ) + (04 + 8θ + 78θ 4 ) sin θ]u + [ 5( cos θ) + ( θ + +8θ 4 ) sin θ]u l+ }}, ()
3 8 TWMS J. PURE APPL. MATH., V., N., 00 b l = 4τ θ {(4 + θ )(F F l+ ) + τ 4 [h (M M l+ ) + 4(u l+ u )]}, c l = 4τ θ[ cos θ + θ( + θ ) sin θ] {[4θ + 6θ(4 + θ ) cos θ + ( 48 θ + θ 4 ) sin θ]f [6θ(4 + θ ) + 4θ cos θ (4 + θ ) sin θ]f l+ + τ {6θ sin θ (8 + θ )m + h[4( cos θ) + θ(4 + θ ) sin θ]m + h[4( cos θ) 4θ sin θ]m l+ + τ[(48 + 6θ ) sin θ](u d l = 6τh [F l+ F + τ (M l+ e l = M )], u l+ )}}, 6θ[ 4( cos θ) + θ( + θ ) sin θ] {[θ cos θ + ( + θ ) sin θ]f + [ 6θ+ and +(6 + θ ) sin θ]f l+ + τ 4 sin θ[6hm + h (M f l = τ 4 sin θ (F i+ F ), + M l+ ) + 6(u u l+ )]}, g l = τ 4 θ[ 4( cos θ) + θ( + θ ) sin θ] {[ 4θ cos θ + (48 4θ + θ 4 ) sin θ ](F F l+ ) 8τ 4 sin θ [6hm + h (M where θ = τh and l =,,..., N. + M l+ ) + 6(u u l+ )]}, (4) From continuity of the first derivative at the point x = x l, we can obtain the following spline relation: m = 4τ θ sin θ( 4θ + θ + 48 sin θ ){ [9 + 8θ + 8(4 + 4θ + θ 4 ) cos θ 8(4 θ + θ 4 ) cos θ ( θ ) cos θ θ(96 0θ + θ 4 ) sin θ 9θ sin θ ]F + [ θ θ 4 + (4 + 6θ + θ 4 ) cos θ + 4( 4 + θ ) cos θ (4 8θ θ 4 ) cos θ 88θ sin θ θ(576 7θ θ 4 ) sin θ 88θ sin θ ]F + [ 576 4θ 4θ ( + θ ) cos θ + 4(4 + θ ) cos θ + 576θ sin θ + θ(88 + 6θ + θ 4 ) sin θ]f l+ ) + 6τ θ sin θ [θ ( 0 + θ ) cos θ 8(sin θ 8 sin θ)]m + 4τ θ sin θ [θ(84 + θ ) cos θ sin θ sin θ]m + τ θ [ 4( cos θ) + θ( + θ ) sin θ]m l τ 4 sin θ [ θ cos θ +
4 J. RASHIDINIA, R. MOHAMMADI: NON-POLYNOMIAL SPLINE APPROXIMATIONS FOR (sin θ + sin θ)]u + 48τ 4 sin θ [θ(48 + θ ) cos θ 4( sin θ + sin θ)]u 4τ 4 [ 4( cos θ) + θ( + θ ) sin θ]u l+ }. (5) Also continuity of the third derivative at the point x = x l, implies that: m = 4τ θ sin θ( θ + sin θ ){ [4 + (4 + 4θ + θ 4 ) cos θ + 4( + θ ) cos θ 8( θ ) cos θ θ(6 7θ ) sin θ 4θ sin θ ]F + [48 + 4θ + (4+ +6θ + θ 4 ) cos θ 48 cos θ (4 8θ θ 4 ) cos θ 4θ sin θ 6θ(8 + θ ) sin θ 4θ sin θ ]F + [ 4 θ ( + θ ) cos θ + 4 cos θ + 4θ sin θ + θ ( + θ ) sin θ]f l+ + τ sin θ [θ( 7θ ) cos θ 4 sin θ + 8θ sin θ]m + +4τ sin θ [ θ(4 + θ ) cos θ + 4 sin θ + θ sin θ]m + τ [ 4( cos θ) + θ( + θ ) sin θ]m l+ + 4τ 4 ( θ + sin θ ) sin θ(u u )}. (6) Finally, from continuity of the fifth derivative at the point x = x l, we have the following spline relation: m = 96τ θ cos θ sin θ {{(4 + 4θ + θ 4 ) cos θ + 8[( + θ ) cos θ θ sin θ ]} F [(4 + 6θ + θ 4 ) cos θ + ( 4 + 8θ + θ 4 ) cos θ 96θ θ cos sin θ ]F + +θ[θ( + θ ) cos θ 4 sin θ ]F l+ 6τ 4 cos θ sin θ [h (M + M )+ where l = ()N. +6(u u )]}, (7) From Eqs. (5) and (6) we have u u + u l+ = +44θ + (96 + 4θ θ 4 ) cos θ]f Also from Eqs. (5) and (7) we have 48τ 4 ( cos θ) {[48 θ + θ 4 (4 + θ ) cos θ]f + [ 96+ +τ θ ( cos θ)(m M θ {[ + θ + cos θ]f [ + ( + θ ) cos θ]f From Eqs. (8) and (9) we obtain τ ( cos θ)(m M + [48 θ + θ 4 (4 + θ ) cos θ]f l+ + + M l+ )}. (8) + [ + θ + cos θ]f l+ + M l+ )} = 0. (9)
5 40 TWMS J. PURE APPL. MATH., V., N., 00 h 4 F = [ + 5θ + ( + θ ) cos θ] {h {[4 θ + θ 4 4 cos θ](m + M l+ )+ +[ 4 + 5θ 4 + θ 4 + ( + θ ) cos θ]m + θ [ θ cos θ](u u + u l+ )}. (0) The elimination of F l s from (9) and (0) gives where h [u i 5 + β (u i + u l+ ) + γ u i +M i+ ) + γm i α = 4( cos θ) + θ ( + θ ) θ ( + θ, + cos θ) β = 48 + θ ( + 5θ ) + (48 + θ + θ 4 ) cos θ 6θ ( + θ, + cos θ) γ = 7 + θ ( + θ ) (6 + θ + 5θ 4 ) cos θ 6θ ( + θ, + cos θ) β = [ 4 + θ + (4 + θ ) cos θ] + θ, + cos θ + u i+ ] = α(m i 5 + M i+ ) + β(m i +, i = ()N, () γ = [ 6 + θ + ( + θ ) cos θ] + θ. () + cos θ When τ 0, that θ 0, then (α, β, γ) 0 (, 56, 46), (β, γ ) (8, 8) and the relations defined by () and () reduce into sextic polynomial spline functions.. Description of methods Now we consider the equation (). At the grid points x i, the proposed differential equation () may be discretized by ɛu i + p i u i = q i, () where p i = p(x i ) and q i = q(x i ). By using moment of spline in () we obtain ɛm i + p i u i = q i. (4) We discretize the given system (4) at the grid points x i, i =,,..., N and using the spline relation (). We obtain (N 4) linear algebraic equation in the (N) unknowns u i, i =,,..., N as ( ɛ + αh p i 5 )u i 5 +( ɛβ + βh p i+ )u i+ + ( ɛβ + βh p i )u i + ( ɛ + αh p i+ )u i+ +βq i+ + αq i+ where p i = p(x i ) and q i = q(x i ), with local truncation error + ( ɛγ + γp i )h 4 u i + + βq i + γq i + = h (αq i 5 ), i = ()N, (5) T i (h) = ɛ( + β + γ )u i + ɛ ( + β + γ )hu i + ɛ 8 (6α + 6β + 8γ 4 0β γ )h u i ɛ 48 (48α + 48β + 4γ 98 6β γ )h u i + ɛ (6α + 480β γ 706 8β γ )h 4 u (4) i ɛ 840 (7840α + 080β + 80γ 88 4β
6 J. RASHIDINIA, R. MOHAMMADI: NON-POLYNOMIAL SPLINE APPROXIMATIONS FOR... 4 γ )h 5 u (5) i + ɛ (8470α β + 0γ β γ )h 6 u (6) i + + ɛ 50 (48476α β + 68γ β γ )h 7 u (7) i (6696α + 650β + 4γ β γ )h 8 u (8) i (87044α β + 88γ β γ )h 9 u (9) i + ɛ 090 ɛ ɛ ( α + 60β + 60γ β γ )h 0 u (0) i + O(h ), where u (k) i = u (k) (ξ i ) x i < ξ i < x i+. (6) By choosing suitable values of parameters β, γ, α, β and γ we obtain the following methods: () If we choose + β + γ = 0 and α + β + γ 4 β = 0 in (5) we obtain second-order schemes for the solution of (). () If we choose (β, γ, α, β, γ) = (4, 0, with truncation error T i = ɛ 80 h6 u (6) (ξ i ). 8, 9 () If we choose (β, γ, α, β, γ) = (4, 0, 0, with truncation error T i = ɛ 50 h8 u (8) (ξ i )., 5) in (5) we obtain a fourth-order scheme 5, 49 0 ) in (5) we obtain a sixth-order scheme (4) If we choose (β, γ, α, β, γ) = ( 8, 8, 465, , ) in (5) we obtain a eight-order scheme with truncation error T i = 79ɛ h0 u (0) (ξ i ), which is the the highest order among the methods for solution of (). 4. Development of boundary formula The relation (5) gives N 4 linear algebraic equations in the n unknowns u i, for i =,,..., n. We need four more equations, two at each end of interval, for the direct computation of u i. By using Taylor series and method of undetermined coefficients the boundary formulas associated with boundary conditions can be determined as follows. We drive the following equations to be associated with fourth-order method as: 0u 0 + 4u u u 5 u 0 u + 0u 4u 5 u 7 = h 48 [M 0 54M 75M M 5 ] h6 u (6) (x 0 )+ +O(h 7 ), i =, (7) = h 48 [0M 0 75M 0M 76M 5 M 7 ] h6 u (6) (x 0 ) + O(h 7 ), i =, (8) and for right boundary we obtain: u N u N + 0u N 4u N 5 u N 7 76M N 5 = h 48 [0M N 75M N 0M N M N 7 ] h6 u (6) (x N ) + O(h 7 ), i = N, (9)
7 4 TWMS J. PURE APPL. MATH., V., N., 00 0u N + 4u N u N u N 5 = h 48 [M N 54M N 75M N M N 5 ] h6 u (6) (x N ) + O(h 7 ), i = N. (0) In same manner for sixth-order method we derive the following boundary formula: 0u 0 + 4u u u 5 u 0 u + 0u 4u 5 u 7 89M 7 = h 440 [ 608M M 4059M + 009M 5 869M M 9 ] h8 u (8) (x 0 ) + O(h 9 ), i =, () = h 0960 [0784M M 6759M 706M M 9 ] h8 u (8) (x 0 ) + O(h 9 ), i =, () u N u N + 0u N 4u N 5 u N M N 706M N 5 89M N [0784M N 57059M N + 406M N 9 ] h8 u (8) (x N )+ +O(h 9 ), i = N, () = h 0u N + 4u N u N u N M N 5 869M N 7 = h 440 [ 608M N 4049M N + 6M N 9 ] h8 u (8) (x N )+ 4059M N + +O(h 9 ), i = N. (4) Finally for eight order method we develop the following boundary equation: 0u 0 + 4u u u 5 = h [ M M M M M M M M ] h0 u (0) (x 0 ) + O(h ), i =, (5) u 0 u + 0u 4u 5 u 7 = h [ M M M M M M M M ] h0 u (0) (x 0 ) + O(h ), i =, (6) u N u N + 0u N 4u N 5 u N 7 = h [ M N M N M N M N M N M N M N M N ]+
8 J. RASHIDINIA, R. MOHAMMADI: NON-POLYNOMIAL SPLINE APPROXIMATIONS FOR h0 u (0) (x N ) + O(h ), i = N, (7) 0u N + 4u N u N u N 5 = h [ M N M N M N M N M N M N M N M N ] h0 u (0) (x N ) + O(h ), i = N. (8) 5. Convergence analysis In this section, we investigate the convergence analysis of our eight order method. The eight order quintic spline solution of () is based on the linear equations given by (5) and boundary formulas (5)-(8). Let Z = (z i ) be N dimensional column vector and A = (a ij ) be an (N N) matrix such that z = (0ɛ h p 0 )η + h [ q q q q q q q q ], z = (ɛ h p 0 )η + h [ q q q q q q q q z i = h 465 (q i q i + 58q i + 688q i+ ], + q i+ ), i = ()N, z N = (ɛ h p N )η + h [ q N q N q N q N q N q N q N q N ], z N = (0ɛ h p N )η + h [ q N q N q N and q N q N q N q N a i,i±k = coefficient of u i k, k = ±, ±, q N ], It is easily seen that the system of N linear algebraic equations given by (5) and boundary formulas (5)-(8) may be written in the form (i) AU = Z + T, (ii) AU = Z. (9)
9 44 TWMS J. PURE APPL. MATH., V., N., 00 The vector T is the local truncation error vector and define as ɛh0 u (0) (x 0 ) + O(h ), i =, t i = ɛh0 u (0) (x 0 ) + O(h ), i =, ɛh0 u (0) (x i ) + O(h ), x i 5 < x i < x i+ i =,..., N, ɛh0 u (0) (x N ) + O(h ), i = N, ɛh0 u (0) (x N ) + O(h ), i = N. (0) The error equation is obtained in the usal manner in the form AE = T, () where E = (e i ), with e i = (u i u i ) and T = (t i ) are the N dimentional vectors. Now the matrix A may be written as where A 0 = A = ɛ(a 0 + A ) + h BP, () , () A = , (4)
10 J. RASHIDINIA, R. MOHAMMADI: NON-POLYNOMIAL SPLINE APPROXIMATIONS FOR b b b b 4 b 5 b 6 b b b b b 4 b 5 b 6 b B = b N b N b N b N 4 b N 5 b N 6 b N b N b N b N b N4 b N5 b N6 b N7 with (b, b, b, b 4, b 5, b 6, b 6 ) = (b N7, b N6, b N5, b N4, b N, b N, b N ) = ( , , , , , , ), (b, b, b, b 4, b 5, b 6, b 7 ) = (b N 7, b N 6, b N 5, b N 4, b N, b N, b N ) = ( , , , , , , ) and P = diag(p i ). (5) Theorem. The symmetric matrix A 0 os irreducible and monotone and A 0 84h 4 [5(b a)4 + 0(b a) h + 9h 4 ]. (6) This theorem proved by Usmani [7]. Theorem. Let u be the solution of singularly perturbed boundary value problem (), and u i be the numerical solution obtained from the difference scheme (9)(ii). Then, we have ɛ E = O(h 8 ) provided p < ( ω) where ω = [5(b a) 4 + 0(b a) h + 9h 4 ] and p = 5ωh 84h 4 Max p(x i ), a p(x i ) b. Proof. From () we have E = A T A T [ɛ(a 0 + A ) + h BP ] T ɛ I + A 0 A + ɛ h A 0 BP ] A 0 T. From (), (4) and (5) we have A + ɛ h A 0 BP ( + ɛ h p ). (7) Also from (0) we have T ɛh0 L 0, where L 0 = Max u (0) (x i ), a x i b.
11 46 TWMS J. PURE APPL. MATH., V., N., 00 By using (I + W ) W E ɛ A 0 T A 0 A +ɛ h A 0 BP ] ɛ A 0 T A 0 A +ɛ h A 0 BP, provided that A 0 A + ɛ h A 0 BP <. Therefore we obtain E provided W <, we get ωh 0 L [ ω( + ɛ h p )] O(h8 ), (8) provided ɛ p < ( ω). (9) 5ωh So that the given method for solving the boundary value problem () for sufficiently small h is a eight-order convergent method. 6. Numerical illustrations In order to test the viability of the proposed methods and to demonstrate their convergence computationally we have solved the following singularly perturbed boundary value problems which their exact solutions are known to us. The maximum absolute errors at the nodal points, max u(x i ) u i are tabulated in Tables. Example.Consider the following singularly perturbed boundary value problem [],[5],[6], []-[7]: with the exact solution ɛu + u = cos πx ɛ cos πx, u(0) = u() = 0 u(x) = [exp((x )/ ɛ) + exp( x/ ɛ)]/[ + exp( / ɛ)] cos πx. This problem has been solved using our methods with different values of N = 6,,, 8, 56 and ɛ = 6,..., 8, computed solutions are compared with the exact solutions in nodal points and the maximum absolute errors in solutions are tabulated in Table, and compared with the methods in Ref. [],[5],[6], []-[7] in Tables -. Our results shows the efficiency and accuracy of our methods.
12 J. RASHIDINIA, R. MOHAMMADI: NON-POLYNOMIAL SPLINE APPROXIMATIONS FOR ɛ N = 6 N = N = N = 8 N = 56 Our eight order method 6 6.9( 9).04( ) 6.86( 4).5( 5).07( 6) 8.48( 9).57( 0).7( 4).44( 5).08( 6) 5.86( 8).6( 9).9( ).44( 5).7( 5) 8 8.5( 7).84( 9) 4.97( ).6( 4).06( 4) Our sixth order method 6.8( 7) 6.74( 0) 4.5( ) 6.5( 4).00( 4) 4.08( 8).8( 0).00( ) 4.( 4).04( 4) 5.84( 7) 4.9( 9).8( ).7( ).68( 4) ( 6) 6.9( 8) 4.4( 0).09( ).0( 4) Our fourth order method 6.8( 4) 4.44( 5) 7.( 6) 9.97( 7).( 7) 5.84( 4).0( 4).07( 5).89( 6).8( 7).50( ).47( 4) 5.69( 5) 8.09( 6).07( 6) 8.7( ) 8.67( 4).50( 4).( 5).00( 6) Sixth order method in [6] 6.9( 7) 6.74( 0) 4.6( ) 6.8( 4).09( 4) 4.08( 8).9( 0).00( ) 4.6( 4) 9.60( 5) 5.84( 7) 4.9( 9).8( ).7( ).6( 4) ( 6) 6.4( 8) 4.4( 0).0( ).84( 4) Fifth order method in [4] 6.( 6) 6.45( 9).40( ).0( ) 8.88( 5).00( 6).68( 8).6( 0).09( ).90( 4) 8.89( 6).6( 7).0( 9).08( ) 9.07( 4) ( 5) 9.98( 7).8( 8).4( 0) 9.9( ) Table. The maximum absolute errors for example. Example. Consider following boundary value problem []: ɛu + u = x, u(0) =, u() = + exp( / ɛ) with the exact solution, u(x) = x + exp( x/ ɛ). This problem has been solved using our method with different values of N = 8, 6,,, 8 and ɛ = 6,..., 5. The maximum absolute errors in solutions are tabulated in Table and compared with the methods in Ref. [].
13 48 TWMS J. PURE APPL. MATH., V., N., 00 ɛ N = 6 N = N = N = 8 N = 56 Fourth order method in [] 6.57( 5) 8.79( 7) 5.( 8).0( 9).05( 0) 8.7( 6) 4.4( 7).6( 8).6( 9).00( 0).84( 5) 8.67( 7) 6.65( 8) 4.9( 9).78( 0) 8.0( 4).6( 6).( 7).54( 8) 9.44( 0) Fourth order method in [] ( 5).5( 6).58( 7) 9.87( 9) 6.7( 0).00( 5).4( 6) 7.74( 8) 4.8( 9).0( 0) 5.45( 5).4( 6).4( 7).4( 8) 8.9( 0) 8.8( 4).( 5) 7.68( 7) 4.8( 8).0( 9) The method in [4] ( ).0( ) 5.08( 4).7( 4).7( 5) 7.( ).79( ) 4.48( 4).( 4).80( 5) 6.58( ).66( ) 4.5( 4).04( 4).60( 5) 8 6.6( ).6( ) 4.0( 4).0( 4).5( 5) The method in [5] 6 4.4( ).0( ).54( 4) 6.5( 5).58( 5).68( ) 9.0( 4) 5.6( 5).40( 5).50( 6).45( ) 8.40( 4).08( 4) 5.0( 5).0( 5) 8.4( ) 8.( 4).0( 4) 5.06( 5).6( 5) The method in [6] 6.0( 4) 7.47( 6) 4.67( 7).90( 8) 4.9( 9).8( 4) 8.00( 6) 5.00( 7).4( 8).99( 9).60( 4).00( 5) 6.6( 7).9( 8).( 9) 8.4( 4).47( 5) 9.( 7) 5.77( 8).7( 9) The method in [7] ( ).77( ) 4.45( 4).( 4).78( 5) 5.68( ).4( ).55( 4) 8.89( 5).( 5) 4.07( ).0( ).54( 4) 6.5( 5).58( 5) ( ).75( ) 4.( 4).08( 4).7( 5) Table. The maximum absolute errors for example. Example. Finally consider the following boundary value problem [7]: ɛu + [ + x( x)]u = + x( x) + [ ɛ x ( x) ( x)] exp[ ]+ ɛ +[ ɛ x( x) ] exp[ x/ ɛ], u(0) = u() = 0,
14 J. RASHIDINIA, R. MOHAMMADI: NON-POLYNOMIAL SPLINE APPROXIMATIONS FOR ɛ N = 8 N = 6 N = N = N = 8 Our eight order method 8.97( 9) 5.0( ).59( 4) 5.45( 5).76( 5) ( 8).7( 0).47( ) 8.( 5).( 5) 8.5( 7).84( 9) 4.96( ) 7.7( 5) 5.99( 5) 9.84( 6) 5.9( 8).6( 0).60( ).8( 4) 8 8.4( 5) 8.( 7).84( 9) 4.96( ) 8.60( 5) ( 4) 9.8( 6) 5.9( 8).6( 0).6( ) 5.0( ) 8.4( 5) 8.( 7).84( 9) 4.96( ) Our sixth order method ( 8) 4.8( 0).( ).70( 4).( 4) 6 6.7( 7) 5.50( 9).08( ).6( ).9( 4) 6.55( 6) 6.45( 8) 4.6( 0).0( ).7( 4) 5.0( 5) 6.70( 7) 5.49( 9).08( ).6( ) 8.6( 4) 6.55( 6) 6.44( 8) 4.6( 0).0( ) 56.55( ) 5.0( 5) 6.70( 7) 5.49( 9).08( ) ( ).6( 4) 6.55( 6) 6.44( 8) 4.6( 0) Our fourth order method 8 7.8( ).7( ) 6.0( 4).6( 4) 4.5( 5) 6.5( ).95( ).( ).04( 4) 7.88( 5).95( ) 7.0( ).( ) 5.84( 4).5( 4).95( ).9( ).9( ).( ).00( 4) 8 4.7( ).94( ) 7.0( ).( ) 5.5( 4) ( ).9( ).0( ).9( ).( ) 5 7.6( ) 4.7( ).95( ) 7.05( ).( ) The method in [] 8.0( 5) 7.0( 7) 4.8( 8).74( 9).7( 0) ( 5).96( 6).85( 7).5( 8) 7.4( 0).78( 4).8( 5) 7.54( 7) 4.67( 8).96( 9) 7.4( 4) 4.74( 5).96( 6).86( 7).6( 8) 8.70( ).78( 4).8( 5) 7.46( 7) 4.67( 8) 56 8.( ) 7.4( 4) 4.74( 5).98( 6).86( 7) 5.05( ).70( ).78( 4).8( 5) 7.46( 7) Table. The maximum absolute errors for example. with the exact solution, u(x) = + (x ) exp[ x/ ɛ] x exp[ ( x)/ ɛ]. This problem has been solved using our method with different values of N = 8, 6,,, 8 and ɛ = 6,..., 56 the maximum absolute errors in solutions are tabulated in Table 4.
15 50 TWMS J. PURE APPL. MATH., V., N., 00 ɛ N = 6 N = N = N = 8 Our eight order method ( 0) 6.00( ) 4.77( 5).( 6) 7.9( 9).( ).6( 4).94( 5).0( 7).74( 0).97( ) 9.0( 5) 8.9( 6) 5.08( 9) 9.8( ).5( 4) Our sixth order method 6.6( 8) 9.( ) 4.4( ) 7.4( 5).47( 7).00( 9) 5.04( ).78( 4).5( 6).05( 8) 6.07( ).75( ) 8.0( 5).05( 7) 7.( 0).57( ) Our fourth order method 6 5.5( 4) 7.79( 5).05( 5).7( 6).0( ).9( 4).70( 5).59( 6).6( ) 4.( 4) 6.9( 5) 9.45( 6) 8 5.7( ).08( ).74( 4).48( 5) The method in [7] 6.6( 8) 9.( ) 4.4( ).4( 5).47( 7).00( 9) 5.04( ).78( 4).5( 6).06( 8) 6.07( ).75( ) 8.0( 5).05( 7) 7.4( 0).57( ) Table 4. The maximum absolute errors for example. References [] Aziz, T., Khan, A., (00), A spline method for second-order singularly perturbed boundary-value problems, J. Comput. Appl. Math., 47(), pp [] Aziz, T., Khan, A., (00), Quintic Spline Approach to the Solution of a Singularly-Perturbed Boundary- Value Problem, J. Optim. Theory Appl., (), pp [] Chin, R.C.Y. and Krasny, R., (98), A hybrid asymptotic-finite element method for stifi two point boundary value problems, SIAM J. Sci. Stat. Comput., 4, pp [4] Jain, M.K. and Aziz, T., (98), Numerical Solution of Stiff and Convection-Diffusion Equations Using Adaptive Spline Function Approximation, Applied Mathematical Modeling, Appl. Math. Model., 7, pp [5] Kadalbajoo, M.K. and Bawa, R.K., (996), Variable-mesh difference scheme for singularly-perturbed boundary-value problem using splines, J. Optim. Theory Applic., 9, pp [6] Khan, A., Khan, I.. and Aziz, T., (006), Sextic spline solution of a singularly perturbed boundary-value problems, Appl. Math. Comput., 8, pp [7] Major, M.R., (985), Prog. Sci. Comput., Birkhauser, Boston, Massachusettes, 5, pp [8] Mohanty, R.K., Navnit Jha, Evans, D.J., (004), Spline in compression method for the numerical solution of singularly perturbed two-point singular boundary-value problems, Intern. J. Computer Math., 8, pp [9] Niijima, K., (980), On a three-point difference scheme for a singular perturbation problem without a first derivative term II, Mem.Numer.Math., 7, pp. -7. [0] Niijima, K., (980), On a three-point difference scheme for a singular perturbation problem without a first derivative term I, Mem.Numer.Math., 7, pp. -0. [] Rashidinia, J., (990), Applications of spline to numerical solution of differential equates, M.phil dissertation A.M.U.India. [] Rashidinia, J., (00), Intern. J. Engin. Sci.,, pp. -.
16 J. RASHIDINIA, R. MOHAMMADI: NON-POLYNOMIAL SPLINE APPROXIMATIONS FOR... 5 [] Rashidinia, J., Ghasemi, M., Mahmoodi, Z., (007), Spline approach to the solution of a singularly-perturbed boundary- value problems, Appl. Math. Comput., 89, pp [4] Surla, K., Herceg, D., and Cvetkovic, L., (99), A finally of exponential spline difference scheme, Review of Research, Faculty of Science, Mathematics Series, University of Novi Sad, 9, pp. -. [5] Surla, K., Stojanovic, M., (988), Solving singularly- perturbed boundary-value problems by splines in tension, J. Comput. Appl. Math., 4, pp [6] Surla, K., Vukoslavcevic, V., (995), A spline difference scheme for boundary value problems with a small parameter, Review of research, Faculty of science, Mathematics series, University of Novi Sad, 5, pp [7] Tirmizi, I. A., Haq, F., Siraj-Ul-Islam, (007), Appl. Math. Comput., to Apear. [8] Usmani, R. A., (99), The use of quartic splines in the numerical solution of a fourth-order boundary-value problem, J. Comput. Appl. Math., 44, pp Jalil Rashidinia - Associate professor, Head Department of Applied Mathematics in School of Mathematics, Iran University of Science and Technology Narmak, Tehran 6844, Iran since 006.He has done Ph.D. In Applied Mathematics (Numerical analysis),in A.M.U India on 994. Research Interests: Applied mathematics and computational sciences Numerical solution of ODE, PDE and IE Spline approximations.
NUMERICAL 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 informationNon-Polynomial Spline Solution of Fourth-Order Obstacle Boundary-Value Problems
Non-Polynomial Spline Solution of Fourth-Order Obstacle Boundary-Value Problems Jalil Rashidinia Reza Jalilian Abstract In this paper we use quintic non-polynomial spline functions to develop numerical
More 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 informationSolution of Fourth Order Boundary Value Problems by Numerical Algorithms Based on Nonpolynomial Quintic Splines
Journal of Numerical Mathematics and Stochastics, 4(1) : 13-25, 2012 http://www.jnmas.org/jnmas4-2.pdf JNM@S Euclidean Press, LLC Online: ISSN 2151-2302 Solution of Fourth Order Boundary Value Problems
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 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 informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationSolution of Fourth Order Obstacle Problems Using Quintic B-Splines
Applied Mathematical Sciences, Vol. 6, 202, no. 94, 465-4662 Solution of Fourth Order Obstacle Problems Using Quintic B-Splines Shahid S. Siddiqi, Ghazala Akram and Kalsoom Arshad Department of Mathematics
More informationNumerical solution of fourth order parabolic partial dierential equation using parametric septic splines
Hacettepe Journal of Mathematics and Statistics Volume 5 20, 07 082 Numerical solution of fourth order parabolic partial dierential equation using parametric septic splines Arshad Khan and Talat Sultana
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 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 solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method
Numerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method Y. Azari Keywords: Local RBF-based finite difference (LRBF-FD), Global RBF collocation, sine-gordon
More 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 informationParametric Spline Method for Solving Bratu s Problem
ISSN 749-3889 print, 749-3897 online International Journal of Nonlinear Science Vol4202 No,pp3-0 Parametric Spline Metod for Solving Bratu s Problem M Zarebnia, Z Sarvari 2,2 Department of Matematics,
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 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 informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
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 informationAMS527: Numerical Analysis II
AMS527: Numerical Analysis II Supplementary Material on Finite Different Methods for Two-Point Boundary-Value Problems Xiangmin Jiao Stony Brook University Xiangmin Jiao Stony Brook University AMS527:
More informationCubic spline Numerov type approach for solution of Helmholtz equation
Journal of Linear and Topological Algebra Vol. 03, No. 01, 2014, 47-54 Cubic spline Numerov type approach for solution of Helmholtz equation J. Rashidinia a, H. S. Shekarabi a and M. Aghamohamadi a a Department
More informationA UNIFORMLY ACCURATE FINITE ELEMENTS METHOD FOR SINGULAR PERTURBATION PROBLEMS. Mirjana Stojanović University of Novi Sad, Serbia and Montenegro
GLASNIK MATEMATIČKI Vol. 38(58(2003, 185 197 A UNIFORMLY ACCURATE FINITE ELEMENTS METHOD FOR SINGULAR PERTURBATION PROBLEMS Mirjana Stojanović University of Novi Sad, Serbia and Montenegro Abstract. We
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 informationThe semi-linear elliptic partial differential equations to be numerically solved, is of the form. y 2 U(x, y) = G(x, y, U(x, y)), (x, y) Ω (1)
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.3(017) No.1, pp.33-45 On the Convergence of Compact Arithmetic Averaging Scheme for Semi-linear d-elliptic Equations
More informationCubic B-spline Collocation Method for Fourth Order Boundary Value Problems. 1 Introduction
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.142012 No.3,pp.336-344 Cubic B-spline Collocation Method for Fourth Order Boundary Value Problems K.N.S. Kasi Viswanadham,
More informationA NEW APPROACH FOR SOLITON SOLUTIONS OF RLW EQUATION AND (1+2)-DIMENSIONAL NONLINEAR SCHRÖDINGER S EQUATION
A NEW APPROACH FOR SOLITON SOLUTIONS OF RLW EQUATION AND (+2-DIMENSIONAL NONLINEAR SCHRÖDINGER S EQUATION ALI FILIZ ABDULLAH SONMEZOGLU MEHMET EKICI and DURGUN DURAN Communicated by Horia Cornean In this
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
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 informationSolution of Non Linear Singular Perturbation Equation. Using Hermite Collocation Method
Applied Mathematical Sciences, Vol. 7, 03, no. 09, 5397-5408 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/ams.03.37409 Solution of Non Linear Singular Perturbation Equation Using Hermite Collocation
More informationA NEW INTEGRATOR FOR SPECIAL THIRD ORDER DIFFERENTIAL EQUATIONS WITH APPLICATION TO THIN FILM FLOW PROBLEM
Indian J. Pure Appl. Math., 491): 151-167, March 218 c Indian National Science Academy DOI: 1.17/s13226-18-259-6 A NEW INTEGRATOR FOR SPECIAL THIRD ORDER DIFFERENTIAL EQUATIONS WITH APPLICATION TO THIN
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 informationThe Homotopy Perturbation Method for Solving the Kuramoto Sivashinsky Equation
IOSR Journal of Engineering (IOSRJEN) e-issn: 2250-3021, p-issn: 2278-8719 Vol. 3, Issue 12 (December. 2013), V3 PP 22-27 The Homotopy Perturbation Method for Solving the Kuramoto Sivashinsky Equation
More informationFusion Higher -Order Parallel Splitting Methods for. Parabolic Partial Differential Equations
International Mathematical Forum, Vol. 7, 0, no. 3, 567 580 Fusion Higher -Order Parallel Splitting Methods for Parabolic Partial Differential Equations M. A. Rehman Department of Mathematics, University
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 informationTHE CONVECTION DIFFUSION EQUATION
3 THE CONVECTION DIFFUSION EQUATION We next consider the convection diffusion equation ɛ 2 u + w u = f, (3.) where ɛ>. This equation arises in numerous models of flows and other physical phenomena. The
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 informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
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 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 informationAdomian Decomposition Method for Solving the Kuramoto Sivashinsky Equation
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn:2319-765x. Volume 1, Issue 1Ver. I. (Jan. 214), PP 8-12 Adomian Decomposition Method for Solving the Kuramoto Sivashinsky Equation Saad A.
More informationComputational Non-Polynomial Spline Function for Solving Fractional Bagely-Torvik Equation
Math. Sci. Lett. 6, No. 1, 83-87 (2017) 83 Mathematical Sciences Letters An International Journal http://dx.doi.org/10.18576/msl/060113 Computational Non-Polynomial Spline Function for Solving Fractional
More informationExact Solutions for a Class of Singular Two-Point Boundary Value Problems Using Adomian Decomposition Method
Applied Mathematical Sciences, Vol 6, 212, no 122, 697-618 Exact Solutions for a Class of Singular Two-Point Boundary Value Problems Using Adomian Decomposition Method Abdelhalim Ebaid 1 and Mona D Aljoufi
More informationAn 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 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 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 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 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 informationAn Exponential High-Order Compact ADI Method for 3D Unsteady Convection Diffusion Problems
An Exponential High-Order Compact ADI Method for 3D Unsteady Convection Diffusion Problems Yongbin Ge, 1 Zhen F. Tian, 2 Jun Zhang 3 1 Institute of Applied Mathematics and Mechanics, Ningxia University,
More informationThe Derivation of Interpolants for Nonlinear Two-Point Boundary Value Problems
European Society of Computational Methods in Sciences and Engineering ESCMSE) Journal of Numerical Analysis, Industrial and Applied Mathematics JNAIAM) vol. 1, no. 1, 2006, pp. 49-58 ISSN 1790 8140 The
More information[21] B. J. McCartin, Theory of exponential splines, Journal of Approximation Theory 66 (1991) 1 23.
15 T. Ak, S. B. G. Karakoc, H. Triki, Numerical simulation for treatment of dispersive sallow water waves wit Rosenau- KdV equation, Te European Pysical Journal Plus 131 (1 (16 356-37. 16 M. Abbas, A.
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationSolution of Seventh Order Boundary Value Problem by Differential Transformation Method
World Applied Sciences Journal 16 (11): 1521-1526, 212 ISSN 1818-4952 IDOSI Publications, 212 Solution of Seventh Order Boundary Value Problem by Differential Transformation Method Shahid S. Siddiqi, Ghazala
More informationComputation Fluid Dynamics
Computation Fluid Dynamics CFD I Jitesh Gajjar Maths Dept Manchester University Computation Fluid Dynamics p.1/189 Garbage In, Garbage Out We will begin with a discussion of errors. Useful to understand
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 semi-lagrangian scheme for advection-diffusion equation
EPiC Series in Engineering Volume 3, 2018, Pages 162 172 HIC 2018. 13th International Conference on Hydroinformatics Engineering A semi-lagrangian scheme for advection-diffusion equation Ersin Bahar 1,
More informationHodograph Transformations in Unsteady MHD Transverse Flows
Applied Mathematical Sciences, Vol. 4, 010, no. 56, 781-795 Hodograph Transformations in Unsteady MHD Transverse Flows Pankaj Mishra Department of Mathematics, Faculty of Science Banaras Hindu University,
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 informationA fourth-order finite difference scheme for the numerical solution of 1D linear hyperbolic equation
203 (203) - Available online at www.ispacs.com/cna Volume 203, Year 203 Article ID cna-0048, Pages doi:0.5899/203/cna-0048 Research Article A fourth-order finite difference scheme for the numerical solution
More informationResearch Article Cubic Spline Iterative Method for Poisson s Equation in Cylindrical Polar Coordinates
International Scolarly Researc Network ISRN Matematical Pysics Volume 202, Article ID 2456, pages doi:0.5402/202/2456 Researc Article Cubic Spline Iterative Metod for Poisson s Equation in Cylindrical
More informationNumerical Solution of an Inverse Diffusion Problem
Applied Mathematical Sciences, Vol. 1, 2007, no. 18, 863-868 Numerical Solution of an Inverse Diffusion Problem Y. Ahmadizadeh Physics Dept., Faculty of Sciences Air University, Postal Code 1384663113,
More informationHydrodynamic Limit with Geometric Correction in Kinetic Equations
Hydrodynamic Limit with Geometric Correction in Kinetic Equations Lei Wu and Yan Guo KI-Net Workshop, CSCAMM University of Maryland, College Park 2015-11-10 1 Simple Model - Neutron Transport Equation
More informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More information1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load
1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load Nader Mohammadi 1, Mehrdad Nasirshoaibi 2 Department of Mechanical
More informationAn exponentially graded mesh for singularly perturbed problems
An exponentially graded mesh for singularly perturbed problems Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with P. Constantinou (UCY), S. Franz (TU Dresden)
More informationThe Trigonometric Cubic B-spline Algorithm for Burgers Equation
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.4(017) No., pp.10-18 The Trigonometric Cubic B-spline Algorithm for Burgers Equation Idris Dag 1, Ozlem Ersoy Hepson,
More informationLecture Notes on Fractals and Multifractals
Lecture Notes on Fractals and Multifractals Topics in Physics of Complex Systems Version: 8th October 2010 1 SELF-AFFINE SURFACES 1 Self-affine surfaces The most common example of a self-affine signal
More informationChapter 1 Numerical approximation of data : interpolation, least squares method
Chapter 1 Numerical approximation of data : interpolation, least squares method I. Motivation 1 Approximation of functions Evaluation of a function Which functions (f : R R) can be effectively evaluated
More informationNumerical Methods for Engineers and Scientists
Numerical Methods for Engineers and Scientists Second Edition Revised and Expanded Joe D. Hoffman Department of Mechanical Engineering Purdue University West Lafayette, Indiana m MARCEL D E К К E R MARCEL
More informationNumerical solution of General Rosenau-RLW Equation using Quintic B-splines Collocation Method
Available online at www.ispacs.com/cna Volume 2012, Year 2012 Article ID cna-00129, 16 pages doi:10.5899/2012/cna-00129 Research Article Numerical solution of General Rosenau-RLW Equation using Quintic
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationJournal of Computational and Applied Mathematics. Multigrid method for solving convection-diffusion problems with dominant convection
Journal of Computational and Applied Mathematics 226 (2009) 77 83 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationInitial value problems for ordinary differential equations
Initial value problems for ordinary differential equations Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 IVP
More informationSOLVING SOLVABLE QUINTICS. D. S. Dummit
D. S. Dummit Abstract. Let f(x) = x 5 + px 3 + qx + rx + s be an irreducible polynomial of degree 5 with rational coefficients. An explicit resolvent sextic is constructed which has a rational root if
More informationError Analyses of Stabilized Pseudo-Derivative MLS Methods
Department of Mathematical Sciences Revision: April 11, 2015 Error Analyses of Stabilized Pseudo-Derivative MLS Methods Don French, M. Osorio and J. Clack Moving least square (MLS) methods are a popular
More informationOn the Solution to Octic Equations
The Mathematics Enthusiast Volume 4 Number 2 Article 6 6-2007 On the Solution to Octic Equations Raghavendra G. Kulkarni Let us know how access to this document benefits you. Follow this and additional
More informationNumerical Simulations on Two Nonlinear Biharmonic Evolution Equations
Numerical Simulations on Two Nonlinear Biharmonic Evolution Equations Ming-Jun Lai, Chun Liu, and Paul Wenston Abstract We numerically simulate the following two nonlinear evolution equations with a fourth
More informationOutline. 1 Boundary Value Problems. 2 Numerical Methods for BVPs. Boundary Value Problems Numerical Methods for BVPs
Boundary Value Problems Numerical Methods for BVPs Outline Boundary Value Problems 2 Numerical Methods for BVPs Michael T. Heath Scientific Computing 2 / 45 Boundary Value Problems Numerical Methods for
More informationFinite 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 informationThe WKB local discontinuous Galerkin method for the simulation of Schrödinger equation in a resonant tunneling diode
The WKB local discontinuous Galerkin method for the simulation of Schrödinger equation in a resonant tunneling diode Wei Wang and Chi-Wang Shu Division of Applied Mathematics Brown University Providence,
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationNUMERICAL METHODS. x n+1 = 2x n x 2 n. In particular: which of them gives faster convergence, and why? [Work to four decimal places.
NUMERICAL METHODS 1. Rearranging the equation x 3 =.5 gives the iterative formula x n+1 = g(x n ), where g(x) = (2x 2 ) 1. (a) Starting with x = 1, compute the x n up to n = 6, and describe what is happening.
More informationList of Publications of Dr. Ghazala Akram
List of Publications List of Publications of Dr. Ghazala Akram Published 1). End Conditions for Interpolatory Septic Ghazala Akram and Shahid S. Siddiqi, International Journal of Computer Mathematics,
More informationSelf-similar solutions for the diffraction of weak shocks
Self-similar solutions for the diffraction of weak shocks Allen M. Tesdall John K. Hunter Abstract. We numerically solve a problem for the unsteady transonic small disturbance equations that describes
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
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 informationA robust uniform B-spline collocation method for solving the generalized PHI-four equation
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 11, Issue 1 (June 2016), pp. 364-376 Applications and Applied Mathematics: An International Journal (AAM) A robust uniform B-spline
More informationChebyshev finite difference method for a Two Point Boundary Value Problems with Applications to Chemical Reactor Theory
Iranian Journal of Mathematical Chemistry, Vol. 3, o., February 22, pp. 7 IJMC Chebyshev finite difference method for a Two Point Boundary Value Problems with Applications to Chemical Reactor Theory ABBAS
More informationBindel, Spring 2012 Intro to Scientific Computing (CS 3220) Week 12: Monday, Apr 16. f(x) dx,
Panel integration Week 12: Monday, Apr 16 Suppose we want to compute the integral b a f(x) dx In estimating a derivative, it makes sense to use a locally accurate approximation to the function around the
More informationHomotopy Perturbation Method for Solving the Second Kind of Non-Linear Integral Equations. 1 Introduction and Preliminary Notes
International Mathematical Forum, 5, 21, no. 23, 1149-1154 Homotopy Perturbation Method for Solving the Second Kind of Non-Linear Integral Equations Seyyed Mahmood Mirzaei Department of Mathematics Faculty
More informationDu Fort Frankel finite difference scheme for Burgers equation
Arab J Math (13) :91 11 DOI 1.17/s465-1-5-1 RESEARCH ARTICLE K. Pandey Lajja Verma Amit K. Verma Du Fort Frankel finite difference scheme for Burgers equation Received: 8 February 1 / Accepted: 1 September
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationMath 578: Assignment 2
Math 578: Assignment 2 13. Determine whether the natural cubic spline that interpolates the table is or is not the x 0 1 2 3 y 1 1 0 10 function 1 + x x 3 x [0, 1] f(x) = 1 2(x 1) 3(x 1) 2 + 4(x 1) 3 x
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationAnalytical Solution to Intra-Phase Mass Transfer in Falling Film Contactors via Homotopy Perturbation Method
International Mathematical Forum, Vol. 6, 2011, no. 67, 3315-3321 Analytical Solution to Intra-Phase Mass Transfer in Falling Film Contactors via Homotopy Perturbation Method Hooman Fatoorehchi 1 and Hossein
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 informationA numerical solution of a Kawahara equation by using Multiquadric radial basis function
Mathematics Scientific Journal Vol. 9, No. 1, (013), 115-15 A numerical solution of a Kawahara equation by using Multiquadric radial basis function M. Zarebnia a, M. Takhti a a Department of Mathematics,
More informationHIGH ACCURACY NUMERICAL METHODS FOR THE SOLUTION OF NON-LINEAR BOUNDARY VALUE PROBLEMS
ABSTRACT Of The Thesis Entitled HIGH ACCURACY NUMERICAL METHODS FOR THE SOLUTION OF NON-LINEAR BOUNDARY VALUE PROBLEMS Submitted To The University of Delhi In Partial Fulfillment For The Award of The Degree
More informationNumerical Solution of Fourth Order Boundary-Value Problems Using Haar Wavelets
Applied Mathematical Sciences, Vol. 5, 20, no. 63, 33-346 Numerical Solution of Fourth Order Boundary-Value Problems Using Haar Wavelets Fazal-i-Haq Department of Maths/Stats/CS Khyber Pukhtoon Khwa Agricultral
More informationConstruction of very high order Residual Distribution Schemes for steady problems
Construction of very high order Residual Distribution Schemes for steady problems Rémi Abgrall, Mario Ricchiuto, Cédric Tavé, Nadège Villedieu and Herman Deconinck Mathématiques Appliquées de Bordeaux,
More informationNew Diagonal-Norm Summation-by-Parts Operators for the First Derivative with Increased Order of Accuracy
AIAA Aviation -6 June 5, Dallas, TX nd AIAA Computational Fluid Dynamics Conference AIAA 5-94 New Diagonal-Norm Summation-by-Parts Operators for the First Derivative with Increased Order of Accuracy David
More information