Two special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
|
|
- Cassandra Johns
- 5 years ago
- Views:
Transcription
1 LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p Boundary value problems in other coordinate system. Laplace s equation in polar coordinates. p Exercises: , , , Two special equations: Bessel s and Legendre s equations The following linear, homogeneous differential equations of the 2nd order play an important role in mathematical physics x 2 y + xy + (x 2 ν 2 )y = 0, ν [0, ), x [0, ) (1) (1 x 2 )y 2xy + n(n + 1)y = 0, n = 1, 2,..., x ( 1, 1) (2) Eqs. (1) are named the Bessel differential equation and the Legendre differential equation, respectively. The solutions of Bessel s equations often appear for problems with rotational symmetry, which are formulated in polar coordinates. Similarly, solutions to Legendre s equation appear in spherical symmetric problems formulated in spherical coordinates. The general solution of (1) and (2) both have the form, cf. Theorem 4.2, p. 130 y(x) = c 1 y 1 (x) + c 2 y 2 (x) (3) y 1 (x) and y 2 (x) are linearly independent solutions to the homogeneous differential equations, and c 1 and c 2 are arbitrary constants. One should compare the Bessel equation (1) with the corresponding Cauchy-Euler equation, see p. 173 x 2 y + xy ν 2 y = 0 (4) The two differential equations are identical save the coefficient of y(x). However, for ν 0 even these factors become identical for x 0. Hence, the solutions of (1) and (4) for ν 0 should become asymptotically identical as x 0. For ν [0, ), ν 0, 1, 2,... (i.e. ν is not integer valued) a fundamental solution system to the Bessel equation can be shown to be 1
2 y 1 (x) = J ν (x) = y 2 (x) = J ν (x) = n=0 n=0 ( 1) n x ) 2n+ν (5) n! Γ(1 + ν + n)( 2 ( 1) n x ) 2n ν (6) n! Γ(1 ν + n)( 2 J ν (x) is known as the Bessel function of the first kind and order ν. Γ(x) signifies the so-called gamma function defined as, see App. 1. Γ(x) = 0 e t t x 1 dt (7) The gamma function fulfills the fundamental relation Γ(x) = (x 1)Γ(x 1), which can be proved by integration by parts of (7). Obvious, Γ(1) = 1. For n = 2, 3,... the following relation may then be derived recursively Γ(n) = (n 1)! (8) The series solutions (5) and (6) have been systematically derived in the textbook. Alternatively, the validity of these series may be proved simple by insertion into (1). Obviously, y 1 (x) and y 2 (x) become identical for ν = 0. Actually, for integer valued ν = n it can be shown that J n (x) and J n (x) become linearly dependent as follows J n (x) = cos(nπ)j n (x) = ( 1) n J n (x) (9) If ν = 0, 1, 2,... we can still use the solution (5). A linear independent solution can in this case be constructed from the limit passing Y n (x) = lim ν n cos(νπ)j ν (x) J ν (x) sin(νπ) (10) Throughout the limit passing the right hand side is a linear combination of the solutions J ν (x) and J ν (x). Hence, the Bessel equation is fulfilled at any time during the limit passing. Both the numerator and the denominator becomes zero in the limit as seen from (9). However, a finite limit value Y n (x) can be shown to exist, which also must be a solution to the Bessel equation. Y n (x) is known as the Bessel function of the second kind and order n. A series solution similar to (5) and (6) can also be derived for these functions. 2
3 The graphs of the lower order Bessel functions J 0 (x), J 1 (x), Y 0 (x), Y 1 (x) have been shown on Fig. 1, along with the asymptotic form of the series solutions valid in the vicinity of x = 0 The following recursive formulas can be derived, see p J ν(x) = ν x J ν(x) J ν+1 (x) Y n(x) = J 0(x) = J 1 (x) Y 0(x) = Y 1 (x) n x Y n(x) Y n+1 (x) (11) Next, x is replaced with λx and dx with λdx in (1). This leads to the so-called parametric Bessel equation of order ν x 2 y + xy + (λ 2 x 2 ν 2 )y = 0 (12) For ν = 0, 1, 2,... the solution becomes y(x) = c 1 J n (λx) + c 2 Y n (λx) (13) (12) can be written in the self-adjoint form d ( xy (x) ) + (λ 2 x n2 ) y(x) = 0 (14) dx x where ν = n has been assumed. (14) is recognized as a differential equation for a Sturm-Liouville problem with r(x) = x, p(x) = x and q(x) = n 2 /x, see Box 1. The differential equation (14) is solved in the interval (0, b). Since, r(0) = 0 the problem is singular at x = 0. Assume the following boundary condition of Sturm-Liouville type is prescribed for solutions of (14) at x = b A 2 y(b) + B 2 y (b) = 0 (15) It then follows from (20) that the eigenfunctions y i (x) of (14), satisfying the boundary condition (15), fulfills the orthogonality condition (21) in Box 1 (with a = 0). Of the two solution J n (λ i x) and Y n (λ i x) only J n (λ i x) is bounded at x = 0. Here and in what follows we shall only consider problems for which the eigenfunctions are of bounded variations, so y i (x) = J n (λ i x). Then, J n (λ i x) are seen to be orthogonal with respect to the function p(x) = x over the interval [0, b] as follows 3
4 0 xj n (λ i x)j n (λ j x)dx = 0, λ i λ j (16) The boundary condition (15) attains the form A 2 J n (λb) + B 2 λj n(λb) = 0 (17) In (17), A 2 and B 2 are given from the physics of the problem. Further, the chain rule has been d used, dx J n(λx) = J n(λx) d dx (λx) = λj n(λx). Here and in what follows J n( ) denotes differentiation with respect to the argument. Eigenvalues λ i are determined as solutions of λ fulfilling (17). For i = j the integral (16) can be evaluated analytically, based on simple manipulations of the differential equation (14), see p With y(x) = J n (λx), the result becomes 0 x Jn(λx) 2 dx = 1 ) ) 2 2 b2( J n(λb) 1 ( + (b 2 n2 ) 2 2 λ 2 J n (λb) (18) a) b) ( x ) 2 J 0 (x) = J 1 (x) = x 2 1 ( x ) Fig.1: a) Bessel functions of the first kind. b) Bessel functions of the second kind. Y 0 (x) = 2 π ln x 2 + Y 1 (x) = 2 π 1 x + Fig. 2: Legendre polynomials. 4
5 Box 1: Sturm-Liouville Problems, p. 449 A Sturm-Liouville problem is defined by the two-point boundary value problem d ( ) ( ) r(x)y + q(x) + λp(x) y(x) = 0, x (a, b) dx A 1 y(a) + B 1 y (a) = 0 A 2 y(b) + B 2 y (b) = 0 (19) where p, q, r and r are continuous functions on the closed interval [a, b], and r(x) > 0 and p(x) > 0 in the open interval (a, b). Obviously, (19) has the trivial solution y(x) 0. Non-trivial solutions y 1 (x), y 2 (x),... exist for discrete values λ 1, λ 2,... of the parameter λ. ( λ i, y i (x) ) represent the eigensolutions of the problem. Based on (19) the following identity may be proved for two different eigensolutions ( λ i, y i (x) ), ( λ j, y j (x) ), cf. p.450 ( ) λi λ j p(x)y i (x)y j (x)dx = a ( r(b) y i (b)y j(b) y j (b)y i(b) ) ( ) r(a) y i (a)y j(a) y j (a)y i(a) (20) The terms within the brackets on the right hand side of (20) vanish due to the boundary conditions in (19). Then, the following orthogonality condition of the eigenfunctions prevail a p(x)y i (x)y j (x)dx = 0 (21) Hence, the eigenfunctions of the Sturm-Liouville problem (19) are orthogonal over the interval of definition weighted with the function p(x). If r(a) > 0 and r(b) > 0 we talk about a regular Sturm-Liouville problem. In this case both the boundary conditions in (19) are necessary in order to prove the orthogonality condition (21). However, if r(a) = 0 or r(b) = 0 it becomes indifferent whether the terms within the brackets on the right hand side of (20) vanish. Hence, if r(a) = 0 or r(b) = 0 the corresponding boundary condition at x = a or x = b needs not be fulfilled of eigenfunctions of the differential equation, in order that the orthogonality condition (21) is valid. Such problems are referred to as singular Sturm-Liouville problems. 5
6 For Legendre s equation with the parameter n one solution can be shown to be an nth degree polynomial y 1 (x) = P n (x) = a n x n + a n 1 x n a 1 x + a 0 (22) The coefficients a n, a n 1,..., a 1, a 0 are determined so Legendre s differential equation is fulfilled. Additional, another solution y 2 (x) in terms of an infinite polynomial series may be derived. The systematic derivation of y 1 (x) and y 2 (x) has been performed in the textbook, but is not that interesting. P n (x) is denoted the Legendre polynomial of order n. The lower order Legendre polynomials are P 0 (x) = 1, P 1 (x) = x P 2 (x) = 1 ( 3x 2 1 ), P 3 (x) = 1 5x 2 2( 3 3x ) P 4 (x) = 1 ( 35x 4 30x ), P 5 (x) = 1 63x 8 8( 5 70x x ) (23) The graph of the functions have been shown on Fig. 2. As seen all polynomials pass through the point (1, 1). Further, Legendre polynomials of even order n are symmetric around the y axis, whereas the polynomials of odd order are all skew symmetric. The validity of the results (23) can be proved by insertion of y(x) = P n (x) into the relevant version of (2). The first few of these are n = 0 : (1 x 2 )P 0 2xP 0 = 0 n = 1 : (1 x 2 )P 1 2xP 1 + 2P 1 = 0 n = 2 : (1 x 2 )P 2 2xP 2 + 6P 2 = 0 n = 3 : (1 x 2 )P 3 2xP P 3 = 0 (24) It is not necessary to know the explicit expression for the Legendre polynomials in order to evaluate function values. These are more efficiently obtained from the recurrence relation, p. 267 P k+1 (x) = 2k + 1 k + 1 xp k(x) k k + 1 P k 1(x), k = 1, 2,... (25) As an example, (25) provides for k = 1: P 2 (x) = 3 2 xp 1(x) 1 2 P 0(x) = 3 2 x x (26) Eq. (2) can be written in the following self-adjoint form 6
7 d ( (1 x 2 )y ) + n(n + 1)y = 0, x ( 1, 1) (27) dx The equivalence of (2) and (27) follows by carrying out the differentiation in (27). (27) is recognized as a differential equation for a Sturm-Liouville problem with r(x) = (1 x 2 ), p(x) = 1, q(x) = 0 and λ 2 = n(n + 1), see Box 1. Since, r( 1) = r(1) = 0 the problem is singular at both the boundaries x = 1 and x = 1. Then, the orthogonality condition (21) with a = 1, b = 1 and p = 1 is fulfilled for any non-trivial solution y n (x) to (27), whatever boundary conditions these may fulfill at x = ±1. Hence, the Legendre polynomials P n (x) are orthogonal with respect to the weight function p(x) = 1 over the interval [ 1, 1], i.e. 1 1 P m (x)p n (x) dx = 0, m n (28) For m = n the integral (28) can be evaluated analytically, p P 2 n(x) dx = 2 2n + 1 (29) Fourier-Bessel and Fourier-Legendre series Bessel functions of the first kind J n (λ i x), and the Legendre polynomials P n (x) have been shown to be eigenfunctions to special, singular Sturm-Liouville eigenvalue problems. Then, it follows that it is possible to expand a function defined on the interval of the said eigenvalue problem into a generalized Fourier series after these eigenfunctions. Fourier-Bessel series We want to expand the function f(x), x [0, b] after the function system {J n (λ 1 x), J n (λ 2 x),...}, where the order n is common for all functions. These deviates by the eigenvalues, which are determined as solutions to (17). The expansion reads f(x) = c j J n (λ j x) (30) j=1 (30) is multiplied with r(x)j n (λ i x) = xj n (λ i x), followed by an integration over the interval [0, b]. From the orthogonality condition (16) follows 7
8 0 c i = x J n (λ i x) f(x) dx = c j x J n (λ i x)j n (λ j x) dx j=1 0 x J n(λ i x) f(x) dx 0 x J = n(λ 2 i x) dx x J n(λ i x) f(x) dx b 2( J n(λ i b) ) ( ) 2 (Jn + b 2 n2 (λ i b) ) 2 λ 2 i (31) where (18) has been used. From (17) follows that J n(λ i b) = h bλ i J n (λ i b), h = ba 2 B 2 (32) Then (31) can be written as c i = 2λ2 i 0 x J n(λ i x) f(x) dx ( )( h 2 n 2 + b 2 λ 2 i Jn (λ i b) ) 2 (33) The boundary condition (32) with the result (33) for the Fourier-Bessel expansion coefficients is referred to as Case II, p Assume that B 2 = 0 in (17), corresponding to the boundary condition J n (λ i b) = 0 (34) Then, from the recurrence relation (11) follows that J n(λ i b) = n λ i b J n(λ i b) J n+1 (λ i b) = J n+1 (λ i b) (35) By the use of (34) and (35), (31) reduces to c i = 2 0 x J n(λ i x) f(x) dx b 2( J n(λ i b) ) 2 = 2 0 x J n(λ i x) f(x) dx b 2( J n+1 (λ i b) ) 2 (36) The boundary condition (34) with the result (36) for the Fourier-Bessel expansion coefficients is referred to as Case I, p Finally, assume that A 2 = 0 and n = 0 in (17), corresponding to the boundary condition J 0(λ i b) = 0 J 1 (λ i b) = 0 (37) 8
9 where (11) has been used. Since, J 1 (0) = 0 it follows that λ 1 = 0 is an eigenvalue to the boundary condition (37). For λ 2, λ 3,..., which are all non-zero, (31) reduces to c i = 2 0 x J 0(λ i x) f(x) dx b 2( J 0 (λ i b) ) 2, i = 2, 3,... (38) For i = 1 the first statement in (31) for n = 0 and λ 1 = 0 provides c 1 = 0 x J 0(0 x) f(x) dx 0 x J 0 2(0 x) dx = 0 x f(x) dx 0 x dx = 2 b 2 0 x f(x) dx (39) where the result J 0 (λ 1 x) = J 0 (0 x) 1 has been used. The boundary condition (37) with the results (38), (39) for the Fourier-Bessel expansion coefficients is referred to as Case III, p In this case the Fourier-Bessel series (30) becomes f(x) = c 1 + c j J 0 (λ j x) (40) j=2 Fourier-Legendre series We want to expand the function f(x), x [ 1, 1] after the function system {P 0 (x), P 1 (x), P 2 (x),...}. The expansion reads f(x) = c n P n (x) (41) n=0 (41) is multiplied with P m (x), followed by an integration over the interval [ 1, 1]. From the orthogonality condition (28) follows 1 1 c m = P m (x) f(x) dx = 1 c n P m (x)p n (x) dx n=0 1 1 P m(x) f(x) dx 1 1 P 2 m(x) dx 1 = 2m P m (x) f(x) dx, m = 0, 1, 2,... (42) where (29) has been used. (41) is a power series quite similar to an truncated Taylor expansion. The difference is that (41) exists even for a non-analytical function, for which a Taylor expansion pr. definition does not exist. 9
10 Fig. 3: Fourier-Legendre approximation to non-analytical functions. Fig. 3 shows a Fourier-Legendre approximation function truncated after the 5th term to a nonanalytical (different Taylor expansions exist for the two branches). Theorem 11.5: Conditions for Convergence If f (x) and f 0 (x) are piecewise continuous on the open intervals (0, b) or ( 1, 1), then the Fourier-Bessel series or the Fourier-Legendre series converge to f (x) at any point, where f (x) is continuous, and to the average 21 [f (x ) + f (x+ )] at a point, where f (x) is discontinuous. 10
11 Box 2: Gradient vector and Laplace operator in polar coordinates Fig. 4: Differentials in polar and Cartesian coordinates. By use of simple projections the following relations between the polar differentials dr, rdθ and the Cartesian differentials dx and dy can be established } dr = cos θdx + sin θdy rdθ = sin θdx + cos θdy r x = cos θ, r y = sin θ θ x = 1 r sin θ, θ y = 1 (43) r cos θ u = u(r, θ) is considered a function of r and θ. By the use of the chain rule of partial differentiation the following result is obtained for the components of the gradient vector in polar coordinates u x = u r r x + u θ θ x u y = u r r y + u θ θ y u = cos θ r 1 u sin θ r θ u = sin θ r + 1 u cos θ r θ Partial differentiation of the gradient components in (44) with respect to x and y and use of the chain rule provides 2 u x 2 = 1 r sin2 θ u r + cos2 θ 2 u r 2 1 r cos θ sin θ 2 u r θ + 2 u cos θ sin θ r2 θ 1 r cos θ sin θ 2 u r θ + 1 r 2 sin2 θ 2 u θ 2 (45) 2 u y 2 = 1 r cos2 θ u r + sin2 θ 2 u r r cos θ sin θ 2 u r θ 2 u cos θ sin θ r2 θ + 1 r cos θ sin θ 2 u r θ + 1 r 2 cos2 θ 2 u θ 2 finally, addition of the components in (45) provides the following result for the Laplace operator in polar coordinates 2 u = 2 u x u y 2 = 2 u r u r r u r 2 θ 2 (46) (44) 11
12 Boundary value problems in other coordinate system. Laplace s equation in polar coordinates Fig. 5: Laplace boundary value problems in Cartesian and polar coordinates. In Section 12.2 we learned how to solve Laplace s differential equation for a rectangular domain by means of the separation method. The method require, that homogeneous boundary conditions (either Dirichlet of Neumann conditions) are prescribed on two opposite sides. In this lecture we shall learn how to solve the same equation within a segment or the entire of a ring shaped domain with inner radius r 0 and outer radius r 1. It turns out to be favorable to solve the problem in polar coordinates. Then we consider u = u(r, θ) as a function of the polar coordinates r and θ. Hence, the first problem is to express the Laplace operator in these coordinates. The derivations is described in Box 2, and the result becomes 2 u = 2 u x u y 2 = 2 u r u r r u r 2 θ 2 = 0 (47) Next, the boundary value problem defined in Fig. 5 can be solved by the separation method. We search for product solutions to (47) of the type u(r, θ) = R(r)Θ(θ) (48) Insertion in (47) provides 2 u = 2 u r u r r u r 2 θ 2 = R (r)θ(θ) + 1 r R (r)θ(θ) + 1 r 2 R(r)Θ (θ) = 0 r 2 R (r) + rr (r) R(r) = Θ (θ) Θ(θ) = λ2 (49) In the last statement of (49) the left hand side of the equation is a function of r, whereas the right hand side is a function of θ. This is only possible if both sides are equal to a constant λ 2, the so-called separation constant, which here at first is assumed to be positive real. Then, (49) is equivalent to the following ordinary differential equations 12
13 r 2 R (r) + rr (r) λ 2 R(r) = 0 (50) Θ (θ) + λ 2 Θ(θ) = 0 (51) (51) has the solution Θ(θ) = c 1 cos(λθ) + c 2 sin(λθ) (52) (50) is a differential equation of Cauchy-Euler s type, see p The solution is given as R(r) = c 3 r λ + c 4 r λ (53) The eigenvalues λ 1, λ 2,... of the parameter λ is determined by (52) by insertion into the homogeneous boundary conditions at θ = 0 and θ = θ 0. By use of the superposition principle the general solution can be written u(r, θ) = n=1 ( ) c 1,n cos(λ n θ) + c 2,n sin(λ n θ) )(c 3,n r λn + c 4,n r λn (54) we shall make a distinction between the following cases: 1 : Normal case : r 0 > 0, r 1 <, θ 0 < 2π a) r 0 = 0 c 4 = 0 b) r 1 = c 3 = 0 2 : Degenerated cases : (r λ, r ) (r λ, r 0) c) θ 0 = 2π λ n = n ( u(r, θ) = u(r, θ + 2π) ) 13
14 Example: Dirichlet boundary values The normal case is considered. u(r, 0) = 0, u(r, θ 0 ) = 0 implies Θ(0) = 0, Θ(θ 0 ) = 0. Insertion of the solution (52) into the indicated boundary conditions provides c 1 cos(λ 0) + c 2 sin(λ 0) = 0 c 1 = 0 c 2 sin(λθ 0 ) = 0 λ n θ 0 = nπ λ n = nπ θ 0 (55) Then, the solution (54) reduces to u(r, θ) = n=1 (A n r nπ θ 0 ) ) + B n r nπ θ 0 sin (nπ θθ0 (56) A n = c 2 c 3 and B n = c 2 c 4 are determined from the non-homogeneous boundary conditions at r = r 0 and r = r 1. 14
LEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationBoundary Value Problems in Cylindrical Coordinates
Boundary Value Problems in Cylindrical Coordinates 29 Outline Differential Operators in Various Coordinate Systems Laplace Equation in Cylindrical Coordinates Systems Bessel Functions Wave Equation the
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationOn an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University
On an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University njrose@math.ncsu.edu 1. INTRODUCTION. The classical eigenvalue problem for the Legendre Polynomials
More informationIntroduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series
CHAPTER 5 Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series We start with some introductory examples. 5.. Cauchy s equation The homogeneous Euler-Cauchy equation (Leonhard
More informationMA22S3 Summary Sheet: Ordinary Differential Equations
MA22S3 Summary Sheet: Ordinary Differential Equations December 14, 2017 Kreyszig s textbook is a suitable guide for this part of the module. Contents 1 Terminology 1 2 First order separable 2 2.1 Separable
More informationLecture IX. Definition 1 A non-singular Sturm 1 -Liouville 2 problem consists of a second order linear differential equation of the form.
Lecture IX Abstract When solving PDEs it is often necessary to represent the solution in terms of a series of orthogonal functions. One way to obtain an orthogonal family of functions is by solving a particular
More informationswapneel/207
Partial differential equations Swapneel Mahajan www.math.iitb.ac.in/ swapneel/207 1 1 Power series For a real number x 0 and a sequence (a n ) of real numbers, consider the expression a n (x x 0 ) n =
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationSeries Solution of Linear Ordinary Differential Equations
Series Solution of Linear Ordinary Differential Equations Department of Mathematics IIT Guwahati Aim: To study methods for determining series expansions for solutions to linear ODE with variable coefficients.
More information23 Elements of analytic ODE theory. Bessel s functions
23 Elements of analytic ODE theory. Bessel s functions Recall I am changing the variables) that we need to solve the so-called Bessel s equation 23. Elements of analytic ODE theory Let x 2 u + xu + x 2
More informationChapter 4. Series Solutions. 4.1 Introduction to Power Series
Series Solutions Chapter 4 In most sciences one generation tears down what another has built and what one has established another undoes. In mathematics alone each generation adds a new story to the old
More information7.3 Singular points and the method of Frobenius
284 CHAPTER 7. POWER SERIES METHODS 7.3 Singular points and the method of Frobenius Note: or.5 lectures, 8.4 and 8.5 in [EP], 5.4 5.7 in [BD] While behaviour of ODEs at singular points is more complicated,
More informationSection 5.2 Series Solution Near Ordinary Point
DE Section 5.2 Series Solution Near Ordinary Point Page 1 of 5 Section 5.2 Series Solution Near Ordinary Point We are interested in second order homogeneous linear differential equations with variable
More informationRelevant sections from AMATH 351 Course Notes (Wainwright): Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls):
Lecture 5 Series solutions to DEs Relevant sections from AMATH 35 Course Notes (Wainwright):.4. Relevant sections from AMATH 35 Course Notes (Poulin and Ingalls): 2.-2.3 As mentioned earlier in this course,
More informationSolving Differential Equations Using Power Series
LECTURE 25 Solving Differential Equations Using Power Series We are now going to employ power series to find solutions to differential equations of the form (25.) y + p(x)y + q(x)y = 0 where the functions
More informationSolving Differential Equations Using Power Series
LECTURE 8 Solving Differential Equations Using Power Series We are now going to employ power series to find solutions to differential equations of the form () y + p(x)y + q(x)y = 0 where the functions
More informationMath 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions
Math 3 Winter 2 Avance Bounary Value Problems I Bessel s Equation an Bessel Functions Department of Mathematical an Statistical Sciences University of Alberta Bessel s Equation an Bessel Functions We use
More informationLinear Operators, Eigenvalues, and Green s Operator
Chapter 10 Linear Operators, Eigenvalues, and Green s Operator We begin with a reminder of facts which should be known from previous courses. 10.1 Inner Product Space A vector space is a collection of
More informationSeries Solutions of Differential Equations
Chapter 6 Series Solutions of Differential Equations In this chapter we consider methods for solving differential equations using power series. Sequences and infinite series are also involved in this treatment.
More informationIntroduction to Sturm-Liouville Theory
Introduction to R. C. Trinity University Partial Differential Equations April 10, 2014 Sturm-Liouville problems Definition: A (second order) Sturm-Liouville (S-L) problem consists of A Sturm-Liouville
More information1 A complete Fourier series solution
Math 128 Notes 13 In this last set of notes I will try to tie up some loose ends. 1 A complete Fourier series solution First here is an example of the full solution of a pde by Fourier series. Consider
More informationSturm-Liouville Theory
More on Ryan C. Trinity University Partial Differential Equations April 19, 2012 Recall: A Sturm-Liouville (S-L) problem consists of A Sturm-Liouville equation on an interval: (p(x)y ) + (q(x) + λr(x))y
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More informationPower Series Solutions to the Legendre Equation
Power Series Solutions to the Legendre Equation Department of Mathematics IIT Guwahati The Legendre equation The equation (1 x 2 )y 2xy + α(α + 1)y = 0, (1) where α is any real constant, is called Legendre
More information(a) Determine the general solution for φ(ρ) near ρ = 0 for arbitary values E. (b) Show that the regular solution at ρ = 0 has the series expansion
Problem 1. Curious Wave Functions The eigenfunctions of a D7 brane in a curved geometry lead to the following eigenvalue equation of the Sturm Liouville type ρ ρ 3 ρ φ n (ρ) = E n w(ρ)φ n (ρ) w(ρ) = where
More informationDifferential Equations
Differential Equations Problem Sheet 1 3 rd November 2011 First-Order Ordinary Differential Equations 1. Find the general solutions of the following separable differential equations. Which equations are
More informationOrdinary Differential Equations II
Ordinary Differential Equations II February 23 2017 Separation of variables Wave eq. (PDE) 2 u t (t, x) = 2 u 2 c2 (t, x), x2 c > 0 constant. Describes small vibrations in a homogeneous string. u(t, x)
More informationConnection to Laplacian in spherical coordinates (Chapter 13)
Connection to Laplacian in spherical coordinates (Chapter 13) We might often encounter the Laplace equation and spherical coordinates might be the most convenient 2 u(r, θ, φ) = 0 We already saw in Chapter
More information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationSolutions to Laplace s Equations- II
Solutions to Laplace s Equations- II Lecture 15: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Laplace s Equation in Spherical Coordinates : In spherical coordinates
More informationPhysics 250 Green s functions for ordinary differential equations
Physics 25 Green s functions for ordinary differential equations Peter Young November 25, 27 Homogeneous Equations We have already discussed second order linear homogeneous differential equations, which
More informationThis ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0
Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More information1 Bessel, Neumann, and Hankel Functions: J ν (x), N ν (x), H (1)
Bessel, Neumann, and Hankel Functions: J ν x), N ν x), H ) ν x), H ) ν x) Bessel functions are solutions of the following differential equation: x y +xy +x ν )y =.) Any two of the following functions are
More informationSeries Solutions Near a Regular Singular Point
Series Solutions Near a Regular Singular Point MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background We will find a power series solution to the equation:
More information2 Series Solutions near a Regular Singular Point
McGill University Math 325A: Differential Equations LECTURE 17: SERIES SOLUTION OF LINEAR DIFFERENTIAL EQUATIONS II 1 Introduction Text: Chap. 8 In this lecture we investigate series solutions for the
More information5.4 Bessel s Equation. Bessel Functions
SEC 54 Bessel s Equation Bessel Functions J (x) 87 # with y dy>dt, etc, constant A, B, C, D, K, and t 5 HYPERGEOMETRIC ODE At B (t t )(t t ), t t, can be reduced to the hypergeometric equation with independent
More informationLecture 4.6: Some special orthogonal functions
Lecture 4.6: Some special orthogonal functions Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4340, Advanced Engineering Mathematics
More informationPhysics 116C Solutions to the Practice Final Exam Fall The method of Frobenius fails because it works only for Fuchsian equations, of the type
Physics 6C Solutions to the Practice Final Exam Fall 2 Consider the differential equation x y +x(x+)y y = () (a) Explain why the method of Frobenius method fails for this problem. The method of Frobenius
More informationQualification Exam: Mathematical Methods
Qualification Exam: Mathematical Methods Name:, QEID#41534189: August, 218 Qualification Exam QEID#41534189 2 1 Mathematical Methods I Problem 1. ID:MM-1-2 Solve the differential equation dy + y = sin
More informationPHYS 404 Lecture 1: Legendre Functions
PHYS 404 Lecture 1: Legendre Functions Dr. Vasileios Lempesis PHYS 404 - LECTURE 1 DR. V. LEMPESIS 1 Legendre Functions physical justification Legendre functions or Legendre polynomials are the solutions
More informationWaves on 2 and 3 dimensional domains
Chapter 14 Waves on 2 and 3 dimensional domains We now turn to the studying the initial boundary value problem for the wave equation in two and three dimensions. In this chapter we focus on the situation
More information12d. Regular Singular Points
October 22, 2012 12d-1 12d. Regular Singular Points We have studied solutions to the linear second order differential equations of the form P (x)y + Q(x)y + R(x)y = 0 (1) in the cases with P, Q, R real
More informationLecture Notes for MAE 3100: Introduction to Applied Mathematics
ecture Notes for MAE 31: Introduction to Applied Mathematics Richard H. Rand Cornell University Ithaca NY 14853 rhr2@cornell.edu http://audiophile.tam.cornell.edu/randdocs/ version 17 Copyright 215 by
More informationMathematics 324 Riemann Zeta Function August 5, 2005
Mathematics 324 Riemann Zeta Function August 5, 25 In this note we give an introduction to the Riemann zeta function, which connects the ideas of real analysis with the arithmetic of the integers. Define
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georgia Tech PHYS 612 Mathematical Methods of Physics I Instructor: Predrag Cvitanović Fall semester 2012 Homework Set #5 due October 2, 2012 == show all your work for maximum credit, == put labels, title,
More informationENGI 9420 Lecture Notes 1 - ODEs Page 1.01
ENGI 940 Lecture Notes - ODEs Page.0. Ordinary Differential Equations An equation involving a function of one independent variable and the derivative(s) of that function is an ordinary differential equation
More informationPower Series Solutions to the Legendre Equation
Department of Mathematics IIT Guwahati The Legendre equation The equation (1 x 2 )y 2xy + α(α + 1)y = 0, (1) where α is any real constant, is called Legendre s equation. When α Z +, the equation has polynomial
More information1 Series Solutions Near Regular Singular Points
1 Series Solutions Near Regular Singular Points All of the work here will be directed toward finding series solutions of a second order linear homogeneous ordinary differential equation: P xy + Qxy + Rxy
More informationBessel s and legendre s equations
Chapter 12 Bessel s and legendre s equations 12.1 Introduction Many linear differential equations having variable coefficients cannot be solved by usual methods and we need to employ series solution method
More informationSeries solutions to a second order linear differential equation with regular singular points
Physics 6C Fall 0 Series solutions to a second order linear differential equation with regular singular points Consider the second-order linear differential equation, d y dx + p(x) dy x dx + q(x) y = 0,
More informationProblem Set 8 Mar 5, 2004 Due Mar 10, 2004 ACM 95b/100b 3pm at Firestone 303 E. Sterl Phinney (2 pts) Include grading section number
Problem Set 8 Mar 5, 24 Due Mar 1, 24 ACM 95b/1b 3pm at Firestone 33 E. Sterl Phinney (2 pts) Include grading section number Useful Readings: For Green s functions, see class notes and refs on PS7 (esp
More information1 Introduction. Green s function notes 2018
Green s function notes 8 Introduction Back in the "formal" notes, we derived the potential in terms of the Green s function. Dirichlet problem: Equation (7) in "formal" notes is Φ () Z ( ) ( ) 3 Z Φ (
More informationElectromagnetism HW 1 math review
Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:
More informationAppendix A Vector Analysis
Appendix A Vector Analysis A.1 Orthogonal Coordinate Systems A.1.1 Cartesian (Rectangular Coordinate System The unit vectors are denoted by x, ŷ, ẑ in the Cartesian system. By convention, ( x, ŷ, ẑ triplet
More informationSPECTRAL METHODS: ORTHOGONAL POLYNOMIALS
SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS 31 October, 2007 1 INTRODUCTION 2 ORTHOGONAL POLYNOMIALS Properties of Orthogonal Polynomials 3 GAUSS INTEGRATION Gauss- Radau Integration Gauss -Lobatto Integration
More informationNotes on Special Functions
Spring 25 1 Notes on Special Functions Francis J. Narcowich Department of Mathematics Texas A&M University College Station, TX 77843-3368 Introduction These notes are for our classes on special functions.
More informationS. Ghorai 1. Lecture XV Bessel s equation, Bessel s function. e t t p 1 dt, p > 0. (1)
S Ghorai 1 1 Gamma function Gamma function is defined by Lecture XV Bessel s equation, Bessel s function Γp) = e t t p 1 dt, p > 1) The integral in 1) is convergent that can be proved easily Some special
More informationPartial Differential Equations
Partial Differential Equations Spring Exam 3 Review Solutions Exercise. We utilize the general solution to the Dirichlet problem in rectangle given in the textbook on page 68. In the notation used there
More information1 Parametric Bessel Equation and Bessel-Fourier Series
1 Parametric Bessel Equation an Bessel-Fourier Series Recall the parametric Bessel equation of orer n: x 2 y + xy + (a 2 x 2 n 2 )y = (1.1) The general solution is given by y = J n (ax) +Y n (ax). If we
More informationSeries Solutions. 8.1 Taylor Polynomials
8 Series Solutions 8.1 Taylor Polynomials Polynomial functions, as we have seen, are well behaved. They are continuous everywhere, and have continuous derivatives of all orders everywhere. It also turns
More informationOrdinary differential equations Notes for FYS3140
Ordinary differential equations Notes for FYS3140 Susanne Viefers, Dept of Physics, University of Oslo April 4, 2018 Abstract Ordinary differential equations show up in many places in physics, and these
More informationEquations with regular-singular points (Sect. 5.5).
Equations with regular-singular points (Sect. 5.5). Equations with regular-singular points. s: Equations with regular-singular points. Method to find solutions. : Method to find solutions. Recall: The
More informationMathematics for Chemistry: Exam Set 1
Mathematics for Chemistry: Exam Set 1 June 18, 017 1 mark Questions 1. The minimum value of the rank of any 5 3 matrix is 0 1 3. The trace of an identity n n matrix is equal to 1-1 0 n 3. A square matrix
More informationMathematics for Chemistry: Exam Set 1
Mathematics for Chemistry: Exam Set 1 March 19, 017 1 mark Questions 1. The maximum value of the rank of any 5 3 matrix is (a (b3 4 5. The determinant of an identity n n matrix is equal to (a 1 (b -1 0
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationKREYSZIG E Advanced Engineering Mathematics (10th ed., Wiley 2011) Chapter 11 - Fourier analysis
KREYSZIG E Advanced Engineering Mathematics (th ed., Wiley ) Chapter - Fourier analysis . CHAPTER Fourier Analysis 474 This chapter on Fourier analysis covers three broad areas: Fourier series in Secs...4,
More information3 Green s functions in 2 and 3D
William J. Parnell: MT34032. Section 3: Green s functions in 2 and 3 57 3 Green s functions in 2 and 3 Unlike the one dimensional case where Green s functions can be found explicitly for a number of different
More informationPhysics 6303 Lecture 11 September 24, LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation
Physics 6303 Lecture September 24, 208 LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation, l l l l l l. Consider problems that are no axisymmetric; i.e., the potential depends
More informationMath Assignment 11
Math 2280 - Assignment 11 Dylan Zwick Fall 2013 Section 8.1-2, 8, 13, 21, 25 Section 8.2-1, 7, 14, 17, 32 Section 8.3-1, 8, 15, 18, 24 1 Section 8.1 - Introduction and Review of Power Series 8.1.2 - Find
More information14 Fourier analysis. Read: Boas Ch. 7.
14 Fourier analysis Read: Boas Ch. 7. 14.1 Function spaces A function can be thought of as an element of a kind of vector space. After all, a function f(x) is merely a set of numbers, one for each point
More informationDegree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m.
Degree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m. Candidates should submit answers to a maximum of four
More informationLegendre s Equation. PHYS Southern Illinois University. October 18, 2016
Legendre s Equation PHYS 500 - Southern Illinois University October 18, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 18, 2016 1 / 11 Legendre s Equation Recall We are trying
More informationPartial Differential Equations for Engineering Math 312, Fall 2012
Partial Differential Equations for Engineering Math 312, Fall 2012 Jens Lorenz July 17, 2012 Contents Department of Mathematics and Statistics, UNM, Albuquerque, NM 87131 1 Second Order ODEs with Constant
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationSeries Solutions of Linear ODEs
Chapter 2 Series Solutions of Linear ODEs This Chapter is concerned with solutions of linear Ordinary Differential Equations (ODE). We will start by reviewing some basic concepts and solution methods for
More informationBackground and Definitions...2. Legendre s Equation, Functions and Polynomials...4 Legendre s Associated Equation and Functions...
Legendre Polynomials and Functions Reading Problems Outline Background and Definitions...2 Definitions...3 Theory...4 Legendre s Equation, Functions and Polynomials...4 Legendre s Associated Equation and
More informationTyn Myint-U Lokenath Debnath. Linear Partial Differential Equations for Scientists and Engineers. Fourth Edition. Birkhauser Boston Basel Berlin
Tyn Myint-U Lokenath Debnath Linear Partial Differential Equations for Scientists and Engineers Fourth Edition Birkhauser Boston Basel Berlin Preface to the Fourth Edition Preface to the Third Edition
More informationReview of Power Series
Review of Power Series MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Introduction In addition to the techniques we have studied so far, we may use power
More informationWelcome to Math 257/316 - Partial Differential Equations
Welcome to Math 257/316 - Partial Differential Equations Instructor: Mona Rahmani email: mrahmani@math.ubc.ca Office: Mathematics Building 110 Office hours: Mondays 2-3 pm, Wednesdays and Fridays 1-2 pm.
More informationMathematical Modeling using Partial Differential Equations (PDE s)
Mathematical Modeling using Partial Differential Equations (PDE s) 145. Physical Models: heat conduction, vibration. 146. Mathematical Models: why build them. The solution to the mathematical model will
More information(L, t) = 0, t > 0. (iii)
. Sturm Liouville Boundary Value Problems 69 where E is Young s modulus and ρ is the mass per unit volume. If the end x = isfixed, then the boundary condition there is u(, t) =, t >. (ii) Suppose that
More informationLAPLACE EQUATION. = 2 is call the Laplacian or del-square operator. In two dimensions, the expression of in rectangular and polar coordinates are
LAPLACE EQUATION If a diffusion or wave problem is stationary (time independent), the pde reduces to the Laplace equation u = u =, an archetype of second order elliptic pde. = 2 is call the Laplacian or
More informationThe method of Fröbenius
Note III.5 1 1 April 008 The method of Fröbenius For the general homogeneous ordinary differential equation y (x) + p(x)y (x) + q(x)y(x) = 0 (1) the series method works, as in the Hermite case, where both
More informationSeries Solutions of ODEs. Special Functions
C05.tex 6/4/0 3: 5 Page 65 Chap. 5 Series Solutions of ODEs. Special Functions We continue our studies of ODEs with Legendre s, Bessel s, and the hypergeometric equations. These ODEs have variable coefficients
More informationBOUNDARY PROBLEMS IN HIGHER DIMENSIONS. kt = X X = λ, and the series solutions have the form (for λ n 0):
BOUNDARY PROBLEMS IN HIGHER DIMENSIONS Time-space separation To solve the wave and diffusion equations u tt = c u or u t = k u in a bounded domain D with one of the three classical BCs (Dirichlet, Neumann,
More informationBefore you begin read these instructions carefully:
NATURAL SCIENCES TRIPOS Part IB & II (General) Tuesday, 29 May, 2012 9:00 am to 12:00 pm MATHEMATICS (1) Before you begin read these instructions carefully: You may submit answers to no more than six questions.
More informationAn Introduction to Partial Differential Equations in the Undergraduate Curriculum
An Introduction to Partial Differential Equations in the Undergraduate Curriculum John C. Polking LECTURE 12 Heat Transfer in the Ball 12.1. Outline of Lecture The problem The problem for radial temperatures
More informationSpherical Coordinates and Legendre Functions
Spherical Coordinates and Legendre Functions Spherical coordinates Let s adopt the notation for spherical coordinates that is standard in physics: φ = longitude or azimuth, θ = colatitude ( π 2 latitude)
More informationSolutions to the Sample Problems for the Final Exam UCLA Math 135, Winter 2015, Professor David Levermore
Solutions to the Sample Problems for the Final Exam UCLA Math 135, Winter 15, Professor David Levermore Every sample problem for the Midterm exam and every problem on the Midterm exam should be considered
More informationChapter 5.8: Bessel s equation
Chapter 5.8: Bessel s equation Bessel s equation of order ν is: x 2 y + xy + (x 2 ν 2 )y = 0. It has a regular singular point at x = 0. When ν = 0,, 2,..., this equation comes up when separating variables
More informationMath221: HW# 7 solutions
Math22: HW# 7 solutions Andy Royston November 7, 25.3.3 let x = e u. Then ln x = u, x2 = e 2u, and dx = e 2u du. Furthermore, when x =, u, and when x =, u =. Hence x 2 ln x) 3 dx = e 2u u 3 e u du) = e
More informationAn Introduction to Bessel Functions
An Introduction to R. C. Trinity University Partial Differential Equations Lecture 17 Bessel s equation Given p 0, the ordinary differential equation x 2 y + xy + (x 2 p 2 )y = 0, x > 0 is known as Bessel
More informationSpotlight on Laplace s Equation
16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian
More informationMath 2930 Worksheet Final Exam Review
Math 293 Worksheet Final Exam Review Week 14 November 3th, 217 Question 1. (* Solve the initial value problem y y = 2xe x, y( = 1 Question 2. (* Consider the differential equation: y = y y 3. (a Find the
More informationFig.1: Non-stationary temperature distribution in a circular plate.
LECTURE Problem 1 13. 5 : Fi.1: Non-stationary temperature distribution in a circular plate. The non-stationary, radial symmetric temperature distribution ur, t in a circular plate of radius c is determined
More informationLecture Notes on. Differential Equations. Emre Sermutlu
Lecture Notes on Differential Equations Emre Sermutlu ISBN: Copyright Notice: To my wife Nurten and my daughters İlayda and Alara Contents Preface ix 1 First Order ODE 1 1.1 Definitions.............................
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More information