Smoothness of solutions with respect to multi-strip integral boundary conditions for nth order ordinary differential equations

Size: px
Start display at page:

Download "Smoothness of solutions with respect to multi-strip integral boundary conditions for nth order ordinary differential equations"

Transcription

1 396 Nonlinear Analysis: Modelling and Control, 2014, Vol. 19, No. 3, ttp://dx.doi.org/ /na Smootness of solutions wit respect to multi-strip integral boundary conditions for nt order ordinary differential equations Jonny Henderson Department of Matematics, Baylor University Waco, Texas, USA Received: 4 December 2013 / Revised: 17 Marc 2014 / Publised online: 30 June 2014 Abstract. Under certain conditions, solutions of te boundary value problem, y n) = fx, y, y,..., y n 1) ), a < x < b, y i 1) x 1) = y i, i = 1,..., n 1, yx m 2) γi yx) dx = y n, a < x 1 < ξ 1 < η 1 < ξ 2 < η 2 < < ξ m < η m < x 2 < b, are differentiated wit respect to te boundary conditions. Keywords: nt order ordinary differential equation, multi-strip integral boundary conditions, smoot dependence, boundary data. 1 Introduction In tis paper, we will be concerned wit differentiating solutions of boundary value problems wit respect to boundary data for te nt order ordinary differential equation, y n) = f x, y, y,..., y n 1)), a < x < b, 1) satisfying te Diriclet and multi-strip integral boundary conditions, y i 1) x 1 ) = y i, i = 1,..., n 1, yx 2 ) yx) dx = y n, 2) were a < x 1 < ξ 1 < η 1 < ξ 2 < η 2 < < ξ m < η m < x 2 < b, R, i = 1,..., m, and y 1,..., y n R, and were we assume: i) fx, s 1, s 2,..., s n ) : a, b) R n R is continuous, ii) f/ s i x, s 1, s 2,..., s n ) : a, b) R n R is continuous, i = 1, 2,..., n, and iii) Solutions of initial value problems for 1) extend to a, b). c Vilnius University, 2014

2 Smootness of solutions wit respect to multi-strip integral BC 397 Condition iii) is not necessary for te results of tis paper, yet, by assuming iii), we avoid continually making statements in terms of solutions maximal intervals of existence. Under uniqueness assumptions on solutions of 1), 2), we will establis analogues of a result tat Hartman 1 attributes to Peano concerning differentiation of solutions of 1) wit respect to initial conditions. For our differentiation wit respect to boundary conditions results, given a solution yx) of 1), we will give muc attention to te variational equation for 1) along yx), wic is defined by z n) = n f s i x, yx), y x),..., y n 1) x) ) z i 1). 3) Interest in nonlocal boundary value problems for differential equations involving integral boundary conditions as been ongoing for several years. To see only few of tese papers, we refer te reader to te papers 2 9. And very recently, Amad and Ntouyas 10 initiated researc regarding multipoint nonlocal integral boundary conditions suc as seen in 2). In describing suc boundary conditions, tey coined te term multistrip integral boundary conditions. Suc boundary conditions can be interpreted in te sense tat a controller at te rigt end of te interval under consideration is influenced by a discrete distribution of finite many sensors or strips) of arbitrary lengt expressed in terms of integral boundary conditions. Subsequent to tat paper, Amad and Ntouyas ave put fort a couple of additional papers devoted to solutions of boundary value problems involving multi-strip integral boundary conditions for bot fractional differential equations and fractional differential inclusions; see 11, 12. It can also be pointed out tat, under suitable measures, te boundary conditions can be considered in te form of Stieltjes integrals; readers can find of interest te papers, and In te same way, tere ave been many papers devoted to smootness of solutions of boundary value problems in regard to smootness of te differential equation s nonlinearity, as well as in regard to te smootness of te boundary conditions. For a view of ow tis work as evolved, involving not only boundary value problems for ordinary differential equations, but also discrete versions, we suggest te manifold results in te classical papers and 8, 9, 29 34, as well as te more current papers 35, 36 and Te teorem for wic we seek an analogue and attributed to Peano by Hartman can be stated in te context of 1) as follows: Teorem 1 Peano). Assume tat wit respect of 1), conditions i) iii) are satisfied. Let x 0 a, b) and yx) := yx, x 0, c 1,..., c n ) denote te solution of 1) satisfying te initial conditions y i 1) x 0 ) = c i, i = 1,..., n. Ten, a) For j = 1,..., n, α j := y/ c j exists on a, b), is te solution of te variational equation 3) along yx), and satisfies te initial conditions, α i 1) j x 0 ) = δ ij, i = 1,..., n. Nonlinear Anal. Model. Control, 2014, Vol. 19, No. 3,

3 398 J. Henderson b) y/ x 0 exists on a, b), and β := y/ x 0 is te solution of te variational equation 3) along yx) satisfying te initial conditions, β i 1) x 0 ) = y i) x 0 ), i = 1,..., n. c) y/ x 0 x) = n j=1 yj) x 0 ) y/ c j x). In addition, our analogue of Teorem 1 depends on uniqueness of solutions of 1), 2), a condition we list as an assumption: iv) Given a < x 1 < ξ 1 < η 1 < ξ 2 < η 2 < < ξ m < η m < x 2 < b, if y i 1) x 1 ) = z i 1) x 1 ), i = 1,..., n 1, and yx 2 ) m γ i yx) dx = zx 2 ) m γ i zx) dx, were yx) and zx) are solutions of 1), ten yx) zx). We will also make extensive use of a similar uniqueness condition on 3) along solutions yx) of 1). v) Given a < x 1 < ξ 1 < η 1 < ξ 2 < η 2 < < ξ m < η m < x 2 < b and a solution yx) of 1), if u i 1) x 1 )=0,,..., n 1, and ux 2 ) m γ i ux) dx=0, were ux) is a solution of 3) along yx), ten ux) 0. Remark 1. We observe tat, if v) is assumed, ten for α n x) in Teorem 1, α n x 2 ) 2 An analogue of Peano s teorem for 1), 2) α n x) dx 0. 4) In tis section, we state and prove our analogue of Teorem 1 for boundary value problem 1), 2). Continuous dependence of solutions on boundary conditions plays a fundamental role for suc a differentiation result. Proof of continuous dependence usually makes application of te Brouwer teorem on invariance of domain. Te spirit of suc arguments can be found in 36 or 38; we state te continuity result ere, but we omit te details. Teorem 2. Assume i) iv) are satisfied wit respect to 1). Let ux) be a solution of 1) on a, b), and let a < c < x 1 < ξ 1 < η 1 < ξ 2 < η 2 < < ξ m < η m < x 2 < d < b be given. Ten, tere exists a δ > 0 suc tat, for x j t j < δ, j = 1, 2, ρ i < δ and η i σ i < δ, i = 1,..., m, u j 1) x 1 ) y j < δ, j = 1,..., n 1, ux 2 ) m γ i ux) dx y n < δ, and γ k ζ k < δ, k = 1,..., m, tere exists a unique solution u δ x) of 1) suc tat u j 1) δ t 1 ) = y j, j = 1,..., n 1, u δ t 2 ) m ζ σi i ρ i u δ x) dx = y n, and {u j) δ x)} converges uniformly to uj) x), as δ 0, on c, d, for j = 0, 1,..., n 1. We now present te result of tis paper.

