Boundary-value problems

Size: px
Start display at page:

Download "Boundary-value problems"

Transcription

1 226

2 Chpter 10 Boundry-vlue problems The initil-vlue problem is chrcterized by the imposition of uxiliry dt t single point: if the eqution is of the nth order, the n otherwise rbitrry constnts in its solution re determined by prescribing the function nd its first n 1 derivtives ll t single vlue of the independent vrible. Boundry-vlue problems provide the uxiliry dt by providing informtion, not t single point, but t two points 1, usully t the endpoints (or boundry) of n intervl. They rise often in the solution of prtil-differentil equtions. Exmple Consider the temperture u (x, t) in wll of thickness L whose interior, t x = 0, is mintined t temperture A nd whose exterior, t x = L, is mintined t temperture B. Initilly, t t = 0, the distribution of temperture is u = u 0 (x). Its subsequent evolution is governed by the prtil-differentil eqution (the het eqution) ρ (x) u t = x ( p (x) u x ). (10.1) Here ρ nd p re known functions depending on the construction of the wll, nd both my be ssumed positive. The problem of finding solution of eqution (10.1) together with the imposed initil dt nd boundry dt is n exmple of n initil-boundry-vlue problem. There is simple, time-independent solution of this eqution stisfying the boundry conditions: x 0 U (x) = A + (B A) p (s) 1 ds L 0 p (s) 1 ds. (10.2) 1 Or possibly more thn two, but the commonly rising cse is tht of two points. 227

3 228 If we write u = U (x) + v in eqution (10.1), we find tht v stisfies exctly the sme eqution (10.1), but its boundry conditions now tke the form v = 0 t x = 0 nd t x = L. Its initil vlue is v = v 0 (x) = u 0 (x) U (x). The method of seprtion of vribles is n ttempt to find solution of eqution (10.1) for v, stisfying the boundry conditions, in the form v (x, t) = X (x) T (t) of product of function of x only with function of t only. Substitution of this expression into eqution (10.1) results in the eqution ρ (x) X (x) T (t) = T (t) (p (x) X (x)). If one further ssumes tht the functions X nd T re not zero except possibly t isolted points, then by dividing either side of this lst eqution by ρxt one obtins T (t) T (t) = (p (x) X (x)) ρ (x) X (x). Since the left-hnd side depends on t only nd the right-hnd side on x only, they must be equl to the sme constnt, λ sy. This requires of these functions tht nd ( d p (x) dx ) λρ (x) X = 0 (10.3) dx dx dt dt = λt, (10.4) where λ is some s yet undetermined constnt. This will provide solution v = X (x) T (t) of the het eqution if eqution (10.3) hs solution stisfying the boundry dt X (0) = 0, X (L) = 0. This combintion of the differentil eqution (10.3) with boundry dt is clled boundry-vlue problem, to distinguish it from the initil-vlue problem. Unlike initil-vlue problems, boundry-vlue problems do not lwys hve solutions, s we shll see below (cf. Exmple ). Even if the boundry-vlue problem consisting of eqution (10.3) does hve solution stisfying the boundry dt, it will not in generl provide solution to the originl initil-boundry-vlue problem, since the initil vlue T (0) X (x) will not in generl gree with v 0 (x). These deficiencies will be delt with below. Eqution (10.3) of the preceding exmple is liner, homogenous nd of the second order. We will mostly confine the discussion to this cse. We begin with certin lgebric properties of liner opertors.

4 Opertors nd their djoints Consider the generl, liner, homogeneous opertor M of the second order defined by the formul Mu f (x) u + g (x) u + h (x) u (10.5) on the intervl (, b). The coefficients f, g, h re ssumed to be C 2 (i.e, continuous together with their first two derivtives). If we multiply M u by n rbitrry C 2 function v, integrte over the intervl, nd use integrtion by prts to move derivtives of u to derivtives of v, we find where vmu dx = { fvu (fv) u + guv } b + (M v) u dx, M v (fv) (gv) + hv = fv + ( 2f g ) v + ( f g + h ) v. (10.6) The opertor M is sid to be djoint to the opertor M. between them cn be summrized s follows: The reltion (vmu um v) dx = { f ( vu uv ) + ( g f ) uv } b. (10.7) The form on the right, which is biliner in the vlues of u, u, v, v t the endpoints, is clled the biliner concomitnt. A cse of specil importnce is tht of the self-djoint opertor, M = M. For this to be true we must hve 2f g = g nd f g + h = h. This in turn holds precisely when g = f ; in this cse the opertor M tkes the form Mu d ( f (x) du ) + h (x) u. (10.8) dx dx This is the form of the opertor ppering on the left-hnd side of eqution (10.3), nd is of very common occurrence in problems of mthemticl physics, which is the origin of its specil importnce. For the self-djoint cse, eqution (10.7) cn be written (vmu umv) dx = { f ( vu uv )} b. (10.9)

5 230 The eqution (10.3) of the preceding exmple hs the form Mu = λρu where the opertor of eqution (10.8) is defined by f = p nd h = 0. The constnt λ tht rose in the course of seprting vribles my pper rbitrry, unrelted to the problem to be solved, but this is not so. It cnnot be ssigned rbitrrily, or the problem of solving eqution (10.3) together with the boundry conditions X (0) = 0 nd X (L) = 0 will hve only the trivil solution X 0. Exmple In the Exmple , tke ρ (x) = p (x) = 1. Then eqution (10.1) tkes the form u t = u xx nd eqution (10.3) tkes the form X = λx. Solutions my be written X = c sin ( λx ) + d cos ( λx ), where we hve ssumed tht λ 0; this will be verified subsequently. Imposition of the boundry condition X (0) = 0 then implies tht d = 0, so imposition of the remining condition t x = L implies c sin ( L λ ) = 0. If c = 0, the solution X 0. To hve nonzero solution we must hve L λ = π, 2π,..., i.e., λ cnnot be rbitry but must tke one or nother of the discrete set of vlues λ n = n 2 π 2 /L 2. There re therefore infinitely mny seprtion solutions of the form v n (x, t) = e n2 π 2 t/l 2 sin (nπx/l). (10.10) Ech stisfies the eqution nd the boundry conditions. Since the eqution nd the boundry conditions re liner nd homogeneous, ny liner combintion with rbitrry coefficients {c n }, v (x, t) = c n v n (x, t), lso stisfies the eqution nd the boundry conditions. If this sum cn be mde to stisfy the initil condition v = v 0 (x), complete solution will be chieved. The ltter condition is written v 0 (x) = c n sin (nπx/l). If the summtion is llowed to run from 1 to, this represents the expnsion of the function v 0 in Fourier sine series. Under pproprite conditions on the function v 0, this is indeed possible by choosing the coefficients {c n } to be the Fourier coefficients of v 0, lthough the convergence of the series for v (x, t) requires further considertion.

