Section 5.3, Exercise 22

Similar documents
Power Series Solutions to the Legendre Equation

LEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.

Legendre s Equation. PHYS Southern Illinois University. October 18, 2016

Power Series Solutions to the Legendre Equation

Background and Definitions...2. Legendre s Equation, Functions and Polynomials...4 Legendre s Associated Equation and Functions...

MAT137 Calculus! Lecture 6

Math 241 Final Exam Spring 2013

Mathematics 22: Lecture 19

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

The polar coordinates

MATH 1372, SECTION 33, MIDTERM 3 REVIEW ANSWERS

Math 214 Spring problem set (a) Consider these two first order equations. (I) dy dx = x + 1 dy

Mathematics 324 Riemann Zeta Function August 5, 2005

d 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.

Electromagnetism HW 1 math review

2.2 The derivative as a Function

Bessel s and legendre s equations

1 The Derivative and Differrentiability

Limits and Continuous Functions. 2.2 Introduction to Limits. We first interpret limits loosely. We write. lim f(x) = L

Review of Power Series

Worksheet 7, Math 10560

Exercises from other sources REAL NUMBERS 2,...,

PHYS 404 Lecture 1: Legendre Functions

Spotlight on the Extended Method of Frobenius

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

h(x) lim H(x) = lim Since h is nondecreasing then h(x) 0 for all x, and if h is discontinuous at a point x then H(x) > 0. Denote

Math Assignment 11

Generating Functions (Revised Edition)

Differentiation - Quick Review From Calculus

2 2 + x =

Definitions & Theorems

Math 4381 / 6378 Symmetry Analysis

LECTURE 10: REVIEW OF POWER SERIES. 1. Motivation

Ma 530 Power Series II

DIFFERENTIATION RULES

1 + lim. n n+1. f(x) = x + 1, x 1. and we check that f is increasing, instead. Using the quotient rule, we easily find that. 1 (x + 1) 1 x (x + 1) 2 =

3.4 Introduction to power series

Differentiation of x n.

This ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0

e x3 dx dy. 0 y x 2, 0 x 1.

14 EE 2402 Engineering Mathematics III Solutions to Tutorial 3 1. For n =0; 1; 2; 3; 4; 5 verify that P n (x) is a solution of Legendre's equation wit

AP Calculus (BC) Chapter 9 Test No Calculator Section Name: Date: Period:

Math 308 Exam I Practice Problems

Special Functions. SMS 2308: Mathematical Methods

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

Chapter 9. Worked-Out Solutions. Check. Chapter 9 Maintaining Mathematical Proficiency (p. 477) y = 2 x 2. y = 1. = (x + 5) 2 4 =?

n f(k) k=1 means to evaluate the function f(k) at k = 1, 2,..., n and add up the results. In other words: n f(k) = f(1) + f(2) f(n). 1 = 2n 2.

Review for Exam 2. Review for Exam 2.

MATH115. Sequences and Infinite Series. Paolo Lorenzo Bautista. June 29, De La Salle University. PLBautista (DLSU) MATH115 June 29, / 16

Power Series and Analytic Function

e x = 1 + x + x2 2! + x3 If the function f(x) can be written as a power series on an interval I, then the power series is of the form

cse547, math547 DISCRETE MATHEMATICS Professor Anita Wasilewska

MA 242 Review Exponential and Log Functions Notes for today s class can be found at

TAYLOR AND MACLAURIN SERIES

Chapter 4 Integration

Notes on Special Functions

AP Calculus Testbank (Chapter 9) (Mr. Surowski)

CHAIN RULE: DAY 2 WITH TRIG FUNCTIONS. Section 2.4A Calculus AP/Dual, Revised /30/2018 1:44 AM 2.4A: Chain Rule Day 2 1

Definition of Derivative

Consider an infinite row of dominoes, labeled by 1, 2, 3,, where each domino is standing up. What should one do to knock over all dominoes?

Problem 1 (Equations with the dependent variable missing) By means of the substitutions. v = dy dt, dv

Math Homework 2

Homework Solutions Week of November = 1(1) + 0(x 1) + 0(x 1) 2 x = 1(1) + 1(x 1) + 0(x 1) 2 x 2 = 1(1) + 2(x 1) + 1(x 1) 2

Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0

