Bessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics

Size: px
Start display at page:

Download "Bessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics"

Transcription

1 Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018

2 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy + (x 2 ν 2 )y = 0 Bessel s equation arises frequently in the mathematical modeling of vibration in polar, cylindrical, or spherical coordinates.

3 Vibrating Circular Membrane (i.e. a drum) Wave equation: u tt = c 2 (u xx + u yy ) for x 2 + y 2 a 2. Boundary condition: u(x, y, t) = 0 for x 2 + y 2 = a 2.

4 Polar Coordinates Since x = r cos θ and y = r sin θ we can re-write the boundary value problem as ( 1 Wave equation: u tt = c 2 r (r u r ) r + 1 ) r 2 u θθ for r a. Boundary condition: u(r, θ, t) = 0 for r = a.

5 Isotropy Assume the solution is independent of θ, i.e. the motion is the same in every direction. Wave equation: u tt = c2 r (r u r ) r for r a. Boundary condition: u(r, t) = 0 for r = a.

6 Isotropy Assume the solution is independent of θ, i.e. the motion is the same in every direction. Wave equation: u tt = c2 r (r u r ) r for r a. Boundary condition: u(r, t) = 0 for r = a. We will assume the solution u(r, t) = R(r)T (t). For PDEs this is called the method of separation of variables.

7 Separation of Variables u tt = c2 r (r u r ) r R(r)T (t) = c2 r (r R (r)) T (t) R(r)T (t) = c2 r (R (r) + r R (r))t (t) 1 T (t) c 2 T (t) = R (r) + r R (r) r R(r)

8 Separation of Variables u tt = c2 r (r u r ) r R(r)T (t) = c2 r (r R (r)) T (t) R(r)T (t) = c2 r (R (r) + r R (r))t (t) 1 T (t) c 2 T (t) = R (r) + r R (r) r R(r) Since the left-hand side is a function of t only and the right-hand side is a function of r only, they must be constant.

9 Constant R (r) + r R (r) r R(r) = λ 2 R (r) + r R (r) = λ 2 r R(r) r R (r) + R (r) + λ 2 r R(r) = 0 r 2 R (r) + r R (r) + λ 2 r 2 R(r) = 0

10 Change of Variable Define x = λr and y(x) = R ( x λ) then R (r) = dy dx λ R (r) = d 2 y dx 2 λ2 and substituting into the ODE yields r 2 R (r) + r R (r) + λ 2 r 2 R(r) = 0 r 2 d 2 y dx 2 λ2 + r dy dx λ + λ2 r 2 y = x 2 y + xy + x 2 y = 0 which is Bessel s equation of order 0.

11 Bessel s Equation of Order Zero x 2 y + xy + x 2 y = 0 has a regular singular point at x = 0. The exponents of singularity are r 1 = r 2 = 0. Two linearly independent solutions are: where H n = y 1 (x) = 1 + n=1 y 2 (x) = y 1 (x) ln x + n k=1 1 k. ( 1) n ( x ) 2n (n!) 2 2 n=1 ( 1) n+1 H ( n x ) 2n (n!) 2 2

12 Remarks Solution y 1 (x) is called the Bessel function of the first kind of order zero and denoted J 0 (x). The Bessel function of the second kind of order zero is denoted Y 0 (x) and defined as Y 0 (x) = 2 π [y 2(x) + (γ ln 2)J 0 (x)] where γ is the Euler Máscheroni constant and γ = lim n (H n ln n) General solution to Bessel s equation of order zero is written y(x) = c 1 J 0 (x) + c 2 Y 0 (x).

13 Graphs 1.0 y x J 0 (x) Y 0 (x)