4 Smootness of solutions wit respect to multi-strip integral BC 399 Teorem 3. Assume conditions i) v) are satisfied. Let ux) be a solution of 1) on a, b). Let a < x 1 < ξ 1 < η 1 < ξ 2 < η 2 < < ξ m < η m < x 2 < b be given, so tat ux) = ux, x 1, x 2, ξ 1, η 1,..., ξ m, η m, u 1,..., u n, γ 1,..., γ m ), were u j 1) x 1 ) = u j, j = 1,..., n 1, and ux 2 ) m ux) dx = u n. Ten, a) For j = 1,..., n 1, r j := u/ u j exists on a, b), is te solution of te variational equation 3) along ux), and satisfies te boundary conditions, r i 1) j x 1 ) = δ ij, i = 1,..., n 1, r j x 2 ) ξ 1 r j x) dx = 0. b) r n := u/ u n exists on a, b), is te solution of 3) along ux), and satisfies te boundary conditions, r i 1) n x 1 ) = 0, i = 1,..., n 1, r n x 2 ) ξ 1 r n x) dx = 1. c) z 1 := u/ x 1 and z 2 := u/ x 2 exist on a, b), are solutions of 3) along ux), and satisfy te respective boundary conditions, z i 1) 1 x 1 ) = u i) x 1 ), i = 1,..., n 1, z 1 x 2 ) z i 1) 2 x 1 ) = 0, i = 1,..., n 1, z 2 x 2 ) z 1 x) dx = 0, z 2 x) dx = u x 2 ). d) For eac j = 1,..., m, w j := u/ exists on a, b), is te solution of 3) along ux), and satisfies te boundary conditions, w i 1) j x 1 ) = 0, i = 1,..., n 1, w j x 2 ) w j x) dx = γ j u ). e) For eac j = 1,..., m, q j := u/ η j exists on a, b), is te solution of 3) along ux), and satisfies te boundary conditions, q i 1) j x 1 ) = 0, i = 1,..., n 1, q j x 2 ) q j x) dx = γ j uη j ). f) For eac j = 1,..., m, p j := u/ γ j exists on a, b), is te solution of 3) along ux), and satisfies te boundary conditions, p i 1) j x 1 ) = 0, i = 1,..., n 1, p j x 2 ) p j x) dx = ux) dx. Nonlinear Anal. Model. Control, 2014, Vol. 19, No. 3,

5 400 J. Henderson g) Te partial derivatives satisfy, n 1 u x) = u j) x 1 ) u x), x 1 u j j=1 u x) = γ j u ) u x), u n u γ j x) = u ux) dx x). u n u x) = u x 2 ) u x), x 2 u n u x) = γ j uη j ) u x), η j u n Proof. Wit ux) = ux, x 1, x 2, ξ 1, η 1,..., ξ m, η m, u 1,..., u n, γ 1,..., γ m ), as given in te statement of te teorem, many of te results will be establised by considering ux) as a solution of an initial value problem for 1). In particular, from te boundary value notation, u i 1) x 1 ) = u i, i = 1,..., n 1, and we will let β n = u n 1) x 1 ). 5) Ten, using te notation of Teorem 1 for solutions of initial value problems for 1) and viewing ux) as a solution of an initial value problem, we will frequently intercange notation by using, ux) = ux, x 1, x 2, ξ 1, η 1,..., ξ m, η m, u 1,..., u n, γ 1,..., γ m ) = yx, x 1, u 1,..., u n 1, β n ). For part a), we fix j = 1,..., n 1, and for notational sortand purposes, we denote ux) by ux,, u j ). Let δ > 0 be as in Teorem 2. Let 0 < < δ be given and define Ten, for every 0, r j x) = 1 ux,, uj + ) ux,, u j ). r j 1) j x 1 ) = 1 u j + u j = 1, and r i 1) j x 1 ) = 1 u i u i = 0, i {1,..., n 1} \ {j}, r j x 2 ) r j x) dx = 1 ux2,, u j + ) ux 2,, u j ) 1 = 1 u n u n = 0. ux,, uj + ) ux,, u j ) dx

6 Smootness of solutions wit respect to multi-strip integral BC 401 Wit β n as defined in 5), let = ) = u n 1) x 1,, u 1 + ) β n. By Teorem 2, = ) 0, as 0. Using te notation of Teorem 1 for solutions of initial value problems for 1) and viewing te solutions u as solutions of initial value problems, we ave r j x) = 1 yx, x1, u 1, u 2,..., u j +,..., u n 1, β n + ) yx, x 1, u 1, u 2,..., u j,..., u n 1, β n ). Ten, by utilizing a telescoping sum, we ave r j x) = 1 { yx, x1, u 1, u 2,..., u j +,..., u n 1, β n + ) yx, x 1, u 1, u 2,..., u j,..., u n 1, β n + ) } + { yx, x 1, u 1, u 2,..., u j,..., u n 1, β n + ) yx, x 1, u 1, u 2,..., u j,..., u n 1, β n ) }. By Teorem 1 and te Mean Value teorem, we obtain r j x) = 1 α j x, yx, x1, u 1,..., u j +,..., u n 1, β n + ) ) u j + u j ) + 1 α n x, yx, x1, u 1,..., u n 1, β n + ) ) β n + β n ), were α j x, y )) is te solution of te variational equation 3) along y ) and satisfies α i 1) j x 1 ) = δ ji, i = 1,..., n, and α n x, y )) is te solution of te variational equation 3) along y ) and satisfies α i 1) n x 1 ) = δ ni, i = 1,..., n. Furtermore, u j + is between u j and u j +, and β n + is between β n and β n +. Now simplifying, r j x) = α j x, yx, x1, u 1,..., u j +,..., u n 1, β n + ) ) + α n x, yx, x1, u 1,..., u n 1, β n + ) ). Tus, to sow lim r j x) exists, it suffices to sow lim / exists. Now, by Remark 1, α n x2, y ) ) ) α n x, y ) dx 0. Nonlinear Anal. Model. Control, 2014, Vol. 19, No. 3,