DIFFERENTIATION. 3.1 Approximate Value and Error (page 151)

Lecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain.

Waves on 2 and 3 dimensional domains

Math 4263 Homework Set 1

Bernoulli Numbers and their Applications

Foundations of Calculus. November 18, 2014

Math Review for Exam Answer each of the following questions as either True or False. Circle the correct answer.

Connection to Laplacian in spherical coordinates (Chapter 13)

1.2. Characteristics of Polynomial Functions. What are the key features of the graphs of polynomial functions?

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

MATH 2413 TEST ON CHAPTER 4 ANSWER ALL QUESTIONS. TIME 1.5 HRS.

Integrals in cylindrical, spherical coordinates (Sect. 15.7)

3! + 4! + Binomial series: if α is a nonnegative integer, the series terminates. Otherwise, the series converges if x < 1 but diverges if x > 1.

Chapter 4 Notes, Calculus I with Precalculus 3e Larson/Edwards

Introduction to Spherical Harmonics

Problem 1 In each of the following problems find the general solution of the given differential

Convergence Concepts of Random Variables and Functions

2017 HSC Mathematics Extension 1 Marking Guidelines

Chapter 4: Partial differentiation

ENGI 9420 Engineering Analysis Solutions to Additional Exercises

AP Calculus BC. Chapter 3: Derivatives 3.3: Rules for Differentiation

Advanced Calculus Math 127B, Winter 2005 Solutions: Final. nx2 1 + n 2 x, g n(x) = n2 x

1. Let A R be a nonempty set that is bounded from above, and let a be the least upper bound of A. Show that there exists a sequence {a n } n N

Midterm Preparation Problems

MATH1231 CALCULUS. Session II Dr John Roberts (based on notes of A./Prof. Bruce Henry) Red Center Room 3065

Chapter 5.3: Series solution near an ordinary point

Copyright 2010 Pearson Education, Inc. Publishing as Prentice Hall.

Chapter 2 Formulas and Definitions:

z = 1 2 x 3 4 y + 3 y dt

Boundary Value Problems in Cylindrical Coordinates

Chapter 8 Integration Techniques and Improper Integrals

Appendix A Vector Analysis

Taylor and Maclaurin Series

Math 651 Introduction to Numerical Analysis I Fall SOLUTIONS: Homework Set 1

TABLE OF CONTENTS 2 CHAPTER 1

ENGI 9420 Lecture Notes 1 - ODEs Page 1.01

Transcription:

The Legendre equation is where α is a constant. Section 5.3, Exercise 22 (1 x 2 ) 2x + α(α + 1) 0 Determine two linearl independent solutions in powers of x for x < 1. Assume (x) a n x n and substitute this into Legendre s equation. 0 (1 x 2 ) n(n 1)a n x n 2 2x na n x n + α(α + 1) a n x n n2 n(n 1)a n x n 2 n2 n1 n(n 1)a n x n n2 (n + 2)(n + 1)a n+2 x n n1 n(n 1)a n x n 2na n x n + α(α + 1) a n x n [(n + 2)(n + 1)a n+2 (n(n + 1) α(α + 1))a n ] x n We obtain the recurrence relation: a n+2 (n(n + 1) α(α + 1))a n. (n + 2)(n + 1) 2na n x n + α(α + 1) a n x n If we let a 0 1 and a 1 0 then we can calculate the following coefficients. α(α + 1) a 2 2! a 4 (2(3) α(α + 1))a 2 (4)(3) α(α + 1)[α2 + α 6] 4! α(α 2)(α + 1)(α + 3) 4! a 6 (4(5) α(α + 1))a 4 (6)(5) α(α 2)(α + 1)(α + 3)(α2 + α 20) 6! α(α 2)(α 4)(α + 1)(α + 3)(α + 5) 6!

