Classroom Tips and Techniques: Eigenvalue Problems for ODEs - Part 1
|
|
- Hillary Goodman
- 5 years ago
- Views:
Transcription
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 of Mathematics and Maple Fellow Maplesoft Some boundary value problems for partial differential equations are amenable to analytic techniques. For example, the constant-coefficient, second-order linear equations called the heat, wave, and potential equations are solved with some type of Fourier series representation obtained from the Sturm-Liouville eigenvalue problem that arises upon separating variables. In this, and the next two articles in this series, we examine the role of Maple in the solution of such boundary value problems. We will show efficient techniques for separating variables, then show how to guide Maple through the solution of the resulting Sturm-Liouville eigenvalue problems. In four previous columns we have demonstrated how to implement the Fourier series calculations in Maple. Together, these articles represent the essence of the typical undergraduate course in boundary value problems. In fact, the author's Advanced Engineering Mathematics electronic book includes a complete course in the subject area. In subsequent articles, we will explore separation of variables for Laplace's equation in cylindrical coordinates, a process that leads to Bessel's equation as the DE in the Sturm-Liouville eigenvalue problem, which will now be singular and not regular as in this article. Finally, we will separate variables for Laplace's equation in spherical coordinates so that the resulting DE in the Sturm- Liouville eigenvalue problem will be Legendre's equation, again a singular problem. The difference between these two problems is that for Bessel's equation, the resulting Bessel functions themselves become eigenfunctions, whereas for Legendre's equation, it is far more difficult to pass from Legendre functions in a general solution to the Legendre polynomials that are the eigenfunctions.
2 Separating Variables in a PDE It is not difficult to implement separation of variables in Maple. For example, the one-dimensional wave equation v 2 u x, t = c2 v2 2 vt vx 2 u x, t v 2 vt 2 u x, t = v 2 c2 vx 2 u x, t separates under the assumption that u x, t can be written as the product (3.1) U := X x T t Substitution leads to U := X x T t (3.2) (3.1) u x, t = U X x d 2 dt 2 T t = d 2 c2 whereupon division by u x, t in its separated form yields dx 2 X x T t (3.3) simplify (3.3) U c 2 d 2 dt 2 T t T t c 2 = d 2 dx 2 X x X x Note that we have also divided through by c 2, the wave speed. The Bernoulli separation constant λ is now introduced, resulting in the two ordinary differential equations (3.4) q 1 d lhs (3.4) =λ q 2 d rhs (3.4) =λ more commonly written as q 1 := q 2 := d 2 dt 2 T t T t c 2 =λ d 2 dx 2 X x X x =λ (3.5) DE 1 d numer normal lhs q 2 Krhs q 2 = 0 DE 2 d numer normal lhs q 1 Krhs q 1 = 0
3 DE 1 := d2 dx 2 X x Kλ X x = 0 DE 2 := d2 dt 2 T t Kλ T t c2 = 0 The first equation governing X x is a formally self-adjoint DE that will become part of the Sturm- Liouville eigenvalue problems considered in the following section. Solving the Sturm-Liouville Eigenvalue Problem Homogeneous Dirichlet Conditions A Sturm-Liouville eigenvalue problem arises if we couple the differential equation (3.6) DE 1 d 2 dx 2 X x Kλ X x = 0 (4.1.1) with the homogeneous Dirichlet boundary conditions X 0 = X π = 0. This is an example of "separated" conditions for which we are guaranteed a complete set of orthogonal eigenfunctions. ne approach to solving this eigenvalue problem might be to invoke Maple's dsolve command, with the result dsolve DE 1, X 0 = 0, X π = 0, X x X x = 0 Unfortunately, Maple returns the zero-solution, and does not even begin to solve the eigenvalue problem. That is why we are writing this series of articles. Maple can actually be of great help in solving such problems, but it must be guided by the user. The characteristic equation for (4.1.1) is obtained by making the exponential guess EG := e m x EG := e m x Substitution and the obvious algebra then give 1 EG collect DE 1 X x = EG, exp m 2 Kλ = 0 It should be obvious that the characteristic roots are m =G λ, and that λ = 0 is the "repeated root" case. Since the form of the general solution of (4.1.1) changes for this case, we consider it first as a special case. The differential equation is then DE 1a d DE 1 λ = 0
4 DE 1a := d2 dx 2 X x = 0 and a solution of this equation with just the left-hand boundary condition is dsolve DE 1a, X 0 = 0, X x X x = _C1 x Applying the right-hand boundary condition to this solution leads to the equation (4.1.2) rhs (4.1.2) = 0 x =π _C1 π = 0 from which it is clear that the multiplicative constant has to be zero. Since this leads to the trivial solution X x = 0, we conclude that λ = 0 is not an eigenvalue. We next consider the case λs 0, and again use the strategy of applying just the left-hand boundary condition. The resulting solution is X 1 d rhs dsolve DE 1, X 0 = 0, X x X 1 :=K_C2 e λ x C_C2 e K λ x Because Maple could return this solution with _C1 as the arbitrary constant, we use the following device to rewrite the solution in terms of the constant c. No matter what name Maple returns for the arbitrary constant, our construction will always be in terms of c. paramdremove has, indets X 1, λ, x 1 : X 2 d X 1 param = c X 2 :=Kc e λ x Cc e K λ x Applying the right-hand boundary condition leads to the algebraic equation X 2 x =π=0 Kc e λ π Cc e K λ π = 0 Without intervention, Maple will provide just the zero solution for λ, as we confirm with (4.1.3) solve (4.1.3),λ 0 The requisite tool is the environment variable _EnvAllSolutions := true _EnvAllSolutions := true which, if set to the value true, enables Maple to return the solution
5 solve (4.1.3),λ K_Z1 2 Maple uses the notation _Zk, k = 1, 2,, to indicate an integer. Unfortunately, each time this solve command is executed, the index k will increment by 1, so again, we cannot reference this symbol in a fail-safe way. Thus, we resort to the following device temp d indets (4.1.4) 1 : (4.1.4) temp = n Kn 2 to write the eigenvalues with a unique notation. The eigenvalues are the negative numbers Kn 2, n = 1, 2,. It is also possible to lead Maple through the steps resulting in this solution for the eigenvalues. We begin by writing (4.1.3) in the form (4.1.4) factor (4.1.3) Kc e λ π Ke K λ π = 0 Since c = 0 would result in the trivial solution X x = 0, we must have cs 0, so that map coeff, (4.1.5), c Ke λ π Ce K λ π = 0 This difference of exponentials suggests rewriting the left-hand side as (4.1.5) (4.1.6) convert (4.1.6), trig K2 sinh λ π = 0 (4.1.7) from which we realize that if λ is real, the only solution of this equation will be λ = 0, again leading to the trivial solution. This inspires the substitution simplify evalc (4.1.7) assuming µ0 λ =Kµ 2 K2 I sin µ π = 0 from which it easily follows that µ is an integer, given in Maple by (4.1.8) solve (4.1.8),µ _Z2~ Homogeneous Neumann Conditions The Sturm-Liouville eigenvalue problem that arises when the homogeneous (and separated) Neumann boundary conditions X# 0 = X# π = 0 are applied to the differential equation
6 DE 1 d 2 dx 2 X x Kλ X x = 0 will have a different set of eigenvalues and eigenfunctions than we obtained with the Dirichlet conditions. Since we already know that λ = 0 must be handled as a special case, we begin by obtaining the general solution of the DE DE 1a d 2 dx 2 X x = 0 We again take precautions to write the solution with unique arbitrary constants, obtaining temp := rhs dsolve DE 1a, X x : X 0 d eval temp, _C1 = b, _C2 = a X 0 := b xca To this general solution, we apply the two boundary conditions, resulting in the algebraic equations X 0 ' = 0 x = 0 b = 0 and X 0 ' = 0 x =π b = 0 Consequently, λ = 0 is an eigenvalue, and the corresponding eigenfunction is X 0 x = a Now, consider the case λs 0. We again find the general solution, and write it with uniquely named arbitrary constants. temp := rhs dsolve DE 1, X x : X 1 d eval temp, _C1 = c 1, _C2 = c 2 X 1 := c 1 e λ x Cc 2 e K λ x To this general solution we apply the boundary conditions, obtaining the homogeneous algebraic equations Eq 1 d X 1 ' = 0 x = 0 Eq 1 := c 1 λ Kc 2 λ = 0 and Eq 2 d X 1 ' = 0 x =π Eq 2 := c 1 λ e λ π Kc 2 λ e K λ π = 0 Maple's solve command cannot determine values of λ for which this system has nontrivial solutions. We guide Maple by writing the system matrix as
7 GenerateMatrix Eq 1, Eq 2, c 1, c 2 1 λ K λ λ e λ π K λ e K λ π (4.2.1) and setting its determinant to zero, the condition under which the system will have nontrivial solutions for c 1 and c 2. Determinant (4.2.1) = 0 Kλ e K λ π Cλ e λ π = 0 Clearly, λ = 0 is a solution, but that belongs to the first case. "Canceling" λ and solving, gives (4.2.2) solve (4.2.2) λ,λ where again, we replace this notation by K_Z3 2 (4.2.3) temp d indets (4.2.3) 1 : (4.2.3) temp = n Kn 2 so that the eigenvalues are more easily seen to be λ =Kn 2, n = 1, 2,. Having determined the eigenvalues, we return to the algebraic equations determining c 1 and c 2, and make the obvious substitution for λ to obtain the equations Eq 3 d simplify Eq 1 λ =Kn 2 Eq 4 d simplify Eq 2 λ =Kn 2 assuming n T posint; assuming n T posint Eq 3 := I c 1 nki c 2 n = 0 The solution of these equations is now Eq 4 := I n K1 n c 1 Kc 2 = 0 solve Eq 3, Eq 4, c 1, c 2 c 2 = c 1, c 1 = c 1 so that the remaining eigenfunctions are temp := subs c 1 = C, c 2 = C,λ =Kn 2, X 1 : simplify evalc temp assuming n T posint 2 C cos n x
8 f course, the constant C can be taken as 1 π, or as 2 2 to obtain orthonormal eigenfunctions. Example: An DE that is Not Self-Adjoint We conclude this article with an example of an eigenvalue problem in which the differential equation is not self-adjoint. Such an example is not covered by the Sturm-Liouville theory, but it does illustrate how Maple can be guided through some of the more difficult calculations that arise in eigenvalue problems. In particular, note the more general separated (but still homogeneous) boundary conditions in the following problem. deq := d2 dx 2 y x C2 d dx y x Cλ y x = 0 bc 1 d 3 D y 0 C5 y 0 = 0 bc 2 d 2 D y π Cy π = 0 deq := d2 dx 2 y x C2 d dx y x Cλ y x = 0 bc 1 := 3 D y 0 C5 y 0 = 0 bc 2 := 2 D y π Cy π = 0 To obtain solutions of the differential equation, we seek a fundamental set of linearly independent members. Making the exponential guess leads to the characteristic equation 1 EG collect deq, exp y x = EG whose solutions are m 2 C2 mcλ = 0 (5.1) solve (5.1), m K1C 1Kλ, K1K 1Kλ From these characteristic roots, we can see that λ = 1 is the repeated-root case for which the solution of the DE assumes the special form (5.2) temp := dsolve deq, y x : λ = 1 Y 1 d unapply eval rhs temp, _C1 = r 1, _C2 = s 1, x : 'Y 1 ' x = Y 1 x Y 1 x = r 1 e Kx Cs 1 e Kx x For λs 1, the general solution of the DE can be written as temp := dsolve deq, y x :
9 Y 2 d unapply eval rhs temp, _C1 = r 2, _C2 = s 2, x : 'Y 2 ' x = Y 2 x Y 2 x = r 2 e K1C KλC1 x Cs 2 e K1K KλC1 x Even though Maple defaults to an exponential form for the members of its fundamental set, we should not assume that λ is characterized by anything other than the condition λs 1. In fact, at this point we don't even know if the eigenvalues are real. Consequently, we avoid separating the calculation into the cases λ! 1 and λ 1 because that would imply that we knew λ was not a complex number. We let the boundary conditions determine λ by considering just the cases λ = 1 and λs 1. In the first case, we get the algebraic equations for k to 2 do q k d bc k y = Y1 end do q 1 := 2 r 1 C3 s 1 = 0 q 2 :=Kr 1 e Kπ Ks 1 e Kπ πc2 s 1 e Kπ = 0 and in the second, for k to 2 do q 2Ck := bc k y = Y 2 end do q 3 := 3 r 2 K1C KλC1 C3 s 2 K1K KλC1 C5 r 2 C5 s 2 = 0 q 4 := 2 r 2 K1C KλC1 e K1C KλC1 π C2 s 2 K1K KλC1 e Cr 2 e K1C KλC1 π Cs 2 e K1K KλC1 π = 0 The determinants of the system matrix in each case are then K1K KλC1 π Q 1 := Determinant GenerateMatrix q 1, q 2, r 1, s 1 1 Q 1 :=K2 e Kπ πc7 e Kπ and (5.3) Q 2 := collect Determinant GenerateMatrix q 3, q 4, r 2, s 2 1, exp Q 2 := 8K7 KλC1 K6 λ e K1C KλC1 π C K8K7 KλC1 (5.4) K1C KλC1 π C6 λ e respectively. Since the determinant in (5.3) is nonzero, we must have r 1 = s 1 = 0, so the corresponding
10 solution is X x h 0, and λ = 1 is therefore not an eigenvalue. There are zeros of the determinant in (5.4), some of which are Analytic Q 2, re =K10..10, im =K , , , , The Analytic command in the RootFinding package finds all real and complex roots in a region of the complex plane. We already know that λ = 1 is not an eigenvalue, so we have the real eigenvalues remove w/is w K1! 0.001, (5.5) , , , While the Analytic command provides a guarantee that all the solutions in a fixed region have been found, it is slow. If we accept (5.5) as evidence that the eigenvalues are positive and real, then we can find them slightly faster by changing the form of the determinant in (5.4) to (5.5) (5.6) simplify evalc Q 2 assuming λ 1 K2 I e Kπ 7 λk1 cos λk1 π K8 sin λk1 π C6 sin λk1 π λ and using the Roots command to find (5.7) (5.7) L := Roots,λ = , numeric = true Kπ I e L := , , , , , For any such eigenvalue, the equations q 3 and q 4 are redundant, so we can solve one of them for s 2, say, in terms of r 2, obtaining (5.8) S := solve q 3, s 2 S := r 2 2C3 K2C3 KλC1 KλC1 For any eigenvalue λ, the corresponding eigenfunction is then given (up to a multiplicative constant) by Yd 1 collect Y r 2 x, 2 s 2 = S r 2, exp Y := e K1C KλC1 x 2C3 C KλC1 K1K KλC1 x e K2C3 KλC1 With λ 1, the exponential terms will have trigonometric equivalents, so we write simplify evalc Y assuming λ 1 K 1 K5C9 λ 2 ekx 9 cos λk1 x K9 cos λk1 x λc6 λk1 sin λk1 x (5.9)
11 and then K4 I sin λk1 x C6 I λk1 cos λk1 x collect (5.9), sin, cos K 2 ekx 6 λk1 K4 I sin λk1 x K5C9 λ (5.10) K 2 ekx 9K9 λc6 I λk1 cos λk1 x K5C9 λ The real and imaginary parts are respectively Ya := evalc R (5.10) assuming λ 1; Yb := evalc I (5.10) assuming λ 1 Ya :=K 12 ekx λk1 sin λk1 x K5C9 λ K 2 ekx 9K9 λ cos λk1 x K5C9 λ Yb := 8 ekx sin λk1 x K5C9 λ K 12 ekx λk1 cos λk1 x K5C9 λ In fact, each of the real and imaginary parts is separately a solution. But they cannot be independent solutions - they must be multiples of each other. This common multiple is f := simplify convert taylor Ya, x, 3, polynom Yb f :=K 3 2 λk1 and we can verify that indeed, the real and imaginary parts of Y x are linearly dependent by showing the linear combination f Y b KY a vanishes. simplify Yb fkya The common multiple can be found with 0 g := Ya x = 0 and a simplified form of the real part is then g :=K 2 9K9 λ K5C9 λ Yc := collect Ya g, exp, cos, sin Yc := cos λk1 x C Using Y c, we obtain the first six eigenfunctions as 6 λk1 sin λk1 x 9K9 λ e Kx
12 unassign 'φ' φ 1 := simplify expand convert Yc for k from 2 to nops L φ k d Yc do λ = L 1, exp assuming x 0; end do λ = L k φ 1 := e K x C e K x φ 2 := cos x K sin x e Kx φ 3 := cos x K sin x e Kx φ 4 := cos x K sin x e Kx φ 5 := cos x K sin x e Kx φ 6 := cos x K sin x e Kx A graph of the first six eigenfunctions is given in Figure 1, where our normalization has all eigenfunctions passing through 0, 1. The eigenfunctions φ k x, k = 1,, 6, are given in dotted black, then solid black, blue, red, green, and magenta. plot seq φ k, k = 1..6, x = 0..π, color = black, black, blue, red, green, magenta, linestyle = 2, 1$5
13 x K0.5 Figure 1 Eigenfunctions φ k, k = 2., 6 are in solid black, blue, red, green, and magenta, respectively, whereas φ 0 is in dotted black Since the differential equation is not self-adjoint, we are not assured of orthogonal eigenfunctions. The inner products φ i,φ j, 1 % i, j% 6, are computed and displayed in Table 1. If these eigenfunctions were orthogonal, the off-diagonal entries would be zero. That they are not is our evidence the eigenfunctions are not orthogonal. Matrix 6, 6, i, j / φ i φ j dx 0 π
14 Table 1 The inner products φ i,φ j, 1 % i, j% 6 We expended some energy separating real and imaginary parts in the solution Y x, more in showing that these parts were multiples of each other. A more direct way to find the simple trigonometric form we used to generate the eigenfunctions is to obtain the general solution of the differential equation under the assumption that λ 1. This leads to temp := dsolve deq, y x assuming λ 1 : Y 3 d unapply eval rhs temp, _C1 = r 3, _C2 = s 3, x : 'Y 3 ' x = Y 3 x Y 3 x = r 3 e Kx sin λk1 x Cs 3 e Kx cos λk1 x and, after applying the boundary conditions, to for k to 2 do q 4Ck := bc k y = Y 3 end do q 5 := 3 r 3 λk1 C2 s 3 = 0 q 6 :=Kr 3 e Kπ sin λk1 π C2 r 3 e Kπ cos λk1 π λk1 Ks 3 e Kπ cos λk1 π K2 s 3 e Kπ sin λk1 π λk1 = 0 If we now express r 3, say, in terms of s 3, as in isolate q 5, r 3 we can write this new general solution as r 3 =K 2 3 s 3 λk1 (5.11) Ydd 1 s 3 collect eval Y 3 x, (5.11), exp, s 3 Yd := K 2 3 sin λk1 x λk1 Ccos λk1 x e Kx
15 To show this is equivalent to our Y c x, we show the difference Y d KY c vanishes. simplify Yd KYc 0 Legal Notice: The copyright for this application is owned by the author(s). Neither Maplesoft nor the author are responsible for any errors contained within and are not liable for any damages resulting from the use of this material. This application is intended for non-commercial, non-profit use only. Contact the author for permission if you wish to use this application in for-profit activities.
Classroom Tips and Techniques: Series Expansions Robert J. Lopez Emeritus Professor of Mathematics and Maple Fellow Maplesoft
Introduction Classroom Tips and Techniques: Series Expansions Robert J. Lopez Emeritus Professor of Mathematics and Maple Fellow Maplesoft Maple has the ability to provide various series expansions and
More informationWithout a Vector Calculus Coordinate System
Classroom Tips and Techniques: Point-and-Click Access to the Differential Operators of Vector Calculus Robert J. Lopez Emeritus Professor of Mathematics and Maple Fellow Maplesoft Introduction The -operator
More informationClassroom Tips and Techniques: Electric Field from Distributed Charge
Classroom Tips and Techniques: Electric Field from Distributed Charge Introduction Robert J. Lopez Emeritus Professor of Mathematics and Maple Fellow Maplesoft This past summer I was asked if Maple could
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 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 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 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 information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More 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 informationOrdinary Differential Equations II
Ordinary Differential Equations II February 23 2017 Separation of variables Wave eq. (PDE) 2 u t (t, x) = 2 u 2 c2 (t, x), x2 c > 0 constant. Describes small vibrations in a homogeneous string. u(t, x)
More informationLEGENDRE 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 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 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 information1 Introduction. Green s function notes 2018
Green s function notes 8 Introduction Back in the "formal" notes, we derived the potential in terms of the Green s function. Dirichlet problem: Equation (7) in "formal" notes is Φ () Z ( ) ( ) 3 Z Φ (
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
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 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 information4.3 - Linear Combinations and Independence of Vectors
- Linear Combinations and Independence of Vectors De nitions, Theorems, and Examples De nition 1 A vector v in a vector space V is called a linear combination of the vectors u 1, u,,u k in V if v can be
More information1 Separation of Variables
Jim ambers ENERGY 281 Spring Quarter 27-8 ecture 2 Notes 1 Separation of Variables In the previous lecture, we learned how to derive a PDE that describes fluid flow. Now, we will learn a number of analytical
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 informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt =
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 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 informationThis is a closed book exam. No notes or calculators are permitted. We will drop your lowest scoring question for you.
Math 54 Fall 2017 Practice Final Exam Exam date: 12/14/17 Time Limit: 170 Minutes Name: Student ID: GSI or Section: This exam contains 9 pages (including this cover page) and 10 problems. Problems are
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 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 informationElements of linear algebra
Elements of linear algebra Elements of linear algebra A vector space S is a set (numbers, vectors, functions) which has addition and scalar multiplication defined, so that the linear combination c 1 v
More informationMaple in Differential Equations
Maple in Differential Equations and Boundary Value Problems by H. Pleym Maple Worksheets Supplementing Edwards and Penney Differential Equations and Boundary Value Problems - Computing and Modeling Preface
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 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 informationMathematics of Imaging: Lecture 3
Mathematics of Imaging: Lecture 3 Linear Operators in Infinite Dimensions Consider the linear operator on the space of continuous functions defined on IR. J (f)(x) = x 0 f(s) ds Every function in the range
More informationwhich arises when we compute the orthogonal projection of a vector y in a subspace with an orthogonal basis. Hence assume that P y = A ij = x j, x i
MODULE 6 Topics: Gram-Schmidt orthogonalization process We begin by observing that if the vectors {x j } N are mutually orthogonal in an inner product space V then they are necessarily linearly independent.
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 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 informationIn this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation.
Eigen Function Expansion and Applications. In this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation. a/ The theory. b/ Example: Solving the Euler equation in two ways.
More informationConsequences of Orthogonality
Consequences of Orthogonality Philippe B. Laval KSU Today Philippe B. Laval (KSU) Consequences of Orthogonality Today 1 / 23 Introduction The three kind of examples we did above involved Dirichlet, Neumann
More informationSAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear Algebra
1.1. Introduction SAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear algebra is a specific branch of mathematics dealing with the study of vectors, vector spaces with functions that
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
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 information1 Review of simple harmonic oscillator
MATHEMATICS 7302 (Analytical Dynamics YEAR 2017 2018, TERM 2 HANDOUT #8: COUPLED OSCILLATIONS AND NORMAL MODES 1 Review of simple harmonic oscillator In MATH 1301/1302 you studied the simple harmonic oscillator:
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 informationHW #3 Solutions: M552 Spring (c)
HW #3 Solutions: M55 Spring 6. (4.-Trefethen & Bau: parts (a), (c) and (e)) Determine the SVDs of the following matrices (by hand calculation): (a) [ 3 (c) (e) [ ANS: In each case we seek A UΣV. The general
More informationMTH Linear Algebra. Study Guide. Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education
MTH 3 Linear Algebra Study Guide Dr. Tony Yee Department of Mathematics and Information Technology The Hong Kong Institute of Education June 3, ii Contents Table of Contents iii Matrix Algebra. Real Life
More informationSpherical Harmonics on S 2
7 August 00 1 Spherical Harmonics on 1 The Laplace-Beltrami Operator In what follows, we describe points on using the parametrization x = cos ϕ sin θ, y = sin ϕ sin θ, z = cos θ, where θ is the colatitude
More informationCentral Difference Approximation of the First Derivative
Differentiation of Continuous Functions Central Difference Approximation of the First Derivative Introduction Ana Catalina Torres, Autar Kaw University of South Florida United States of America kaw@eng.usf.edu
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 informationarxiv: v2 [math.gm] 8 Apr 2007
published elec. at www.maplesoft.com arxiv:0704.07v [math.gm] 8 Apr 007 1 Introduction A PROCEDURE FOR FINDING THE k TH POWER OF A MATRIX Branko Malešević Faculty of Electrical Engineering, University
More informationSecond Order Linear Equations
Second Order Linear Equations Linear Equations The most general linear ordinary differential equation of order two has the form, a t y t b t y t c t y t f t. 1 We call this a linear equation because the
More informationRelevant sections from AMATH 351 Course Notes (Wainwright): 1.3 Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls): 1.1.
Lecture 8 Qualitative Behaviour of Solutions to ODEs Relevant sections from AMATH 351 Course Notes (Wainwright): 1.3 Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls): 1.1.1 The last few
More information18.06 Problem Set 8 Solution Due Wednesday, 22 April 2009 at 4 pm in Total: 160 points.
86 Problem Set 8 Solution Due Wednesday, April 9 at 4 pm in -6 Total: 6 points Problem : If A is real-symmetric, it has real eigenvalues What can you say about the eigenvalues if A is real and anti-symmetric
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 information2.4 Eigenvalue problems
2.4 Eigenvalue problems Associated with the boundary problem (2.1) (Poisson eq.), we call λ an eigenvalue if Lu = λu (2.36) for a nonzero function u C 2 0 ((0, 1)). Recall Lu = u. Then u is called an eigenfunction.
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 information1.6 Computing and Existence
1.6 Computing and Existence 57 1.6 Computing and Existence The initial value problem (1) y = f(x,y), y(x 0 ) = y 0 is studied here from a computational viewpoint. Answered are some basic questions about
More informationEXAM MATHEMATICAL METHODS OF PHYSICS. TRACK ANALYSIS (Chapters I-V). Thursday, June 7th,
EXAM MATHEMATICAL METHODS OF PHYSICS TRACK ANALYSIS (Chapters I-V) Thursday, June 7th, 1-13 Students who are entitled to a lighter version of the exam may skip problems 1, 8-11 and 16 Consider the differential
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 informationPhysics 6303 Lecture 11 September 24, LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation
Physics 6303 Lecture September 24, 208 LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation, l l l l l l. Consider problems that are no axisymmetric; i.e., the potential depends
More informationMATH-UA 263 Partial Differential Equations Recitation Summary
MATH-UA 263 Partial Differential Equations Recitation Summary Yuanxun (Bill) Bao Office Hour: Wednesday 2-4pm, WWH 1003 Email: yxb201@nyu.edu 1 February 2, 2018 Topics: verifying solution to a PDE, dispersion
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 informationMATHEMATICS (MATH) Calendar
MATHEMATICS (MATH) This is a list of the Mathematics (MATH) courses available at KPU. For information about transfer of credit amongst institutions in B.C. and to see how individual courses transfer, go
More informationReview and problem list for Applied Math I
Review and problem list for Applied Math I (This is a first version of a serious review sheet; it may contain errors and it certainly omits a number of topic which were covered in the course. Let me know
More informationNOTES ON LINEAR ALGEBRA CLASS HANDOUT
NOTES ON LINEAR ALGEBRA CLASS HANDOUT ANTHONY S. MAIDA CONTENTS 1. Introduction 2 2. Basis Vectors 2 3. Linear Transformations 2 3.1. Example: Rotation Transformation 3 4. Matrix Multiplication and Function
More informationv(x, 0) = g(x) where g(x) = f(x) U(x). The solution is where b n = 2 g(x) sin(nπx) dx. (c) As t, we have v(x, t) 0 and u(x, t) U(x).
Problem set 4: Solutions Math 27B, Winter216 1. The following nonhomogeneous IBVP describes heat flow in a rod whose ends are held at temperatures u, u 1 : u t = u xx < x < 1, t > u(, t) = u, u(1, t) =
More information5.1 Second order linear differential equations, and vector space theory connections.
Math 2250-004 Wed Mar 1 5.1 Second order linear differential equations, and vector space theory connections. Definition: A vector space is a collection of objects together with and "addition" operation
More informationLinear Algebra using Dirac Notation: Pt. 2
Linear Algebra using Dirac Notation: Pt. 2 PHYS 476Q - Southern Illinois University February 6, 2018 PHYS 476Q - Southern Illinois University Linear Algebra using Dirac Notation: Pt. 2 February 6, 2018
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 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 informationSolutions to the Sample Problems for the Final Exam UCLA Math 135, Winter 2015, Professor David Levermore
Solutions to the Sample Problems for the Final Exam UCLA Math 135, Winter 15, Professor David Levermore Every sample problem for the Midterm exam and every problem on the Midterm exam should be considered
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 ecture Notes 8 - PDEs Page 8.0 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More information(L, t) = 0, t > 0. (iii)
. Sturm Liouville Boundary Value Problems 69 where E is Young s modulus and ρ is the mass per unit volume. If the end x = isfixed, then the boundary condition there is u(, t) =, t >. (ii) Suppose that
More informationChapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form
Chapter 5. Linear Algebra A linear (algebraic) equation in n unknowns, x 1, x 2,..., x n, is an equation of the form a 1 x 1 + a 2 x 2 + + a n x n = b where a 1, a 2,..., a n and b are real numbers. 1
More informationSecond-Order Linear ODEs (Textbook, Chap 2)
Second-Order Linear ODEs (Textbook, Chap ) Motivation Recall from notes, pp. 58-59, the second example of a DE that we introduced there. d φ 1 1 φ = φ 0 dx λ λ Q w ' (a1) This equation represents conservation
More informationSection 8.2 : Homogeneous Linear Systems
Section 8.2 : Homogeneous Linear Systems Review: Eigenvalues and Eigenvectors Let A be an n n matrix with constant real components a ij. An eigenvector of A is a nonzero n 1 column vector v such that Av
More informationThe Gram-Schmidt Process 1
The Gram-Schmidt Process In this section all vector spaces will be subspaces of some R m. Definition.. Let S = {v...v n } R m. The set S is said to be orthogonal if v v j = whenever i j. If in addition
More informationUNDETERMINED COEFFICIENTS SUPERPOSITION APPROACH *
4.4 UNDETERMINED COEFFICIENTS SUPERPOSITION APPROACH 19 Discussion Problems 59. Two roots of a cubic auxiliary equation with real coeffi cients are m 1 1 and m i. What is the corresponding homogeneous
More information2.2 Separable Equations
82 2.2 Separable Equations An equation y = f(x, y) is called separable provided algebraic operations, usually multiplication, division and factorization, allow it to be written in a separable form y =
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More informationLet's begin by evaluating several limits. sqrt x 2 C4 K2. limit. x 2, x = 0 ; O limit. 1K cos t t. , t = 0 ; 1C 1 n. limit. 6 x 3 C5.
Dr. Straight's Maple Examples Example II: Calculus perations and Functions Functions Covered: diff, dsolve, implicitdiff, implicitplot, int, limit, odeplot, piecewise, plot, solve, subs, sum, taylor Let's
More informationEIGENVALUES AND EIGENVECTORS
EIGENVALUES AND EIGENVECTORS Diagonalizable linear transformations and matrices Recall, a matrix, D, is diagonal if it is square and the only non-zero entries are on the diagonal This is equivalent to
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 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 information,, rectilinear,, spherical,, cylindrical. (6.1)
Lecture 6 Review of Vectors Physics in more than one dimension (See Chapter 3 in Boas, but we try to take a more general approach and in a slightly different order) Recall that in the previous two lectures
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt =
More information6 Non-homogeneous Heat Problems
6 Non-homogeneous Heat Problems Up to this point all the problems we have considered for the heat or wave equation we what we call homogeneous problems. This means that for an interval < x < l the problems
More informationPartial Differential Equations with MATLAB
CHAPMAN & HALL/CRC APPLIED MATHEMATICS AND NONLINEAR SCIENCE SERIES An Introduction to Partial Differential Equations with MATLAB Second Edition Matthew P. Coleman CHAPMAN & HALL/CRC APPLIED MATHEMATICS
More informationCourse Outline. Date Lecture Topic Reading
Course Outline Date Lecture Topic Reading Graduate Mathematical Physics Tue 24 Aug Linear Algebra: Theory 744 756 Vectors, bases and components Linear maps and dual vectors Inner products and adjoint operators
More informationMath 224, Fall 2007 Exam 3 Thursday, December 6, 2007
Math 224, Fall 2007 Exam 3 Thursday, December 6, 2007 You have 1 hour and 20 minutes. No notes, books, or other references. You are permitted to use Maple during this exam, but you must start with a blank
More informationCHAPTER 3. Matrix Eigenvalue Problems
A SERIES OF CLASS NOTES FOR 2005-2006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 3 A COLLECTION OF HANDOUTS ON SYSTEMS OF ORDINARY DIFFERENTIAL
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 informationRemark 1 By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.
Sec 5 Eigenvectors and Eigenvalues In this chapter, vector means column vector Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called
More informationSAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to Linear Algebra
SAMPLE OF THE STUDY MATERIAL PART OF CHAPTER 1 Introduction to 1.1. Introduction Linear algebra is a specific branch of mathematics dealing with the study of vectors, vector spaces with functions that
More informationLinear Operators, Eigenvalues, and Green s Operator
Chapter 10 Linear Operators, Eigenvalues, and Green s Operator We begin with a reminder of facts which should be known from previous courses. 10.1 Inner Product Space A vector space is a collection of
More informationLINEAR SYSTEMS AND MATRICES
CHAPTER 3 LINEAR SYSTEMS AND MATRICES SECTION 3. INTRODUCTION TO LINEAR SYSTEMS This initial section takes account of the fact that some students remember only hazily the method of elimination for and
More informationCh 5.4: Regular Singular Points
Ch 5.4: Regular Singular Points! Recall that the point x 0 is an ordinary point of the equation if p(x) = Q(x)/P(x) and q(x)= R(x)/P(x) are analytic at at x 0. Otherwise x 0 is a singular point.! Thus,
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 informationMathematical Modeling using Partial Differential Equations (PDE s)
Mathematical Modeling using Partial Differential Equations (PDE s) 145. Physical Models: heat conduction, vibration. 146. Mathematical Models: why build them. The solution to the mathematical model will
More informationReduction to the associated homogeneous system via a particular solution
June PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 5) Linear Algebra This study guide describes briefly the course materials to be covered in MA 5. In order to be qualified for the credit, one
More informationMatrices and Matrix Algebra.
Matrices and Matrix Algebra 3.1. Operations on Matrices Matrix Notation and Terminology Matrix: a rectangular array of numbers, called entries. A matrix with m rows and n columns m n A n n matrix : a square
More informationLecture 10: Fourier Sine Series
Introductory lecture notes on Partial Differential Equations - c Anthony Peirce. Not to be copied, used, or revised without eplicit written permission from the copyright owner. ecture : Fourier Sine Series
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More information