7 402 J. Henderson However, we derived above tat r j x 2 ) m obtain + γ i r j x) dx = 0, from wic we = α j x 2, yx, x 1,, u j +,, β n + )) α n x 2, yx, x 1, u 1,, β n + )) m γ i α n x, yx, x 1, u 1,, β n + )) dx m γ i α j x, yx, x 1,, u j +,, β n + )) dx α n x 2, yx, x 1, u 1,, β n + )) m α n x, yx, x 1, u 1,, β n + )) dx. As a consequence of continuous dependence, we can let 0, so tat lim = α jx 2, yx, x 1, u 1,, β n )) + m α n x 2, yx, x 1, u 1,, β n )) m γ i = α jx 2, ux)) + m α n x 2, ux)) m γ i γ i α j x, ux)) dx := A. α n x, ux)) dx Let r j x) = lim r j x), and note by construction of r j x), Furtermore, α j x, yx, x 1, u 1,, β n )) dx α n x, yx, x 1, u 1,, β n )) dx r j x) = u u j x, x 1, x 2, ξ 1, η 1,..., ξ m, η m, u 1,..., u n, γ 1,..., γ m ). r j x) = lim r j x) = α j x, yx, x1, u 1,, β n ) ) + Aα 2 x, yx, x1, u 1,, β n ) ) = α 1 x, ux) ) + Aα2 x, ux) ), wic is a solution of te variational equation 3) along ux). In addition because of te boundary conditions satisfied by r j x), we also ave, r i 1) j x 1 ) = δ ji, i = 1,..., n 1, and r j x 2 ) Tis completes te argument for u/ u j. r j x) dx = 0. For part b), tere are similarities wit te previous argument, yet tere are significant enoug differences for us to include te details concerning te caracterization of u/ u n. For tis consideration, we denote ux) by ux,, u n ). Again, let δ > 0 be as in Teorem 2. Let 0 < < δ be given and define Tis time, for 0, r n x) = 1 ux,, un + ) ux,, u n ). r i 1) n x 1 ) = 1 u i u i = 0, i = 1,..., n 1,

8 Smootness of solutions wit respect to multi-strip integral BC 403 and r n x 2 ) r n x) dx = 1 ux2,, u n + ) ux 2,, u n ) 1 Again wit β n defined in 5), let = 1 u n + u n = 1. ux,, un + ) ux,, u n ) dx = ) = u n 1) x 1,, u n + ) β n. As before, = ) 0, as 0. Employing te notation of Teorem 1 for solutions of initial value problems for 1) and viewing te solutions u as solutions of initial value problems, we ave r n x) = 1 yx, x1, u 1,..., u n 1, β n + ) yx, x 1, u 1,..., u n 1, β n ). By Teorem 1 and te Mean Value teorem, we obtain r n x) = 1 α n x, yx, x1, u 1,..., u n 1, β n + ) ) β n + β n ) = α n x, yx, x1, u 1,..., u n 1, β n + ) ), were β n + is between β n and β n +, and α n x, y )) is te solution of te variational equation 3) along y ) and satisfies, α i 1) n x 1 ) = δ ni, i = 1,..., n. Tus, to sow lim r n x) exists, it suffices to sow lim / exists. By Remark 1, α n x2, y ) ) ) α n x, y ) dx 0, and we also ave above tat r n x 2 ) m γ i r n x) dx = 1, from wic we obtain = 1 α n x 2, yx, x 1, u 1,, β n + )) m γ i α n x, yx, x 1, u 1,, β n + )) dx. Nonlinear Anal. Model. Control, 2014, Vol. 19, No. 3,

9 404 J. Henderson By continuous dependence, we can let 0, so tat lim = 1 α n x 2, yx, x 1, u 1,, β n )) m γ i α n x, yx, x 1, u 1,, β 2 )) dx 1 = α n x 2, ux)) m α n x, ux)) dx = B. Let r n x) = lim r n x), and ten by construction of r n x), Moreover, r n x) = u u n x, x 1, x 2, ξ 1, η 1,..., ξ m, η m, u 1,..., u n, γ 1,..., γ m ). r n x) = lim r n x) = Bα n x, yx, x1, u 1,, u n 1, β n ) ) = Bα n x, ux) ), wic is a solution of te variational equation 3) along ux). Because of te boundary conditions satisfied by r n x), we also ave, r i 1) n x 1 ) = 0, i = 1,..., n 1, and r n x 2 ) And tis completes te argument for u/ u n. r n x) dx = 1. For part c) of te teorem, we will produce te details for u/ x 2, wit te arguments for u/ x 1 being somewat along te same lines. For tis case, we denote ux) by ux, x 2, ). Wit δ > 0 as in Teorem 2, let 0 < < δ be given, and define z 2 x) = 1 ux, x2 +, ) ux, x 2, ). First, we consider boundary conditions. We ave z i 1) 2 x 1 ) = 1 u i u i = 0, i = 1,..., n 1. Next, by employing te Mean Value teorem, z 2 x 2 ) = 1 1 z 2 x) dx ux 2, x 2 +, ) ux 2, x 2, ) ux, x 2 +, ) dx ux, x 2, ) dx

