AMS 212A Applied Mathematical Methods I Appendices of Lecture 06 Copyright by Hongyun Wang, UCSC. ( ) cos2
|
|
- Tabitha Dorsey
- 5 years ago
- Views:
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 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 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 informationAMS 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 informationLecture 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 informationSturm-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 informationMidterm 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 informationSturm-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 informationMath 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 informationIntroduction 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 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 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 informationMath 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 informationJuly 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 informationAbstract 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 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 informationPDE 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 informationMATH 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 informationAPPENDIX 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 informationMATH 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 information3 Green s functions in 2 and 3D
William J. Parnell: MT34032. Section 3: Green s functions in 2 and 3 57 3 Green s functions in 2 and 3 Unlike the one dimensional case where Green s functions can be found explicitly for a number of different
More 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 informationHouston 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 informationAMS 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 informationMath 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 informationThe 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 informationLecture 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 informationSection 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 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 informationEconomics 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 informationExtreme 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 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 informationMAT137 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 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 informationIntroduction 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 informationSturm-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 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 informationUNIQUENESS 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 informationPartial 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 informationAMS 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 informationPhysics 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 informationMath 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 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 informationLECTURE 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 informationPutnam 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 informationMath123 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 informationj=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 informationName: 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 informationMAT137 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 informationReflected 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 informationThrough 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 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 information1. 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 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 informationDS-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 informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationAnalysis 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 informationSome 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 informationMath 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 informationGreedy 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 information1 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 information23 Elements of analytic ODE theory. Bessel s functions
23 Elements of analytic ODE theory. Bessel s functions Recall I am changing the variables) that we need to solve the so-called Bessel s equation 23. Elements of analytic ODE theory Let x 2 u + xu + x 2
More informationMAT137 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 informationCompletion 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 informationLecture 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 informationChapter 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 information1 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 informationEigenvalues 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 informationNumerical 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 informationEigenvalue 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 informationIV. 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 informationDefinition 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 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 informationLecture 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 informationFunctional 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 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 informationMATH 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 informationSecond 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 informationEigenvalues 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 informationPhysical 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 information100 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 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 informationIntroduction 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 informationUnit 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
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 informationMATH 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 informationPower 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 informationDS-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 informationAn 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 informationNumerical 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 informationTheorem (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 informationSpectral 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 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 informationLecture 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 informationEscape 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 informationLECTURE 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 informationTEST 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 informationGeneralized 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 informationAMS 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 informationMath 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