6 231 In the preceding exmple, the vlues of the prmeter λ for which nontrivil solutions exist re clled eigenvlues, nd the corresponding functions sin ( λx ) re clled eigenfunctions. Eigenvlues nd their eigenfunctions ply centrl role in boundry-vlue problems. We therefore consider them in greter generlity Eigenvlue problems Consider the self-djoint differentil opertor on the intervl [, b], which we now write in the form Mu ( p (x) u ) + q (x) u = λρ (x) u, (10.11) where the coefficients stisfy the following conditions: p C 1 nd p > 0 on [, b]; q, ρ C on [, b] nd ρ > 0 on (, b). (10.12) This form of the opertor is generl enough to pply to problems rising in diverse pplictions. Tht the coefficient p is continuous nd positive on the closed intervl [, b] implies tht it is bounded wy from zero (i.e., tht p (x) δ for some δ > 0 for ech x [, b]), nd therefore the equtions Mu λρu = 0 nd Mu λρu = r cn be put into the stndrd form to which the theory of Chpter 2 pplies. The minus sign hs been introduced for lter convenience : pplied to Exmple , the new sign convention of eqution (10.11) would mke the eigenvlues turn out positive rther thn negtive. The integrl reltion (10.9) now tkes the form (vmu umv) dx = { p ( uv vu )} b. (10.13) The biliner concomitnt hs been reduced to the product of the Wronskin of the two functions multiplied by the coefficient p. The eqution (10.11) must be supplied with boundry conditions. These we will tke to be homogeneous boundry conditions generlizing those of Exmples nd The generl structure of homogeneous boundry conditions is: 11 u () + 12 u () + b 11 u (b) + b 12 u (b) = 0, nd 21 u () + 22 u () + b 21 u (b) + b 22 u (b) = 0. (10.14) Here the numbers { ij, b ij }, i = 1, 2, j = 1, 2 re prescribed: they determine the boundry dt. Rrely is one clled on to consider boundry dt

7 232 of such generlity. In the exmples bove, 11 nd b 21 re nonzero, nd ll the others vnish. But before pssing on to more specific cses, we mke couple of observtions tht follow from the liner nd homogeneous nture of both the eqution nd the boundry conditions. One is tht u 0 is lwys solution of eqution (10.11) together with the uxiliry dt (10.14). Another is tht if u is solution so lso is cu for ny rbitrry constnt c nd, more generlly, if u nd v re solutions so lso is ny liner combintion cu + dv with rbitrry constnts c nd d. Given n opertor M, we know how to construct the djoint opertor M, s in eqution (10.6) bove. We my consider, long with the boundryvlue problem Mu = λρu, < x < b nd conditions (10.15) the relted boundry-vlue problem M v = µρv, < x < b nd conditions 10.17, (10.16) where the boundry conditions for this relted problem, 11v () + 12v () + b 11v (b) + b 12v (b) = 0, nd 21v () + 22v () + b 21v (b) + b 22v (b) = 0, (10.17) re similr in structure { to } equtions (10.14) but with possibly different vlues for the coefficients ij. Suppose now tht we cn choose these coefficients such tht, whenever v stisfies the conditions (10.17) while u stisfies the conditions (10.14), the biliner concomitnt vnishes, i.e., by eqution (10.7) bove, vmu dx = um v dx, (10.18) for ll C 2 functions u, v stisfying conditions (10.14), (10.17), respectively. In this cse, the eigenvlue problem (10.16) will be referred to s the djoint problem to the eigenvlue problem (10.14). Tht the djoint boundry conditions cn in fct be found such tht the biliner comcomitnt vnishes my be seen s follows. Let the four-component vector U be defined by the equtions U = (u(), u (), u(b), u (b)) = (u 1, u 2, u 3, u 4 ),

8 233 nd similrly for V : V = (v(), v (), v(b), v (b)) = (v 1, v 2, v 3, v 4 ). The pir of boundry conditions (10.14) my be written A 1 U = 0 nd A 2 U = 0, (10.19) where A 1 nd A 2 re ppropritely chosen four-component vectors. These must be linerly independent if the two equtions (10.14) re to represent two distinct boundry conditions, nd we my therefore, without loss of generlity, suppose they re orthogonl: A 1 A 2 = 0. We cn extend these to n orthogonl bsis A 1, A 2, A 3, A 4 for R 4, nd observe tht U = αa 3 + βa 4 by virtue of the boundry conditions (10.19). The biliner concomitnt, the right-hnd side of eqution (10.7), is biliner in U nd V. If we cll it C, we my write C = V KU where K is four-by-four mtrix whose entries cn be red off eqution (10.7). Therefore C = α (V KA 3 ) + β (V KA 4 ). The boundry conditions on u re stisfied for ny choices of α nd β so the vnishing of the biliner concomitnt C requires tht V KA 3 = 0 nd V KA 4 = 0. This implies there re t lest two linerly independent choices of the vector V tht result in the vnishing of the biliner comcomitnt, proving the clim. If the mtrix K is nonsingulr, there re precisely two, since the vectors KA 3 nd KA 4 re then linerly independent: reding K from eqution (10.7) shows tht it is indeed nonsingulr, so there re precisely two. We shll, beginning in the next section, restrict considertion to the cse when one of the boundry conditions refers to the left-hnd endpoint only, the other to the right-hnd endpoint only, but we provide here simple exmples of boundry conditions tht ech involve both endpoints. Exmple Boundry-vlue problems re by no mens confined to equtions of the second order. A simple exmple is the first-order eqution u = λu on the intervl [0, 2π] with the single endpoint condition u(0) = 2u(2π). If the opertor L is defined by Lu = u then the djoint opertor is

9 234 defined by L v = v. It is esy to see tht the djoint boundry conditions is v(0) = (1/2)v(2π). The solution is u = u 0 exp (λx) nd the eigenvlue reltion is exp (2πλ) = 1/2. The permissible vlues of λ re therefore (ln 2/2π)+ik, where k tkes on integer vlues. Exmple Consider the eqution u + λu = 0 on the intervl [0, 2π], with the boundry vlues u (0) u (2π) = 0 nd u (0) u (2π) = 0. This problem is self-djoint nd the djoint boundry conditions re the sme s those bove for u. The following ssertions cn be verified explicitly: The solution is identiclly zero unless λ is equl to one or nother of the eigenvlues 0, 1, 4,... k 2,.... If so, the corresponding solutions (the eigenfunctions) re 1; cos x, sin x; cos 2x, sin 2x;... ; cos kx, sin kx;..., where to ech eigenvlue k 2 except for k = 0 there correspond the pir of eigenfunctions cos (kx), sin (kx). By n eigenfunction we men nontrivil solution of homogeneous boundry-vlue problem. This does not exist for generl vlues of the prmeter λ: it requires tht λ hve one of specil set of vlues, i.e. n eigenvlue. The problem of finding these eigenvlues nd eigenfunctions we refer to s n eigenvlue problem. The definitions of the terms djoint nd self-djoint tht we hve used re used in functionl nlysis in mnner tht is relted but not identicl. There the notion of djointness is introduced in connection with tht of liner functionl on liner spce. A typicl such functionl my be written < v, u > where v lbels the functionl nd u rnges over the liner spce B. If M : B 0 B is liner opertor on the spce B 0, then < v, Mw > for w B 0 is gin liner functionl on B 0, which is denoted by < M v, w >. This defines the djoint mpping M : B B0: < v, Mu >=< M v, u >. The connection of this with our discussion bove is s follows. Suppose tht the functions w(x) we wish to consider re in the spce of C 2 functions on [, b] which stisfy the boundry conditions

