1 Solutions in cylindrical coordinates: Bessel functions
|
|
- Malcolm Blake
- 6 years ago
- Views:
Transcription
1 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 is: 1 Φ Φ + Φ Separate variables: Let Φ () () () Then we find: The last term is a function of only, while the sum of the first two terms is a function of and only. Thus we take each part to be a constant called Then and the solutions are ± This is the appropriate solution outside of a charge distribution, say above a plane, (Φ as ± ) or inside a cylinder with grounded walls and non-zero potential on one end. The remaining equation is: 1 Now multiply through by : Here the last term is a function of only and the first two terms are functions of only. Again we often want a solution that is periodic with period so we choose a negative separation constant: ± Finally we have the equation for the function of : + 1
2 To see that this equation is of Sturm-Liouville form, divide through by : + (1) Now we have a Sturm-Liouville equation with () () eigenvalue and weighting function () Equation (1) is Bessel s equation. The solutions are orthogonal functions. Since () we do not need to specify any boundary condition at if our range is as is frequently the case. We do need a boundary condition at It is simpler and more elegant to solve Bessel s equation if we change to the dimensionless variable Then: + + The equation has a singular point at So we look for a series solution of the Frobenius type (cf Chapter ): X ( + ) + 1 ( + ) + 1 Then the equation becomes: ( + ) X The indicial equation is given by the coefficient of 1 : ± + 1 Thus one of the solutions (with ) is analytic at and one (with ) is not. To find the recursion relation, look at the + 1 power of : ( + ) + and so ( + ) ( ± )
3 Let s look first at the solution with + Wecanstepdowntofind each If we start the series with then will always be even, and 1 1 ( +) ( ) ( +) 4 ( 1) 3 3 ( 1) ( ) 3 ( + )( + 1) ( + ) 6 ( 1)! The usual convention is to take Then 1 ( + )( + 1) ( +1) 1 Γ ( +1) 1 ( 1) 1 Γ ( +1)! ( + )( + 1) ( +1) ( 1)!Γ ( + +1) + (3) and the solution is the Bessel function: () () ( 1) + (4)!Γ ( + +1)³ The function () has only even powers if is an even integer and only odd powers if is an odd integer. The series converges for all values of Let s see what the second solution looks like. With the recursion relation is: (5) ( ) where again Now if is an integer we will not be able to determine because the recursion relation blows up. One solution to this dilemma is to start the series with Then we can find the succeeding : (+) ( 1)+ ( + ) ( )+ 4 ( + )( + 1) ( 1) ( 1) Γ ( +1)!Γ ( + +1) which is the same recursion relation we had before. (Compare the equation above with equation (3). Thus we do not get a linearly independent solution this way. (This dilemma does not arise if the separation constant is taken to be with non-integer. In that case the second recursion relation provides aseries () that is linearly independent of the first.) Indeed we find: ( 1) Γ ( +1) ()!Γ ( + +1) (+) ( 1) Γ ( +1) ³ +!Γ ( + +1) 3
4 and if we choose ( 1) Γ ( +1) then () ( 1) () (6) With this choice () is a continuous function of (Notice that we can also express the series using equation (3) for the coefficients, with and + and where for ) We still have to determine the second, linearly independent solution of the Bessel equation. We can find it by taking the limit as of a linear combination of and knownastheneumannfunction () : () ()cos () lim () lim sin + ()cos( + ) (+) () lim sin ( + ) + ()(cos cos sin sin ) (+) () lim sin cos +cos sin + ()( 1) cos (+) () lim ( 1) sin Now we expand the functions to first order in We use a Taylor series for the Bessel functions. Note that appears in the index, not the argument, so we have to differentiate with respect to. + ()( 1) (+) () () lim ( 1) ( 1) lim ( 1) + + Using relation (6), we have: () 1 ( 1) The derivative has a logarithmic term: Ã! X ( 1)!Γ ( + +1)³ and ( 1) +!Γ ( + +1)³ ln ln ln ln 4 ( 1)!Γ ( + +1)³
5 and so has a term containing ln This term diverges as provided that () is not zero, i.e. for The function () also diverges as for 6 because it contains negative powers of (The series for starts with a term ) is finite as because goes to zero sufficiently fast. Two additional functions called Hankel functions are defined as linear combinations of and : (1) () ()+ () (7) and () () () () (8) Compare the relation between sine, cosine, and exponential: ± cos ± sin 1. Properties of the functions The Bessel functions (s) are well behaved both at the origin and as They have infinitely many zeroes. All of them, except for are zero at The first few functions are shown in the figure. J x The first three Bessel functions. 1 (red) and For small values of the argument, we may approximate the function with the first term in the series: 1 () for 1 (9) Γ ( +1)³ The Neumann functions are not well behaved at has a logarithmic singularity, and for diverges as an inverse power of : () ln for 1 ( 1)! () for 1 (1) 5
6 For large values of the argument, both and oscillate: they are like damped cosine or sine functions: r ³ () cos for À 1 (11) 4 r ³ () sin for À 1 (1) 4 and thus the Hankel functions are like complex exponentials: r h ³ (1) exp ± i for À 1 (13) 4 Notice that if 1 the large argument expansions apply for À rather than the usual À Relations between the functions As we found with the Legendre functions, we can determine a set of recursion relations that relate successive () For example: () +1 () (14) which is valid for In particular, with we obtain: 1 () () (15) ( ()) 1 () (16) (17) and similarly, by subtracting instead of adding, we can find +1 1 Thesamerelationsholdforthes andthes. (18) 1.4 Orthogonality of the Since the Bessel equation is of Sturm-Liouville form, the Bessel functions are orthogonal if we demand that they satisfy boundary conditions of the form (??). In particular, suppose the region of interest is to and the boundary conditions are () We do not need a boundary condition at because the function () is zero there. Then the eigenvalues are 6
7 where is the th zero of (The zeros are tabulated in standard references such as Abramowitz and Stegun. Also programs such as Mathematica and Maple can compute them.) Then Z 1.5 Solving a potential problem. [ ( )] [ ( )] (19) Example. A cylinder of radius and height has its curved surface and its bottom grounded. The top surface has potential What is the potential inside the cylinder? The potential has no dependence on and so only eigenfunctions with contribute. The potential is zero at sothesolutionweneedis ()with eigenvalues chosen to make () Thus the eigenvalues are given by where are the zeros of the function The remaining function of must be zero at so we choose the hyperbolic sine. Thus the potential is: Φ ( ) ( )sinh( ) Nowweevaluatethisat : 1 Φ ( ) ( )sinh( ) 1 Next we make use of the orthogonality of the Bessel functions. Multiply both sides by ( )and integrate from to (Note here that the weight function () This is the first time we have a weight function that is not 1.) Only one term in the sum, with survives the integration. Z ( ) Z ( ) ( )sinh( ) 1 Z [ ( )] sinh ( ) [ ( )] sinh ( ) To evaluate the left hand side, we use equation (16) with : Z Z () 1 ( 1 ()) 1 1 () 1 () 7
8 So 1 ( ) [ ( )] sinh ( ) 1 ( )sinh( ) where we used the result from equation (15) that 1 Finally our solution is: ( ) sinh ( ) Φ 1 ( ) sinh ( ) 1 The first two zeros of are: 1.448, 5.51, and thus the first two terms in the potential are: Ã 448 sinh 448 Φ (448) sinh sinh! (551) sinh Modified Bessel functions Suppose we change the potential problem so that the top and bottom of the cylinder are grounded but the outer wall at has potential ( ) Then we would need to choose a negative separation constant so that the solutions of the -equation are trigonometric functions: sin + cos At () so we need the sine, and therefore set We also need () so we choose the eigenvalue. This change in sign of the separation constant also affects the equation for the function ()because the sign of the term changes. or, changing variables to : () which is called the modified Bessel equation. Thesolutionstothisequation are () It is usual to define the modified Bessel function () by the relation: () 1 () (1) 8
9 so that the function is always real (whether or not is an integer). Using equation (4) we can write a series expansion for : () 1 ( 1)!Γ ( + +1) ()!Γ ( + +1)³ As with the s, if is an integer, is not independent of : in fact: () () ( 1) () ( 1) () () (3) The second independent solution is usually chosen to be: Then these functions have the limiting forms: and () () +1 (1) () (4) 1 for 1 (5) Γ ( +1)³ () 577 ln for 1 (6) () Γ () for 1 (7) At large À 1, the asymptotic forms are: () 1 (8) and r () (9) (cf Lea Chapter 3 Example 3.9) These functions, like the real exponentials, do not have multiple zeros and are not orthogonal functions. Note that the s are well behaved at the origin but diverge at infinity. For the s the reverse is true. They diverge at the origin but are well behaved at infinity (See Figure below). 9
10 I,K Modified Bessel functions. solid dashed Black Blue 1 Red x The recursion relations satisifed by the modified Bessel functions are similar to, but not identical to, the relations satisfied by the s. For the s, again we can start with the series: Now let 1: !Γ ( + +1)³ 1!Γ ( + +1)³ 1 +1!Γ ( + +)³ 1 ++1!Γ ( + +)³ +1 () (3) and similarly ( ) 1 (31) Expanding out and combining, we get: (3) For the s, the relations are: ( ) 1 ; 1 () +1 ()
11 and consequently: 1.7 Combining functions (33) When solving a physics problem, we start with a partial differential equation and a set of boundary conditions. Separation of variables produces a set of coupled ordinary differential equations in the various coordinates. The standard solution method ("Orthogonal" notes ) requires that we choose the separation constants by fitting the zero boundary conditions first. In a standard 3-dimensional problem, once we have chosen the two separation constants we have no more freedom and the third function is determined. When solving Laplace s equation in cylindrical coordinates, the functions couple as follows: Zero boundary conditions in : ³ The eigenfunctions are of the form: ± ³ sinh + cosh The set of functions ± form a complete orthogonal set on the surfaces constant that bound the region. Zero boundary conditions in : The eigenfunctions are of the form: ³ ³ ³ ³ + sin ± The set of functions sin ± form a complete orthogonal set on the boundary surface constant. Thus in solutions of Laplace s equation, s in always couple with the hyperbolic sines and hyperbolic cosines (or real exponentials) in while the s and s in always couple with the sines and cosines (or complex exponentials) in Example: Suppose the potential on the curved wall of the cylinder is ( ) sin The top and bottom are grounded The solution is of the form + ³ ³ ³ ³ Φ ( ) sin + 1 The solution must be finite on axis at so Now we evaluate the potential at Φ ( ) + sin 1 ³ 11 ³ sin
12 We make use of orthogonality by mutltiplying by sin and integrating from to + 1 Z sin ³ ³ sin + ³ The result is non-zero only for odd Next we make use of the orthogonality of the + Z ³ Z even sin Z sin sin cos sin [( 1) 1] sin To do the integral on the RHS, express the sine in terms of exponentials. Z Z sin ( 1 1) Only one term with survives the integration on the LHS. Thus Φ ( ) 4 1 odd But 1 1 (eqn 3), so ³ ( 1 1) 1 ( 1 1) 4 sin Φ ( ) 4 ³ 1 " 1 odd sin 1 ³ 1 sin # odd The dimensions are correct, and the 1 gives us reasonable convergence. The ratio 1 1 for 1 The potential at is zero, as expected. The 1"source" (the potenial on the surface) generates an 1response. Again we can look at the first few 1
13 terms (here 1). Let s take and 4 (red) and Φ (black) phi/v rho/a 1.8 Continuous set of eigenvalues: the Fourier Bessel Transform In Chapter 7 we approached the Fourier transform by letting the length of the domain in a Fourier series problem become infinite. The orthogonality relation for the exponential functions: Z 1 ³ exp ³ exp becomes Z 1 ( ) That is, the Kronecker delta becomes a delta function, and the countable set of eigenvalues becomes a continuous set of values. The same thing happens with Bessel functions. With a finite domain in say we can determine a countable set of eigenvalues from the set of zeros of the Bessel functions If our domain in becomes infinite, then we cannot determine the eigenvalues, and instead we have a continuous set. The orthogonality relation: Z ³ ³ [ ( )] 13
14 becomes: Z () ( ) ( ) (34) (The proof of this relation is in Lea Appendix 7.) Then the solution to the physics problem is determined as an integral over For example, a solution of Laplace s equation may be written as: Φ ( ) X Z () () () where () depends on the boundary conditions in It will be a combination of the exponentials and + ³ Example: Suppose the potential is () sin on a plane at and we want to find the potential for Then the appropriate function of is, chosen so that Φ as ( a long way from the plane). We also need the which remain finite at Then the solution is of the form: Φ ( ) X Z () () Evaluating Φ on the plane at we get: sin Φ ( ) X Z () () Now we can make use of the orthogonality of the Multiply both sides by and integrate over the range to On the LHS, only the term with survives the integration, and on the RHS only the term with survives. Z sin Z X Z () () sin Z () () a Fourier Bessel transform. Next 1 multiply both sides by ( ) integrate from to in and use equation (34) to get: Z sin Z Z ( ) () () ( ) Z () ( ) ( ) which determines the coefficient ( ) in terms of the known potential on the plane. Only is non-zero, as expected from the azimuthal symmetry. 1 We also drop the primes on the for convenience. 14
15 On the RHS, let We get (dropping the primes on ) Z sin () () This integral is GR 6.671#7. So ( 1 () () So finally we have Φ ( ) Z 1 () q1 () Z 1 1 For we have Z 1 Φ ( ) 1 1 sin which gives us back the potential on the plane. Thechargedensityontheplaneisgivenby ˆ GR 6.554# So we need Φ Z 1 Z ( 1) ( 1) [!] Let sin cos Then the integral is Thus Z sin (+1) cos cos Z [!] ³ ³ Z 1 sin (+1) (+1) 1 ( +1)!! ( +)!! GR 3.61#3, LeaChP9d Φ 4 ( 1) [!] ( 1) (!) 15 ³ ( +1)!! +1 ( +1)! ³ ( +1)!! ( +1)!
16 and so the charge density on the plane is 4 ( 1) (!) Theplotshowstheseriesupto 4 8 ( +1)!! ( +1)! sigma a/epsv rho/a Convince yourself that this agrees with the field line diagram. (The dashed line is the potential.) 16
1 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 information1 Expansion in orthogonal functions
Notes "orthogonal" January 9 Expansion in orthogonal functions To obtain a more useful form of the Green s function, we ll want to expand in orthogonal functions that are (relatively) easy to integrate.
More informationPhysics 6303 Lecture 15 October 10, Reminder of general solution in 3-dimensional cylindrical coordinates. sinh. sin
Physics 6303 Lecture 15 October 10, 2018 LAST TIME: Spherical harmonics and Bessel functions Reminder of general solution in 3-dimensional cylindrical coordinates,, sin sinh cos cosh, sin sin cos cos,
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 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 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 informationSolutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and
More informationPhysics 6303 Lecture 8 September 25, 2017
Physics 6303 Lecture 8 September 25, 2017 LAST TIME: Finished tensors, vectors, and matrices At the beginning of the course, I wrote several partial differential equations (PDEs) that are used in many
More informationPARTIAL DIFFERENTIAL EQUATIONS and BOUNDARY VALUE PROBLEMS
PARTIAL DIFFERENTIAL EQUATIONS and BOUNDARY VALUE PROBLEMS NAKHLE H. ASMAR University of Missouri PRENTICE HALL, Upper Saddle River, New Jersey 07458 Contents Preface vii A Preview of Applications and
More informationSurvival Guide to Bessel Functions
Survival Guide to Bessel Functions December, 13 1 The Problem (Original by Mike Herman; edits and additions by Paul A.) For cylindrical boundary conditions, Laplace s equation is: [ 1 s Φ ] + 1 Φ s s s
More informationFOURIER SERIES, TRANSFORMS, AND BOUNDARY VALUE PROBLEMS
fc FOURIER SERIES, TRANSFORMS, AND BOUNDARY VALUE PROBLEMS Second Edition J. RAY HANNA Professor Emeritus University of Wyoming Laramie, Wyoming JOHN H. ROWLAND Department of Mathematics and Department
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 informationIntroductions to ExpIntegralEi
Introductions to ExpIntegralEi Introduction to the exponential integrals General The exponential-type integrals have a long history. After the early developments of differential calculus, mathematicians
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 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 informationPhysics 6303 Lecture 9 September 17, ct' 2. ct' ct'
Physics 6303 Lecture 9 September 17, 018 LAST TIME: Finished tensors, vectors, 4-vectors, and 4-tensors One last point is worth mentioning although it is not commonly in use. It does, however, build on
More informationBessel function - Wikipedia, the free encyclopedia
Bessel function - Wikipedia, the free encyclopedia Bessel function Page 1 of 9 From Wikipedia, the free encyclopedia In mathematics, Bessel functions, first defined by the mathematician Daniel Bernoulli
More informationADVANCED ENGINEERING MATHEMATICS
ADVANCED ENGINEERING MATHEMATICS DENNIS G. ZILL Loyola Marymount University MICHAEL R. CULLEN Loyola Marymount University PWS-KENT O I^7 3 PUBLISHING COMPANY E 9 U Boston CONTENTS Preface xiii Parti ORDINARY
More informationTwo 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 informationLecture 4b. Bessel functions. Introduction. Generalized factorial function. 4b.1. Using integration by parts it is easy to show that
4b. Lecture 4b Using integration by parts it is easy to show that Bessel functions Introduction In the previous lecture the separation of variables method led to Bessel's equation y' ' y ' 2 y= () 2 Here
More informationMATHEMATICAL FORMULAS AND INTEGRALS
MATHEMATICAL FORMULAS AND INTEGRALS ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom Academic Press San Diego New York Boston London
More informationADVANCED ENGINEERING MATHEMATICS MATLAB
ADVANCED ENGINEERING MATHEMATICS WITH MATLAB THIRD EDITION Dean G. Duffy Contents Dedication Contents Acknowledgments Author Introduction List of Definitions Chapter 1: Complex Variables 1.1 Complex Numbers
More informationModified Bessel functions : Iα, Kα
Modified Bessel functions : Iα, Kα The Bessel functions are valid even for complex arguments x, and an important special case is that of a purely imaginary argument. In this case, the solutions to the
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationIntroduction to Mathematical Physics
Introduction to Mathematical Physics Methods and Concepts Second Edition Chun Wa Wong Department of Physics and Astronomy University of California Los Angeles OXFORD UNIVERSITY PRESS Contents 1 Vectors
More information8.5 Taylor Polynomials and Taylor Series
8.5. TAYLOR POLYNOMIALS AND TAYLOR SERIES 50 8.5 Taylor Polynomials and Taylor Series Motivating Questions In this section, we strive to understand the ideas generated by the following important questions:
More informationSpecial Functions of Mathematical Physics
Arnold F. Nikiforov Vasilii B. Uvarov Special Functions of Mathematical Physics A Unified Introduction with Applications Translated from the Russian by Ralph P. Boas 1988 Birkhäuser Basel Boston Table
More information1 Current Flow Problems
Physics 704 Notes Sp 08 Current Flow Problems The current density satisfies the charge conservation equation (notes eqn 7) thusinasteadystate, is solenoidal: + =0 () =0 () In a conducting medium, we may
More informationMATHEMATICAL FORMULAS AND INTEGRALS
HANDBOOK OF MATHEMATICAL FORMULAS AND INTEGRALS Second Edition ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom ACADEMIC PRESS A Harcourt
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 informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationStrauss PDEs 2e: Section Exercise 4 Page 1 of 6
Strauss PDEs 2e: Section 5.3 - Exercise 4 Page of 6 Exercise 4 Consider the problem u t = ku xx for < x < l, with the boundary conditions u(, t) = U, u x (l, t) =, and the initial condition u(x, ) =, where
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 informationChapter 4 Sequences and Series
Chapter 4 Sequences and Series 4.1 Sequence Review Sequence: a set of elements (numbers or letters or a combination of both). The elements of the set all follow the same rule (logical progression). The
More informationBESSEL FUNCTIONS APPENDIX D
APPENDIX D BESSEL FUNCTIONS D.1 INTRODUCTION Bessel functions are not classified as one of the elementary functions in mathematics; however, Bessel functions appear in the solution of many physical problems
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
More informationBessel functions. The drum problem: Consider the wave equation (with c = 1) in a disc with homogeneous Dirichlet boundary conditions:
Bessel functions The drum problem: Consider the wave equation (with c = 1) in a disc with homogeneous Dirichlet boundary conditions: u = u t, u(r 0,θ,t) = 0, u u(r,θ,0) = f(r,θ), (r,θ,0) = g(r,θ). t (Note
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 informationAPPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems
APPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems Fourth Edition Richard Haberman Department of Mathematics Southern Methodist University PEARSON Prentice Hall PEARSON
More informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
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 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 informationLast Update: March 1 2, 201 0
M ath 2 0 1 E S 1 W inter 2 0 1 0 Last Update: March 1 2, 201 0 S eries S olutions of Differential Equations Disclaimer: This lecture note tries to provide an alternative approach to the material in Sections
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 information1 Fundamentals. 1.1 Overview. 1.2 Units: Physics 704 Spring 2018
Physics 704 Spring 2018 1 Fundamentals 1.1 Overview The objective of this course is: to determine and fields in various physical systems and the forces and/or torques resulting from them. The domain of
More informationThe Sommerfeld Polynomial Method: Harmonic Oscillator Example
Chemistry 460 Fall 2017 Dr. Jean M. Standard October 2, 2017 The Sommerfeld Polynomial Method: Harmonic Oscillator Example Scaling the Harmonic Oscillator Equation Recall the basic definitions of the harmonic
More informationAND NONLINEAR SCIENCE SERIES. Partial Differential. Equations with MATLAB. Matthew P. Coleman. CRC Press J Taylor & Francis Croup
CHAPMAN & HALL/CRC APPLIED MATHEMATICS AND NONLINEAR SCIENCE SERIES An Introduction to Partial Differential Equations with MATLAB Second Edition Matthew P Coleman Fairfield University Connecticut, USA»C)
More informationElectromagnetic waves in free space
Waveguide notes 018 Electromagnetic waves in free space We start with Maxwell s equations for an LIH medum in the case that the source terms are both zero. = =0 =0 = = Take the curl of Faraday s law, then
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 informationFourier Analysis Fourier Series C H A P T E R 1 1
C H A P T E R Fourier Analysis 474 This chapter on Fourier analysis covers three broad areas: Fourier series in Secs...4, more general orthonormal series called Sturm iouville epansions in Secs..5 and.6
More informationLaplace 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 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 informationBoundary. DIFFERENTIAL EQUATIONS with Fourier Series and. Value Problems APPLIED PARTIAL. Fifth Edition. Richard Haberman PEARSON
APPLIED PARTIAL DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems Fifth Edition Richard Haberman Southern Methodist University PEARSON Boston Columbus Indianapolis New York San Francisco
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 informationMath 425 Fall All About Zero
Math 425 Fall 2005 All About Zero These notes supplement the discussion of zeros of analytic functions presented in 2.4 of our text, pp. 127 128. Throughout: Unless stated otherwise, f is a function analytic
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 informationOn Exponential Decay and the Riemann Hypothesis
On Exponential Decay and the Riemann Hypothesis JEFFREY N. COOK ABSTRACT. A Riemann operator is constructed in which sequential elements are removed from a decaying set by means of prime factorization,
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 informationIndex. Excerpt from "Calculus" 2013 AoPS Inc. Copyrighted Material INDEX
Index #, 2 \, 5, 4 [, 4 - definition, 38 ;, 2 indeterminate form, 2 22, 24 indeterminate form, 22 23, see Euler s Constant 2, 2, see infinity, 33 \, 6, 3, 2 3-6-9 triangle, 2 4-dimensional sphere, 85 45-45-9
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 informationMTH5201 Mathematical Methods in Science and Engineering 1 Fall 2014 Syllabus
MTH5201 Mathematical Methods in Science and Engineering 1 Fall 2014 Syllabus Instructor: Dr. Aaron Welters; O ce: Crawford Bldg., Room 319; Phone: (321) 674-7202; Email: awelters@fit.edu O ce hours: Mon.
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 informationSteady and unsteady diffusion
Chapter 5 Steady and unsteady diffusion In this chapter, we solve the diffusion and forced convection equations, in which it is necessary to evaluate the temperature or concentration fields when the velocity
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 informationSUMMATION TECHNIQUES
SUMMATION TECHNIQUES MATH 53, SECTION 55 (VIPUL NAIK) Corresponding material in the book: Scattered around, but the most cutting-edge parts are in Sections 2.8 and 2.9. What students should definitely
More informationwhich implies that we can take solutions which are simultaneous eigen functions of
Module 1 : Quantum Mechanics Chapter 6 : Quantum mechanics in 3-D Quantum mechanics in 3-D For most physical systems, the dynamics is in 3-D. The solutions to the general 3-d problem are quite complicated,
More informationElectrodynamics I Midterm - Part A - Closed Book KSU 2005/10/17 Electro Dynamic
Electrodynamics I Midterm - Part A - Closed Book KSU 5//7 Name Electro Dynamic. () Write Gauss Law in differential form. E( r) =ρ( r)/ɛ, or D = ρ, E= electricfield,ρ=volume charge density, ɛ =permittivity
More informationBOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES
1 BOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES 1.1 Separable Partial Differential Equations 1. Classical PDEs and Boundary-Value Problems 1.3 Heat Equation 1.4 Wave Equation 1.5 Laplace s Equation
More informationz b k P k p k (z), (z a) f (n 1) (a) 2 (n 1)! (z a)n 1 +f n (z)(z a) n, where f n (z) = 1 C
. Representations of Meromorphic Functions There are two natural ways to represent a rational function. One is to express it as a quotient of two polynomials, the other is to use partial fractions. The
More informationAdvanced. Engineering Mathematics
Advanced Engineering Mathematics A new edition of Further Engineering Mathematics K. A. Stroud Formerly Principal Lecturer Department of Mathematics, Coventry University with additions by Dexter j. Booth
More informationSolutions to Laplace s Equations
Solutions to Laplace s Equations Lecture 14: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay We have seen that in the region of space where there are no sources of charge,
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 informationFourier Analysis Partial Differential Equations (PDEs)
PART C Fourier Analysis. Partial Differential Equations (PDEs) CHAPTER CHAPTER Fourier Analysis Partial Differential Equations (PDEs) Chapter and Chapter are directly related to each other in that Fourier
More informationClassroom Tips and Techniques: Eigenvalue Problems for ODEs - Part 1
Classroom Tips and Techniques: Eigenvalue Problems for DEs - Part 1 Initializations restart with LinearAlgebra : with Student Calculus1 : with RootFinding : Introduction Robert J. Lopez Emeritus Professor
More informationDifferential Equations with Boundary Value Problems
Differential Equations with Boundary Value Problems John Polking Rice University Albert Boggess Texas A&M University David Arnold College of the Redwoods Pearson Education, Inc. Upper Saddle River, New
More informationGreen s function for the wave equation
Green s function for the wave equation Non-relativistic case January 2018 1 The wave equations In the Lorentz Gauge, the wave equations for the potentials are (Notes 1 eqns 44 and 43): 1 2 2 2 2 0 (1)
More informationLast Update: April 7, 201 0
M ath E S W inter Last Update: April 7, Introduction to Partial Differential Equations Disclaimer: his lecture note tries to provide an alternative approach to the material in Sections.. 5 in the textbook.
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 4884 NOVEMBER 9, 7 Summary This is an introduction to ordinary differential equations We
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 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 informationBoundary value problems
1 Introduction Boundary value problems Lecture 5 We have found that the electric potential is a solution of the partial differential equation; 2 V = ρ/ǫ 0 The above is Poisson s equation where ρ is the
More informationComputation of the scattering amplitude in the spheroidal coordinates
Computation of the scattering amplitude in the spheroidal coordinates Takuya MINE Kyoto Institute of Technology 12 October 2015 Lab Seminar at Kochi University of Technology Takuya MINE (KIT) Spheroidal
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 informationTAYLOR AND MACLAURIN SERIES
TAYLOR AND MACLAURIN SERIES. Introduction Last time, we were able to represent a certain restricted class of functions as power series. This leads us to the question: can we represent more general functions
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 informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
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 informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationOn Exponential Decay and the Riemann Hypothesis
On Exponential Decay and the Riemann Hypothesis JEFFREY N. COOK ABSTRACT. A Riemann operator is constructed in which sequential elements are removed from a decaying set by means of prime factorization,
More informationAs shown in the text, we can write an arbitrary azimuthally-symmetric solution to Laplace s equation in spherical coordinates as:
Remarks: In dealing with spherical coordinates in general and with Legendre polynomials in particular it is convenient to make the sustitution c = cosθ. For example, this allows use of the following simplification
More informationQuantum Orbits. Quantum Theory for the Computer Age Unit 9. Diving orbit. Caustic. for KE/PE =R=-3/8. for KE/PE =R=-3/8. p"...
W.G. Harter Coulomb Obits 6-1 Quantum Theory for the Computer Age Unit 9 Caustic for KE/PE =R=-3/8 F p' p g r p"... P F' F P Diving orbit T" T T' Contact Pt. for KE/PE =R=-3/8 Quantum Orbits W.G. Harter
More informationSecond-order ordinary differential equations Special functions, Sturm-Liouville theory and transforms R.S. Johnson
Second-order ordinary differential equations Special functions, Sturm-Liouville theory and transforms R.S. Johnson R.S. Johnson Second-order ordinary differential equations Special functions, Sturm-Liouville
More informationSome functions and their derivatives
Chapter Some functions an their erivatives. Derivative of x n for integer n Recall, from eqn (.6), for y = f (x), Also recall that, for integer n, Hence, if y = x n then y x = lim δx 0 (a + b) n = a n
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 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 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 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 informationLinear Partial Differential Equations for Scientists and Engineers
Tyn Myint-U Lokenath Debnath Linear Partial Differential Equations for Scientists and Engineers Fourth Edition Birkhäuser Boston Basel Berlin Tyn Myint-U 5 Sue Terrace Westport, CT 06880 USA Lokenath Debnath
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 informationSOLVED PROBLEMS ON TAYLOR AND MACLAURIN SERIES
SOLVED PROBLEMS ON TAYLOR AND MACLAURIN SERIES TAYLOR AND MACLAURIN SERIES Taylor Series of a function f at x = a is ( f k )( a) ( x a) k k! It is a Power Series centered at a. Maclaurin Series of a function
More information