14 Asymptotic Properties (1 of 2) when x. 0 = x 2 y + xy + (x 2 ν 2 )y = y + 1 ) x y + (1 ν2 y y + y x 2

15 Asymptotic Properties (1 of 2) when x. 0 = x 2 y + xy + (x 2 ν 2 )y = y + 1 ) x y + (1 ν2 y y + y Since cos(x c) and sin(x c) solve y + y = 0 we should expect J 0 (x) and Y 0 (x) to be oscillatory. In fact x 2 as x. J 0 (x) Y 0 (x) ( 2 πx ( 2 πx ) 1/2 ( cos x π ) 4 ) 1/2 ( sin x π ) 4

16 Asymptotic Properties (2 of 2) y J 0 (x) 2 π x cos x - π x

17 Bessel s Equation of Order 1/2 ( x 2 y + x y + x 2 1 ) y = 0 4 has a regular singular point at x = 0. Indicial equation: ( x ) p 0 = lim x x 0 x 2 q 0 = lim x 0 x 2 = 1 ) ( x x 2 = = r(r 1) + r 1 4 = r Exponents of singularity: r 1 = 1 2 and r 2 = 1 2.

18 Solution (1 of 8) Assuming y(x) = a n x r+n, differentiating, and substituting into the ODE we obtain 0 = x 2 (r + n)(r + n 1)a n x r+n 2 + x (r + n)a n x r+n 1 = = ( + x 2 1 ) a n x r+n 4 (r + n)(r + n 1)a n x r+n + (r + n)a n x r+n 1 4 a n x r+n + a n x r+n+2 [ (r + n)(r + n 1) + (r + n) 1 ] a n x r+n + 4 a n x r+n+2

19 Solution (2 of 8) 0 = = = = [ (r + n)(r + n 1) + (r + n) 1 ] a n x r+n + 4 [ (r + n) 2 1 ] a n x r+n + 4 ( r a n x r+n+2 [ (r + n) 2 1 ] a n x r+n + 4 n=2 ) ( a 0 x r + (r + 1) n=2 [ (r + n) 2 1 ] a n x r+n + 4 a n 2 x r+n ) a 1 x r+1 a n 2 x r+n n=2 a n x r+n+2

20 Solution (3 of 8) 0 = = ( r 2 1 ) ( a 0 x r + (r + 1) 2 1 ) a 1 x r n=2 ( r n=2 [ (r + n) 2 1 ] a n x r+n + 4 ) ( a 0 x r + (r + 1) n=2 a n 2 x r+n ) a 1 x r+1 ([ (r + n) 2 1 ] ) a n + a n 2 x r+n 4 Recurrence relation: a n (r) = a n 2(r) (r + n) 2 1 4

21 Solution (4 of 8) When r = r 1 = 1 2 then the expression ( (r + 1) 2 1 ) a 1 x r+1 = 0 4 when a 1 = 0.

22 Solution (4 of 8) When r = r 1 = 1 2 then the expression ( (r + 1) 2 1 ) a 1 x r+1 = 0 4 when a 1 = 0. Thus we will choose a 0 = 1 and a 1 = 0. a 1 = 0 implies a 2n+1 = 0 for all n N.

23 Solution (5 of 8) a 2 (1/2) = a 0 (2)(3) = 1 3! a 4 (1/2) = a 2 (4)(5) = 1 5! a 6 (1/2) = a 4 (6)(7) = 1 7!. a 2n (1/2) = ( 1) n (2n + 1)! y 1 (x) = x 1/2 [1 + = x 1/2 sin x. n=1 ] ( 1) n x 2n = x 1/2 ( 1) n x 2n+1 (2n + 1)! (2n + 1)!

24 Solution (6 of 8) Note: r 1 r 2 = 1 N but when r = r 2 = 1 2 then the expression ( (r + 1) 2 1 ) a 1 x r+1 = 0 4 for any value of a 1.

25 Solution (6 of 8) Note: r 1 r 2 = 1 N but when r = r 2 = 1 2 then the expression ( (r + 1) 2 1 ) a 1 x r+1 = 0 4 for any value of a 1. Thus we will again choose a 0 = 1 and a 1 = 0. a 1 = 0 implies a 2n+1 = 0 for all n N.

26 Solution (7 of 8) y 2 (x) = x 1/2 [1 + a 2 ( 1/2) = a 0 (2)(1) = 1 2! a 4 ( 1/2) = a 2 (4)(3) = 1 4! a 6 ( 1/2) = a 4 (6)(5) = 1 6! a 2n ( 1/2) = ( 1)n (2n)! = x 1/2 cos x. n=1. ] ( 1) n x 2n = x 1/2 ( 1) n x 2n (2n)! (2n)!

27 Solution (8 of 8) Bessel function of the first kind of order one-half J 1/2 (x) = ( ) 2 1/2 ( ) 2 1/2 y 1 (x) = sin x. π πx Bessel function of the second kind of order one-half J 1/2 (x) = ( ) 2 1/2 ( ) 2 1/2 y 2 (x) = cos x. π πx

28 Graphs y x J 1 (x) 2 Y 1 (x)

29 Bessel s Equation of Order One x 2 y + xy + (x 2 1)y = 0 The value x = 0 is a regular singular point since ( x ) lim x x 0 x 2 = 1 = p 0 ( x lim x 2 2 ) 1 = 1 = q 0. x 0 The exponents of singularity are the solutions of so r 1 = 1 and r 2 = 1. x 2 0 = F(r) = r(r 1) + r 1 = r 2 1

30 Solution (1 of 12) Assuming y(x) = a n x r+n, differentiating, and substituting into the ODE we obtain 0 = x 2 (r + n)(r + n 1)a n x r+n 2 + x (r + n)a n x r+n 1 = = + (x 2 1) a n x r+n (r + n)(r + n 1)a n x r+n + (r + n)a n x r+n a n x r+n + a n x r+n+2 [(r + n)(r + n 1) + (r + n) 1] a n x r+n + a n x r+n+2

31 Solution (2 of 12) 0 = = = [(r + n)(r + n 1) + (r + n) 1] a n x r+n + a n x r+n+2 [ ] (r + n) 2 1 a n x r+n + [ ] (r + n) 2 1 a n x r+n + = (r 2 1)a 0 x r + + n=2 ( (r + 1) 2 1 a n x r+n+2 a n 2 x r+n n=2 [ ] (r + n) 2 1 a n x r+n + ) a 1 x r+1 a n 2 x r+n n=2

32 Solution (3 of 12) 0 = = ( ) ( ) r 2 1 a 0 x r + (r + 1) 2 1 a 1 x r+1 + n=2 [ ] (r + n) 2 1 a n x r+n + a n 2 x r+n n=2 ( ) ( ) r 2 1 a 0 x r + (r + 1) 2 1 a 1 x r+1 ([ ] ) + (r + n) 2 1 a n + a n 2 x r+n n=2 Recurrence relation: a n (r) = a n 2(r) (r + n) 2 1 for n 2.

33 Solution (4 of 12) When r = r 1 = 1 then the expression ( ) (r + 1) 2 1 a 1 x r+1 = 0 when a 1 = 0.

34 Solution (4 of 12) When r = r 1 = 1 then the expression ( ) (r + 1) 2 1 a 1 x r+1 = 0 when a 1 = 0. Thus we will choose a 0 = 1 and a 1 = 0. a 1 = 0 implies a 2n+1 = 0 for all n N.

35 Solution (5 of 12) a 2 (1) = a 0 (2)(4) = 1 4(2!)(1!) a 4 (1) = a 2 (4)(6) = (3!)(2!) a 6 (1) = a 4 (6)(8) = (4!)(3!). a 2n (1) = y 1 (x) = x J 1 (x) = ( 1) n 4 n (n + 1)!(n!) [ 1 + n=1 ] ( 1) n x 2n 4 n (n + 1)!(n!) ( 1) n ( x ) 2n+1. n!(n + 1)! 2

36 Solution (6 of 12) Since r 1 r 2 = 2 the second solution will be of the form where y 2 (x) = ay 1 (x) ln x + x 1 [1 + ] c n ( 1)x n n=1 a = lim r 1 [(r ( 1))a 2(r)] = lim r 1 = lim r 1 = lim r 1 (r + 1)a 0 (r + 2) 2 1 (r + 1) r 2 + 4r + 3 (r + 1) (r + 3)(r + 1) 1 = lim r 1 r + 3 = 1 2

37 Solution (7 of 12) Recall: the coefficients of the power series c n (r 2 ) = d dr [(r r 2)a n (r)] r=r2, a 1 (±1) = 0 which implies a 1 ( 1) = 0, this in turn implies c 1 ( 1) = 0, from the recurrence relation then c 2n+1 ( 1) = 0 for n N.

38 Solution (8 of 12) From the recurrence relation a n (r) = a n 2(r) (r + n) 2 1 for n 2, and the fact that a 1 ( 1) = 0 then a 2n+1 ( 1) = 0 and a 2n+1 ( 1) = 0 for n N. a 2 (r) = 1 (r + 1)(r + 3) a 4 (r) = 1 (r + 1)(r + 3)(r + 3)(r + 5). a 2n (r) = ( 1) n (r + 1) (r + 2n 1)(r + 3) (r + 2n + 1)

39 Solution (9 of 12) Thus (r + 1)a 2 (r) = 1 r (r + 1)a 4 (r) = (r + 3)(r + 3)(r + 5) (r + 1)a 2n (r) =. ( 1) n (r + 3) (r + 2n 1)(r + 3) (r + 2n + 1)

40 Solution (10 of 12) Lemma If f (x) = (x α 1 ) β 1(x α 2 ) β 2 (x α n ) βn and x / {α 1, α 2,..., α n } then This implies f (x) f (x) = β 1 + β β n. x α 1 x α 2 x α n f (x) = f (x) [ β1 x α 1 + β 2 x α β n x α n ].

41 Solution (11 of 12) Applying the previous lemma we get [ d 2 dr [(r + 1)a 2n(r)] r= 1 = ( 1) n (r + 1)a 2n (r) r r + 5 ] r + 2n r + 2n + 1 ( 1) n+1 [ 2 = 2 2n 1 (n 1)!(n!) n ] 2n = c 2n ( 1) = ( 1) n+1 4 n (n 1)!(n!) r= 1 [ n n ( 1) n+1 4 n (n 1)!(n!) (H n 1 + H n ) ]

42 Solution (12 of 12) Finally, the second solution takes the form [ ] y 2 (x) = J 1 (x) ln x + 1 ( 1) n (H n 1 + H n ) ( x ) 2n 1. x n!(n 1)! 2 n=1 The Bessel function of the second kind of order one is defined to be Y 1 (x) = 2 π ( y 2(x) + (γ ln 2)J 1 (x)).

43 Graphs y x J 1 (x) Y 1 (x)

44 Useful Bessel Function Identities d dx [x n J n (x)] = x n J n 1 (x) d [ x n J n (x) ] dx = x n J n+1 (x) J n+1 (x) = 2n x J n(x) J n 1 (x)

45 Homework Read Section 5.7 Exercises: 5 8

Series Solutions Near a Regular Singular Point

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

Review for Exam 2. Review for Exam 2.

Review for Exam 2. Review for Exam 2. Review for Exam 2. 5 or 6 problems. No multiple choice questions. No notes, no books, no calculators. Problems similar to homeworks. Exam covers: Regular-singular points (5.5). Euler differential equation

More information

Chapter 5.8: Bessel s equation

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

An Introduction to Bessel Functions

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

Boundary Value Problems in Cylindrical Coordinates

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

Homogeneous Equations with Constant Coefficients

Homogeneous Equations with Constant Coefficients Homogeneous Equations with Constant Coefficients MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 General Second Order ODE Second order ODEs have the form

More information

ODE Homework Series Solutions Near an Ordinary Point, Part I 1. Seek power series solution of the equation. n(n 1)a n x n 2 = n=0

ODE Homework Series Solutions Near an Ordinary Point, Part I 1. Seek power series solution of the equation. n(n 1)a n x n 2 = n=0 ODE Homework 6 5.2. Series Solutions Near an Ordinary Point, Part I 1. Seek power series solution of the equation y + k 2 x 2 y = 0, k a constant about the the point x 0 = 0. Find the recurrence relation;

More information

Equations with regular-singular points (Sect. 5.5).

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

1 Series Solutions Near Regular Singular Points

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

MATH 312 Section 6.2: Series Solutions about Singular Points

MATH 312 Section 6.2: Series Solutions about Singular Points MATH 312 Section 6.2: Series Solutions about Singular Points Prof. Jonathan Duncan Walla Walla University Spring Quarter, 2008 Outline 1 Classifying Singular Points 2 The Method of Frobenius 3 Conclusions

More information

Math 2930 Worksheet Final Exam Review

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

Problem 1 Kaplan, p. 436: 2c,i

Problem 1 Kaplan, p. 436: 2c,i Mathematical Methods I - Fall 03 Homework Solutions Page Problem Kaplan, p 36: c,i Find the first three nonzero terms of the following Taylor series: c) ln + x) about x = 0 i) arctanh x about x = 0 The