10 235 (10.14). This is liner spce; cll it L. Now, for fixed function v(x) s yet unspecified, consider f v [u] = v(x)u(x)dx = (v, u) where u L, the spce of continuous functions on [, b]. This is esily seen to be liner functionl on L, for ny fixed choice of v. Then, noting tht M : L L for the opertor of eqution (10.5), put u = Mw If we further restrict v to be in C 2 nd to stisfy the djoint boundry conditions (10.17), then, with our definitions of M nd M, we derive (v, Mu) = (M v, u). Among the differences in the two pproches is one of generlity: we do not consider ll liner functionls on L, but only those defined by integrtion ginst certin specil clss of functions v The Sturm-Liouville problem We shll emphsize the fmily of boundry-vlue problems in which the endpoint conditions re seprted: i.e., there is one condition t ech endpoint, so the conditions (10.14) tke the form αu () + α u () = 0 nd βu (b) + β u (b) = 0, (10.20) where it is ssumed tht the coefficients α nd α re not both zero, nd likewise for β, β. The ssumption (10.20) excludes certin cses for which the endpoint conditions re not seprted but re mixed, s in the more generl conditions (10.14), or Exmples nd bove. The Sturm-Liouville problem is the eigenvlue problem consisting of the self-djoint eqution (10.11), under conditions on the coefficients given by eqution (10.12), together with the seprted endpoint conditions (10.20). It is the bsic eigenvlue problem of ordinry differentil equtions, nd tkes its nme from two of the mthemticins who developed its theory. We shll study the theory in greter depth in Chpter 11. In this section we estblish some of its bsic lgebric properties. A simple but useful observtion is the following. Lemm Suppose two functions u nd v both stisfy the endpoint condition (10.20) t x =. Then their Wronskin vnishes there.

11 236 Proof: The ssumption mens tht αu () + α u () = 0, αv () + α v () = 0. This is system of two liner equtions for α nd α. The determinnt of this system is the Wronskin of the functions evluted t x =. If it were nonzero, we would rrive t the contrdiction tht α nd α both vnish. Remrks: A similr conclusion holds for the other endpoint x = b. If the functions u nd v re llowed to be complex-vlued functions of the rel vrible x, this conclusion is unchnged. We now esily find: Theorem For the eigenvlue problem (10.11), (10.20) with seprted endpoint conditions, there is t most one linerly independent eigenfunction corresponding to single vlue of λ. Proof: Suppose there were two linerly independent solutions u nd v. They would both hve to stisfy the endpoint condition t x =, so their Wronskin would hve to vnish there. But they re both solutions of the eqution (10.11) nd if they were linerly independent their Wronskin could not vnish nywhere, hence they cnnot be linerly independent. Remrks: Exmple bove, for which the endpoint conditions re not seprted, demonstrtes tht there my be two linerly independent eigenfunctions corresponding to single eigenvlue in tht cse. We infer tht for n eigenvlue problem with seprted endpoint conditions, ny eigenfunction is determined uniquely up to multipliction by n rbitrry constnt. Exmple shows tht for some boundry-vlue problems, even though the coefficients in the eqution nd in the boundry conditions re rel, the eigenvlues nd eigenfunctions my be complex. We shll show first tht, for the Sturm-Liouville boundry-vlue problem, the eigenvlues re rel. To estblish this we must llow for the possibility tht they be complex. We continue to restrict considertion to rel

12 237 vlues of the independent vrible x, nd to rel vlues of the coefficient functions p, q, ρ nd of the prmeters α, α nd β, β of the boundry conditions. If λ is complex, the solution u of eqution (10.11) is necessrily lso complex: u (x) = u 1 (x) + iu 2 (x) where u 1 nd u 2 re rel. If we write u = u 1 iu 2 for the complex conjugte, we find Mu = λρu, wheres if u stisfies the boundry conditions (10.20) so lso does u. We therefore hve the two equtions Mu = λρu, Mu = λρu where both u nd u stisfy the boundry conditions (10.20). Applying the integrl condition (10.13) with v = u we find (umu umu) dx = ( λ λ ) ρ u 2 dx = { p ( uu uu )} b. The lst term vnishes in view of Lemm nd we therefore find ( ) λ λ ρ u 2 dx = 0. If the integrl vnishes, then u 0 on [, b]: the integrnd is non-negtive nd continuous, so if it should be nonzero t ny point, there would be neighborhood of tht point throughout which it is positive, nd this neighborhood would mke positive contribution to the integrl, which would therefore be nonzero. Since u is n eigenfunction, it cnnot vnish identiclly nd we conclude tht the integrl is nonzero, nd therefore λ = λ, i.e., λ is rel. We hve proved Theorem Any eigenvlue of the Sturm-Liouville problem is rel. In investigting eigenvlues, we cn now restrict our ttention to rel vlues of λ. We cn then lso restrict ttention to rel-vlued solutions u, since ny complex-vlued solution is obtined by multiplying rel-vlued solution by complex constnt. The integrl uv dx of the product of the functions u nd v represents n nlog of the inner product (or sclr product, or dot product) of vectors. It is n nlogy tht cn be pushed quite fr. It is therefore nturl to refer to two functions for which this integrl vnishes s orthogonl. Definition The functions u nd v re orthogonl on (, b) with respect to the weight function ρ if u (x) v (x) ρ (x) dx = 0. (10.21)

13 238 Theorem Eigenfunctions of the Sturm-Liouville problem belonging to different eigenvlues re orthogonl with respect to the weight function ρ. Proof: Let u nd v be eigenfunctions belonging to different eigenvlues λ nd µ respectively. Then (vmu umv) dx = (λ µ) uvρ dx = { p ( uv vu )} b = 0. The conclusion now follows since λ µ. In Exmple bove, n importnt element ws the bility to represent n essentilly rbitrry function s n infinite sum of sine functions. More generlly it is likewise importnt to be ble to represent n essentilly rbitrry function f on (, b) in the form f (x) = c n u n (x) n=0 where {u n } is sequence of eigenfunctions of Sturm-Liouville problem. If we ssume such development is possible nd tht we cn interchnge summtion nd integrtion, then the orthogonlity of the eigenfunctions provides the coefficients in this sum: c n = f (x) u n (x) ρ (x) dx u n (x) 2 ρ (x) dx. PROBLEM SET Refer to Exmple Verify tht the function U (x) given there is indeed solution of eqution (10.1), tht if u = U + v then v indeed stisfies (10.1) if u does nd tht X (x) T (t) is indeed solution of (10.1) if X nd T re solutions of equtions (10.3) nd (10.4) respectively. 2. Work out (M ). 3. For the opertor M given by the formul (10.5), find function I (x) such tht IM is self-djoint. Wht do you need to ssume bout the coefficients of M?