10 Smootness of solutions wit respect to multi-strip integral BC 405 = 1 1 ux 2 +, x 2 +, ) ux 2, x 2, ) ux, x 2, ) dx ux, x 2 +, ) dx 1 ux2 +, x 2 +, ) ux 2, x 2 +, ) = 1 u n u n 1 u ν, x 2 +, ) = u ν, x 2 +, ), were ν is between x 2 and x 2 +. In passing to te limit, we ave lim { z 2 x 2 ) z 2 x) dx } = u x 2, x 2, ) = u x 2 ). Next, we deal wit te existence of lim z 2 x). Wit β n as defined in 5), tis time we let = ) = u n 1) x 1, x 2 +, ) β n, and by Teorem 2, = ) 0, as 0. As in parts a) and b), we use te notation of Teorem 1 for solutions of initial value problems for 1) and viewing te solutions u as solutions of initial value problems, we ave z 2 x) = 1 ux, x2 +, ) ux, x 2, ) = 1 yx, x1, u 1,..., u n 1, β n + ) yx, x 1, u 1,..., u n 1, β n ) = 1 α n x, yx, x1, u 1,..., u n 1, β n + ) ) = α n x, yx, x1, u 1,..., u n 1, β n + ) ), were α n x, y )) is te solution of 3) along y ) and satisfies α i 1) n x 1 ) = δ ni, i = 1,..., n, and β n + lies between β n and β n +. As before, to sow lim z 2 x) exists, it suffices to sow lim / exists. Now, recalling from above tat z 2 x 2 ) m γ i z 2 x) dx = u ν, x 2 +, ), it follows from Remark 1 tat, = u ν, x 2 +, ) α n x 2, yx, x 1, u 1,, β n + )) m γ i α n x, yx, x 1, u 1,, β n + )) dx. Nonlinear Anal. Model. Control, 2014, Vol. 19, No. 3,

11 406 J. Henderson And passing to te limit due to continuous dependence, we ave, lim From above, = u x 2 ) α n x 2, ux)) m α n x, ux)) dx z 2 x) = α n x, yx, x1, u 1,..., β n + ) ), := C. from wic we can evaluate te limit as 0, and if we let z 2 x) = lim z 2 x), we ave z 2 x) = u/ x 2. Tat is, we obtain z 2 x) = u x 2 x, x 1, x 2, ξ 1, η 1,..., ξ m, η m, u 1,..., u n, γ 1,..., γ m ) = lim z 2 x) = Cα n x, ux) ), wic is a solution of 3) along ux). In addition, from above computations, z 2 x) satisfies te boundary conditions, z 2 x 2 ) z i 1) 2 x 1 ) = lim z i 1) 2 x 1 ) = 0, i = 1,..., n 1, z 2 x) dx = lim Tis completes te proof of part c). For part d), fix j {1,..., m}, and define z 2 x 2 ) J := {1,..., m} \ {j}. z 2 x) dx = u x 2 ). In dealing wit caracterization of u/, we denote ux) by ux,, ). Let δ > 0 as in Teorem 2 and let 0 < < δ be given. Define w j x) = 1 ux,, ξj + ) ux,, ). We first look at boundary conditions satisfied by w j x). To begin wit, w i 1) j x 1 ) = 1 u i 1) x 1,, + ) u i 1) x 1,, ) = 1 u i u i = 0, i = 1,..., n 1.

12 Smootness of solutions wit respect to multi-strip integral BC 407 Next, by employing te Mean Value teorem for integrals, w j x 2 ) = 1 = 1 γ j γ j + γ j ξ j+ w j x) dx ux,, + ) dx + ux,, + ) dx + } { { ux 2,, + ) i J ux 2,, ) ux,, + ) dx + u n u n = 1 γ juc,, + ) = γ j uc,, + ), ux,, + ) dx ux,, ) dx for some c inclusively between and +. By Teorem 2, we can compute te limit, lim w j x 2 ) w j x) dx = γ j u ). Now, we deal wit te existence of lim w j x). Let β n be as in 5), and set = u n 1) x 1,, + ) β n. By Teorem 2, 0, as 0, and upon employing initial value solutions notation, w j x) = 1 yx, x1, u 1,..., u n 1, β n + ) yx, x 1, u 1,..., u n 1, β n ) = α n x, yx, x1, u 1,..., u n 1, β n + ) ), were β n + lies between β n and β n +, and α n x, y )) is as in te cases above. From Remark 1, we can solve for = w j x 2 ) m γ i α n x 2, y )) m from wic, using te above limit, we can calculate, lim = γ j u ) α n x 2, ux)) m w j x) dx α n x, y )) dx, α n x, ux)) dx := D. } Nonlinear Anal. Model. Control, 2014, Vol. 19, No. 3,

13 408 J. Henderson As as consequence, lim w j x) exists, and we define w j x) := lim w j x). In particular, w j x) = u x, x 1, x 2, ξ 1, η 1,..., ξ m, η m, u 1,..., u n, γ 1,..., γ m ) = lim w j x) = Dα n x, ux) ), wic is a solution of 3) along ux). In addition, w j x) satisfies te boundary conditions, w j x 2 ) w i 1) j x 1 ) = lim w j x) dx = lim Te proof of part d) is complete. j x 1 ) = 0, i = 1,..., n 1, η i w j x 2 ) w j x) dx = γ j u ). w i 1) For part e), te arguments are completely analogous to tose just given for part d). For part f), fix j {1,..., m}, and let J be as defined in te proof of part d). In caracterizing u/ γ j, we will designate ux) by ux,, γ j ). Wit δ > 0 and 0 < < δ given as usual, define p j x) = 1 ux,, γj + ) ux,, γ j ). As in te previous cases, we first consider boundary conditions satisfied by p j x). To begin wit, and p j x 2 ) = 1 = 1 i J p i 1) j x 1 ) = 1 u i 1) x 1,, γ j + ) u i 1) x 1,, γ j ) = 1 ui u i = 0, i = 1,..., n 1, p j x) dx ux,, γ j + ) dx + ux 2,, γ j + ) γ j + ) ux,, γ j + ) dx { ux,, γ j + ) dx + u n u n ux 2,, γ j ) = ux,, γ j + ) dx, ux,, γ j + ) dx ux,, γ j ) dx }

14 Smootness of solutions wit respect to multi-strip integral BC 409 from wic we can take te limit, lim p j x 2 ) p j x) dx = ux) dx. Finally, in considering te existence of lim p j x), let β n be as in 5), and again, define = u n 1) x 1,, γ j + ) β n, for wic by continuous dependence 0, as 0. Using initial value solutions notation, we ave p j x) = 1 yx, x1, u 1,..., u n 1, β n + ) yx, x 1, u 1,..., u n 1, β n ) = α 2 x, yx, x1, u 1,..., u n 1, β n + ) ), wit β n + between β n and β n +, and α n x, y )) as usual. In view of Remark 1, we can solve for = p j x 2 ) m γ i p j x) dx α n x 2, y )) η m. i α n x, y )) dx Using te last above limit, we can calculate, lim = ηj ux) dx α n x 2, ux)) m It follows tat p j x) := lim p j x) exists, and α n x, ux)) dx := E. p j x) = u γ j x, x 1, x 2, ξ 1, η 1,..., ξ m, η m, u 1,..., u n, γ 1,..., γ m ) = Eα n x, ux) ), is a solution of 3) along ux). And p j x) satisfies te boundary conditions, p i 1) j x 1 ) = 0, i = 1,..., n 1, and p j x 2 ) Tis completes te proof of part f). p j x) dx = ux) dx. Part g) of te teorem is immediate by verifying tat eac side of te respective equations are solutions of 3) along ux) and satisfy te same boundary conditions, and ten assumption v) establises te equalities. Te proof is complete. Nonlinear Anal. Model. Control, 2014, Vol. 19, No. 3,

