Boundary Value Problems (Sect. 10.1). Two-point Boundary Value Problem.

Similar documents
7.1. Simple Eigenfunction Problems. The main idea of this chapter is to solve the heat equation. main ideas to solve that equation.

Solving the Heat Equation (Sect. 10.5).

Find the Fourier series of the odd-periodic extension of the function f (x) = 1 for x ( 1, 0). Solution: The Fourier series is.

Lecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems

Essential Ordinary Differential Equations

Lecture IX. Definition 1 A non-singular Sturm 1 -Liouville 2 problem consists of a second order linear differential equation of the form.

Physics 250 Green s functions for ordinary differential equations

Boundary-value Problems in Rectangular Coordinates

Section 12.6: Non-homogeneous Problems

Review: Power series define functions. Functions define power series. Taylor series of a function. Taylor polynomials of a function.

1 Wave Equation on Finite Interval

x 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

Ch 10.1: Two Point Boundary Value Problems

Chapter 4: Higher-Order Differential Equations Part 1

Math Assignment 14

Review for Exam 2. Review for Exam 2.

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

MB4018 Differential equations

dx n a 1(x) dy

The One-Dimensional Heat Equation

CHAPTER 10 NOTES DAVID SEAL

CHALLENGE! (0) = 5. Construct a polynomial with the following behavior at x = 0:

CHAPTER 1. Theory of Second Order. Linear ODE s

California State University Northridge MATH 280: Applied Differential Equations Midterm Exam 1

PDE and Boundary-Value Problems Winter Term 2014/2015

Calculus I Exam 1 Review Fall 2016

Method of Separation of Variables

Assignment # 7, Math 370, Fall 2018 SOLUTIONS: y = e x. y = e 3x + 4xe 3x. y = e x cosx.

Overview of Fourier Series (Sect. 6.2). Origins of the Fourier Series.

Linear Algebra Practice Problems

SECOND-ORDER LINEAR DIFFERENTIAL EQUATIONS

MA Chapter 10 practice

ODEs Cathal Ormond 1

Math 5440 Problem Set 7 Solutions

7.1. Simple Eigenfunction Problems. the. The main idea of this chapter is to solve the heat equation. main ideas to solve that equation.

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

1 Solution to Homework 4

Ordinary Differential Equations II

Partial Differential Equations

Consequences of Orthogonality

Test #2 Math 2250 Summer 2003

Topics from Algebra and Pre-Calculus. (Key contains solved problems)

Linear Variable coefficient equations (Sect. 2.1) Review: Linear constant coefficient equations

Polytechnic Institute of NYU MA 2132 Final Practice Answers Fall 2012

Math 4B Notes. Written by Victoria Kala SH 6432u Office Hours: T 12:45 1:45pm Last updated 7/24/2016

The Laplace Transform and the IVP (Sect. 6.2).

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

Ordinary Differential Equation Theory

Partial Differential Equations for Engineering Math 312, Fall 2012

Next, we ll use all of the tools we ve covered in our study of trigonometry to solve some equations.

ENGI 9420 Lecture Notes 8 - PDEs Page 8.01

Boundary-Value Problems

The lecture of 1/23/2013: WHY STURM-LIOUVILLE?

Partial Differential Equations

Linear Variable coefficient equations (Sect. 1.2) Review: Linear constant coefficient equations

Section 8.0 Introduction to Boundary Value Problems. How do initial value problems (IVPs) and boundary value problems (BVPs) differ?

Ordinary Differential Equations

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

Degree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m.

0 3 x < x < 5. By continuing in this fashion, and drawing a graph, it can be seen that T = 2.

Power Series Solutions to the Legendre Equation

Section x7 +

LECTURE 19: SEPARATION OF VARIABLES, HEAT CONDUCTION IN A ROD

systems of linear di erential If the homogeneous linear di erential system is diagonalizable,

Math 23 Practice Quiz 2018 Spring

June 2011 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations

Theorem The BVP for the one-space dimensional heat equation, infinitely. where k > 0, L > 0 are constants, has infinitely. #, c n R.

6 Non-homogeneous Heat Problems

The Laplace Transform (Sect. 4.1). The Laplace Transform (Sect. 4.1).

Lecture Notes for Math 251: ODE and PDE. Lecture 12: 3.3 Complex Roots of the Characteristic Equation

Linear transformations

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

Section 9.8 Higher Order Linear Equations

