AMS 212A Applied Mathematical Methods I Appendices of Lecture 06 Copyright by Hongyun Wang, UCSC. ( ) cos2

Size: px
Start display at page:

Download "AMS 212A Applied Mathematical Methods I Appendices of Lecture 06 Copyright by Hongyun Wang, UCSC. ( ) cos2"

Transcription

1 AMS 22A Applied Mathematical Methods I Appendices of Lecture 06 Copyright by Hongyun Wang UCSC Appendix A: Proof of Lemma Lemma : Let (x ) be the solution of x ( r( x)+ q( x) )sin 2 + ( a) 0 < cos2 where p(x) > 0 and r(x) > 0 for x in [a b] (from the Sturm-Liouville problem). Then (x ) has the 3 properties below. For x > a a) (x ) is a strictly increasing function of. b) lim ( x )+ + c) lim ( x ) 0 Proof of a): Consider 2 >. We need to show (x 2 ) (x ) > 0 for x > a. We prove it in a few steps. Step A: Suppose (x 0 2 ) (x 0 ) 0 and sin( 0 ) 0 at some x 0. Then (x 2 ) (x ) is a strictly increasing function of x in a neighborhood of x 0. Proof of Step A: ( x ) ( r( x)+ q( x) )sin > 0 x cos2 > x ( x 0 2 ) x ( x 0 ) ( 2 )r( x 0 )sin 2 Taylor expansion of (x 2 ) (x ) at x 0 : ( x 0 + s 2 ) ( x 0 + s ) ( 2 )r( x 0 )sin 2 0 s + O s 2 Step B: The conclusion of Step A is still true if sin( 0 ) 0. Suppose (x 0 2 ) (x 0 ) 0 and sin( 0 ) 0 at some x 0. Then (x 2 ) (x ) is a strictly increasing function of x in a neighborhood of x

2 AMS 22A Applied Mathematical Methods I Proof of Step B: x ( x 0 2 ) x ( x 0 ) 0 Let us calculate xx (x ). xx ( x ) ( r( x)+ q( x) )sin ( r( x)+ q( x) )2sin cos x > xx ( x 0 2 ) xx ( x 0 ) 0 cos 2 2cos sin x ext we calculate xxx (x ). We will only look at its value at x 0. At x 0 we have (x 0 2 ) (x 0 ) 0 where sin( 0 ) 0. It follows that x xx0 0 independent of sin 0 independent of xx 0 These results significantly simplify the calculation of xxx (x 0 ). xxx xx0 cos 2 + 2( r( x)+ q( x) )cos 2 x xx p( x) cos2 ( x ) 2 xx 0 2 xx 0 > xxx ( x 0 2 ) xxx ( x 0 ) 2( 2 )r( x 0 ) Taylor expansion of (x 2 ) (x ) at x 0 : 0 2 > 0 Step C: ( x 0 + s 2 ) ( x 0 + s ) ( 3 2 )r( x 0 ) 0 2 s 3 + O( s 4 ) (x 2 ) and (x ) are the same at x a : (a). Results of Steps A and B tell us that (x 2 ) (x ) > 0 for x slightly above a (it does not matter whether or not sin 0)