15 410 J. Henderson References 1. P. Hartman, Ordinary Differential Equations, Wiley, New York, B. Amad, S.K. Ntouyas, A study of iger-order nonlinear ordinary differential equations wit four-point nonlocal integral boundary conditions, J. Appl. Mat. Comput., 39:97 108, 2012). 3. M. Bencora, F. Beroun, J.J. Nieto, Existence results for impulsive boundary value problem wit integral boundary conditions, Dyn. Syst. Appl., 19: , M. Bencora, S. Hamani, J.J. Nieto, Te metod of upper and lower solutions for second order differential inclusions wit integral boundary conditions, Rocky Mt. J. Mat., 40:13 26, T. Jankowski, Differential equations wit integral boundary conditions, J. Comput. Appl. Mat., 147:1 8, S. Liu, M. Jia, Y. Tian, Existence of positive solutions for boundary-value problems wit integral boundary conditions and sign canging nonlinearities, Electron. J. Differ. Equ., ):1 12, S. Pečiulytė, O. Štikonienė, A. Štikonas, Sturm Liouville problem for stationary differential operator wit nonlocal integral boundary condition, Mat. Model. Anal., 10: , M. Sapagovas, R. Čiupaila, Ž. Jokšienė, T. Meškauskas, Computational experiment for stability analysis of difference scemes wit nonlocal conditions, Informatica, 24: , M. Sapagovas, K. Jakubėlienė, Alternating direction metod for two-dimensional parabolic equation wit nonlocal integral condition, Nonlinear Anal. Model. Control, 17:91 98, B. Amad, S.K. Ntouyas, Nonlinear fractional differential equations and inclusions of arbitrary order and multi-strip boundary condtions, Electron. J. Differ. Equ., ):1 22, B. Amad, S.K. Ntouyas, A study of nonlinear fractional differential equations of arbitrary order wit Riemann Liouville type multi-strip boundary conditions, Mat. Probl. Eng., 2013, Article ID , 9 pp., B. Amad, S.K. Ntouyas, Existence results for iger order fractional differential inclusions wit multi-strip fractional integral boundary conditions, Electron. J. Qual. Teory Differ. Equ., ):1 19, G. Infante, J.R.L. Webb, Nonlinear nonlocal boundary value problems and perturbed Hammerstein integral equations, Proc. Edinb. Mat. Soc., 49: , G.L. Karakostas, P.C. Tsamatos, Existence of multiple positive solutions for a nonlocal boundary value problem, Topol. Metods Nonlinear Anal., 19: , G.L. Karakostas, P.C. Tsamatos, Multiple positive solutions of some Fredolm integral equations arisen from nonlocal boundary value problem, Electron. J. Differ. Equ., ):1 17, J.R.L. Webb, G. Infante, Positive solutions of nonlocal boundary value problems: A unified approac, J. Lond. Mat. Soc., 74: ,

16 Smootness of solutions wit respect to multi-strip integral BC J.R.L. Webb, Existence of positive solutions for a termostat model, Nonlinear Anal., Real World Appl., 13: , W.M. Wyburn, Differential equations wit general boundary conditions, Bull. Am. Mat. Soc., 48: , A. Datta, Differences wit respect to boundary points for rigt focal boundary conditions, J. Differ. Equ. Appl., 4: , J. Eme, Differentiation of solutions of boundary value problems wit respect to nonlinear boundary conditions, J. Differ. Equations, 101: , J. Eme, P.W. Eloe, J. Henderson, Differentiability wit respect to boundary conditions and deviating argument for functional-differential systems, Differ. Equ. Dyn. Syst., 1:59 71, J. Eme, J. Henderson, Differentiation of solutions of boundary value problems wit respect to boundary conditions, Appl. Anal., 46: , J. Eme, J. Henderson, Functional boundary value problems and smootness of solutions, Nonlinear Anal., Teory Metods Appl., 26: , J. Eme, B. Lawrence, Linearized problems and continuous dependence for finite difference equations, Panam. Mat. J., 10:13 24, J. Henderson, Rigt focal point boundary value problems for ordinary differential equation and variational equations, J. Mat. Anal. Appl., 98: , J. Henderson, Disconjugacy, disfocality and differentiation wit respect to boundary conditions, J. Mat. Anal. Appl., 121:1 9, J. Henderson, B. Lawrence, Smoot dependence on boundary matrices, J. Differ. Equ. Appl., 2: , B. Lawrence, A variety of differentiability results for a multi-point boundary value problem, J. Comput. Appl. Mat., 141: , A.C. Peterson, Comparison teorems and existence teorems for ordinary differential equations, J. Mat. Anal. Appl., 55: , A.C. Peterson, Existence-uniqueness for ordinary differential equations, J. Mat. Anal. Appl., 64: , A.C. Peterson, Existence-uniqueness for focal point boundary value problems, SIAM J. Mat. Anal., 12: , A.C. Peterson, Existence and uniqueness teorems for nonlinear difference equations, J. Mat. Anal. Appl., 125: , J. Spencer, Relations between boundary value functions for a nonlinear differential equation and its variational equation, Can. Mat. Bull., 18: , D. Sukup, On te existence of solutions to multipoint boundary value problems, Rocky Mt. J. Mat., 6: , J. Erke, J. Henderson, C. Kunkel, Q. Seng, Boundary data smootness for solutions of nonlocal boundary value problems for second order differential equations, J. Mat. Anal. Appl., 333: , Nonlinear Anal. Model. Control, 2014, Vol. 19, No. 3,

