Assignment Number 3 Solutions

Similar documents
Chapter 4. Fourier Series

PROBLEMS AND SOLUTIONS 2

Solutions to Final Exam Review Problems

SOLUTION SET VI FOR FALL [(n + 2)(n + 1)a n+2 a n 1 ]x n = 0,

Løsningsførslag i 4M

Math 220B Final Exam Solutions March 18, 2002

Assignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.

Ma 530 Introduction to Power Series

Question 1: The magnetic case

Math 475, Problem Set #12: Answers

Ma 4121: Introduction to Lebesgue Integration Solutions to Homework Assignment 5

SOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.

j=1 dz Res(f, z j ) = 1 d k 1 dz k 1 (z c)k f(z) Res(f, c) = lim z c (k 1)! Res g, c = f(c) g (c)

Sequences and Limits

Sequences and Series of Functions

MATH2007* Partial Answers to Review Exercises Fall 2004

MIDTERM 3 CALCULUS 2. Monday, December 3, :15 PM to 6:45 PM. Name PRACTICE EXAM SOLUTIONS

UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014

Seunghee Ye Ma 8: Week 5 Oct 28

MATH 6101 Fall 2008 Newton and Differential Equations

Math 113, Calculus II Winter 2007 Final Exam Solutions

( ) (( ) ) ANSWERS TO EXERCISES IN APPENDIX B. Section B.1 VECTORS AND SETS. Exercise B.1-1: Convex sets. are convex, , hence. and. (a) Let.

Fourier Series and the Wave Equation

6.3 Testing Series With Positive Terms

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

Fourier Series and their Applications

MATH 304: MIDTERM EXAM SOLUTIONS

Strauss PDEs 2e: Section Exercise 4 Page 1 of 5. u tt = c 2 u xx ru t for 0 < x < l u = 0 at both ends u(x, 0) = φ(x) u t (x, 0) = ψ(x),

PROBLEM SET 5 SOLUTIONS 126 = , 37 = , 15 = , 7 = 7 1.

TMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods

Problem. Consider the sequence a j for j N defined by the recurrence a j+1 = 2a j + j for j > 0

MATH 10550, EXAM 3 SOLUTIONS

Math 110 Assignment #6 Due: Monday, February 10

PLEASE MARK YOUR ANSWERS WITH AN X, not a circle! 1. (a) (b) (c) (d) (e) 3. (a) (b) (c) (d) (e) 5. (a) (b) (c) (d) (e) 7. (a) (b) (c) (d) (e)

Math 132, Fall 2009 Exam 2: Solutions

Inverse Nodal Problems for Differential Equation on the Half-line

Section 11.8: Power Series

Math 12 Final Exam, May 11, 2011 ANSWER KEY. 2sinh(2x) = lim. 1 x. lim e. x ln. = e. (x+1)(1) x(1) (x+1) 2. (2secθ) 5 2sec2 θ dθ.

MATH 31B: MIDTERM 2 REVIEW

Math 113 Exam 3 Practice

Solutions for Math 411 Assignment #8 1

Quiz. Use either the RATIO or ROOT TEST to determine whether the series is convergent or not.

Chapter 6 Infinite Series

Enumerative & Asymptotic Combinatorics

REVIEW 1, MATH n=1 is convergent. (b) Determine whether a n is convergent.

Zeros of Polynomials

Math 128A: Homework 1 Solutions

PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES

MATH Exam 1 Solutions February 24, 2016

Solutions to home assignments (sketches)

NICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =

MH1101 AY1617 Sem 2. Question 1. NOT TESTED THIS TIME

LECTURE 21. DISCUSSION OF MIDTERM EXAM. θ [0, 2π). f(θ) = π θ 2

Orthogonal Functions

MAT1026 Calculus II Basic Convergence Tests for Series

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

Mathematics 170B Selected HW Solutions.

1+x 1 + α+x. x = 2(α x2 ) 1+x

1 Heat Equation Dirichlet Boundary Conditions

: Transforms and Partial Differential Equations

Where do eigenvalues/eigenvectors/eigenfunctions come from, and why are they important anyway?

Chapter 10 Partial Differential Equations and Fourier Series

6. Uniform distribution mod 1

MTH 142 Exam 3 Spr 2011 Practice Problem Solutions 1

Appendix: The Laplace Transform

Sturm-Liouville Eigenvalue Problem

Math 155 (Lecture 3)

CHAPTER 10 INFINITE SEQUENCES AND SERIES

Recurrence Relations

Problem 4: Evaluate ( k ) by negating (actually un-negating) its upper index. Binomial coefficient

The Riemann Zeta Function

Review for Test 3 Math 1552, Integral Calculus Sections 8.8,

SCORE. Exam 2. MA 114 Exam 2 Fall 2017

x !1! + 1!2!

MAT 271 Project: Partial Fractions for certain rational functions

b i u x i U a i j u x i u x j

Quantum Theory Assignment 3

SECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES

Matsubara-Green s Functions

1 Lecture 2: Sequence, Series and power series (8/14/2012)

PHY4905: Nearly-Free Electron Model (NFE)