More information

MATH 241 Practice Second Midterm Exam - Fall 2012

MATH 241 Practice Second Midterm Exam - Fall 2012 MATH 41 Practice Second Midterm Exam - Fall 1 1. Let f(x = { 1 x for x 1 for 1 x (a Compute the Fourier sine series of f(x. The Fourier sine series is b n sin where b n = f(x sin dx = 1 = (1 x cos = 4

More information

Chapter 4. Series Solutions. 4.1 Introduction to Power Series

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

Series Solution of Linear Ordinary Differential Equations

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

Series solutions of second order linear differential equations

Series solutions of second order linear differential equations Series solutions of second order linear differential equations We start with Definition 1. A function f of a complex variable z is called analytic at z = z 0 if there exists a convergent Taylor series

More information

1 A complete Fourier series solution

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

Review of Power Series

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

1 Solution to Homework 4

1 Solution to Homework 4 Solution to Homework Section. 5. The characteristic equation is r r + = (r )(r ) = 0 r = or r =. y(t) = c e t + c e t y = c e t + c e t. y(0) =, y (0) = c + c =, c + c = c =, c =. To find the maximum value

More information

Math Assignment 11

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

7.3 Singular points and the method of Frobenius

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

Power Series and Analytic Function

Power Series and Analytic Function Dr Mansoor Alshehri King Saud University MATH204-Differential Equations Center of Excellence in Learning and Teaching 1 / 21 Some Reviews of Power Series Differentiation and Integration of a Power Series

More information

The method of Fröbenius

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

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

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

Series Solutions Near an Ordinary Point

Series Solutions Near an Ordinary Point Series Solutions Near an Ordinary Point MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 Ordinary Points (1 of 2) Consider the second order linear homogeneous

More information

Two special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p

Two special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.

More information

Math 334 A1 Homework 3 (Due Nov. 5 5pm)

Math 334 A1 Homework 3 (Due Nov. 5 5pm) Math 334 A1 Homework 3 Due Nov. 5 5pm No Advanced or Challenge problems will appear in homeworks. Basic Problems Problem 1. 4.1 11 Verify that the given functions are solutions of the differential equation,

More information

Math 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions

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

Relevant sections from AMATH 351 Course Notes (Wainwright): Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls):

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

y + α x s y + β x t y = 0,

y + α x s y + β x t y = 0, 80 Chapter 5. Series Solutions of Second Order Linear Equations. Consider the differential equation y + α s y + β t y = 0, (i) where α = 0andβ = 0 are real numbers, and s and t are positive integers that

More information

Math 4263 Homework Set 1

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

Taylor and Maclaurin Series

Taylor and Maclaurin Series Taylor and Maclaurin Series MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Background We have seen that some power series converge. When they do, we can think of them as

More information

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

Heat Equation on Unbounded Intervals

Heat Equation on Unbounded Intervals Heat Equation on Unbounded Intervals MATH 467 Partial Differential Equations J. Robert Buchanan Department of Mathematics Fall 28 Objectives In this lesson we will learn about: the fundamental solution

More information

Separable Differential Equations

Separable Differential Equations Separable Differential Equations MATH 6 Calculus I J. Robert Buchanan Department of Mathematics Fall 207 Background We have previously solved differential equations of the forms: y (t) = k y(t) (exponential

More information

The Laplacian in Polar Coordinates

The Laplacian in Polar Coordinates The Laplacian in Polar Coordinates R. C. Trinity University Partial Differential Equations March 17, 15 To solve boundary value problems on circular regions, it is convenient to switch from rectangular

More information

Partial Differential Equations

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

The Method of Frobenius

The Method of Frobenius The Method of Frobenius Department of Mathematics IIT Guwahati If either p(x) or q(x) in y + p(x)y + q(x)y = 0 is not analytic near x 0, power series solutions valid near x 0 may or may not exist. If either

More information

Introduction to Differential Equations

Introduction to Differential Equations Chapter 1 Introduction to Differential Equations 1.1 Basic Terminology Most of the phenomena studied in the sciences and engineering involve processes that change with time. For example, it is well known

More information

Chapter 4: Partial differentiation

Chapter 4: Partial differentiation Chapter 4: Partial differentiation It is generally the case that derivatives are introduced in terms of functions of a single variable. For example, y = f (x), then dy dx = df dx = f. However, most of

More information

Autumn 2015 Practice Final. Time Limit: 1 hour, 50 minutes

Autumn 2015 Practice Final. Time Limit: 1 hour, 50 minutes Math 309 Autumn 2015 Practice Final December 2015 Time Limit: 1 hour, 50 minutes Name (Print): ID Number: This exam contains 9 pages (including this cover page) and 8 problems. Check to see if any pages

More information

Calculus and Parametric Equations

Calculus and Parametric Equations Calculus and Parametric Equations MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Introduction Given a pair a parametric equations x = f (t) y = g(t) for a t b we know how

More information

Ch 5.7: Series Solutions Near a Regular Singular Point, Part II

Ch 5.7: Series Solutions Near a Regular Singular Point, Part II Ch 5.7: Series Solutions Near a Regular Singular Point, Part II! Recall from Section 5.6 (Part I): The point x 0 = 0 is a regular singular point of with and corresponding Euler Equation! We assume solutions

More information

Lecture Notes for MAE 3100: Introduction to Applied Mathematics

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

Partial Differential Equations for Engineering Math 312, Fall 2012

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

Additional material: Linear Differential Equations

Additional material: Linear Differential Equations Chapter 5 Additional material: Linear Differential Equations 5.1 Introduction The material in this chapter is not formally part of the LTCC course. It is included for completeness as it contains proofs

More information

Power Series Solutions to the Bessel Equation

Power Series Solutions to the Bessel Equation Power Series Solutions to the Bessel Equation Department of Mathematics IIT Guwahati The Bessel equation The equation x 2 y + xy + (x 2 α 2 )y = 0, (1) where α is a non-negative constant, i.e α 0, is called

More information

Lecture 21 Power Series Method at Singular Points Frobenius Theory

Lecture 21 Power Series Method at Singular Points Frobenius Theory Lecture 1 Power Series Method at Singular Points Frobenius Theory 10/8/011 Review. The usual power series method, that is setting y = a n 0 ) n, breaks down if 0 is a singular point. Here breaks down means

More information

An Introduction to Bessel Functions

An Introduction to Bessel Functions An Introduction to R. C. Trinity University Partial Differential Equations March 25, 2014 Bessel s equation Given p 0, the ordinary differential equation x 2 y +xy +(x 2 p 2 )y = 0, x > 0 (1) is known

More information

ODE. Philippe Rukimbira. Department of Mathematics Florida International University PR (FIU) MAP / 92

ODE. Philippe Rukimbira. Department of Mathematics Florida International University PR (FIU) MAP / 92 ODE Philippe Rukimbira Department of Mathematics Florida International University PR (FIU) MAP 2302 1 / 92 4.4 The method of Variation of parameters 1. Second order differential equations (Normalized,

More information

Problem (p.613) Determine all solutions, if any, to the boundary value problem. y + 9y = 0; 0 < x < π, y(0) = 0, y (π) = 6,

Problem (p.613) Determine all solutions, if any, to the boundary value problem. y + 9y = 0; 0 < x < π, y(0) = 0, y (π) = 6, Problem 10.2.4 (p.613) Determine all solutions, if any, to the boundary value problem y + 9y = 0; 0 < x < π, y(0) = 0, y (π) = 6, by first finding a general solution to the differential equation. Solution.

More information

Elementary Differential Equations, Section 2 Prof. Loftin: Practice Test Problems for Test Find the radius of convergence of the power series

Elementary Differential Equations, Section 2 Prof. Loftin: Practice Test Problems for Test Find the radius of convergence of the power series Elementary Differential Equations, Section 2 Prof. Loftin: Practice Test Problems for Test 2 SOLUTIONS 1. Find the radius of convergence of the power series Show your work. x + x2 2 + x3 3 + x4 4 + + xn

More information

Linear Independence. MATH 322, Linear Algebra I. J. Robert Buchanan. Spring Department of Mathematics

Linear Independence. MATH 322, Linear Algebra I. J. Robert Buchanan. Spring Department of Mathematics Linear Independence MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Introduction Given a set of vectors {v 1, v 2,..., v r } and another vector v span{v 1, v 2,...,

More information

12d. Regular Singular Points

12d. 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 information

Vibrating-string problem

Vibrating-string problem EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u

More information

In what follows, we examine the two-dimensional wave equation, since it leads to some interesting and quite visualizable solutions.

In what follows, we examine the two-dimensional wave equation, since it leads to some interesting and quite visualizable solutions. ecture 22 igher-dimensional PDEs Relevant section of text: Chapter 7 We now examine some PDEs in higher dimensions, i.e., R 2 and R 3. In general, the heat and wave equations in higher dimensions are given

More information

Extrema of Functions of Several Variables

Extrema of Functions of Several Variables Extrema of Functions of Several Variables MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Background (1 of 3) In single-variable calculus there are three important results

More information

The Method of Frobenius

The Method of Frobenius The Method of Frobenius R. C. Trinity University Partial Differential Equations April 7, 2015 Motivating example Failure of the power series method Consider the ODE 2xy +y +y = 0. In standard form this

More information

Department of Mathematics

Department of Mathematics INDIAN INSTITUTE OF TECHNOLOGY, BOMBAY Department of Mathematics MA 04 - Complex Analysis & PDE s Solutions to Tutorial No.13 Q. 1 (T) Assuming that term-wise differentiation is permissible, show that

More information

swapneel/207

swapneel/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 information

FINAL EXAM, MATH 353 SUMMER I 2015

FINAL EXAM, MATH 353 SUMMER I 2015 FINAL EXAM, MATH 353 SUMMER I 25 9:am-2:pm, Thursday, June 25 I have neither given nor received any unauthorized help on this exam and I have conducted myself within the guidelines of the Duke Community

More information

MATH 391 Test 1 Fall, (1) (12 points each)compute the general solution of each of the following differential equations: = 4x 2y.

MATH 391 Test 1 Fall, (1) (12 points each)compute the general solution of each of the following differential equations: = 4x 2y. MATH 391 Test 1 Fall, 2018 (1) (12 points each)compute the general solution of each of the following differential equations: (a) (b) x dy dx + xy = x2 + y. (x + y) dy dx = 4x 2y. (c) yy + (y ) 2 = 0 (y

More information

Absolute Convergence and the Ratio Test

Absolute Convergence and the Ratio Test Absolute Convergence and the Ratio Test MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Bacground Remar: All previously covered tests for convergence/divergence apply only

More information

Special Instructions:

Special Instructions: Be sure that this examination has 20 pages including this cover The University of British Columbia Sessional Examinations - December 2016 Mathematics 257/316 Partial Differential Equations Closed book

More information

Spotlight on Laplace s Equation

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

Lecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018

Lecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018 Lecture 10 Partial derivatives Dan Nichols nichols@math.umass.edu MATH 233, Spring 2018 University of Massachusetts February 27, 2018 Last time: functions of two variables f(x, y) x and y are the independent

More information

Chapter 3 Second Order Linear Equations

Chapter 3 Second Order Linear Equations Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,

More information

Series Solutions of ODEs. Special Functions

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

Method of Frobenius. General Considerations. L. Nielsen, Ph.D. Dierential Equations, Fall Department of Mathematics, Creighton University

Method of Frobenius. General Considerations. L. Nielsen, Ph.D. Dierential Equations, Fall Department of Mathematics, Creighton University Method of Frobenius General Considerations L. Nielsen, Ph.D. Department of Mathematics, Creighton University Dierential Equations, Fall 2008 Outline 1 The Dierential Equation and Assumptions 2 3 Main Theorem

More information

Final Exam May 4, 2016

Final Exam May 4, 2016 1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.

More information

2 Series Solutions near a Regular Singular Point

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

Mathematics of Physics and Engineering II: Homework problems

Mathematics of Physics and Engineering II: Homework problems Mathematics of Physics and Engineering II: Homework problems Homework. Problem. Consider four points in R 3 : P (,, ), Q(,, 2), R(,, ), S( + a,, 2a), where a is a real number. () Compute the coordinates

More information

Section 5.2 Series Solution Near Ordinary Point

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

Math 21B - Homework Set 8

Math 21B - Homework Set 8 Math B - Homework Set 8 Section 8.:. t cos t dt Let u t, du t dt and v sin t, dv cos t dt Let u t, du dt and v cos t, dv sin t dt t cos t dt u v v du t sin t t sin t dt [ t sin t u v ] v du [ ] t sin t

More information

Math 180, Final Exam, Fall 2012 Problem 1 Solution

Math 180, Final Exam, Fall 2012 Problem 1 Solution Math 80, Final Exam, Fall 0 Problem Solution. Find the derivatives of the following functions: (a) ln(ln(x)) (b) x 6 + sin(x) e x (c) tan(x ) + cot(x ) (a) We evaluate the derivative using the Chain Rule.

More information

Circular Membranes. Farlow, Lesson 30. November 21, Department of Mathematical Sciences Florida Atlantic University Boca Raton, FL 33431

Circular Membranes. Farlow, Lesson 30. November 21, Department of Mathematical Sciences Florida Atlantic University Boca Raton, FL 33431 The Problem Polar coordinates Solving the Problem by Separation of Variables Circular Membranes Farlow, Lesson 30 Department of Mathematical Sciences Florida Atlantic University Boca Raton, FL 33431 November

More information

This practice exam is intended to help you prepare for the final exam for MTH 142 Calculus II.

This practice exam is intended to help you prepare for the final exam for MTH 142 Calculus II. MTH 142 Practice Exam Chapters 9-11 Calculus II With Analytic Geometry Fall 2011 - University of Rhode Island This practice exam is intended to help you prepare for the final exam for MTH 142 Calculus

More information

Ma 221 Final Exam Solutions 5/14/13

Ma 221 Final Exam Solutions 5/14/13 Ma 221 Final Exam Solutions 5/14/13 1. Solve (a) (8 pts) Solution: The equation is separable. dy dx exy y 1 y0 0 y 1e y dy e x dx y 1e y dy e x dx ye y e y dy e x dx ye y e y e y e x c The last step comes

More information

UNIVERSITY OF MANITOBA

UNIVERSITY OF MANITOBA Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic

More information

M412 Assignment 5 Solutions

M412 Assignment 5 Solutions M41 Assignment 5 Solutions 1. [1 pts] Haberman.5.1 (a). Solution. Setting u(x, y) =X(x)Y(y), we find that X and Y satisfy X + λx = Y λy =, X () = X () = Y() =. In this case, we find from the X equation

More information

Welcome to Math 257/316 - Partial Differential Equations

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

Differential Equations

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

Stokes Theorem. MATH 311, Calculus III. J. Robert Buchanan. Summer Department of Mathematics. J. Robert Buchanan Stokes Theorem

Stokes Theorem. MATH 311, Calculus III. J. Robert Buchanan. Summer Department of Mathematics. J. Robert Buchanan Stokes Theorem tokes Theorem MATH 311, alculus III J. Robert Buchanan Department of Mathematics ummer 2011 Background (1 of 2) Recall: Green s Theorem, M(x, y) dx + N(x, y) dy = R ( N x M ) da y where is a piecewise

More information

MATH 425, FINAL EXAM SOLUTIONS

MATH 425, FINAL EXAM SOLUTIONS MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u

More information

Fall Math 3410 Name (Print): Solution KEY Practice Exam 2 - November 4 Time Limit: 50 Minutes

Fall Math 3410 Name (Print): Solution KEY Practice Exam 2 - November 4 Time Limit: 50 Minutes Fall 206 - Math 340 Name (Print): Solution KEY Practice Exam 2 - November 4 Time Limit: 50 Minutes This exam contains pages (including this cover page) and 5 problems. Check to see if any pages are missing.

More information

Chain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics

Chain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics 3.33pt Chain Rule MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Spring 2019 Single Variable Chain Rule Suppose y = g(x) and z = f (y) then dz dx = d (f (g(x))) dx = f (g(x))g (x)

More information

Examples of the Fourier Theorem (Sect. 10.3). The Fourier Theorem: Continuous case.

Examples of the Fourier Theorem (Sect. 10.3). The Fourier Theorem: Continuous case. s of the Fourier Theorem (Sect. 1.3. The Fourier Theorem: Continuous case. : Using the Fourier Theorem. The Fourier Theorem: Piecewise continuous case. : Using the Fourier Theorem. The Fourier Theorem:

More information

ENGI 9420 Lecture Notes 1 - ODEs Page 1.01

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

Power Series Solutions We use power series to solve second order differential equations

Power Series Solutions We use power series to solve second order differential equations Objectives Power Series Solutions We use power series to solve second order differential equations We use power series expansions to find solutions to second order, linear, variable coefficient equations

More information

23 Elements of analytic ODE theory. Bessel s functions

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

5.4 Bessel s Equation. Bessel Functions

5.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 information

Differential Equations Practice: Euler Equations & Regular Singular Points Page 1

Differential Equations Practice: Euler Equations & Regular Singular Points Page 1 Differential Equations Practice: Euler Equations & Regular Singular Points Page 1 Questions Eample (5.4.1) Determine the solution to the differential equation y + 4y + y = 0 that is valid in any interval

More information

Problem 1. Use a line integral to find the plane area enclosed by the curve C: r = a cos 3 t i + b sin 3 t j (0 t 2π). Solution: We assume a > b > 0.

Problem 1. Use a line integral to find the plane area enclosed by the curve C: r = a cos 3 t i + b sin 3 t j (0 t 2π). Solution: We assume a > b > 0. MATH 64: FINAL EXAM olutions Problem 1. Use a line integral to find the plane area enclosed by the curve C: r = a cos 3 t i + b sin 3 t j ( t π). olution: We assume a > b >. A = 1 π (xy yx )dt = 3ab π

More information

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (I)

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (I) Laplace s Equation in Cylindrical Coordinates and Bessel s Equation I) 1 Solution by separation of variables Laplace s equation is a key equation in Mathematical Physics. Several phenomena involving scalar

More information

Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series

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

Partial Derivatives October 2013

Partial Derivatives October 2013 Partial Derivatives 14.3 02 October 2013 Derivative in one variable. Recall for a function of one variable, f (a) = lim h 0 f (a + h) f (a) h slope f (a + h) f (a) h a a + h Partial derivatives. For a

More information

ODE Homework Solutions of Linear Homogeneous Equations; the Wronskian

ODE Homework Solutions of Linear Homogeneous Equations; the Wronskian ODE Homework 3 3.. Solutions of Linear Homogeneous Equations; the Wronskian 1. Verify that the functions y 1 (t = e t and y (t = te t are solutions of the differential equation y y + y = 0 Do they constitute

More information

4 Power Series Solutions: Frobenius Method

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

Series solutions to a second order linear differential equation with regular singular points

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

Analyse 3 NA, FINAL EXAM. * Monday, January 8, 2018, *

Analyse 3 NA, FINAL EXAM. * Monday, January 8, 2018, * Analyse 3 NA, FINAL EXAM * Monday, January 8, 08, 4.00 7.00 * Motivate each answer with a computation or explanation. The maximum amount of points for this exam is 00. No calculators!. (Holomorphic functions)

More information