17 412 J. Henderson 36. J. Henderson, B. Karna, C.C. Tisdell, Uniqueness implies existence for multipoint boundary value problems second order equations, Proc. Am. Mat. Soc., 133: , J. Henderson, J. Lyons, Caracterization of partial derivatives wit respect to boundary conditions for nonlocal boundary value problems for n-t order differential equations, Int. J. Pure Appl. Mat., 56: , J. Henderson, C.C. Tisdell, Boundary data smootness for solutions of tree point boundary value problems for second order ordinary differential equations, Z. Anal. Anwend., 44: , B. Hopkins, E. Kim, J. Lyons, B. Speer, Boundary data smootness for solutions of nonlocal boundary value problems for second order difference equations, Commun. Appl. Nonlinear Anal., 16:1 12,

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra

More information

Function Composition and Chain Rules

Function Composition and Chain Rules Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function

More information

Math 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0

Math 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0 3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,

More information

Poisson Equation in Sobolev Spaces

Poisson Equation in Sobolev Spaces Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on

More information

Click here to see an animation of the derivative

Click here to see an animation of the derivative Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,

More information

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx. Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions

More information

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example, NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing

More information

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations

Parameter 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 information

Overlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions

Overlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions Proc. Indian Acad. Sci. (Mat. Sci.) Vol. 121, No. 4, November 2011, pp. 481 493. c Indian Academy of Sciences Overlapping domain decomposition metods for elliptic quasi-variational inequalities related

More information

1. Introduction. We consider the model problem: seeking an unknown function u satisfying

1. Introduction. We consider the model problem: seeking an unknown function u satisfying A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced

More information

arxiv: v1 [math.dg] 4 Feb 2015

arxiv: v1 [math.dg] 4 Feb 2015 CENTROID OF TRIANGLES ASSOCIATED WITH A CURVE arxiv:1502.01205v1 [mat.dg] 4 Feb 2015 Dong-Soo Kim and Dong Seo Kim Abstract. Arcimedes sowed tat te area between a parabola and any cord AB on te parabola

More information

ch (for some fixed positive number c) reaching c

ch (for some fixed positive number c) reaching c GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan

More information

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems 5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we

More information

Key words. Finite element method; convection-diffusion-reaction; nonnegativity; boundedness

Key words. Finite element method; convection-diffusion-reaction; nonnegativity; boundedness PRESERVING NONNEGATIVITY OF AN AFFINE FINITE ELEMENT APPROXIMATION FOR A CONVECTION-DIFFUSION-REACTION PROBLEM JAVIER RUIZ-RAMÍREZ Abstract. An affine finite element sceme approximation of a time dependent

More information

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225 THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:

More information

lecture 26: Richardson extrapolation

lecture 26: Richardson extrapolation 43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)

More information

THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein

THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul

More information

Influence of the Stepsize on Hyers Ulam Stability of First-Order Homogeneous Linear Difference Equations

Influence of the Stepsize on Hyers Ulam Stability of First-Order Homogeneous Linear Difference Equations International Journal of Difference Equations ISSN 0973-6069, Volume 12, Number 2, pp. 281 302 (2017) ttp://campus.mst.edu/ijde Influence of te Stepsize on Hyers Ulam Stability of First-Order Homogeneous

More information

Reflection Symmetries of q-bernoulli Polynomials

Reflection Symmetries of q-bernoulli Polynomials Journal of Nonlinear Matematical Pysics Volume 1, Supplement 1 005, 41 4 Birtday Issue Reflection Symmetries of q-bernoulli Polynomials Boris A KUPERSHMIDT Te University of Tennessee Space Institute Tullaoma,

More information

1 Lecture 13: The derivative as a function.

1 Lecture 13: The derivative as a function. 1 Lecture 13: Te erivative as a function. 1.1 Outline Definition of te erivative as a function. efinitions of ifferentiability. Power rule, erivative te exponential function Derivative of a sum an a multiple

More information

More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms

More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,

More information

Order of Accuracy. ũ h u Ch p, (1)

Order of Accuracy. ũ h u Ch p, (1) Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical

More information

Dedicated to the 70th birthday of Professor Lin Qun

Dedicated to the 70th birthday of Professor Lin Qun Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,

More information

Topics in Generalized Differentiation

Topics in Generalized Differentiation Topics in Generalized Differentiation J. Marsall As Abstract Te course will be built around tree topics: ) Prove te almost everywere equivalence of te L p n-t symmetric quantum derivative and te L p Peano

More information

Math Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim

Math Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z

More information

A Finite Element Primer

A Finite Element Primer A Finite Element Primer David J. Silvester Scool of Matematics, University of Mancester d.silvester@mancester.ac.uk. Version.3 updated 4 October Contents A Model Diffusion Problem.................... x.

More information

ERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*

ERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER* EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T

More information

Differentiation in higher dimensions

Differentiation in higher dimensions Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends

More information

Symmetry Labeling of Molecular Energies

Symmetry Labeling of Molecular Energies Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry

More information

Section 15.6 Directional Derivatives and the Gradient Vector

Section 15.6 Directional Derivatives and the Gradient Vector Section 15.6 Directional Derivatives and te Gradient Vector Finding rates of cange in different directions Recall tat wen we first started considering derivatives of functions of more tan one variable,

More information

Polynomial Interpolation

Polynomial Interpolation Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc

More information

Numerical Solution of One Dimensional Nonlinear Longitudinal Oscillations in a Class of Generalized Functions

Numerical Solution of One Dimensional Nonlinear Longitudinal Oscillations in a Class of Generalized Functions Proc. of te 8t WSEAS Int. Conf. on Matematical Metods and Computational Tecniques in Electrical Engineering, Bucarest, October 16-17, 2006 219 Numerical Solution of One Dimensional Nonlinear Longitudinal

More information

Some Results on the Growth Analysis of Entire Functions Using their Maximum Terms and Relative L -orders

Some Results on the Growth Analysis of Entire Functions Using their Maximum Terms and Relative L -orders Journal of Matematical Extension Vol. 10, No. 2, (2016), 59-73 ISSN: 1735-8299 URL: ttp://www.ijmex.com Some Results on te Growt Analysis of Entire Functions Using teir Maximum Terms and Relative L -orders

More information

An approximation method using approximate approximations

An approximation method using approximate approximations Applicable Analysis: An International Journal Vol. 00, No. 00, September 2005, 1 13 An approximation metod using approximate approximations FRANK MÜLLER and WERNER VARNHORN, University of Kassel, Germany,

More information

How to Find the Derivative of a Function: Calculus 1

How to Find the Derivative of a Function: Calculus 1 Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te

More information

LECTURE 14 NUMERICAL INTEGRATION. Find

LECTURE 14 NUMERICAL INTEGRATION. Find LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use