Unit - 7 Vibration of Continuous System

Exam Basics. midterm. 1 There will be 9 questions. 2 The first 3 are on pre-midterm material. 3 The next 1 is a mix of old and new material.

Advanced Mathematics Unit 2 Limits and Continuity

Advanced Mathematics Unit 2 Limits and Continuity

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

(L, t) = 0, t > 0. (iii)

Vector Spaces ปร ภ ม เวกเตอร

Linear Algebra- Final Exam Review

Math 480 The Vector Space of Differentiable Functions

Math Assignment 6

MT804 Analysis Homework II

Assignment # 8, Math 370, Fall 2018 SOLUTIONS:

swapneel/207

SERIES SOLUTION OF DIFFERENTIAL EQUATIONS

Department of Mathematics IIT Guwahati

This is a closed book exam. No notes or calculators are permitted. We will drop your lowest scoring question for you.

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.

2.12: Derivatives of Exp/Log (cont d) and 2.15: Antiderivatives and Initial Value Problems

A matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and

x 2 = 1 Clearly, this equation is not true for all real values of x. Nevertheless, we can solve it by taking careful steps:

Chapter P: Preliminaries

Review session Midterm 1

A matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and

Ma 221. The material below was covered during the lecture given last Wed. (1/29/14). Homogeneous Linear Equations with Constant Coefficients

Math 333 Qualitative Results: Forced Harmonic Oscillators

Chapter 2. Limits and Continuity 2.6 Limits Involving Infinity; Asymptotes of Graphs

A = 3 1. We conclude that the algebraic multiplicity of the eigenvalues are both one, that is,

Transcription:

Boundary Value Problems (Sect. 10.1). Two-point BVP. from physics. Comparison: IVP vs BVP. Existence, uniqueness of solutions to BVP. Particular case of BVP: Eigenvalue-eigenfunction problem. Two-point Boundary Value Problem. Definition A two-point BVP is the following: Given functions p, q, g, and constants x 1 < x 2, y 1, y 2, b 1, b 2, b1, b 2, find a function y solution of the differential equation y + p(x) y + q(x) y = g(x), together with the extra, boundary conditions, Remarks: b 1 y(x 1 ) + b 2 y (x 1 ) = y 1, b 1 y(x 2 ) + b 2 y (x 2 ) = y 2. Both y and y might appear in the boundary condition, evaluated at the same point. In this notes we only study the case of constant coefficients, y + a 1 y + a 0 y = g(x).

Two-point Boundary Value Problem. s of BVP. Assume x 1 x 2. (1) Find y solution of y + a 1 y + a 0 y = g(x), y(x 1 ) = y 1, y(x 2 ) = y 2. (2) Find y solution of y + a 1 y + a 0 y = g(x), y (x 1 ) = y 1, y (x 2 ) = y 2. (3) Find y solution of y + a 1 y + a 0 y = g(x), y(x 1 ) = y 1, y (x 2 ) = y 2. Boundary Value Problems (Sect. 10.1). Two-point BVP. from physics. Comparison: IVP vs BVP. Existence, uniqueness of solutions to BVP. Particular case of BVP: Eigenvalue-eigenfunction problem.

from physics. Problem: The equilibrium (time independent) temperature of a bar of length L with insulated horizontal sides and the bar vertical extremes kept at fixed temperatures T 0, T L is the solution of the BVP: T (x) = 0, x (0, L), T (0) = T 0, T (L) = T L, y insulation T 0 T L 0 x L x z insulation Boundary Value Problems (Sect. 10.1). Two-point BVP. from physics. Comparison: IVP vs BVP. Existence, uniqueness of solutions to BVP. Particular case of BVP: Eigenvalue-eigenfunction problem.

Comparison: IVP vs BVP. Review: IVP: Find the function values y(t) solutions of the differential equation y + a 1 y + a 0 y = g(t), together with the initial conditions y(t 0 ) = y 1, y (t 0 ) = y 2. Remark: In physics: y(t): Position at time t. Initial conditions: Position and velocity at the initial time t 0. Comparison: IVP vs BVP. Review: BVP: Find the function values y(x) solutions of the differential equation y + a 1 y + a 0 y = g(x), together with the initial conditions y(x 1 ) = y 1, y(x 2 ) = y 2. Remark: In physics: y(x): A physical quantity (temperature) at a position x. Boundary conditions: Conditions at the boundary of the object under study, where x 1 x 2.