Hence we have one solution. a 2m () mα (α 2m + 2)(α + 1) (α + 2m 1) (2m)! α(α + 1) 1 (x) 1 x 2 α(α 2)(α + 1)(α + 3) + x 4 2! 4! α (α 2m + 2)(α + 1) (α + 2m 1) + x 2m. (2m)! m3 Using the Ratio Test for absolute convergence we see that lim a 2(m+1) x 2(m+1) m a 2m x 2m lim m lim m x 2. α (α 2m)(α+1) (α+2m+1) x 2m+2 (2m+2)! α (α 2m+2)(α+1) (α+2m) (2m)! (α 2m)(α + 2m + 1) x 2 (2m + 2)(2m + 1) Thus the power series converges for x < 1. If we let a 0 0 and a 1 1 then we can calculate the following coefficients. 1(2) α(α + 1) a 3 3! (α 1)(α + 2) 3! a 5 (3(4) α(α + 1))a 3 (5)(4) (α 1)(α + 2)(α2 + α 12) 5! (α 1)(α 3)(α + 2)(α + 4) 5! a 7 (5(6) α(α + 1))a 5 (7)(6) (α 1)(α 3)(α + 2)(α + 4)(α2 + α 30) 7! (α 1)(α 3)(α 5)(α + 2)(α + 4)(α + 6) 7! x 2m. a 2m+1 () m(α 1) (α 2m + 1)(α + 2) (α + 2m) (2m + 1)!

Hence we have another solution (α 1)(α + 2) 2 (x) x x 3 (α 1)(α 3)(α + 2)(α + 4) + x 5 3! 5! + () m(α 1) (α 2m + 1)(α + 2) (α + 2m) x 2m+1. (2m + 1)! m3 Using the Ratio Test for absolute convergence we see that lim a 2(m+1)+1 x 2(m+1)+1 (α) (α 2m)(α+2) (α+2m+2) x 2m+3 (2m+3)! m a 2m+1 x 2m+1 lim m (α) (α 2m+1)(α+2) (α+2m) x (2m+1)! 2m+1 lim (α 2m 1)(α + 2m + 2) m x 2 (2m + 3)(2m + 2) x 2. Thus the power series converges for x < 1. We ma verif the solutions are linearl independent using the Wronskian. Section 5.3, Exercise 23 W( 1, 2 )(0) 1 (0) 2(0) 1(0) 2 (0) 1 0 Show that if α is zero or a positive even integer 2n, the series solution 1 reduces to a polnomial of degree 2n containing onl even powers of x. Find the polnomials corresponding to α 0, 2, and 4. Show that if α is a positive odd integer 2n+1, the series solution 2 reduces to a polnomial of degree 2n + 1 containing onl odd powers of x. Find the polnomials corresponding to α 1, 3, and 5. If α 0 then 1 (x) 1 which is a polnomial of degree zero. If α 2n for some n N then according to the recurrence relation a 2n+2k 0 for all k N which implies that α(α + 1) 1 (x) 1 x 2 + 2! + 2n m3 α(α 2)(α + 1)(α + 3) x 4 4! α (α 2m + 2)(α + 1) (α + 2m 1) x 2m, (2m)! a polnomial of degree 2n containing onl even powers of x. α 0: (x) 1 α 2: (x) 1 3x 2 α 4: (x) 1 10x 2 + 35 3 x4

If α 1 then 2 (x) x which is a polnomial of degree 1 containing onl odd powers of x. If α 2n + 1 for some n N then according to the recurrence relation a 2n+1+2k 0 for all k N which implies that (α 1)(α + 2) 2 (x) x x 3 (α 1)(α 3)(α + 2)(α + 4) + x 5 3! 5! 2n+1 + () m(α 1) (α 2m + 1)(α + 2) (α + 2m) x 2m+1, (2m + 1)! m3 a polnomial of degree 2n + 1 containing onl odd powers of x. α 1: α 3: Section 5.3, Exercise 24 (x) x (x) x 5 3 x3 α 5: (x) x 14 3 x3 + 21 5 x5 The Legendre polnomial P n (x) is defined as the polnomial solution of the Legendre equation with α n that also satisfies the condition P n (1) 1. (a) Using the results of Problem 23, find the Legendre polnomials P 0 (x),..., P 5 (x). We can easil see that P 0 (x) 1 and P 1 (x) x. Since then a /2 and we have P 2 (x) a(1 3x 2 ) P 2 (1) a(1 3) 1 P 2 (x) 1 2 (1 3x2 ). In the same wa we see that P 3 (x) a (x 53 ) x3 ( P 3 (1) a 1 5 ) 1 3 a 3 2 P 4 (x) a (1 10x 2 + 353 ) x4 ( P 4 (1) a 1 10 + 35 ) 1 3