14 In Exmple there ppers to be n implicit ssumption tht λ < 0. Assume insted tht λ = µ 2, where µ is rel number. Cn you chieve nontrivil solution of the boundry-vlue problem of tht exmple? 5. Prove the conclusions sserted in Exmple Consider the mixed-endpoint boundry-vlue problem d ( p (t) du ) + q (t) u = λu, u (0) = u (T ), u (0) = u (T ), dt dt (10.22) where the coefficients p nd q re defined nd continuous for ll rel t nd periodic with period T : p (t + T ) = p (t) nd q (t + T ) = q (t); in ddition ssume tht p is positive nd C 1. Show tht ny solution stisfying the given dt is likewise periodic with period T. 7. Suppose λ = λ 1 + iλ 2 nd u (x) = u 1 (x) + iu 2 (x) represent complex eigenvlue nd eigenfunction, respectively, of the Sturm-Liouville problem M u = λρu together with the boundry conditions (10.20). Write this problem s pir of rel differentil equtions for u 1 nd u 2 nd obtin Theorem in this wy. 8. Prove the sttement in the text following Theorem tht if the eigenvlue λ is rel then ny complex eigenfunction is complex multiple of rel eigenfunction. 9. The eigenfunctions of Exmple re u n (x) = sin (nπx/l). Show by explicit integrtion tht these re orthogonl on (0, L) with weight function ρ(x) Consider the fourth-order opertor Mu p 0 u (4) + p 1 u + p 2 u + p 3 u + p 4 u. Assume tht the coefficients p i, i = 0,..., 4, re sufficiently differentible nd tht p 0 does not vnish. Find the djoint opertor M by successive integrtions by prts. Find the conditions on the coefficients for M to be self-djoint, nd re-express the self-djoint opertor in the form ( ) Mu = d2 dx 2 P d2 u dx 2 + d ( Q du ) + Ru, dx dx

15 240 identifying the coefficients P, Q, R in terms of the originl coefficients {p i }. Show tht if the boundry dt re such tht the corresponding biliner concomitnt vnishes, i.e., if (vmu umv) dx = 0, then the eigenvlues of ny such self-djoint opertor re rel. 11. For Exmple the opertor is M = d/dx nd M = M (verify this); M is sid to be skew-djoint. For the third-order liner opertor Mu p 0 u + p 1 u + p 2 u + p 3 u form the djoint opertor by successive integrtions by prts, nd show tht the opertor M cnnot be self-djoint, but cn be skewdjoint under pproprite conditions on the coefficients. Show further tht if the boundry dt re such tht the corresponding biliner concomitnt vnishes, i.e., if (vmu + umv) dx = 0, then the eigenvlues of ny such skew djoint opertor re pure-imginry. 12. Let {P n (x)} 0 be sequence of polynomils on [, b], with P n of degree n, nd suppose they re orthogonl with respect to positive weight function ρ. Show tht P n hs n zeros in [, b] Green s Functions Consider the generl problem of solving the nth-order, inhomogeneous, differentil eqution Lu = r, < x < b (10.23) where Lu p 0 (x) u (n) + p 1 (x) u (n 1) + + p n (x) u nd the coefficients r nd {p i } re given, continuous functions on [, b] nd p 0 does not vnish there. Now, however, we suppose the dditionl conditions 2 This describes the Legendre polynomils, which indeed stisfy differentil eqution of self-djoint form, but the proof should be independent of the theory of differentil equtions

16 241 on the solution u nd its derivtives to be given in the form of homogeneous boundry dt like those of eqution (10.14) rther thn initil dt. In prticulr, we suppose the boundry conditions to hve the form n ij u (j 1) () + b ij u (j 1) (b) = 0, i = 1, 2,..., n. (10.24) j=1 We shll find tht this boundry-vlue problem my fil to hve solution. However, if it does hve solution, it cn be expressed in the form u (x) = G (x, ξ) r (ξ) dξ, (10.25) where the function G (x, ξ) ppering in the integrnd depends on the opertor L nd the boundry dt, but not on the function r. The function G, clled the Green s function for this boundry-vlue problem, therefore plys similr role to tht plyed by the influence function for initil-vlue problems (cf. eqution 2.26). We outline the principl result in It tkes its simplest form in the (importnt) cse of the Sturm-Liouville boundry-vlue problem, so we consider this in detil The Sturm-Liouville Cse Consider now the inhomogeneous eqution Mu λρu = r (10.26) with the sme definition of the opertor M nd the sme seprted boundry conditions s in 10.3; in prticulr, M is given by eqution (10.11). If λ is n eigenvlue for the homogeneous problem nd v the corresponding eigenfunction, then vr dx = v (Mu λρu) dx = { p ( uv vu )} b + u (Mv λρv) dx = 0, since the integrl vnishes becuse of the ssumptions on v nd the integrted prt vnishes becuse u nd v both stisfy the boundry conditions. We conclude tht eqution (10.26) cn fil to hve solution: if λ is n eigenvlue nd r is chosen so tht vr dx 0 (e.g., r = v), there cn be no solution. We then inquire whether there is lwys, i.e., for every continuous function r on [, b], solution of this boundry-vlue problem provided λ is

17 242 not n eigenvlue of the homogeneous problem. Here the nswer will be ffirmtive, nd we cn express the solution in the form u (x) = G (x, ξ) r (ξ) dξ. (10.27) The function G (x, ξ) is clled the Green s function, nd is the counterprt for boundry-vlue problems to the influence function nd the vrition-ofprmeters formul obtined erlier (Chpter 2) in the discussion of liner initil-vlue problems. In fct, our derivtion of the Green s function, to which we now proceed, will be mde to depend on tht formul. Let u 1 (x, λ) nd u 2 (x, λ) be ny bsis of solutions of the homogenous differentil eqution Mu λρu = 0. (10.28) Then there is lwys prticulr integrl of eqution (10.26) which my be written x ( ) u2 (ξ) u 1 (x) u 2 (x) u 1 (ξ) U (x) = r (ξ) dξ, (10.29) p (ξ) W (ξ) which is esily checked on writing the eqution (10.26) in the stndrd form u + p p u + λρ q p = r p. The solutions u 1 nd u 2 depend on the prmeter λ but since λ will not chnge in wht follows, we suppress this dependence in eqution (10.29) nd in other reltions below. Since the most generl solution of eqution (10.26) cn be written in the form u (x) = U (x) + c 1 u 1 (x) + c 2 u 2 (x) (10.30) for some choice of constnts c 1 nd c 2, this must lso be true for the solution of the boundry-vlue problem (10.26) with the boundry conditions (10.20). It is just mtter of clculting c 1 nd c 2. The clcultion is simplified somewht if we choose the bsis u 1, u 2 ppropritely. In prticulr, we ll choose u 1 to stisfy the left-hnd boundry condition nd u 2 to stisfy the right-hnd boundry condition. This cn lwys be done: for exmple, choose u 1 to solve the initil-vlue problem consisting of the differentil eqution (10.28) together with the initil dt u 1 () = α, u 1 () = α, nd similrly for u 2: it cn be chosen to be solution of eqution (10.28) with initil dt u 2 (b) = β, u 2 (b) = β. These