Boundary Value Problems (Sect. 10.1). Two-point BVP. from physics. Comparison: IVP vs BVP. Existence, uniqueness of solutions to BVP. Particular case of BVP: Eigenvalue-eigenfunction problem. Existence, uniqueness of solutions to BVP. Review: The initial value problem. Theorem (IVP) Consider the homogeneous initial value problem: y + a 1 y + a 0 y = 0, y(t 0 ) = y 0, y (t 0 ) = y 1, and let r ± be the roots of the characteristic polynomial p(r) = r 2 + a 1 r + a 0. If r + r, real or complex, then for every choice of y 0, y 1, there exists a unique solution y to the initial value problem above. Summary: The IVP above always has a unique solution, no matter what y 0 and y 1 we choose.

Existence, uniqueness of solutions to BVP. Theorem (BVP) Consider the homogeneous boundary value problem: y + a 1 y + a 0 y = 0, y(0) = y 0, y(l) = y 1, and let r ± be the roots of the characteristic polynomial p(r) = r 2 + a 1 r + a 0. (A) If r + r, real, then for every choice of L 0 and y 0, y 1, there exists a unique solution y to the BVP above. (B) If r ± = α ± iβ, with β 0, and α, β R, then the solutions to the BVP above belong to one of these possibilities: (1) There exists a unique solution. (2) There exists no solution. (3) There exist infinitely many solutions. Existence, uniqueness of solutions to BVP. Proof of IVP: We study the case r + r -. The general solution is y(t) = c 1 e r - t + e r + t, c 1, R. The initial conditions determine c 1 and as follows: y 0 = y(t 0 ) = c 1 e r - t 0 + e r + t 0 y 1 = y (t 0 ) = c 1 r - e r - t 0 + r + e r + t 0 Using matrix notation, e r - t 0 e r + t 0 c1 r - e r -t 0 r + e r + t 0 = The linear system above has a unique solution c 1 and for every constants y 0 and y 1 iff the det(z) 0, where e r - t 0 e Z = r + t 0 c1 y0 r - e r - t 0 r + e r + t 0 Z =. y0 y 1. y 1

Existence, uniqueness of solutions to BVP. Proof of IVP: e r - t 0 e Recall: Z = r + t 0 r - e r - t 0 r + e r + t 0 Z c1 = y0 y 1. A simple calculation shows det(z) = ( r + r - ) e (r + +r - ) t0 0 r + r -. Since r + r -, the matrix Z is invertible and so c1 = Z 1 y0. y 1 We conclude that for every choice of y 0 and y 1, there exist a unique value of c 1 and, so the IVP above has a unique solution. Existence, uniqueness of solutions to BVP. Proof of BVP: The general solution is y(x) = c 1 e r - x + e r + x, c 1, R. The boundary conditions determine c 1 and as follows: y 0 = y(0) = c 1 +. y 1 = y(l) = c 1 e r - L + e r + L Using matrix notation, 1 1 c1 e r - L e r + L = The linear system above has a unique solution c 1 and for every constants y 0 and y 1 iff the det(z) 0, where 1 1 c1 y0 Z = Z =. e r - L e r + L y0 y 1. y 1

Existence, uniqueness of solutions to BVP. 1 1 Proof of IVP: Recall: Z = A simple calculation shows e r - L e r + L Z c1 det(z) = e r + L e r - L 0 e r + L e r - L. (A) If r + r - and real-valued, then det(z) 0. = We conclude: For every choice of y 0 and y 1, there exist a unique value of c 1 and, so the BVP in (A) above has a unique solution. (B) If r ± = α ± iβ, with α, β R and β 0, then det(z) = e αl( e iβl e iβl) y0 y 1. det(z) = 2i e αl sin(βl). Since det(z) = 0 iff βl = nπ, with n integer, (1) If βl nπ, then BVP has a unique solution. (2) If βl = nπ then BVP either has no solutions or it has infinitely many solutions. Existence, uniqueness of solutions to BVP. Find y solution of the BVP y + y = 0, y(0) = 1, y(π) = 1. Solution: The characteristic polynomial is The general solution is The boundary conditions are p(r) = r 2 + 1 r ± = ±i. y(x) = c 1 cos(x) + sin(x). 1 = y(0) = c 1, 1 = y(π) = c 1 c 1 = 1, free. We conclude: y(x) = cos(x) + sin(x), with R. The BVP has infinitely many solutions.