To summarize a 3 8( P 5 (x) a x 14 3 x3 + 21 ) 5 x5 ( P 5 (1) a 1 14 3 + 21 ) 1 5 a 15 8. Legendre Polnomials P 0 (x) 1 P 1 (x) x P 2 (x) 1 2 (3x2 1) P 3 (x) 1 2 (5x3 3x) P 4 (x) 1 8 (35x4 30x 2 + 3) P 5 (x) 1 8 (63x5 70x 3 + 15x) (b) Plot the graphs of P 0 (x),..., P 5 (x) for x 1. P 0 (x) 1 2.0 1.5 0.5 0.5 0.5 x

P 1 (x) x 0.5 0.5 0.5 x 0.5 P 2 (x) 1 2 (3x2 1) 0.8 0.6 0.4 0.2 0.5 0.5 x 0.2 0.4 P 3 (x) 1 2 (5x3 3x) 0.5 0.5 0.5 x 0.5

P 4 (x) 1 8 (35x4 30x 2 + 3) 0.8 0.6 0.4 0.2 0.5 0.5 x 0.2 0.4 P 5 (x) 1 8 (63x5 70x 3 + 15x) 0.5 0.5 0.5 x 0.5 (c) Find the zeros of P 0 (x),..., P 5 (x). P 0 (x) 1 has no zeros. P 1 (x) x has a zero at x 0. P 2 (x) 1 2 (3x2 1) has zeros at x ± 1. 3 3 P 3 (x) 1 2 (5x3 3x) has zeros at x 0 and at x ± 5. P 4 (x) 1 8 (35x4 30x 2 15 ± 30 + 3) has zeros at x ±. 35 1 8 (63x5 70x 3 + 15x) has zeros at x 0 and at x ± 35 ± 2 70. 63

Section 5.3, Exercise 25 It can be shown that the general formula for P n (x) is P n (x) 1 2 n [n/2] k0 () k (2n 2k)! k!(n k)!(n 2k)! xn 2k, where [n/2] denotes the greatest integer less than or equal to n/2. B observing the form of P n (x) for n even and n odd, show that P n () () n. B exercise 23 we know P n (x) is a polnomial of even (odd) degree when n is even (odd) and contains onl even (odd) powers of x. Thus P n (x) is an even (odd) function when n is even (odd). B exercise 24, P n (1) 1, thus when n is even P n () 1 () n and when n is odd P n () P n (1) () n. Section 5.3, Exercise 26 The Legendre polnomials pla an important role in mathematical phsics. For example in solving Laplace s equation (the potential equation) in spherical coordinates, we encounter the equation d 2 F(ϕ) dϕ 2 + cotϕ df(ϕ) dϕ + n(n + 1)F(ϕ) 0, 0 < ϕ < π, where n is a positive integer. Show that the change of variable x cosϕ leads to the Legendre equation with α n for f(x) F(arccosx). Using the chain rule for derivatives we see that df(ϕ) dϕ df(ϕ) dx dx dϕ df(arccosx) ( sin ϕ) df dx dx ( sin ϕ) sin ϕ. Differentiating once more we have d 2 F(ϕ) dϕ 2 d dϕ [ sin ϕ ] cos ϕ + sin 2 ϕ. Substituting into the equation above we have 0 0 cosϕ + sin 2 ϕ + cot ϕ ( sin ϕ ) + n(n + 1) sin 2 ϕ 2 cosϕ + n(n + 1) (1 cos 2 ϕ) 2 cosϕ + n(n + 1) (1 x 2 ) 2x + n(n + 1).