18 243 solutions re linerly independent: if they were linerly dependent we would hve u 1 (x) = Cu 2 (x) where C is constnt, nd then u 1 would stisfy both sets of boundry conditions nd hence be n eigenfunction, contrdicting the ssumption tht λ is not n eigenvlue. It s convenient to introduce shorthnd nottion for the boundry conditions. For ny C 1 function v on [, b] define A[v] = αv () + α v (), B[v] = βv (b) + β v (b). (10.31) In terms of these liner opertors A nd B the bsis functions u 1, u 2 stisfy A[u 1 ] = 0, B[u 2 ] = 0; A[u 2 ] 0, B[u 1 ] 0. Note lso tht our choice of prticulr integrl U is such tht U nd U both vnish t x = nd, s result A[U] = 0. Now pply the boundry conditions A[u] = 0 nd B[u] = 0 to the generl solution (10.30) bove. We find A[u] = c 2 A[u 2 ] = 0 nd B[u] = B[U] + c 1 B[u 1 ] = 0. This implies tht c 2 = 0 nd c 1 = B[U]/B[u 1 ]. But it is strightforwrd clcultion using the explicit form of U in eqution (10.29) tht u 2 (ξ) B[U] = B[u 1 ] r (ξ) dξ. (10.32) p (ξ) W (ξ) This gives for the solution u the formul x u (x) = identifying G (x, ξ) = u 1 (ξ) u 2 (x) b r (ξ) dξ p (ξ) W (ξ) x { u2 (x) u 1 (ξ) /p (ξ) W (ξ) if ξ < x u 2 (ξ) u 1 (x) /p (ξ) W (ξ) if ξ > x u 1 (x) u 2 (ξ) r (ξ) dξ, (10.33) p (ξ) W (ξ) }. (10.34) One cn now verify tht the formul (10.33) stisfies both the eqution (10.26) nd the boundry conditions (10.20). We now hve the following Theorem Suppose λ is not n eigenvlue of the Sturm-Liouville problem (10.11), (10.20). Then, for ny continuous function r on [, b], there is unique solution of equtions (10.26), (10.20).