Existence, uniqueness of solutions to BVP. Find y solution of the BVP y + y = 0, y(0) = 1, y(π) = 0. Solution: The characteristic polynomial is The general solution is The boundary conditions are p(r) = r 2 + 1 r ± = ±i. y(x) = c 1 cos(x) + sin(x). 1 = y(0) = c 1, 0 = y(π) = c 1 The BVP has no solution. Existence, uniqueness of solutions to BVP. Find y solution of the BVP y + y = 0, y(0) = 1, y(π/2) = 1. Solution: The characteristic polynomial is The general solution is The boundary conditions are p(r) = r 2 + 1 r ± = ±i. y(x) = c 1 cos(x) + sin(x). 1 = y(0) = c 1, 1 = y(π/2) = c 1 = = 1. We conclude: y(x) = cos(x) + sin(x). The BVP has a unique solution.

Boundary Value Problems (Sect. 10.1). Two-point BVP. from physics. Comparison: IVP vs BVP. Existence, uniqueness of solutions to BVP. Particular case of BVP: Eigenvalue-eigenfunction problem. Particular case of BVP: Eigenvalue-eigenfunction problem. Problem: Find a number λ and a non-zero function y solutions to the boundary value problem Remark: This problem is similar to the eigenvalue-eigenvector problem in Linear Algebra: Given an n n matrix A, find λ and a non-zero n-vector v solutions of Differences: A Av λ v = 0. { } computing a second derivative and v {a function y}. applying the boundary conditions.

Particular case of BVP: Eigenvalue-eigenfunction problem. Find every λ R and non-zero functions y solutions of the BVP Remarks: We will show that: (1) If λ 0, then the BVP has no solution. (2) If λ > 0, then there exist infinitely many eigenvalues λ n and eigenfunctions y n, with n any positive integer, given by ( nπ ) 2, ( nπx ) λ n = yn (x) = sin, L L (3) Analogous results can be proven for the same equation but with different types of boundary conditions. For example, for y(0) = 0, y (L) = 0; or for y (0) = 0, y (L) = 0. Particular case of BVP: Eigenvalue-eigenfunction problem. Find every λ R and non-zero functions y solutions of the BVP Solution: Case λ = 0. The equation is y = 0 y(x) = c 1 + x. The boundary conditions imply 0 = y(0) = c 1, 0 = c 1 + L c 1 = = 0. Since y = 0, there are NO non-zero solutions for λ = 0.

Particular case of BVP: Eigenvalue-eigenfunction problem. Find every λ R and non-zero functions y solutions of the BVP Solution: Case λ < 0. Introduce the notation λ = µ 2. The characteristic equation is The general solution is The boundary condition are p(r) = r 2 µ 2 = 0 r ± = ±µ. y(x) = c 1 e µx + e µx. 0 = y(0) = c 1 +, 0 = y(l) = c 1 e µl + e µl. Particular case of BVP: Eigenvalue-eigenfunction problem. Find every λ R and non-zero functions y solutions of the BVP Solution: Recall: y(x) = c 1 e µx + e µx and c 1 + = 0, c 1 e µl + e µl = 0. We need to solve the linear system 1 1 c1 0 = Z 0 e µl e µl c1 = 0 1 1, Z = 0 e µl e µl Since det(z) = e µl e µl 0 for L 0, matrix Z is invertible, so the linear system above has a unique solution c 1 = 0 and = 0. Since y = 0, there are NO non-zero solutions for λ < 0.

Particular case of BVP: Eigenvalue-eigenfunction problem. Find every λ R and non-zero functions y solutions of the BVP Solution: Case λ > 0. Introduce the notation λ = µ 2. The characteristic equation is The general solution is The boundary condition are p(r) = r 2 + µ 2 = 0 r ± = ±µi. y(x) = c 1 cos(µx) + sin(µx). 0 = y(0) = c 1, y(x) = sin(µx). 0 = y(l) = sin(µl), 0 sin(µl) = 0. Particular case of BVP: Eigenvalue-eigenfunction problem. Find every λ R and non-zero functions y solutions of the BVP Solution: Recall: c 1 = 0, 0, and sin(µl) = 0. The non-zero solution condition is the reason for 0. Hence sin(µl) = 0 µ n L = nπ µ n = nπ L. Recalling that λ n = µ 2 n, and choosing = 1, we conclude ( nπ λ n = L ) 2, ( nπx ) yn (x) = sin. L