More information

Stationary Gaussian Markov Processes As Limits of Stationary Autoregressive Time Series

Stationary Gaussian Markov Processes As Limits of Stationary Autoregressive Time Series Stationary Gaussian Markov Processes As Limits of Stationary Autoregressive Time Series Lawrence D. Brown, Pilip A. Ernst, Larry Sepp, and Robert Wolpert August 27, 2015 Abstract We consider te class,

More information

arxiv:math/ v1 [math.ca] 1 Oct 2003

arxiv:math/ v1 [math.ca] 1 Oct 2003 arxiv:mat/0310017v1 [mat.ca] 1 Oct 2003 Cange of Variable for Multi-dimensional Integral 4 Marc 2003 Isidore Fleiscer Abstract Te cange of variable teorem is proved under te sole ypotesis of differentiability

More information

Research Article New Results on Multiple Solutions for Nth-Order Fuzzy Differential Equations under Generalized Differentiability

Research Article New Results on Multiple Solutions for Nth-Order Fuzzy Differential Equations under Generalized Differentiability Hindawi Publising Corporation Boundary Value Problems Volume 009, Article ID 395714, 13 pages doi:10.1155/009/395714 Researc Article New Results on Multiple Solutions for Nt-Order Fuzzy Differential Equations

More information

CELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION SCHEME WITHOUT FICTITIOUS DOMAINS

CELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION SCHEME WITHOUT FICTITIOUS DOMAINS INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 7, Number 1, Pages 1 9 c 010 Institute for Scientific Computing and Information CELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION

More information

(4.2) -Richardson Extrapolation

(4.2) -Richardson Extrapolation (.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as

More information

On convexity of polynomial paths and generalized majorizations

On convexity of polynomial paths and generalized majorizations On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal

More information

, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim

, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim Mat 132 Differentiation Formulas Stewart 2.3 So far, we ave seen ow various real-world problems rate of cange and geometric problems tangent lines lead to derivatives. In tis section, we will see ow to

More information

Finite Difference Method

Finite Difference Method Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According

More information

CONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL SINGULARLY PERTURBED DIFFERENTIAL EQUATIONS ON EQUIDISTRIBUTING MESHES

CONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL SINGULARLY PERTURBED DIFFERENTIAL EQUATIONS ON EQUIDISTRIBUTING MESHES INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Number, Pages 57 74 c 5 Institute for Scientific Computing and Information CONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL

More information

Function Composition and Chain Rules

Function Composition and Chain Rules Function Composition an Cain Rules James K. Peterson Department of Biological Sciences an Department of Matematical Sciences Clemson University November 2, 2018 Outline Function Composition an Continuity

More information

Global Existence of Classical Solutions for a Class Nonlinear Parabolic Equations

Global Existence of Classical Solutions for a Class Nonlinear Parabolic Equations Global Journal of Science Frontier Researc Matematics and Decision Sciences Volume 12 Issue 8 Version 1.0 Type : Double Blind Peer Reviewed International Researc Journal Publiser: Global Journals Inc.

More information

Finding and Using Derivative The shortcuts

Finding and Using Derivative The shortcuts Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex

More information

Boundary value problem with integral condition for a Blasius type equation

Boundary value problem with integral condition for a Blasius type equation ISSN 1392-5113 Nonlinear Analysis: Modelling and Control, 216, Vol. 21, No. 1, 114 12 http://dx.doi.org/1.15388/na.216.1.8 Boundary value problem with integral condition for a Blasius type equation Sergey

More information

AN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS

AN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES

More information

MATH745 Fall MATH745 Fall

MATH745 Fall MATH745 Fall MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext

More information

Polynomial Interpolation

Polynomial Interpolation Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x

More information

LIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT

LIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as

More information

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note

More information

University Mathematics 2

University Mathematics 2 University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at

More information

Kernel Density Estimation

Kernel Density Estimation Kernel Density Estimation Univariate Density Estimation Suppose tat we ave a random sample of data X 1,..., X n from an unknown continuous distribution wit probability density function (pdf) f(x) and cumulative

More information

New Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems

New Fourth Order Quartic Spline Method for Solving Second Order Boundary Value Problems MATEMATIKA, 2015, Volume 31, Number 2, 149 157 c UTM Centre for Industrial Applied Matematics New Fourt Order Quartic Spline Metod for Solving Second Order Boundary Value Problems 1 Osama Ala yed, 2 Te

More information

Math 312 Lecture Notes Modeling

Math 312 Lecture Notes Modeling Mat 3 Lecture Notes Modeling Warren Weckesser Department of Matematics Colgate University 5 7 January 006 Classifying Matematical Models An Example We consider te following scenario. During a storm, a

More information

A method for solving first order fuzzy differential equation

A method for solving first order fuzzy differential equation Available online at ttp://ijim.srbiau.ac.ir/ Int. J. Industrial Matematics (ISSN 2008-5621) Vol. 5, No. 3, 2013 Article ID IJIM-00250, 7 pages Researc Article A metod for solving first order fuzzy differential

More information

POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY

POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat

More information

Bob Brown Math 251 Calculus 1 Chapter 3, Section 1 Completed 1 CCBC Dundalk

Bob Brown Math 251 Calculus 1 Chapter 3, Section 1 Completed 1 CCBC Dundalk Bob Brown Mat 251 Calculus 1 Capter 3, Section 1 Completed 1 Te Tangent Line Problem Te idea of a tangent line first arises in geometry in te context of a circle. But before we jump into a discussion of

More information

Analytic Functions. Differentiable Functions of a Complex Variable

Analytic Functions. Differentiable Functions of a Complex Variable Analytic Functions Differentiable Functions of a Complex Variable In tis capter, we sall generalize te ideas for polynomials power series of a complex variable we developed in te previous capter to general

More information

AMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.

AMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error. Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error

More information

Brazilian Journal of Physics, vol. 29, no. 1, March, Ensemble and their Parameter Dierentiation. A. K. Rajagopal. Naval Research Laboratory,

Brazilian Journal of Physics, vol. 29, no. 1, March, Ensemble and their Parameter Dierentiation. A. K. Rajagopal. Naval Research Laboratory, Brazilian Journal of Pysics, vol. 29, no. 1, Marc, 1999 61 Fractional Powers of Operators of sallis Ensemble and teir Parameter Dierentiation A. K. Rajagopal Naval Researc Laboratory, Wasington D. C. 2375-532,

More information

Mass Lumping for Constant Density Acoustics

Mass Lumping for Constant Density Acoustics Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial

More information

The Laplace equation, cylindrically or spherically symmetric case

The Laplace equation, cylindrically or spherically symmetric case Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,

More information

The Sinc-Collocation Method for Solving the Telegraph Equation

The Sinc-Collocation Method for Solving the Telegraph Equation Te Sinc-Collocation Metod for Solving te Telegrap Equation E. Hesameddini *1, E. Asadolaifard Department of Matematics, Faculty of Basic Sciences, Siraz University of Tecnology, Siraz, Iran *1 esameddini@sutec.ac.ir;

More information

An L p di erentiable non-di erentiable function

An L p di erentiable non-di erentiable function An L di erentiable non-di erentiable function J. Marsall As Abstract. Tere is a a set E of ositive Lebesgue measure and a function nowere di erentiable on E wic is di erentible in te L sense for every

More information

Numerical Differentiation

Numerical Differentiation Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function

More information

3.4 Worksheet: Proof of the Chain Rule NAME

3.4 Worksheet: Proof of the Chain Rule NAME Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are

More information

Stability properties of a family of chock capturing methods for hyperbolic conservation laws

Stability properties of a family of chock capturing methods for hyperbolic conservation laws Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation

More information

Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS

Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Opuscula Matematica Vol. 26 No. 3 26 Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Abstract. In tis work a new numerical metod is constructed for time-integrating

More information

MANY scientific and engineering problems can be

MANY scientific and engineering problems can be A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial

More information

Section 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is

Section 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit

More information

On a continued fraction formula of Wall

On a continued fraction formula of Wall c Te Ramanujan Journal,, 7 () Kluwer Academic Publisers, Boston. Manufactured in Te Neterlands. On a continued fraction formula of Wall DONGSU KIM * Department of Matematics, KAIST, Taejon 305-70, Korea

More information

A finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations

A finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations ANZIAM J. 54 (CTAC2012) pp.c681 C698, 2013 C681 A finite element approximation for te quasi-static Maxwell Landau Lifsitz Gilbert equations Kim-Ngan Le 1 T. Tran 2 (Received 31 October 2012; revised 29

More information

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator. Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative

More information

Parametric Spline Method for Solving Bratu s Problem

Parametric 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 information

FINITE DIFFERENCE APPROXIMATIONS FOR NONLINEAR FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS

FINITE DIFFERENCE APPROXIMATIONS FOR NONLINEAR FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS UNIVERSITATIS IAGELLONICAE ACTA MATHEMATICA, FASCICULUS XL 2002 FINITE DIFFERENCE APPROXIMATIONS FOR NONLINEAR FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS by Anna Baranowska Zdzis law Kamont Abstract. Classical

More information

Preface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.

Preface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed. Preface Here are my online notes for my course tat I teac ere at Lamar University. Despite te fact tat tese are my class notes, tey sould be accessible to anyone wanting to learn or needing a refreser

More information

Combining functions: algebraic methods

Combining functions: algebraic methods Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)

More information

Math 161 (33) - Final exam

Math 161 (33) - Final exam Name: Id #: Mat 161 (33) - Final exam Fall Quarter 2015 Wednesday December 9, 2015-10:30am to 12:30am Instructions: Prob. Points Score possible 1 25 2 25 3 25 4 25 TOTAL 75 (BEST 3) Read eac problem carefully.

More information

THE IMPLICIT FUNCTION THEOREM

THE IMPLICIT FUNCTION THEOREM THE IMPLICIT FUNCTION THEOREM ALEXANDRU ALEMAN 1. Motivation and statement We want to understand a general situation wic occurs in almost any area wic uses matematics. Suppose we are given number of equations

More information

Continuity. Example 1

Continuity. Example 1 Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)