19 244 Proof: One solution u is provided by the formul (10.33). To see tht it is unique, suppose there were second v. Then the difference u v would stisfy equtions (10.11), (10.20) of the Sturm-Liouville problem. Unless u v is the zero solution, the ltter hs n eigenfunction, i.e., λ is n eigenvlue. This is contrdiction. There is nother chrcteriztion of the Green s function for this problem tht is often useful. Fix ξ in the intervl (, b). Then G is the unique, continuous function stisfying the boundry conditions t ech endpoint, stisfying the differentil eqution MG = 0 t ech point of the intervl except t x = ξ, nd stisfying lso the jump-discontinuity condition [ ] G = 1/p (ξ). x Here the brcket nottion mens, for function f (x), ξ [f] ξ = lim x ξ + f (x) lim x ξ f (x), i.e., the difference of the limiting vlues pproched from the right nd from the left. It is not difficult to check tht the explicit expression (10.34) stisfies these conditions The Generl Cse The ide of the Green s function cn be generlized beyond the Sturm- Liouville boundry-vlue problem, nd indeed beyond the theory of ordinry differentil equtions. In this brief section we merely outline the generliztion to the boundry-vlue problem consisting of the differentil eqution (10.23) together with the liner, homogeneous boundry dt (10.24). Note tht the sign convention for the opertor L of eqution (10.23) differs from tht of the preceding subsection. We suppose tht the coefficient p 0 does not vnish on the intervl [, b]. We suppose further tht the homogeneous problem Lu = 0 together with the boundry dt (10.24) hs only the trivil solution u 0 on [, b]. Under these conditions, one cn define function G (x, ξ) s follows: Definition For fixed ξ in (, b), G stisfies the boundry dt (10.24), it stisfies the differentil eqution Lu = 0 except t x = ξ nd, t x = ξ [ k ] { G 0 if k = 0, 1,..., n 2 x k = 1/p 0 (ξ) if k = n 1. ξ

20 245 This function is the Green s function for the boundry-vlue problem, i.e., the (unique) solution of the ltter is given by the formul u (x) = PROBLEM SET G (x, ξ) r (ξ) dξ. (10.35) 1. Consider the differentil eqution Mu λρu = 0 on [, b], but with the inhomogeneous boundry dt αu () + α u () = A, βu (b) + β u (b) = B. Convert this to n inhomogeneous differentil eqution of the form (10.26) with homogeneous boundry dt. 2. Verify eqution (10.32). 3. Verify tht the formul (10.33) provides solution to equtions (10.26) nd (10.20). 4. Show tht the denomintor p (ξ) W (ξ) in eqution (10.34) is constnt (i.e., independent of ξ). Under wht conditions does it vnish? 5. Find the Green s function for the eqution u + λu = r (x) on the intervl [0, π] with boundry conditions u (0) = 0 nd u (π) = 0. For wht vlues of λ does it become singulr? 6. Show tht the Green s function (10.34) is symmetric: G (x, ξ) = G (ξ, x). 7. Suppose tht λ is n eigenvlue with eigenfunction u 1 (x), nd suppose the condition u 1r dx = 0 is stisfied. Show tht if in eqution (10.33) the function u 1 is this eigenfunction nd u 2 is ny linerly independent solution of the differentil eqution, the inhomogeneous boundry-vlue problem (10.26), (10.20) is stisfied. 8. Verify tht the expression (10.35) solves the boundry-vlue problem (10.23), (10.24); show tht the solution is unique. 9. Consider the liner opertor Lu = u nd the single boundry condition αu() + βu(b) = 0. Under wht conditions on α nd β is it true tht the boundry-vlue problem Lu = 0 hs only the trivil solution? Under this condition, work out the Green s function for the solution of the inhomogeneous problem (10.23) with the sme boundry condition.

21 The sme s problem (9) but with Lu = u + u.

22 314

23 Bibliogrphy [1] G. Birkhoff nd G.-C. Rot. Ordinry Differentil Equtions. John Wiley nd Sons, Inc., New York, [2] P.F. Byrd nd M. D. Friedmn. Hndbook of elliptic integrls for engineers nd physicists. Springer, Berlin, [3] Jck Crr. Applictions of Centre Mnifold Theory. Springer-Verlg, New York, [4] Erl A. Coddington nd Normn Levinson. Theory of Ordinry Differentil Equtions. New York: McGrw-Hill, [5] J. Guckenheimer nd P. Holmes. Nonliner Oscilltions, Dynmicl Systems, nd Bifurction of Vector Fields. Springer-Verlg, New York, [6] Einr Hille. Lectures on Ordinry Differentil Equtions. London: Addison-Wesley Publishing Compny, [7] M. Hirsch nd S. Smle. Differentil Equtions, Dynmicl Systems, nd Liner Algebr. Acdemic Press, New York, [8] E.L. Ince. Ordinry Differentil Equtions. Dover Publictions, New York, [9] S. McLne nd G. Birkhoff. Algebr. New York: Mcmilln, [10] Jerrold E. Mrsden, Anthony J. Tromb, nd Aln Weinstein. Bsic Multivrible Clculus. New York: Springer-Verlg:W.H. Freemn, [11] I.G. Petrovski. Ordinry Differentil Equtions. Prentice-Hll Inc., Englewood Cliffs,

24 316 [12] Wlter Rudin. Principles of Mthemticl Anlysis. McGrw-Hill, New York, 1964.

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner

More information

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1 MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further

More information

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS Throughout, we let [, b] be bounded intervl in R. C 2 ([, b]) denotes the spce of functions with derivtives of second order continuous up to the endpoints. Cc 2

More information

Theoretical foundations of Gaussian quadrature

Theoretical foundations of Gaussian quadrature Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of

More information

1.9 C 2 inner variations

1.9 C 2 inner variations 46 CHAPTER 1. INDIRECT METHODS 1.9 C 2 inner vritions So fr, we hve restricted ttention to liner vritions. These re vritions of the form vx; ǫ = ux + ǫφx where φ is in some liner perturbtion clss P, for

More information

1 1D heat and wave equations on a finite interval

1 1D heat and wave equations on a finite interval 1 1D het nd wve equtions on finite intervl In this section we consider generl method of seprtion of vribles nd its pplictions to solving het eqution nd wve eqution on finite intervl ( 1, 2. Since by trnsltion

More information

The Regulated and Riemann Integrals

The Regulated and Riemann Integrals Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue

More information

Lecture 1. Functional series. Pointwise and uniform convergence.

Lecture 1. Functional series. Pointwise and uniform convergence. 1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is

More information

Sturm-Liouville Eigenvalue problem: Let p(x) > 0, q(x) 0, r(x) 0 in I = (a, b). Here we assume b > a. Let X C 2 1

Sturm-Liouville Eigenvalue problem: Let p(x) > 0, q(x) 0, r(x) 0 in I = (a, b). Here we assume b > a. Let X C 2 1 Ch.4. INTEGRAL EQUATIONS AND GREEN S FUNCTIONS Ronld B Guenther nd John W Lee, Prtil Differentil Equtions of Mthemticl Physics nd Integrl Equtions. Hildebrnd, Methods of Applied Mthemtics, second edition

More information

Variational Techniques for Sturm-Liouville Eigenvalue Problems

Variational Techniques for Sturm-Liouville Eigenvalue Problems Vritionl Techniques for Sturm-Liouville Eigenvlue Problems Vlerie Cormni Deprtment of Mthemtics nd Sttistics University of Nebrsk, Lincoln Lincoln, NE 68588 Emil: vcormni@mth.unl.edu Rolf Ryhm Deprtment

More information

Review of Calculus, cont d

Review of Calculus, cont d Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

Abstract inner product spaces

Abstract inner product spaces WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

A BRIEF INTRODUCTION TO UNIFORM CONVERGENCE. In the study of Fourier series, several questions arise naturally, such as: c n e int

A BRIEF INTRODUCTION TO UNIFORM CONVERGENCE. In the study of Fourier series, several questions arise naturally, such as: c n e int A BRIEF INTRODUCTION TO UNIFORM CONVERGENCE HANS RINGSTRÖM. Questions nd exmples In the study of Fourier series, severl questions rise nturlly, such s: () (2) re there conditions on c n, n Z, which ensure

More information

Linearity, linear operators, and self adjoint eigenvalue problems

Linearity, linear operators, and self adjoint eigenvalue problems Linerity, liner opertors, nd self djoint eigenvlue problems 1 Elements of liner lgebr The study of liner prtil differentil equtions utilizes, unsurprisingly, mny concepts from liner lgebr nd liner ordinry

More information

Lecture 14: Quadrature

Lecture 14: Quadrature Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

u t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx

u t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx Chpter 9: Green s functions for time-independent problems Introductory emples One-dimensionl het eqution Consider the one-dimensionl het eqution with boundry conditions nd initil condition We lredy know

More information

7.2 The Definite Integral

7.2 The Definite Integral 7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties

More information

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct

More information

Math 1B, lecture 4: Error bounds for numerical methods

Math 1B, lecture 4: Error bounds for numerical methods Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the

More information

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004 Advnced Clculus: MATH 410 Notes on Integrls nd Integrbility Professor Dvid Levermore 17 October 2004 1. Definite Integrls In this section we revisit the definite integrl tht you were introduced to when

More information

Consequently, the temperature must be the same at each point in the cross section at x. Let:

Consequently, the temperature must be the same at each point in the cross section at x. Let: HW 2 Comments: L1-3. Derive the het eqution for n inhomogeneous rod where the therml coefficients used in the derivtion of the het eqution for homogeneous rod now become functions of position x in the

More information

Review of basic calculus

Review of basic calculus Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below

More information

ODE: Existence and Uniqueness of a Solution

ODE: Existence and Uniqueness of a Solution Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =

More information

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph

More information

UNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3

UNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3 UNIFORM CONVERGENCE Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3 Suppose f n : Ω R or f n : Ω C is sequence of rel or complex functions, nd f n f s n in some sense. Furthermore,

More information

Orthogonal Polynomials

Orthogonal Polynomials Mth 4401 Gussin Qudrture Pge 1 Orthogonl Polynomils Orthogonl polynomils rise from series solutions to differentil equtions, lthough they cn be rrived t in vriety of different mnners. Orthogonl polynomils

More information

Applied Partial Differential Equations with Fourier Series and Boundary Value Problems 5th Edition Richard Haberman

Applied Partial Differential Equations with Fourier Series and Boundary Value Problems 5th Edition Richard Haberman Applied Prtil Differentil Equtions with Fourier Series nd Boundry Vlue Problems 5th Edition Richrd Hbermn Person Eduction Limited Edinburgh Gte Hrlow Essex CM20 2JE Englnd nd Associted Compnies throughout

More information

MATH 174A: PROBLEM SET 5. Suggested Solution

MATH 174A: PROBLEM SET 5. Suggested Solution MATH 174A: PROBLEM SET 5 Suggested Solution Problem 1. Suppose tht I [, b] is n intervl. Let f 1 b f() d for f C(I; R) (i.e. f is continuous rel-vlued function on I), nd let L 1 (I) denote the completion

More information

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008 Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of

More information

21.6 Green Functions for First Order Equations

21.6 Green Functions for First Order Equations 21.6 Green Functions for First Order Equtions Consider the first order inhomogeneous eqution subject to homogeneous initil condition, B[y] y() = 0. The Green function G( ξ) is defined s the solution to

More information

Advanced Calculus: MATH 410 Uniform Convergence of Functions Professor David Levermore 11 December 2015

Advanced Calculus: MATH 410 Uniform Convergence of Functions Professor David Levermore 11 December 2015 Advnced Clculus: MATH 410 Uniform Convergence of Functions Professor Dvid Levermore 11 December 2015 12. Sequences of Functions We now explore two notions of wht it mens for sequence of functions {f n

More information

STURM-LIOUVILLE PROBLEMS

STURM-LIOUVILLE PROBLEMS STURM-LIOUVILLE PROBLEMS Mrch 8, 24 We hve seen tht in the process of solving certin liner evolution equtions such s the het or wve equtions we re led in very nturl wy to n eigenvlue problem for second

More information

1 2-D Second Order Equations: Separation of Variables

1 2-D Second Order Equations: Separation of Variables Chpter 12 PDEs in Rectngles 1 2-D Second Order Equtions: Seprtion of Vribles 1. A second order liner prtil differentil eqution in two vribles x nd y is A 2 u x + B 2 u 2 x y + C 2 u y + D u 2 x + E u +

More information

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007 A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus

More information

Conservation Law. Chapter Goal. 5.2 Theory

Conservation Law. Chapter Goal. 5.2 Theory Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very

More information

Chapter 28. Fourier Series An Eigenvalue Problem.

Chapter 28. Fourier Series An Eigenvalue Problem. Chpter 28 Fourier Series Every time I close my eyes The noise inside me mplifies I cn t escpe I relive every moment of the dy Every misstep I hve mde Finds wy it cn invde My every thought And this is why

More information

UNIFORM CONVERGENCE MA 403: REAL ANALYSIS, INSTRUCTOR: B. V. LIMAYE

UNIFORM CONVERGENCE MA 403: REAL ANALYSIS, INSTRUCTOR: B. V. LIMAYE UNIFORM CONVERGENCE MA 403: REAL ANALYSIS, INSTRUCTOR: B. V. LIMAYE 1. Pointwise Convergence of Sequence Let E be set nd Y be metric spce. Consider functions f n : E Y for n = 1, 2,.... We sy tht the sequence

More information

(9) P (x)u + Q(x)u + R(x)u =0

(9) P (x)u + Q(x)u + R(x)u =0 STURM-LIOUVILLE THEORY 7 2. Second order liner ordinry differentil equtions 2.1. Recll some sic results. A second order liner ordinry differentil eqution (ODE) hs the form (9) P (x)u + Q(x)u + R(x)u =0

More information

Review of Riemann Integral

Review of Riemann Integral 1 Review of Riemnn Integrl In this chpter we review the definition of Riemnn integrl of bounded function f : [, b] R, nd point out its limittions so s to be convinced of the necessity of more generl integrl.

More information

Section 6.1 INTRO to LAPLACE TRANSFORMS

Section 6.1 INTRO to LAPLACE TRANSFORMS Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform

More information

Best Approximation in the 2-norm

Best Approximation in the 2-norm Jim Lmbers MAT 77 Fll Semester 1-11 Lecture 1 Notes These notes correspond to Sections 9. nd 9.3 in the text. Best Approximtion in the -norm Suppose tht we wish to obtin function f n (x) tht is liner combintion

More information

MAA 4212 Improper Integrals

MAA 4212 Improper Integrals Notes by Dvid Groisser, Copyright c 1995; revised 2002, 2009, 2014 MAA 4212 Improper Integrls The Riemnn integrl, while perfectly well-defined, is too restrictive for mny purposes; there re functions which

More information

Calculus of Variations

Calculus of Variations Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function

More information

Lecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature

Lecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 6: Singulr Integrls, Open Qudrture rules, nd Guss Qudrture (Compiled 6 August 7) In this lecture we discuss the

More information

The Riemann-Lebesgue Lemma

The Riemann-Lebesgue Lemma Physics 215 Winter 218 The Riemnn-Lebesgue Lemm The Riemnn Lebesgue Lemm is one of the most importnt results of Fourier nlysis nd symptotic nlysis. It hs mny physics pplictions, especilly in studies of

More information

MA 124 January 18, Derivatives are. Integrals are.

MA 124 January 18, Derivatives are. Integrals are. MA 124 Jnury 18, 2018 Prof PB s one-minute introduction to clculus Derivtives re. Integrls re. In Clculus 1, we lern limits, derivtives, some pplictions of derivtives, indefinite integrls, definite integrls,

More information

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1 Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution

More information

4. Calculus of Variations

4. Calculus of Variations 4. Clculus of Vritions Introduction - Typicl Problems The clculus of vritions generlises the theory of mxim nd minim. Exmple (): Shortest distnce between two points. On given surfce (e.g. plne), nd the

More information

Math Theory of Partial Differential Equations Lecture 2-9: Sturm-Liouville eigenvalue problems (continued).

Math Theory of Partial Differential Equations Lecture 2-9: Sturm-Liouville eigenvalue problems (continued). Mth 412-501 Theory of Prtil Differentil Equtions Lecture 2-9: Sturm-Liouville eigenvlue problems (continued). Regulr Sturm-Liouville eigenvlue problem: d ( p dφ ) + qφ + λσφ = 0 ( < x < b), dx dx β 1 φ()

More information

Lecture 19: Continuous Least Squares Approximation

Lecture 19: Continuous Least Squares Approximation Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for

More information

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0. STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA STEPHEN SCHECTER. The unit step function nd piecewise continuous functions The Heviside unit step function u(t) is given by if t

More information

STURM-LIOUVILLE THEORY, VARIATIONAL APPROACH

STURM-LIOUVILLE THEORY, VARIATIONAL APPROACH STURM-LIOUVILLE THEORY, VARIATIONAL APPROACH XIAO-BIAO LIN. Qudrtic functionl nd the Euler-Jcobi Eqution The purpose of this note is to study the Sturm-Liouville problem. We use the vritionl problem s

More information

Hilbert Spaces. Chapter Inner product spaces

Hilbert Spaces. Chapter Inner product spaces Chpter 4 Hilbert Spces 4.1 Inner product spces In the following we will discuss both complex nd rel vector spces. With L denoting either R or C we recll tht vector spce over L is set E equipped with ddition,

More information

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

1.2. Linear Variable Coefficient Equations. y + b ! = a y + b  Remark: The case b = 0 and a non-constant can be solved with the same idea as above. 1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 1 Sturm-Liouville Theory In the two preceing lectures I emonstrte the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series re just the tip of the iceerg of the theory

More information

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl

More information

Taylor Polynomial Inequalities

Taylor Polynomial Inequalities Tylor Polynomil Inequlities Ben Glin September 17, 24 Abstrct There re instnces where we my wish to pproximte the vlue of complicted function round given point by constructing simpler function such s polynomil

More information

Green s functions. f(t) =

Green s functions. f(t) = Consider the 2nd order liner inhomogeneous ODE Green s functions d 2 u 2 + k(t)du + p(t)u(t) = f(t). Of course, in prctice we ll only del with the two prticulr types of 2nd order ODEs we discussed lst

More information

Math 360: A primitive integral and elementary functions

Math 360: A primitive integral and elementary functions Mth 360: A primitive integrl nd elementry functions D. DeTurck University of Pennsylvni October 16, 2017 D. DeTurck Mth 360 001 2017C: Integrl/functions 1 / 32 Setup for the integrl prtitions Definition:

More information

Math 113 Fall Final Exam Review. 2. Applications of Integration Chapter 6 including sections and section 6.8

Math 113 Fall Final Exam Review. 2. Applications of Integration Chapter 6 including sections and section 6.8 Mth 3 Fll 0 The scope of the finl exm will include: Finl Exm Review. Integrls Chpter 5 including sections 5. 5.7, 5.0. Applictions of Integrtion Chpter 6 including sections 6. 6.5 nd section 6.8 3. Infinite

More information

Improper Integrals, and Differential Equations

Improper Integrals, and Differential Equations Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted

More information

c n φ n (x), 0 < x < L, (1) n=1

c n φ n (x), 0 < x < L, (1) n=1 SECTION : Fourier Series. MATH4. In section 4, we will study method clled Seprtion of Vribles for finding exct solutions to certin clss of prtil differentil equtions (PDEs. To do this, it will be necessry

More information

Math 270A: Numerical Linear Algebra

Math 270A: Numerical Linear Algebra Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner

More information

1 The Lagrange interpolation formula

1 The Lagrange interpolation formula Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x

More information

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes Jim Lmbers MAT 169 Fll Semester 2009-10 Lecture 4 Notes These notes correspond to Section 8.2 in the text. Series Wht is Series? An infinte series, usully referred to simply s series, is n sum of ll of

More information

Recitation 3: More Applications of the Derivative

Recitation 3: More Applications of the Derivative Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech

More information

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)

More information

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a). The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples

More information

DEFINITION The inner product of two functions f 1 and f 2 on an interval [a, b] is the number. ( f 1, f 2 ) b DEFINITION 11.1.

DEFINITION The inner product of two functions f 1 and f 2 on an interval [a, b] is the number. ( f 1, f 2 ) b DEFINITION 11.1. 398 CHAPTER 11 ORTHOGONAL FUNCTIONS AND FOURIER SERIES 11.1 ORTHOGONAL FUNCTIONS REVIEW MATERIAL The notions of generlized vectors nd vector spces cn e found in ny liner lger text. INTRODUCTION The concepts

More information

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019 ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS MATH00030 SEMESTER 208/209 DR. ANTHONY BROWN 7.. Introduction to Integrtion. 7. Integrl Clculus As ws the cse with the chpter on differentil

More information

AMATH 731: Applied Functional Analysis Fall Additional notes on Fréchet derivatives

AMATH 731: Applied Functional Analysis Fall Additional notes on Fréchet derivatives AMATH 731: Applied Functionl Anlysis Fll 214 Additionl notes on Fréchet derivtives (To ccompny Section 3.1 of the AMATH 731 Course Notes) Let X,Y be normed liner spces. The Fréchet derivtive of n opertor

More information

Continuous Random Variables

Continuous Random Variables STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht

More information

1 E3102: a study guide and review, Version 1.0

1 E3102: a study guide and review, Version 1.0 1 E3102: study guide nd review, Version 1.0 Here is list of subjects tht I think we ve covered in clss (your milege my vry). If you understnd nd cn do the bsic problems in this guide you should be in very

More information

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy

More information

SOME PROPERTIES OF CHEBYSHEV SYSTEMS

SOME PROPERTIES OF CHEBYSHEV SYSTEMS SOME PROPERTIES OF CHEBYSHEV SYSTEMS RICHARD A. ZALIK Abstrct. We study Chebyshev systems defined on n intervl, whose constituent functions re either complex or rel vlued, nd focus on problems tht my hve

More information

Green function and Eigenfunctions

Green function and Eigenfunctions Green function nd Eigenfunctions Let L e regulr Sturm-Liouville opertor on n intervl (, ) together with regulr oundry conditions. We denote y, φ ( n, x ) the eigenvlues nd corresponding normlized eigenfunctions

More information

8 Laplace s Method and Local Limit Theorems

8 Laplace s Method and Local Limit Theorems 8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved

More information

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite

More information

Here we study square linear systems and properties of their coefficient matrices as they relate to the solution set of the linear system.

Here we study square linear systems and properties of their coefficient matrices as they relate to the solution set of the linear system. Section 24 Nonsingulr Liner Systems Here we study squre liner systems nd properties of their coefficient mtrices s they relte to the solution set of the liner system Let A be n n Then we know from previous

More information

f(x)dx . Show that there 1, 0 < x 1 does not exist a differentiable function g : [ 1, 1] R such that g (x) = f(x) for all

f(x)dx . Show that there 1, 0 < x 1 does not exist a differentiable function g : [ 1, 1] R such that g (x) = f(x) for all 3 Definite Integrl 3.1 Introduction In school one comes cross the definition of the integrl of rel vlued function defined on closed nd bounded intervl [, b] between the limits nd b, i.e., f(x)dx s the

More information

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz

More information

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response

More information

(4.1) D r v(t) ω(t, v(t))

(4.1) D r v(t) ω(t, v(t)) 1.4. Differentil inequlities. Let D r denote the right hnd derivtive of function. If ω(t, u) is sclr function of the sclrs t, u in some open connected set Ω, we sy tht function v(t), t < b, is solution

More information

Introduction to the Calculus of Variations

Introduction to the Calculus of Variations Introduction to the Clculus of Vritions Jim Fischer Mrch 20, 1999 Abstrct This is self-contined pper which introduces fundmentl problem in the clculus of vritions, the problem of finding extreme vlues

More information

Math Advanced Calculus II

Math Advanced Calculus II Mth 452 - Advnced Clculus II Line Integrls nd Green s Theorem The min gol of this chpter is to prove Stoke s theorem, which is the multivrible version of the fundmentl theorem of clculus. We will be focused

More information

Chapter 4. Additional Variational Concepts

Chapter 4. Additional Variational Concepts Chpter 4 Additionl Vritionl Concepts 137 In the previous chpter we considered clculus o vrition problems which hd ixed boundry conditions. Tht is, in one dimension the end point conditions were speciied.

More information

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1 Tor Kjellsson Stockholm University Chpter 5 5. Strting with the following informtion: R = m r + m r m + m, r = r r we wnt to derive: µ = m m m + m r = R + µ m r, r = R µ m r 3 = µ m R + r, = µ m R r. 4

More information

Summary: Method of Separation of Variables

Summary: Method of Separation of Variables Physics 246 Electricity nd Mgnetism I, Fll 26, Lecture 22 1 Summry: Method of Seprtion of Vribles 1. Seprtion of Vribles in Crtesin Coordintes 2. Fourier Series Suggested Reding: Griffiths: Chpter 3, Section

More information

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by. NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with

More information

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1 3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =

More information

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS. THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem

More information

1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ),

1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ), 1. Guss-Jcobi qudrture nd Legendre polynomils Simpson s rule for evluting n integrl f(t)dt gives the correct nswer with error of bout O(n 4 ) (with constnt tht depends on f, in prticulr, it depends on

More information

Matrices, Moments and Quadrature, cont d

Matrices, Moments and Quadrature, cont d Jim Lmbers MAT 285 Summer Session 2015-16 Lecture 2 Notes Mtrices, Moments nd Qudrture, cont d We hve described how Jcobi mtrices cn be used to compute nodes nd weights for Gussin qudrture rules for generl

More information

Sections 5.2: The Definite Integral

Sections 5.2: The Definite Integral Sections 5.2: The Definite Integrl In this section we shll formlize the ides from the lst section to functions in generl. We strt with forml definition.. The Definite Integrl Definition.. Suppose f(x)

More information

Chapter 0. What is the Lebesgue integral about?

Chapter 0. What is the Lebesgue integral about? Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous

More information

ENGI 9420 Lecture Notes 7 - Fourier Series Page 7.01

ENGI 9420 Lecture Notes 7 - Fourier Series Page 7.01 ENGI 940 ecture Notes 7 - Fourier Series Pge 7.0 7. Fourier Series nd Fourier Trnsforms Fourier series hve multiple purposes, including the provision of series solutions to some liner prtil differentil

More information