COMM 602: Digital Signal Processing

Math 113 Exam 3 Practice

Section 1.4. Power Series

Complex Analysis Spring 2001 Homework I Solution

6.003 Homework #3 Solutions

M17 MAT25-21 HOMEWORK 5 SOLUTIONS

Mathematics 116 HWK 21 Solutions 8.2 p580

Homework 3 Solutions

Solutions to quizzes Math Spring 2007

INFINITE SERIES PROBLEMS-SOLUTIONS. 3 n and 1. converges by the Comparison Test. and. ( 8 ) 2 n. 4 n + 2. n n = 4 lim 1

HOMEWORK #10 SOLUTIONS

Castiel, Supernatural, Season 6, Episode 18

Iterative Techniques for Solving Ax b -(3.8). Assume that the system has a unique solution. Let x be the solution. Then x A 1 b.

Solutions to Homework 1

Indian Statistical Institute, Bangalore Centre Solution set of M.Math II Year, End-Sem Examination 2015 Fourier Analysis

COMPUTING FOURIER SERIES

1 6 = 1 6 = + Factorials and Euler s Gamma function

Calculus 2 - D. Yuen Final Exam Review (Version 11/22/2017. Please report any possible typos.)

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

Transcription:

Math 4354, Assigmet Number 3 Solutios 1. u t (x, t) = u xx (x, t), < x <, t > (1) u(, t) =, u(, t) = u(x, ) = x ( 1) +1 u(x, t) = e t si(x). () =1 Solutio: Look for simple solutios i the form u(x, t) = X(x)T (t). We obtai X λx =, The solutio of the secod equatio is T λt =. T (t) = Ce kλt (3) where C is a arbitrary costat. Furthermore, the boudary coditios give X()T (t) =, X()T (t) = for all t. Sice T (t) is ot idetically zero we obtai the desired eigevalue problem X (x) λx(x) =, X() =, X() =. (4) (a) If λ = the X(x) = ax + b so applyig the boudary coditios we get Zero is ot a eigevalue. (b) If λ = µ > the = X() = b, = X() = a a = b =. X(x) = a cosh(µx) + b sih(µx). Applyig the boudary coditios we have = X() = a a = = X() = b sih(µ) b =. Therefore, there are o positive eigevalues. Cosider the followig alterative argumet: If X (x) = λx(x) the multiplyig by X we have X(x)X (x) = λx(x). Itegrate this expressio from x = to x =. We have λ X(x) dx = X(x)X (x) dx = X (x) dx + X(x)X (x). 1

Sice X() = X() = we coclude λ = X (x) dx X(x) dx ad we see that λ must be less tha or equal to zero. (c) So, fially, cosider λ = µ so that X(x) = a cos(µx) + b si(µx). Applyig the boudary coditios we have = X() = a a = = X() = b si(µ). From this we coclude si(µ) = which implies ad therefore µ = λ = µ =, X (x) = si(x), = 1,,.. (5) From (14) we also have the associated fuctios T (t) = e t. The oly problem remaiig is to somehow pick the costats a so that the iitial coditio u(x, ) = ϕ(x) is satisfied. To do this we cosider what we leared from Fourier series. I particular we look for u as a ifiite sum u(x, t) = b e t si (x) ad we try to fid {b } satisfyig =1 ϕ(x) = u(x, ) = b si (x). But this othig more tha a Sie expasio of the fuctio ϕ o the iterval (, ). b = = = = [ x = ( 1)+1 ϕ(x) si (x) dx =1 x si (x) dx/ ( ) cos (x) x dx ( ) cos (x) u(x, t) = ( ) cos (x) x dx] ( 1) +1 e t si(x). (6) =1

. u t (x, t) = u xx (x, t), < x <, t > (7) u x (, t) =, u x (, t) = u(x, ) = x Solutio: Look for simple solutios i the form We obtai u(x, t) = X(x)T (t). X λx =, The solutio of the secod equatio is T λt =. T (t) = Ce λt (8) where C is a arbitrary costat. Furthermore, the boudary coditios give X ()T (t) =, X ()T (t) = for all t. Sice T (t) is ot idetically zero we obtai the desired eigevalue problem There are, i geeral, three cases: X (x) λx(x) =, X () =, X () =. (9) (a) If λ = the X(x) = ax + b so applyig the boudary coditios we get = X () = a, = X () = a a =. Notice that b is still a arbitrary costat. We coclude that λ = is a eigevalue with eigefuctio ϕ (x) = 1. (b) If λ = µ > the ad X(x) = a cosh(µx) + b sih(µx) X (x) = aµ sih(µx) + bµ cosh(µx). Applyig the boudary coditios we have = X () = bµ b = = X () = aµ sih(µ) a =. Therefore, there are o positive eigevalues. Cosider the followig alterative argumet: If X (x) = λx(x) the multiplyig by X we have X(x)X (x) = λx(x). Itegrate this expressio from x = to x =. We have λ X(x) dx = X(x)X (x) dx = X (x) dx + X(x)X (x). Sice X () = X () = we coclude λ = X (x) dx X(x) dx ad we see that λ must be less tha or equal to zero ( zero oly if X = ). 3