Section 5.3, Exercise 27 Show that for n 0, 1, 2, 3, the corresponding Legendre polnomial is given b P n (x) 1 d n 1) n. 2 n n! dx n(x2 This formula, known as Rodrigues formula, is true for all positive integers n. For the case n 0 P 0 (x) 1 d 0 1) 0 1. 2 0 0! dx 0(x2 For the case n 1 P 1 (x) 1 d 2 1 1! dx (x2 1) 1 (2x) x. 2 For the case n 2 P 2 (x) 1 d 2 1) 2 1 d 2 2x 2 + 1) 1 d 2 2 2! dx 2(x2 8 dx 2(x4 8 dx (4x3 4x) 1 8 (12x2 4) 3 2 x2 1 2. For the case n 3 P 3 (x) 1 d 3 1) 3 1 d 3 3x 4 + 3x 2 1) 1 d 2 12x 3 + 6x) 2 3 3! dx 3(x2 48 dx 3(x6 48 dx 2(6x5 1 d 48 dx (30x4 36x 2 + 6) 1 48 (120x3 72x) 5 2 x3 3 2 x. Section 5.3, Exercise 28 Show that the Legendre equation can also be written as Then it follows that [(1 x 2 ) ] α(α + 1). [(1 x 2 )P n (x)] n(n + 1)P n (x) and [(1 x 2 )P m (x)] m(m + 1)P m (x). B multipling the first equation b P m (x) and the second equation b P n (x), integrating b parts, and then subtracting one equation from the other, show that P n (x)p m (x) dx 0 if n m. This propert of the Legendre polnomials is known as the orthogonalit propert. If m n, it can be show that the value of the preceding integral is 2/(2n + 1).

Starting with the Legendre equation 0 (1 x 2 ) 2x + α(α + 1) (1 x 2 )( ) + (1 x 2 ) + α(α + 1) α(α + 1) [(1 x 2 ) ] Thus when (x) P n (x) we have (product rule). [(1 x 2 )P n (x)] n(n + 1)P n (x). Suppose we multipl this equation b P m (x) where m n, then [(1 x 2 )P n (x)] P m (x) n(n + 1)P n (x)p m (x) [(1 x 2 )P n(x)] P m (x) dx n(n + 1) P n (x)p m (x) dx We can appl integration b parts to the integral on the left-hand side of this equation. Thus u P m (x) du P m (x) v (1 x 2 )P n(x) dv [(1 x 2 )P n (x)] dx [(1 x 2 )P n (x)] P m (x) dx P m (x)[(1 x 2 )P n (x)] 1 Similarl when (x) P m (x) we have (1 x 2 )P n(x)p m(x) dx. [(1 x 2 )P m(x)] m(m + 1)P m (x). Suppose we multipl this equation b P n (x), then [(1 x 2 )P m(x)] P n (x) m(m + 1)P m (x)p n (x) [(1 x 2 )P m (x)] P n (x) dx m(m + 1) (1 x 2 )P n (x)p m (x) dx P m (x)p n (x) dx We can appl integration b parts to the integral on the left-hand side of this equation. Thus u P n (x) du P n(x) v (1 x 2 )P m (x) dv [(1 x 2 )P m(x)] dx [(1 x 2 )P m(x)] P n (x) dx P n (x)[(1 x 2 )P m(x)] 1 (1 x 2 )P m(x)p n(x) dx (1 x 2 )P m (x)p n (x) dx.

Combining the equations we see that m(m + 1) [m(m + 1) n(n + 1)] (1 x 2 )P m (x)p n (x) dx (1 x 2 )P n (x)p m (x) dx P m (x)p n (x) dx n(n + 1) P m (x)p n (x) dx 0. Thus either P m (x)p n (x) dx 0 which is the orthogonalit propert, or 0 [m(m + 1) n(n + 1)] (m n)(m + n + 1) P n (x)p m (x) dx Since m, n N and we have assumed m n then this equation is not satisfied. Section 5.3, Exercise 29 Given a polnomial f of degree n, it is possible to express f as a linear combination of P 0, P 1, P 2,..., P n : n f(x) a k P k (x). Using the result of Problem 28, show that then Suppose 2m + 1 2 a k 2k + 1 2 f(x) k0 f(x)p k (x) dx. n a k P k (x). k0 f(x)p m (x) f(x)p m (x) dx n a k P k (x)p m (x) k0 n a k P k (x)p m (x) dx k0 f(x)p m (x) dx a m 2 2m + 1 a m