More information

(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?

(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a? Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is

More information

1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist

1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter

More information

Robotic manipulation project

Robotic manipulation project Robotic manipulation project Bin Nguyen December 5, 2006 Abstract Tis is te draft report for Robotic Manipulation s class project. Te cosen project aims to understand and implement Kevin Egan s non-convex

More information

Math 31A Discussion Notes Week 4 October 20 and October 22, 2015

Math 31A Discussion Notes Week 4 October 20 and October 22, 2015 Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes

More information

Error estimates for a semi-implicit fully discrete finite element scheme for the mean curvature flow of graphs

Error estimates for a semi-implicit fully discrete finite element scheme for the mean curvature flow of graphs Interfaces and Free Boundaries 2, 2000 34 359 Error estimates for a semi-implicit fully discrete finite element sceme for te mean curvature flow of graps KLAUS DECKELNICK Scool of Matematical Sciences,

More information

FEM solution of the ψ-ω equations with explicit viscous diffusion 1

FEM solution of the ψ-ω equations with explicit viscous diffusion 1 FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible

More information

= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)

= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c) Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''

More information

The Column and Row Hilbert Operator Spaces

The Column and Row Hilbert Operator Spaces Te Column and Row Hilbert Operator Spaces Roy M Araiza Department of Matematics Purdue University Abstract Given a Hilbert space H we present te construction and some properties of te column and row Hilbert

More information

Recent Progress in the Integration of Poisson Systems via the Mid Point Rule and Runge Kutta Algorithm

Recent Progress in the Integration of Poisson Systems via the Mid Point Rule and Runge Kutta Algorithm Recent Progress in te Integration of Poisson Systems via te Mid Point Rule and Runge Kutta Algoritm Klaus Bucner, Mircea Craioveanu and Mircea Puta Abstract Some recent progress in te integration of Poisson

More information

4.2 - Richardson Extrapolation

4.2 - Richardson Extrapolation . - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence

More information

MA455 Manifolds Solutions 1 May 2008

MA455 Manifolds Solutions 1 May 2008 MA455 Manifolds Solutions 1 May 2008 1. (i) Given real numbers a < b, find a diffeomorpism (a, b) R. Solution: For example first map (a, b) to (0, π/2) and ten map (0, π/2) diffeomorpically to R using

More information

A = h w (1) Error Analysis Physics 141

A = h w (1) Error Analysis Physics 141 Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.

More information

Taylor Series and the Mean Value Theorem of Derivatives

Taylor Series and the Mean Value Theorem of Derivatives 1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential

More information

REVIEW LAB ANSWER KEY

REVIEW LAB ANSWER KEY REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g

More information

Math 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006

Math 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2

More information