(c) So, fially, cosider λ = µ so that X(x) = a cos(µx) + b si(µx) ad X (x) = aµ si(µx) + bµ cos(µx). Applyig the boudary coditios we have = X () = bµ b = = X () = aµ si(µ). From this we coclude si(µ) = which implies µ = ad therefore λ = µ = (), X (x) = cos(x), = 1,,.. (1) From (8) we also have the associated fuctios T (t) = e t. Our work o Fourier series showed us that a = 1 ϕ(x) dx, a = ϕ(x) cos (x) dx. (11) As a explicit example for the iitial coditio cosider = 1 ad ϕ(x) = x(1 x). I this case (11) becomes We have ad a = 1 a = ϕ(x) dx, a = a = 1 ϕ(x) dx = 1 = 1 [ ] x =. ϕ(x) cos (x) dx = ϕ(x) cos (x) dx. x dx x cos (x) dx = = [ si (x) x = ( ) si (x) x dx ( ) cos (x) = (( 1) 1) = si (x) ] dx 4, odd., eve 4

I order to elimiate the odd terms i the expasio we itroduce a ew idex, k by = k where k = 1,,. So fially we arrive at the solutio u(x, t) = 4 k=1 1 (k 1) e (k 1)t cos((k 1)x). (1) 3. u t (x, t) = u xx (x, t), < x <, t > (13) u(, t) =, u x (, t) = u(x, ) = x Solutio: Look for simple solutios i the form u(x, t) = X(x)T (t). We obtai X λx =, T λt =. The solutio of the secod equatio is T (t) = Ce λt (14) where C is a arbitrary costat. Furthermore, the boudary coditios give X()T (t) =, X ()T (t) = for all t. Sice T (t) is ot idetically zero we obtai the desired eigevalue problem X (x) λx(x) =, X() =, X () =. (15) There are, i geeral, three cases: (a) If λ = the X(x) = ax + b so applyig the boudary coditios we get = X() = b, = X () = a a = b =. Zero is ot a eigevalue. (b) If λ = µ > the X(x) = a cosh(µx) + b sih(µx) ad X (x) = aµ sih(µx) + bµ cosh(µx). Applyig the boudary coditios we have = X () = aµ a = = X () = bµ cosh(µ) b =. Therefore, there are o positive eigevalues. 5

Cosider the followig alterative argumet: If X (x) = λx(x) the multiplyig by X we have X(x)X (x) = λx(x). Itegrate this expressio from x = to x =. We have λ X(x) dx = X(x)X (x) dx = X (x) dx + X(x)X (x). Sice X() = X () = we coclude λ = X (x) dx X(x) dx ad we see that λ must be less tha or equal to zero. (c) So, fially, cosider λ = µ so that ad X(x) = a cos(µx) + b si(µx) X (x) = aµ si(µx) + bµ cos(µx). Applyig the boudary coditios we have = X() = aµ a = = X () = bµ cos(µ). From this we coclude cos(µ) = which implies ad therefore µ = ( 1) ( ) ( 1) λ = µ =, X (x) = si(µ x), = 1,,.. (16) From (14) we also have the associated fuctios T (t) = e λt. We look for u as a ifiite sum u(x, t) = ad we try to fid {b } satisfyig ( ) ( 1)x b e λt si =1 ϕ(x) = u(x, ) = ( ) ( 1)x b si. =1 But this othig more tha a Sie type expasio of the fuctio ϕ o the iterval (, ). Usig ( ) ( 1)x X (x) = si 6

we have ϕ(x) = b k X k (x). k=1 We proceed as usual by multiplyig both sides by X (x) ad itegratig from to ad usig the orthogoality (described below i (18), () ). X (x)ϕ(x) dx = b k X (x)x k (x) dx k=1 which implies orthogoality: To see this recall that X j λ λ k ad therefore λ X (x)x k (x) dx = Therefore For = k we have b = ϕ(x)x (x) dx. (17), = k X (x)x k (x) dx =., k = λ j X j ad X j () =, X j() =. First cosider k so = = X (x)x k (x) dx (18) X (x)x k(x) dx + [X (x)x k (x)] X (x)x k (x) dx + [X (x)x k (x) X (x)x k(x)] = λ k X (x)x k (x) dx. (19) (λ λ k ) X (x)x k (x) dx = X (x) dx = X (x)x k (x) dx =. ( ) ( 1)x si dx () = 1 (1 cos (( 1)x)) dx = 1 ( ) 1 si (( 1)x) 1 =. (1) 7

( ) 1 For the iitial coditio ϕ(x) = x. Recall that µ = b = = = [ ϕ(x)x (x) dx = ( x cos (µ ) x) dx µ x cos(µ x) µ + = [ cos(µ ) + si(µ ] x) µ = µ [ cos(( 1)/) + µ = 8( 1)+1 ( 1) x si (µ x) dx ] cos(µ x) dx µ si(( 1)/) (( 1)/) ] u(x, t) = 8 =1 ( ) ( 1) +1 ( 1) ( 1) e [( 1)/]t si x. () 8