3 AMS 22A Applied Mathematical Methods I Suppose x 0 is the location where (x 2 ) (x ) transitions from positive to non-positive. Mathematically that is ( x 0 2 ) ( x 0 ) 0 ( x 2 ) ( x ) 0 for x slightly above x 0. But this contradicts with the results of Steps A and B. Therefore we conclude that (x 2 ) (x ) > 0 for x > a. Proof of b): We first calculate an integral. sin 2 + cos 2 d d tan 2 2 tan dx x ds 2 s as + Since r(x) and p(x) are positive on [a b] there exist r 0 q 0 p 0 and such that r( x) r 0 > 0 q( x) q 0 for x [ a b] p( x) p 0 > 0 for x [ a b] Let ( c) q 0 / r 0. For sufficiently large we have r( x)+ q( x) ( ( c) )r( x)+ ( c) r( x)+ q x ( ( c) for x [ a b] The differential equation for becomes x ( r( x)+ q( x) )sin 2 + ( c) ( sin 2 + cos 2 cos2 for x [ a b] Thus for sufficiently large (x ) is an increasing function of x. Let x n be the location where (x ) hits n. We re-write the differential equation for as d ( ( c) sin 2 + cos 2 dx - 3 -

4 AMS 22A Applied Mathematical Methods I Integrate from x n to x n+ we have ( n+) d ( ( c) ( x sin 2 + cos 2 n+ x n ) Using the analytical expression of the integral we obtain d ( x n+ x n ) ( ( c) sin 2 + cos as + ( ( c) > ( x a) x n x n n + ( x 0 a) 0 as + ( ( c) (Recall that we showed that (x ) is an increasing function of x when is large.) It follows that for any x > a and > 0 we can make (x ) > when is large enough Therefore we conclude lim ( x )+ + ote: The proof shows that when is large (b ) hits at O( 2 ) which means O( 2 ) This is consistent with the solution of u u. u( 0) 0 u( L) 0 Proof of c): We first show that (x ) 0 for x > a. At 0 the evolution equation of tells us r( x)+ q x x 0 ( ) sin 2 + p( x) cos2 > (x ) will never go below 0 for x > a. 0 p( x) > 0 Since r(x) and p(x) are positive on [a b] there exist r 0 q 0 p 0 and such that - 4 -

5 AMS 22A Applied Mathematical Methods I r( x) r 0 > 0 q( x) q 0 for x [ a b] p( x) p 0 > 0 for x [ a b] Let ( c) q 0 / r 0. For negative and large we have r( x)+ q( x) ( + ( c) )r( x) ( c) r( x)+ q x ( + ( c) for x [ a b] > x ( r( x)+ q( x) )sin 2 + cos2 ( + ( c) sin 2 + for x [ a b] p 0 For any small > 0 we select M r 0 sin 2 + p 0. When ( + ( c) )< M the two statements below are true. As long as (x ) satisfies sin 2 sin 2 we have x That is (x ) decreases with large negative slope -/ until sin 2 sin 2. Once (x ) it will never go above. That is (x ) will stay below. The condition for (x ) at x a is ( a) 0 < In the case of 0 (x ) will stay below when is negative and large. In the case of > > 0 we pick small > 0 satisfying sin 2 < sin 2. When is negative and large (x ) decreases with large negative slope -/ until and stays below thereafter. Therefore we conclude lim ( x ) 0 Appendix B: Proof of Rayleigh principle Rayleigh principle: - 5 -

6 Let W span{ v 0 v } W 2 uuc BC AMS 22A Applied Mathematical Methods I { and u W } Consider the minimization of R(u) over W. Let u + arg min Ru uw u0 Claim 2: u + is an eigenfunction corresponding to eigenvalue + and + Proof of Claim 2: Step : We first prove that min Ru uw u0 r( x) Lu W This is verified directly. for u W Let v W. By definition v is a linear combination of {v 0 v v }. v c n r( x) Lu v u r x L( v) u r( x) L c n u ( n )c n 0 Ste We consider a function of a scalar variable s gs Ru ( + + su) where u W By definition of u + we have g( 0) 0 for any u W The derivative of R(u + + s u) has the expression d ds Ru ( + su) + s0 2 u u + u + r( x) Lu ( +) + Ru ( + ) u u

7 AMS 22A Applied Mathematical Methods I 2 u + u + u r( x) Lu ( +)+ Ru ( + )u + The derivative of R(u + + s u) is zero for any u W. d ds Ru ( + + su) 0 for any u W s0 > u r( x) Lu ( +)+ Ru ( + )u + 0 for any u W From Step we know r( x) Lu ( +)+ Ru + r( x) Lu ( +)+ Ru ( + )u + 0 > Lu ( + ) R( u + )r( x)u + u + W. We conclude > R(u + ) is an eigenvalue and u + is a corresponding eigenfunction. By definition of u + we have Ru ( + ) min Ru uw > Ru ( + ) Rv ( + ) + (v + is an eigenfunction corresponding to + ). On the other hand + is the smallest eigenvalue whose eigenfunctions are orthogonal to W. R(u + ) is an eigenvalue whose eigenfunctions are orthogonal to W. > + Ru ( + ) Combining these two results we conclude + Ru ( + ) min Ru uw Appendix C: Proof of the completeness of eigenfunctions For f( x)c 2 BC we need to show that for any > 0 there exist and { c n } such that f c n < where the norm is based on the inner product - 7 -

8 AMS 22A Applied Mathematical Methods I Proof: Let u u u u x b a ( ) 2 r x dx b n f v n and h f b n. We only need to show that there exits such that h h < 2 We prove it through a sequence of steps. Step: h W Ste: It is straightforward to verify that h 0 0 n. Rayleigh principle tells us + min Ru. uw u0 > + Rh ( ) > + h r( x) L ( h ) h h > h h + h r( x) L ( h ) Step 3: Using the fact that is the eigenfunction corresponding to n we have f r( x) L ( ) f n n b n r( x) L ( ) n n Step 4: We use the result of Step 3 to calculate h r( x) L ( h ). h r( x) L ( h ) f b n r x L f b n - 8 -

9 AMS 22A Applied Mathematical Methods I f r( x) L( f) 2 f r( x) L b n + b n r x L b n Step 5: f r( x) L( f) + 2 n b n Combining the results of Ste and Step 4 we obtain h h f + r( x) L( f) + b n n 2 f + r( x) L( f) K 0 2 n b n where (K 0 + ) is the number of negative eigenvalues. Since K 0 is fixed and lim + + we conclude + h h 0 as

AMS 212A Applied Mathematical Methods I Lecture 06 Copyright by Hongyun Wang, UCSC. ( ), v (that is, 1 ( ) L i

AMS 212A Applied Mathematical Methods I Lecture 06 Copyright by Hongyun Wang, UCSC. ( ), v (that is, 1 ( ) L i AMS A Applied Mthemticl Methods I Lecture 6 Copyright y Hongyun Wng, UCSC Recp of Lecture 5 Clssifiction of oundry conditions Dirichlet eumnn Mixed Adjoint opertor, self-djoint opertor Sturm-Liouville

More information

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

Lecture 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 information

AMS 212A Applied Mathematical Methods I Lecture 14 Copyright by Hongyun Wang, UCSC. with initial condition

AMS 212A Applied Mathematical Methods I Lecture 14 Copyright by Hongyun Wang, UCSC. with initial condition Lecture 14 Copyright by Hongyun Wang, UCSC Recap of Lecture 13 Semi-linear PDE a( x, t) u t + b( x, t)u = c ( x, t, u ) x Characteristics: dt d = a x, t dx d = b x, t with initial condition t ( 0)= t 0

More information

Lecture 4 Orthonormal vectors and QR factorization

Lecture 4 Orthonormal vectors and QR factorization Orthonormal vectors and QR factorization 4 1 Lecture 4 Orthonormal vectors and QR factorization EE263 Autumn 2004 orthonormal vectors Gram-Schmidt procedure, QR factorization orthogonal decomposition induced

More information

Sturm-Liouville Theory

Sturm-Liouville Theory More on Ryan C. Trinity University Partial Differential Equations April 19, 2012 Recall: A Sturm-Liouville (S-L) problem consists of A Sturm-Liouville equation on an interval: (p(x)y ) + (q(x) + λr(x))y

More information

Midterm Solution

Midterm Solution 18303 Midterm Solution Problem 1: SLP with mixed boundary conditions Consider the following regular) Sturm-Liouville eigenvalue problem consisting in finding scalars λ and functions v : [0, b] R b > 0),

More information

Sturm-Liouville operators have form (given p(x) > 0, q(x)) + q(x), (notation means Lf = (pf ) + qf ) dx

Sturm-Liouville operators have form (given p(x) > 0, q(x)) + q(x), (notation means Lf = (pf ) + qf ) dx Sturm-Liouville operators Sturm-Liouville operators have form (given p(x) > 0, q(x)) L = d dx ( p(x) d ) + q(x), (notation means Lf = (pf ) + qf ) dx Sturm-Liouville operators Sturm-Liouville operators

More information

Math 3363 Examination I Solutions

Math 3363 Examination I Solutions Math 3363 Examination I Solutions Spring 218 Please use a pencil and do the problems in the order in which they are listed. No books, notes, calculators, cell phones, smart watches, or other electronics.

More information

Introduction to Sturm-Liouville Theory

Introduction to Sturm-Liouville Theory Introduction to R. C. Trinity University Partial Differential Equations April 10, 2014 Sturm-Liouville problems Definition: A (second order) Sturm-Liouville (S-L) problem consists of A Sturm-Liouville

More information

Ordinary Differential Equations II

Ordinary 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 information

In this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation.

In 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 information

Consequences of Orthogonality

Consequences 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 information

Math 250B Final Exam Review Session Spring 2015 SOLUTIONS

Math 250B Final Exam Review Session Spring 2015 SOLUTIONS Math 5B Final Exam Review Session Spring 5 SOLUTIONS Problem Solve x x + y + 54te 3t and y x + 4y + 9e 3t λ SOLUTION: We have det(a λi) if and only if if and 4 λ only if λ 3λ This means that the eigenvalues

More information

July 21 Math 2254 sec 001 Summer 2015

July 21 Math 2254 sec 001 Summer 2015 July 21 Math 2254 sec 001 Summer 2015 Section 8.8: Power Series Theorem: Let a n (x c) n have positive radius of convergence R, and let the function f be defined by this power series f (x) = a n (x c)

More information

Abstract Vector Spaces and Concrete Examples

Abstract Vector Spaces and Concrete Examples LECTURE 18 Abstract Vector Spaces and Concrete Examples Our discussion of linear algebra so far has been devoted to discussing the relations between systems of linear equations, matrices, and vectors.

More information

Physics 250 Green s functions for ordinary differential equations

Physics 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 information

PDE and Boundary-Value Problems Winter Term 2014/2015

PDE and Boundary-Value Problems Winter Term 2014/2015 PDE and Boundary-Value Problems Winter Term 2014/2015 Lecture 6 Saarland University 17. November 2014 c Daria Apushkinskaya (UdS) PDE and BVP lecture 6 17. November 2014 1 / 40 Purpose of Lesson To show

More information

MATH 412 Fourier Series and PDE- Spring 2010 SOLUTIONS to HOMEWORK 5

MATH 412 Fourier Series and PDE- Spring 2010 SOLUTIONS to HOMEWORK 5 MATH 4 Fourier Series PDE- Spring SOLUTIONS to HOMEWORK 5 Problem (a: Solve the following Sturm-Liouville problem { (xu + λ x u = < x < e u( = u (e = (b: Show directly that the eigenfunctions are orthogonal

More information

APPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.

APPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2. APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product

More information

MATH 590: Meshfree Methods

MATH 590: Meshfree Methods MATH 590: Meshfree Methods The Connection to Green s Kernels Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu MATH 590 1 Outline 1 Introduction

More information

3 Green s functions in 2 and 3D

3 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 information

Waves on 2 and 3 dimensional domains

Waves 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 information

Houston Journal of Mathematics. c 2015 University of Houston Volume 41, No. 4, 2015

Houston Journal of Mathematics. c 2015 University of Houston Volume 41, No. 4, 2015 Houston Journal of Mathematics c 25 University of Houston Volume 4, No. 4, 25 AN INCREASING FUNCTION WITH INFINITELY CHANGING CONVEXITY TONG TANG, YIFEI PAN, AND MEI WANG Communicated by Min Ru Abstract.

More information

AMS 147 Computational Methods and Applications Lecture 17 Copyright by Hongyun Wang, UCSC

AMS 147 Computational Methods and Applications Lecture 17 Copyright by Hongyun Wang, UCSC Lecture 17 Copyright by Hongyun Wang, UCSC Recap: Solving linear system A x = b Suppose we are given the decomposition, A = L U. We solve (LU) x = b in 2 steps: *) Solve L y = b using the forward substitution

More information

Math Theory of Partial Differential Equations Lecture 3-2: Spectral properties of the Laplacian. Bessel functions.

Math Theory of Partial Differential Equations Lecture 3-2: Spectral properties of the Laplacian. Bessel functions. Math 412-501 Theory of Partial ifferential Equations Lecture 3-2: Spectral properties of the Laplacian. Bessel functions. Eigenvalue problem: 2 φ + λφ = 0 in, ( αφ + β φ ) n = 0, where α, β are piecewise

More information

The total error in numerical differentiation

The total error in numerical differentiation AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and

More information

Lecture 16: Bessel s Inequality, Parseval s Theorem, Energy convergence

Lecture 16: Bessel s Inequality, Parseval s Theorem, Energy convergence Introductory lecture notes on Partial Differential Equations - c Anthony Peirce. ot to be copied, used, or revised without explicit written permission from the copyright owner. ecture 6: Bessel s Inequality,

More information

Section 12.6: Non-homogeneous Problems

Section 12.6: Non-homogeneous Problems Section 12.6: Non-homogeneous Problems 1 Introduction Up to this point all the problems we have considered are we what we call homogeneous problems. This means that for an interval < x < l the problems

More information

Physics 6303 Lecture 8 September 25, 2017

Physics 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 information

Economics 204 Fall 2011 Problem Set 2 Suggested Solutions

Economics 204 Fall 2011 Problem Set 2 Suggested Solutions Economics 24 Fall 211 Problem Set 2 Suggested Solutions 1. Determine whether the following sets are open, closed, both or neither under the topology induced by the usual metric. (Hint: think about limit

More information

Extreme Values and Positive/ Negative Definite Matrix Conditions

Extreme Values and Positive/ Negative Definite Matrix Conditions Extreme Values and Positive/ Negative Definite Matrix Conditions James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 016 Outline 1

More information

Elements of linear algebra

Elements 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 information

MAT137 Calculus! Lecture 45

MAT137 Calculus! Lecture 45 official website http://uoft.me/mat137 MAT137 Calculus! Lecture 45 Today: Taylor Polynomials Taylor Series Next: Taylor Series Power Series Definition (Power Series) A power series is a series of the form

More information

1 A complete Fourier series solution

1 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 information

Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series

Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series CHAPTER 5 Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series We start with some introductory examples. 5.. Cauchy s equation The homogeneous Euler-Cauchy equation (Leonhard

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information

2.4 Eigenvalue problems

2.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 information

UNIQUENESS OF POSITIVE SOLUTION TO SOME COUPLED COOPERATIVE VARIATIONAL ELLIPTIC SYSTEMS

UNIQUENESS OF POSITIVE SOLUTION TO SOME COUPLED COOPERATIVE VARIATIONAL ELLIPTIC SYSTEMS TRANSACTIONS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9947(XX)0000-0 UNIQUENESS OF POSITIVE SOLUTION TO SOME COUPLED COOPERATIVE VARIATIONAL ELLIPTIC SYSTEMS YULIAN

More information

Partial Differential Equations

Partial Differential Equations Partial Differential Equations Spring Exam 3 Review Solutions Exercise. We utilize the general solution to the Dirichlet problem in rectangle given in the textbook on page 68. In the notation used there

More information

AMS 147 Computational Methods and Applications Lecture 13 Copyright by Hongyun Wang, UCSC

AMS 147 Computational Methods and Applications Lecture 13 Copyright by Hongyun Wang, UCSC Lecture 13 Copyright y Hongyun Wang, UCSC Recap: Fitting to exact data *) Data: ( x j, y j ), j = 1,,, N y j = f x j *) Polynomial fitting Gis phenomenon *) Cuic spline Convergence of cuic spline *) Application

More information

Physics 137A Quantum Mechanics Fall 2012 Midterm II - Solutions

Physics 137A Quantum Mechanics Fall 2012 Midterm II - Solutions Physics 37A Quantum Mechanics Fall 0 Midterm II - Solutions These are the solutions to the exam given to Lecture Problem [5 points] Consider a particle with mass m charge q in a simple harmonic oscillator

More information

Math 4263 Homework Set 1

Math 4263 Homework Set 1 Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that

More information

EIGENVALUES AND EIGENVECTORS

EIGENVALUES 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 information

LECTURE 10: REVIEW OF POWER SERIES. 1. Motivation

LECTURE 10: REVIEW OF POWER SERIES. 1. Motivation LECTURE 10: REVIEW OF POWER SERIES By definition, a power series centered at x 0 is a series of the form where a 0, a 1,... and x 0 are constants. For convenience, we shall mostly be concerned with the

More information

Putnam f(x) g(x) + h(x) = 3x +2 if 1 x 0 2x +2 if x > 0

Putnam f(x) g(x) + h(x) = 3x +2 if 1 x 0 2x +2 if x > 0 Putnam 999 A. Find polynomials f(x),g(x), andh(x), iftheyexist, suchthatforallx, if x< f(x) g(x) + h(x) 3x + if x 0 x + if x > 0 Solution. Let h(x) ax + b. Note that there are corners in the graph of this

More information

Math123 Lecture 1. Dr. Robert C. Busby. Lecturer: Office: Korman 266 Phone :

Math123 Lecture 1. Dr. Robert C. Busby. Lecturer: Office: Korman 266 Phone : Lecturer: Math1 Lecture 1 Dr. Robert C. Busby Office: Korman 66 Phone : 15-895-1957 Email: rbusby@mcs.drexel.edu Course Web Site: http://www.mcs.drexel.edu/classes/calculus/math1_spring0/ (Links are case

More information

j=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent.

j=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent. Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. Let u = [u

More information

Name: Math Homework Set # 5. March 12, 2010

Name: Math Homework Set # 5. March 12, 2010 Name: Math 4567. Homework Set # 5 March 12, 2010 Chapter 3 (page 79, problems 1,2), (page 82, problems 1,2), (page 86, problems 2,3), Chapter 4 (page 93, problems 2,3), (page 98, problems 1,2), (page 102,

More information

MAT137 Calculus! Lecture 48

MAT137 Calculus! Lecture 48 official website http://uoft.me/mat137 MAT137 Calculus! Lecture 48 Today: Taylor Series Applications Next: Final Exams Important Taylor Series and their Radii of Convergence 1 1 x = e x = n=0 n=0 x n n!

More information

Reflected Brownian Motion

Reflected Brownian Motion Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide

More information

Through level curves and contour lines give some information about the function

Through level curves and contour lines give some information about the function Module 10 : Scaler fields, Limit and Continuity Lecture 29 : Limit of scaler fields [Section 29.1] Objectives In this section you will learn the following : The notion of limits for scalar fields. 29.1

More information

Physics 6303 Lecture 9 September 17, ct' 2. ct' ct'

Physics 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 information

1. What is the determinant of the following matrix? a 1 a 2 4a 3 2a 2 b 1 b 2 4b 3 2b c 1. = 4, then det

1. What is the determinant of the following matrix? a 1 a 2 4a 3 2a 2 b 1 b 2 4b 3 2b c 1. = 4, then det What is the determinant of the following matrix? 3 4 3 4 3 4 4 3 A 0 B 8 C 55 D 0 E 60 If det a a a 3 b b b 3 c c c 3 = 4, then det a a 4a 3 a b b 4b 3 b c c c 3 c = A 8 B 6 C 4 D E 3 Let A be an n n matrix

More information

Physics 6303 Lecture 15 October 10, Reminder of general solution in 3-dimensional cylindrical coordinates. sinh. sin

Physics 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 information

DS-GA 1002 Lecture notes 10 November 23, Linear models

DS-GA 1002 Lecture notes 10 November 23, Linear models DS-GA 2 Lecture notes November 23, 2 Linear functions Linear models A linear model encodes the assumption that two quantities are linearly related. Mathematically, this is characterized using linear functions.

More information

Vectors in Function Spaces

Vectors 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 information

Analysis II - few selective results

Analysis II - few selective results Analysis II - few selective results Michael Ruzhansky December 15, 2008 1 Analysis on the real line 1.1 Chapter: Functions continuous on a closed interval 1.1.1 Intermediate Value Theorem (IVT) Theorem

More information

Some definitions. Math 1080: Numerical Linear Algebra Chapter 5, Solving Ax = b by Optimization. A-inner product. Important facts

Some definitions. Math 1080: Numerical Linear Algebra Chapter 5, Solving Ax = b by Optimization. A-inner product. Important facts Some definitions Math 1080: Numerical Linear Algebra Chapter 5, Solving Ax = b by Optimization M. M. Sussman sussmanm@math.pitt.edu Office Hours: MW 1:45PM-2:45PM, Thack 622 A matrix A is SPD (Symmetric

More information

Math 255 Honors: Gram-Schmidt Orthogonalization on the Space of Polynomials

Math 255 Honors: Gram-Schmidt Orthogonalization on the Space of Polynomials Math 55 Honors: Gram-Schmidt Orthogonalization on the Space of Polynomials David Moore May, 03 Abstract Gram-Schmidt Orthogonalization is a process to construct orthogonal vectors from some basis for a

More information

Greedy Algorithms. CSE 101: Design and Analysis of Algorithms Lecture 10

Greedy Algorithms. CSE 101: Design and Analysis of Algorithms Lecture 10 Greedy Algorithms CSE 101: Design and Analysis of Algorithms Lecture 10 CSE 101: Design and analysis of algorithms Greedy algorithms Reading: Kleinberg and Tardos, sections 4.1, 4.2, and 4.3 Homework 4

More information

1 Solutions to selected problems

1 Solutions to selected problems Solutions to selected problems Section., #a,c,d. a. p x = n for i = n : 0 p x = xp x + i end b. z = x, y = x for i = : n y = y + x i z = zy end c. y = (t x ), p t = a for i = : n y = y(t x i ) p t = p

More information

23 Elements of analytic ODE theory. Bessel s functions

23 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 information

MAT137 Calculus! Lecture 6

MAT137 Calculus! Lecture 6 MAT137 Calculus! Lecture 6 Today: 3.2 Differentiation Rules; 3.3 Derivatives of higher order. 3.4 Related rates 3.5 Chain Rule 3.6 Derivative of Trig. Functions Next: 3.7 Implicit Differentiation 4.10

More information

Completion Date: Monday February 11, 2008

Completion Date: Monday February 11, 2008 MATH 4 (R) Winter 8 Intermediate Calculus I Solutions to Problem Set #4 Completion Date: Monday February, 8 Department of Mathematical and Statistical Sciences University of Alberta Question. [Sec..9,

More information

Lecture No. 1 Introduction to Method of Weighted Residuals. Solve the differential equation L (u) = p(x) in V where L is a differential operator

Lecture No. 1 Introduction to Method of Weighted Residuals. Solve the differential equation L (u) = p(x) in V where L is a differential operator Lecture No. 1 Introduction to Method of Weighted Residuals Solve the differential equation L (u) = p(x) in V where L is a differential operator with boundary conditions S(u) = g(x) on Γ where S is a differential

More information

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12 Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider

More information

1 Overview. 2 A Characterization of Convex Functions. 2.1 First-order Taylor approximation. AM 221: Advanced Optimization Spring 2016

1 Overview. 2 A Characterization of Convex Functions. 2.1 First-order Taylor approximation. AM 221: Advanced Optimization Spring 2016 AM 221: Advanced Optimization Spring 2016 Prof. Yaron Singer Lecture 8 February 22nd 1 Overview In the previous lecture we saw characterizations of optimality in linear optimization, and we reviewed the

More information

Eigenvalues and Eigenvectors

Eigenvalues and Eigenvectors LECTURE 3 Eigenvalues and Eigenvectors Definition 3.. Let A be an n n matrix. The eigenvalue-eigenvector problem for A is the problem of finding numbers λ and vectors v R 3 such that Av = λv. If λ, v are

More information

Numerical Optimization

Numerical Optimization Numerical Optimization Unit 2: Multivariable optimization problems Che-Rung Lee Scribe: February 28, 2011 (UNIT 2) Numerical Optimization February 28, 2011 1 / 17 Partial derivative of a two variable function

More information

Eigenvalue comparisons in graph theory

Eigenvalue comparisons in graph theory Eigenvalue comparisons in graph theory Gregory T. Quenell July 1994 1 Introduction A standard technique for estimating the eigenvalues of the Laplacian on a compact Riemannian manifold M with bounded curvature

More information

IV. Transport Phenomena Lecture 18: Forced Convection in Fuel Cells II

IV. Transport Phenomena Lecture 18: Forced Convection in Fuel Cells II IV. Transport Phenomena Lecture 18: Forced Convection in Fuel Cells II MIT Student (and MZB) As discussed in the previous lecture, we are interested in forcing fluid to flow in a fuel cell in order to

More information

Definition and basic properties of heat kernels I, An introduction

Definition and basic properties of heat kernels I, An introduction Definition and basic properties of heat kernels I, An introduction Zhiqin Lu, Department of Mathematics, UC Irvine, Irvine CA 92697 April 23, 2010 In this lecture, we will answer the following questions:

More information

v(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).

v(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 information

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.

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. 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. For example f(x) = 1 1 x = 1 + x + x2 + x 3 + = ln(1 + x) = x x2 2

More information

Functional Analysis HW #5

Functional Analysis HW #5 Functional Analysis HW #5 Sangchul Lee October 29, 2015 Contents 1 Solutions........................................ 1 1 Solutions Exercise 3.4. Show that C([0, 1]) is not a Hilbert space, that is, there

More information

LEGENDRE 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 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 information

MATH The Chain Rule Fall 2016 A vector function of a vector variable is a function F: R n R m. In practice, if x 1, x n is the input,

MATH The Chain Rule Fall 2016 A vector function of a vector variable is a function F: R n R m. In practice, if x 1, x n is the input, MATH 20550 The Chain Rule Fall 2016 A vector function of a vector variable is a function F: R n R m. In practice, if x 1, x n is the input, F(x 1,, x n ) F 1 (x 1,, x n ),, F m (x 1,, x n ) where each

More information

Second Order Elliptic PDE

Second Order Elliptic PDE Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic

More information

Eigenvalues and Eigenfunctions of the Laplacian

Eigenvalues and Eigenfunctions of the Laplacian The Waterloo Mathematics Review 23 Eigenvalues and Eigenfunctions of the Laplacian Mihai Nica University of Waterloo mcnica@uwaterloo.ca Abstract: The problem of determining the eigenvalues and eigenvectors

More information

Physical Metaphors for Graphs

Physical Metaphors for Graphs Graphs and Networks Lecture 3 Physical Metaphors for Graphs Daniel A. Spielman October 9, 203 3. Overview We will examine physical metaphors for graphs. We begin with a spring model and then discuss resistor

More information

100 CHAPTER 4. SYSTEMS AND ADAPTIVE STEP SIZE METHODS APPENDIX

100 CHAPTER 4. SYSTEMS AND ADAPTIVE STEP SIZE METHODS APPENDIX 100 CHAPTER 4. SYSTEMS AND ADAPTIVE STEP SIZE METHODS APPENDIX.1 Norms If we have an approximate solution at a given point and we want to calculate the absolute error, then we simply take the magnitude

More information

Welcome to Math 257/316 - Partial Differential Equations

Welcome 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 information

Introduction to Hilbert spaces. I.

Introduction to Hilbert spaces. I. Teoretisk Fysik, KTH September 5, 26 5A135 Fysikens matematiska metoder, ht6. Introduction to Hilbert spaces. I. My aim here is to give a down-to-earth introduction to the theory of Hilbert spaces which

More information

Unit 2: Solving Scalar Equations. Notes prepared by: Amos Ron, Yunpeng Li, Mark Cowlishaw, Steve Wright Instructor: Steve Wright

Unit 2: Solving Scalar Equations. Notes prepared by: Amos Ron, Yunpeng Li, Mark Cowlishaw, Steve Wright Instructor: Steve Wright cs416: introduction to scientific computing 01/9/07 Unit : Solving Scalar Equations Notes prepared by: Amos Ron, Yunpeng Li, Mark Cowlishaw, Steve Wright Instructor: Steve Wright 1 Introduction We now

More information

( ), λ. ( ) =u z. ( )+ iv z

( ), λ. ( ) =u z. ( )+ iv z AMS 212B Perturbation Metho Lecture 18 Copyright by Hongyun Wang, UCSC Method of steepest descent Consider the integral I = f exp( λ g )d, g C =u + iv, λ We consider the situation where ƒ() and g() = u()

More information

MATH 425, HOMEWORK 3 SOLUTIONS

MATH 425, HOMEWORK 3 SOLUTIONS MATH 425, HOMEWORK 3 SOLUTIONS Exercise. (The differentiation property of the heat equation In this exercise, we will use the fact that the derivative of a solution to the heat equation again solves the

More information

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

Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0 Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can

More information

DS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.

DS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra. DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1

More information

An Inverse Problem for the Matrix Schrödinger Equation

An Inverse Problem for the Matrix Schrödinger Equation Journal of Mathematical Analysis and Applications 267, 564 575 (22) doi:1.16/jmaa.21.7792, available online at http://www.idealibrary.com on An Inverse Problem for the Matrix Schrödinger Equation Robert

More information

Numerical Methods. King Saud University

Numerical Methods. King Saud University Numerical Methods King Saud University Aims In this lecture, we will... find the approximate solutions of derivative (first- and second-order) and antiderivative (definite integral only). Numerical Differentiation

More information

Theorem (4.11). Let M be a closed subspace of a Hilbert space H. For any x, y E, the parallelogram law applied to x/2 and y/2 gives.

Theorem (4.11). Let M be a closed subspace of a Hilbert space H. For any x, y E, the parallelogram law applied to x/2 and y/2 gives. Math 641 Lecture #19 4.11, 4.12, 4.13 Inner Products and Linear Functionals (Continued) Theorem (4.10). Every closed, convex set E in a Hilbert space H contains a unique element of smallest norm, i.e.,

More information

Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem

Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem L u x f x BC u x g x with the weak problem find u V such that B u,v

More information

1 Solutions in cylindrical coordinates: Bessel functions

1 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 information

Lecture 8: Basic convex analysis

Lecture 8: Basic convex analysis Lecture 8: Basic convex analysis 1 Convex sets Both convex sets and functions have general importance in economic theory, not only in optimization. Given two points x; y 2 R n and 2 [0; 1]; the weighted

More information

Escape Problems & Stochastic Resonance

Escape Problems & Stochastic Resonance Escape Problems & Stochastic Resonance 18.S995 - L08 & 09 dunkel@mit.edu Examples Illustration by J.P. Cartailler. Copyright 2007, Symmation LLC. transport & evolution: stochastic escape problems http://evolutionarysystemsbiology.org/intro/

More information

LECTURE 33: NONHOMOGENEOUS HEAT CONDUCTION PROBLEM

LECTURE 33: NONHOMOGENEOUS HEAT CONDUCTION PROBLEM LECTURE 33: NONHOMOGENEOUS HEAT CONDUCTION PROBLEM 1. General Solving Procedure The general nonhomogeneous 1-dimensional heat conduction problem takes the form Eq : [p(x)u x ] x q(x)u + F (x, t) = r(x)u

More information

TEST CODE: MIII (Objective type) 2010 SYLLABUS

TEST CODE: MIII (Objective type) 2010 SYLLABUS TEST CODE: MIII (Objective type) 200 SYLLABUS Algebra Permutations and combinations. Binomial theorem. Theory of equations. Inequalities. Complex numbers and De Moivre s theorem. Elementary set theory.

More information

Generalized trigonometric functions

Generalized trigonometric functions 7/211, Santiago de Compostela (Etwas allgemein machen heisst, es denken. G.W.F. Hegel) Based on: * David E. Edmunds, Petr Gurka, J.L., Properties of generalized Trigonometric functions, Preprint * David

More information

AMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.

AMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error. Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error

More information

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

Math Review for Exam Answer each of the following questions as either True or False. Circle the correct answer. Math 22 - Review for Exam 3. Answer each of the following questions as either True or False. Circle the correct answer. (a) True/False: If a n > 0 and a n 0, the series a n converges. Soln: False: Let

More information