Bifurcation Theory. Solutions to Problem Sheet 2. u(0) = u(1) = 0

Size: px
Start display at page:

Download "Bifurcation Theory. Solutions to Problem Sheet 2. u(0) = u(1) = 0"

Transcription

1 Karlsruhe Institute of Technology (KIT) Institute for Analysis Dr. Rainer Mandel M.Sc. Janina Gärtner, M.Sc. Dominic Scheider Summer Term Bifurcation Theory Solutions to Problem Sheet Problem 4 (On the Energy Method) Find such continuous functions g : R R R with g(, µ) = for all µ R that { u = g(u, µ) in (, 1), u() = u(1) = (1) has the asserted properties. (a) There is such µ R that (1) does not admit a nontrivial solution for µ µ. (b) There exist µ R and sequences (u n, µ n ) n N, (ũ n, µ n ) n N C ([, 1]) R of solutions of (1) which bifurcate from the trivial family at (, µ ) and satisfy u n, ũ n, µ n < µ < µ n for all n N. (c) We have g C (R R, R) and there exist α >, j N, λ R with g( z, µ ) = g(z, µ ), g(z, µ )z > for < z < α ; (G(α, µ ) G(z, µ )) 1 < for < α < α ; g z (, µ ) = π (j + 1). According to Theorem II.3 from the lecture, the transversality condition g zµ (, µ ) implies that j-nodal solutions of (1) bifurcate from (, µ ) in C ([, 1]). Demonstrate that bifurcation need not occur if the transversality condition is violated. Solution (a) For arbitrary µ R, consider g(z, µ) := (µ µ)z. Then a direct calculation shows that, for µ µ, the differential equation in (1) has a solution of the form u(t) = α cosh( µ µ t) + β sinh( µ µ t) with coefficients α, β R. The homogeneous boundary conditions now give = u() = α, = u(1) = α cosh( µ µ ) + β sinh( µ µ ) = β sinh( µ µ ), hence α = β = and u. Bifurcation Theory.5.17

2 (b) Consider, as a slight modification of the pendulum equation discussed in the lecture, g(z, µ) = (λ j + (µ λ j ) ) sin(z) where λ j = ((j + 1)π) for some j N, and consider µ := λ j. We know from the lecture that, for u = λ sin(u), u() = u(1) =, j-nodal solutions bifurcate at the point (, λ j ), parametrized as (α, λ j (α)) <α<α with λ j (α) λ j as α. With g chosen as above, we find two parameters λ corresponding to each value of α via λ j (α) = λ j + (µ λ j ) µ = λ j ± λ j (α) λ j, which yields to distinct families of bifurcating j-nodal solutions of (1) parametrized as ( ) ( ) α, λ j + λ j (α) λ j, α, λ j λ j (α) λ j. <α<α <α<α In the context of this problem, we choose α n, < α n < α and let µ n := λ j λ j (α n ) λ j, µ n := λ j + λ j (α n ) λ j, u n = ũ n with the suitable j-nodal solution of u = λ j (α n ) sin(u), u() = u(1) = with u n = α n. (c) Again, we modify the pendulum equation. We introduce g(z, µ) = (λ j (µ λ j ) ) sin(z) where λ j = ((j + 1)π) for some j N, and consider µ := λ j. However, using the notation from part (b), the equation λ j (α) = λ j (µ λ j ) does not have solutions, and the bifurcation diagram for u = λ sin(u), u() = u(1) = from the lecture reveals that there is no bifurcation of (1) at (, λ j ). As in the lecture, we see that all conditions of Theorem II.3 are satisfied, in particular but for the transversality condition g z (, λ j ) = λ j cos() = π (j + 1), g zµ (, λ j ) = (λ j λ j ) cos(1) =. Bifurcation Theory.5.17

3 Problem 5 (An application of the Energy Method) Consider the boundary value problem { u + λ(u u 3 ) = u() = u(π) =. in (, π), () Prove that, for k N, nontrivial (k 1)-nodal solutions of () bifurcate from the trivial branch at the point λ k = k. Solution We intend to apply Theorem II.3 and define g : R R R, g(z, λ) := λ(z z 3 ). Then g C (R R, R), g(, λ) = for every λ R. Fixing λ >, we check the assumptions of Theorem II.3: (i) For all z R, we have g( z, λ) = g(z, λ). (ii) g(z, λ) z = λ(z z 4 ) = λz (1 z ) > holds for z R with < z < 1. (iii) We let G(z, λ) := z g(s, λ)ds = 1 4 λ(z z 4 ) be the primitive of g(, λ). For α (, 1), we calculate (G(α, λ) G(z, λ)) 1 α = λ (α α 4 z + z 4 ) 1 = ((1 z ) (1 α ) ) 1. λ In order to prove finiteness of this integral, we use the mean value theorem to find ξ z (z, α) with (1 z ) (1 α ) d = (z α) dξ (1 ξ ) = (α z) 4ξ z (1 ξz) ξ=ξz and then estimate as follows: (α z) 4z(1 α ), ((1 z ) (1 α ) ) 1 ((α z) 4z(1 α )) 1 = = α + (α z) z 1 α α + α z 1 α α 1 α z = α α α z z 1 α α z α 1 α α 1 α α <. Bifurcation Theory.5.17

4 As g C 1 (R R, R), Theorem II.3 (with T = π and α = 1) gives the following necessary condition for bifurcation of (k 1)-nodal solutions from the trivial family at a point (, λ k ): z g(, λ k ) = ( ) πk, ( ) T and equivalently λ k = k. Further, since the transversality condition (ii) of Theorem II.3, λ z g(, λ k ) = 1 holds, the condition ( ) is sufficient for bifurcation, and hence the proof is complete. Bifurcation Theory.5.17

5 Problem 6 (A problem without periodic Solutions) Let g C(R R, R) with g(, λ) = (λ R). Consider b C(R R, R) with b(x, λ) for all x, λ R. Prove that the differential equation u + b(x, λ) u = g(u, λ) in R, (3) does not admit a periodic solution in C (R) unless it is constant. Solution First, we note that by assumption, b is continuous and does not have any zero, so b is either negative or positive on all of R R. We assume that u C (R) is a periodic solution of (3). We introduce E : R R, E(t) := 1 [u (t)] + G(u(t), λ) where G(z, λ) := z g(s, λ) ds for λ, z R. Then, E C1 (R), and for t R E (t) = u (t) (u (t) + g(u(t), λ)) = b(t, λ) [u (t)] ( ) where we have inserted the differential equation in the last step. Since b does not change sign, this implies that E is a monotone function. As u is periodic, so is E, and we conclude that E is constant. Hence, E, and by ( ), u. So u is constant, and the assertion is proved. Bifurcation Theory.5.17

6 Problem 7 (Variational Methods) Let n N. We consider the n-dimensional bifurcation problem x I(x, λ) = (4) where I : R n R R, I(x, λ) := 1 x Ax λ sin( x p ) for some p (1, ) and a symmetric, positive definite matrix A R n n. (a) Show that, for λ R, I(, λ) is differentiable with respect to x and calculate x I(, λ). Verify that x I(, λ) = for all λ R. (b) Prove that, for λ R, I(, λ) admits a minimizer x λ R n. (c) Prove that (, ) is a bifurcation point for problem (4). Solution (a) We fix λ R and i {1,..., n}. For x R n, we have I(x, λ) = 1 n a jk x j x k λ sin( x p ). j,k=1 For x, the chain rule gives partial differentiability and, using the symmetry of A, I xi (x, λ) = 1 = n a jk (δ ij x k + x j δ ik ) λ cos( x p ) px i x p j,k=1 For x =, we compute directly n a ij x j λ cos( x p ) px i x p = j I(he i, λ) I(, λ) I xi (, λ) = lim h h lim I x i (h, λ) = h [ ] Ax λ cos( x p ) p x p x. i 1 = lim a iih λ sin( h p ) h h since p > 1. We conclude that I is indeed continuously differentiable with { Ax λ cos( x p ) p x p x, λ R, x, x I(x, λ) =, λ R, x = =, and and in particular, we have the trivial family of solutions x I(, λ) =, λ R. Bifurcation Theory.5.17

7 (b) For x R n, λ > we estimate I(x, λ) α x λ sin( x p ) λ, ( ) hence j λ := inf x R n I(x, λ) λ and we find a minimizing sequence (x (n) λ ) n N, i.e. lim n I(x(n) λ, λ) = j λ. Taking another look at the first estimate in ( ), we see that (x (n) λ ) n N is necessarily bounded and hence, working in finite dimension, the Bolzano-Weierstrass theorem asserts that we can assume without loss of generality (possibly after passing to a subsequence) that (x (n) λ ) n N converges to some element x λ R n. It then follows from the continuity of I that j λ = lim n I(x (n) λ, λ) = I(x λ, λ), and thus x λ is a minimizer for I(, λ). (c) To find non-trivial solutions of (4), we consider the minimizers x λ of I(, λ), which are stationary points of I(, λ) and hence provide solutions of (4). We will prove the following: Assertion 1 : For every λ >, x λ. Assertion : We have x λ as λ. This finally yields that (, ) is a bifurcation point of problem (4) since it is a point where the non-trivial family {(x λ, λ) : λ > } meets the trivial branch of solutions. We now prove the assertions 1 and. First of all, since A is positive definite, we find β α > with α x x Ax β x for x R n. ( ) Proof of Assertion 1 : If x λ = held true, then we would find j λ = and therefore I(, λ). We aim to show, however, that I(, λ) attains negative values. To this end, we recall that sin(ξ) [ π ξ for ξ, π ]. For x R n with x p π, we estimate I(x, λ) β x λ π x p = β ( x p x p 4λ ). βπ { ( Hence, choosing x R n with < x < min π ) 1 p, ( 4λ βπ ) } 1 p, we have I(x, λ) <, and the assertion is proved. Note that such a choice is possible because p <. Bifurcation Theory.5.17

8 Proof of Assertion : Since x λ, we have for all λ > x I(x λ, λ) Ax λ = λ cos( x λ p ) p x λ p x λ. Multiplying from the left with x λ, we find in view of ( ): ( <) α λ cos( x λ p ) p x λ p β. Letting λ, the estimate from below can only hold if cos( x λ p ) x λ p. As the cosine is bounded and p <, we infer x λ as λ. Bifurcation Theory

Ch 10.1: Two Point Boundary Value Problems

Ch 10.1: Two Point Boundary Value Problems Ch 10.1: Two Point Boundary Value Problems In many important physical problems there are two or more independent variables, so the corresponding mathematical models involve partial differential equations.

More information

MATH 220 solution to homework 4

MATH 220 solution to homework 4 MATH 22 solution to homework 4 Problem. Define v(t, x) = u(t, x + bt), then v t (t, x) a(x + u bt) 2 (t, x) =, t >, x R, x2 v(, x) = f(x). It suffices to show that v(t, x) F = max y R f(y). We consider

More information

Chapter 1: Systems of Linear Equations

Chapter 1: Systems of Linear Equations Chapter : Systems of Linear Equations February, 9 Systems of linear equations Linear systems Lecture A linear equation in variables x, x,, x n is an equation of the form a x + a x + + a n x n = b, where

More information

Convex Optimization M2

Convex Optimization M2 Convex Optimization M2 Lecture 3 A. d Aspremont. Convex Optimization M2. 1/49 Duality A. d Aspremont. Convex Optimization M2. 2/49 DMs DM par email: dm.daspremont@gmail.com A. d Aspremont. Convex Optimization

More information

Hyperbolic Functions and the Twin Paradox

Hyperbolic Functions and the Twin Paradox Hyperbolic Functions and the Twin Paradox Basics of hyperbolic functions The basic hyperbolic functions are defined as cosh x 1 2 (ex + e x ), sinh 1 2 (ex e x ). The notation ch, sh is also used (especially

More information

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential

More information

Proof. We indicate by α, β (finite or not) the end-points of I and call

Proof. We indicate by α, β (finite or not) the end-points of I and call C.6 Continuous functions Pag. 111 Proof of Corollary 4.25 Corollary 4.25 Let f be continuous on the interval I and suppose it admits non-zero its (finite or infinite) that are different in sign for x tending

More information

Linear Differential Equations. Problems

Linear Differential Equations. Problems Chapter 1 Linear Differential Equations. Problems 1.1 Introduction 1.1.1 Show that the function ϕ : R R, given by the expression ϕ(t) = 2e 3t for all t R, is a solution of the Initial Value Problem x =

More information

Lecture: Duality.

Lecture: Duality. Lecture: Duality http://bicmr.pku.edu.cn/~wenzw/opt-2016-fall.html Acknowledgement: this slides is based on Prof. Lieven Vandenberghe s lecture notes Introduction 2/35 Lagrange dual problem weak and strong

More information

1. Solve the boundary-value problems or else show that no solutions exist. y (x) = c 1 e 2x + c 2 e 3x. (3)

1. Solve the boundary-value problems or else show that no solutions exist. y (x) = c 1 e 2x + c 2 e 3x. (3) Diff. Eqns. Problem Set 6 Solutions. Solve the boundary-value problems or else show that no solutions exist. a y + y 6y, y, y 4 b y + 9y x + e x, y, yπ a Assuming y e rx is a solution, we get the characteristic

More information

1 + lim. n n+1. f(x) = x + 1, x 1. and we check that f is increasing, instead. Using the quotient rule, we easily find that. 1 (x + 1) 1 x (x + 1) 2 =

1 + lim. n n+1. f(x) = x + 1, x 1. and we check that f is increasing, instead. Using the quotient rule, we easily find that. 1 (x + 1) 1 x (x + 1) 2 = Chapter 5 Sequences and series 5. Sequences Definition 5. (Sequence). A sequence is a function which is defined on the set N of natural numbers. Since such a function is uniquely determined by its values

More information

INDIAN INSTITUTE OF TECHNOLOGY BOMBAY MA205 Complex Analysis Autumn 2012

INDIAN INSTITUTE OF TECHNOLOGY BOMBAY MA205 Complex Analysis Autumn 2012 INDIAN INSTITUTE OF TECHNOLOGY BOMBAY MA205 Complex Analysis Autumn 2012 September 5, 2012 Mapping Properties Lecture 13 We shall once again return to the study of general behaviour of holomorphic functions

More information

Analysis III Solutions - Serie 12

Analysis III Solutions - Serie 12 .. Necessary condition Let us consider the following problem for < x, y < π, u =, for < x, y < π, u y (x, π) = x a, for < x < π, u y (x, ) = a x, for < x < π, u x (, y) = u x (π, y) =, for < y < π. Find

More information

Ordinary Differential Equations (ODEs)

Ordinary Differential Equations (ODEs) Chapter 13 Ordinary Differential Equations (ODEs) We briefly review how to solve some of the most standard ODEs. 13.1 First Order Equations 13.1.1 Separable Equations A first-order ordinary differential

More information

The Implicit and Inverse Function Theorems Notes to supplement Chapter 13.

The Implicit and Inverse Function Theorems Notes to supplement Chapter 13. The Implicit and Inverse Function Theorems Notes to supplement Chapter 13. Remark: These notes are still in draft form. Examples will be added to Section 5. If you see any errors, please let me know. 1.

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

Solving Quadratic Equations

Solving Quadratic Equations Solving Quadratic Equations MATH 101 College Algebra J. Robert Buchanan Department of Mathematics Summer 2012 Objectives In this lesson we will learn to: solve quadratic equations by factoring, solve quadratic

More information

Convex Optimization Boyd & Vandenberghe. 5. Duality

Convex Optimization Boyd & Vandenberghe. 5. Duality 5. Duality Convex Optimization Boyd & Vandenberghe Lagrange dual problem weak and strong duality geometric interpretation optimality conditions perturbation and sensitivity analysis examples generalized

More information

5. Duality. Lagrangian

5. Duality. Lagrangian 5. Duality Convex Optimization Boyd & Vandenberghe Lagrange dual problem weak and strong duality geometric interpretation optimality conditions perturbation and sensitivity analysis examples generalized

More information

A Variational Analysis of a Gauged Nonlinear Schrödinger Equation

A Variational Analysis of a Gauged Nonlinear Schrödinger Equation A Variational Analysis of a Gauged Nonlinear Schrödinger Equation Alessio Pomponio, joint work with David Ruiz Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Variational and Topological

More information

Lecture Notes on PDEs

Lecture Notes on PDEs Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential

More information

Mathematical Methods and its Applications (Solution of assignment-12) Solution 1 From the definition of Fourier transforms, we have.

Mathematical Methods and its Applications (Solution of assignment-12) Solution 1 From the definition of Fourier transforms, we have. For 2 weeks course only Mathematical Methods and its Applications (Solution of assignment-2 Solution From the definition of Fourier transforms, we have F e at2 e at2 e it dt e at2 +(it/a dt ( setting (

More information

2. Dual space is essential for the concept of gradient which, in turn, leads to the variational analysis of Lagrange multipliers.

2. Dual space is essential for the concept of gradient which, in turn, leads to the variational analysis of Lagrange multipliers. Chapter 3 Duality in Banach Space Modern optimization theory largely centers around the interplay of a normed vector space and its corresponding dual. The notion of duality is important for the following

More information

Math 60. Rumbos Spring Solutions to Assignment #17

Math 60. Rumbos Spring Solutions to Assignment #17 Math 60. Rumbos Spring 2009 1 Solutions to Assignment #17 a b 1. Prove that if ad bc 0 then the matrix A = is invertible and c d compute A 1. a b Solution: Let A = and assume that ad bc 0. c d First consider

More information

UNIVERSITY OF DUBLIN

UNIVERSITY OF DUBLIN UNIVERSITY OF DUBLIN TRINITY COLLEGE JS & SS Mathematics SS Theoretical Physics SS TSM Mathematics Faculty of Engineering, Mathematics and Science school of mathematics Trinity Term 2015 Module MA3429

More information

Econ Lecture 14. Outline

Econ Lecture 14. Outline Econ 204 2010 Lecture 14 Outline 1. Differential Equations and Solutions 2. Existence and Uniqueness of Solutions 3. Autonomous Differential Equations 4. Complex Exponentials 5. Linear Differential Equations

More information

In this chapter we study several functions that are useful in calculus and other areas of mathematics.

In this chapter we study several functions that are useful in calculus and other areas of mathematics. Calculus 5 7 Special functions In this chapter we study several functions that are useful in calculus and other areas of mathematics. 7. Hyperbolic trigonometric functions The functions we study in this

More information

SPECTRAL PROPERTIES AND NODAL SOLUTIONS FOR SECOND-ORDER, m-point, BOUNDARY VALUE PROBLEMS

SPECTRAL PROPERTIES AND NODAL SOLUTIONS FOR SECOND-ORDER, m-point, BOUNDARY VALUE PROBLEMS SPECTRAL PROPERTIES AND NODAL SOLUTIONS FOR SECOND-ORDER, m-point, BOUNDARY VALUE PROBLEMS BRYAN P. RYNNE Abstract. We consider the m-point boundary value problem consisting of the equation u = f(u), on

More information

Determine whether the following system has a trivial solution or non-trivial solution:

Determine whether the following system has a trivial solution or non-trivial solution: Practice Questions Lecture # 7 and 8 Question # Determine whether the following system has a trivial solution or non-trivial solution: x x + x x x x x The coefficient matrix is / R, R R R+ R The corresponding

More information

Applied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.

Applied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm. Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the

More information

Solutions Manual for Homework Sets Math 401. Dr Vignon S. Oussa

Solutions Manual for Homework Sets Math 401. Dr Vignon S. Oussa 1 Solutions Manual for Homework Sets Math 401 Dr Vignon S. Oussa Solutions Homework Set 0 Math 401 Fall 2015 1. (Direct Proof) Assume that x and y are odd integers. Then there exist integers u and v such

More information

Math 592 Spring Applying the identity sin(u) + sin(v) = 2 sin ( ) ( sin(f 1 2πt + ϕ) + sin(f 2 2πt + µ) = f1 + f 2 2 sin 2.

Math 592 Spring Applying the identity sin(u) + sin(v) = 2 sin ( ) ( sin(f 1 2πt + ϕ) + sin(f 2 2πt + µ) = f1 + f 2 2 sin 2. Math 59 Spring 015 Homework Drew Armstrong 1. Beats Again. Show that the phenomenon of beats is independent of phase shifts. [Hint: Consider the superposition sin(f 1 πt + ϕ) + sin(f πt + µ).] Applying

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

Monday, 6 th October 2008

Monday, 6 th October 2008 MA211 Lecture 9: 2nd order differential eqns Monday, 6 th October 2008 MA211 Lecture 9: 2nd order differential eqns 1/19 Class test next week... MA211 Lecture 9: 2nd order differential eqns 2/19 This morning

More information

SMSTC (2007/08) Probability.

SMSTC (2007/08) Probability. SMSTC (27/8) Probability www.smstc.ac.uk Contents 12 Markov chains in continuous time 12 1 12.1 Markov property and the Kolmogorov equations.................... 12 2 12.1.1 Finite state space.................................

More information

Periodic solutions of weakly coupled superlinear systems

Periodic solutions of weakly coupled superlinear systems Periodic solutions of weakly coupled superlinear systems Alessandro Fonda and Andrea Sfecci Abstract By the use of a higher dimensional version of the Poincaré Birkhoff theorem, we are able to generalize

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

Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:

Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant

More information

Note that in the example in Lecture 1, the state Home is recurrent (and even absorbing), but all other states are transient. f ii (n) f ii = n=1 < +

Note that in the example in Lecture 1, the state Home is recurrent (and even absorbing), but all other states are transient. f ii (n) f ii = n=1 < + Random Walks: WEEK 2 Recurrence and transience Consider the event {X n = i for some n > 0} by which we mean {X = i}or{x 2 = i,x i}or{x 3 = i,x 2 i,x i},. Definition.. A state i S is recurrent if P(X n

More information

4. Functions of one variable

4. Functions of one variable 4. Functions of one variable These lecture notes present my interpretation of Ruth Lawrence s lecture notes (in Hebrew) 1 In this chapter we are going to meet one of the most important concepts in mathematics:

More information

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

Georgia Tech PHYS 6124 Mathematical Methods of Physics I Georgia Tech PHYS 612 Mathematical Methods of Physics I Instructor: Predrag Cvitanović Fall semester 2012 Homework Set #5 due October 2, 2012 == show all your work for maximum credit, == put labels, title,

More information

Second-Order Linear ODEs

Second-Order Linear ODEs Chap. 2 Second-Order Linear ODEs Sec. 2.1 Homogeneous Linear ODEs of Second Order On pp. 45-46 we extend concepts defined in Chap. 1, notably solution and homogeneous and nonhomogeneous, to second-order

More information

EE/AA 578, Univ of Washington, Fall Duality

EE/AA 578, Univ of Washington, Fall Duality 7. Duality EE/AA 578, Univ of Washington, Fall 2016 Lagrange dual problem weak and strong duality geometric interpretation optimality conditions perturbation and sensitivity analysis examples generalized

More information

Name: Answer Key David Arnold. Math 50B Integral Calculus May 13, Final Exam

Name: Answer Key David Arnold. Math 50B Integral Calculus May 13, Final Exam Math 5B Integral Calculus May 3, 7 Final Exam Name: Answer Key David Arnold Instructions. (9 points) Follow the directions exactly! Whatever you are asked to do, you must do to receive full credit for

More information

Solving a linear equation in a set of integers II

Solving a linear equation in a set of integers II ACTA ARITHMETICA LXXII.4 (1995) Solving a linear equation in a set of integers II by Imre Z. Ruzsa (Budapest) 1. Introduction. We continue the study of linear equations started in Part I of this paper.

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

Linear Algebra 1 Exam 2 Solutions 7/14/3

Linear Algebra 1 Exam 2 Solutions 7/14/3 Linear Algebra 1 Exam Solutions 7/14/3 Question 1 The line L has the symmetric equation: x 1 = y + 3 The line M has the parametric equation: = z 4. [x, y, z] = [ 4, 10, 5] + s[10, 7, ]. The line N is perpendicular

More information

BIHARMONIC WAVE MAPS INTO SPHERES

BIHARMONIC WAVE MAPS INTO SPHERES BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.

More information

CS100: DISCRETE STRUCTURES. Lecture 3 Matrices Ch 3 Pages:

CS100: DISCRETE STRUCTURES. Lecture 3 Matrices Ch 3 Pages: CS100: DISCRETE STRUCTURES Lecture 3 Matrices Ch 3 Pages: 246-262 Matrices 2 Introduction DEFINITION 1: A matrix is a rectangular array of numbers. A matrix with m rows and n columns is called an m x n

More information

h(x) lim H(x) = lim Since h is nondecreasing then h(x) 0 for all x, and if h is discontinuous at a point x then H(x) > 0. Denote

h(x) lim H(x) = lim Since h is nondecreasing then h(x) 0 for all x, and if h is discontinuous at a point x then H(x) > 0. Denote Real Variables, Fall 4 Problem set 4 Solution suggestions Exercise. Let f be of bounded variation on [a, b]. Show that for each c (a, b), lim x c f(x) and lim x c f(x) exist. Prove that a monotone function

More information

Math113: Linear Algebra. Beifang Chen

Math113: Linear Algebra. Beifang Chen Math3: Linear Algebra Beifang Chen Spring 26 Contents Systems of Linear Equations 3 Systems of Linear Equations 3 Linear Systems 3 2 Geometric Interpretation 3 3 Matrices of Linear Systems 4 4 Elementary

More information

Lecture 7. Econ August 18

Lecture 7. Econ August 18 Lecture 7 Econ 2001 2015 August 18 Lecture 7 Outline First, the theorem of the maximum, an amazing result about continuity in optimization problems. Then, we start linear algebra, mostly looking at familiar

More information

Applied/Numerical Analysis Qualifying Exam

Applied/Numerical Analysis Qualifying Exam Applied/Numerical Analysis Qualifying Exam Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless, the exam

More information

EXERCISE SET 5.1. = (kx + kx + k, ky + ky + k ) = (kx + kx + 1, ky + ky + 1) = ((k + )x + 1, (k + )y + 1)

EXERCISE SET 5.1. = (kx + kx + k, ky + ky + k ) = (kx + kx + 1, ky + ky + 1) = ((k + )x + 1, (k + )y + 1) EXERCISE SET 5. 6. The pair (, 2) is in the set but the pair ( )(, 2) = (, 2) is not because the first component is negative; hence Axiom 6 fails. Axiom 5 also fails. 8. Axioms, 2, 3, 6, 9, and are easily

More information

Chapter 9. Analytic Continuation. 9.1 Analytic Continuation. For every complex problem, there is a solution that is simple, neat, and wrong.

Chapter 9. Analytic Continuation. 9.1 Analytic Continuation. For every complex problem, there is a solution that is simple, neat, and wrong. Chapter 9 Analytic Continuation For every complex problem, there is a solution that is simple, neat, and wrong. - H. L. Mencken 9.1 Analytic Continuation Suppose there is a function, f 1 (z) that is analytic

More information

Part IB Geometry. Theorems. Based on lectures by A. G. Kovalev Notes taken by Dexter Chua. Lent 2016

Part IB Geometry. Theorems. Based on lectures by A. G. Kovalev Notes taken by Dexter Chua. Lent 2016 Part IB Geometry Theorems Based on lectures by A. G. Kovalev Notes taken by Dexter Chua Lent 2016 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after lectures.

More information

NORMS ON SPACE OF MATRICES

NORMS ON SPACE OF MATRICES NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system

More information

13 Lecture 13 L Hospital s Rule and Taylor s Theorem

13 Lecture 13 L Hospital s Rule and Taylor s Theorem 3 Lecture 3 L Hospital s Rule and Taylor s Theorem 3 L Hospital s Rule Theorem 3 L Hospital s Rule) Let a R or a = and functions f,g : a,b) R satisfy the following properties ) f,g are differentiable on

More information

Some lecture notes for Math 6050E: PDEs, Fall 2016

Some lecture notes for Math 6050E: PDEs, Fall 2016 Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.

More information

Chapter 2 Convex Analysis

Chapter 2 Convex Analysis Chapter 2 Convex Analysis The theory of nonsmooth analysis is based on convex analysis. Thus, we start this chapter by giving basic concepts and results of convexity (for further readings see also [202,

More information

3.4 Elementary Matrices and Matrix Inverse

3.4 Elementary Matrices and Matrix Inverse Math 220: Summer 2015 3.4 Elementary Matrices and Matrix Inverse A n n elementary matrix is a matrix which is obtained from the n n identity matrix I n n by a single elementary row operation. Elementary

More information

Brownian Motion and Conditional Probability

Brownian Motion and Conditional Probability Math 561: Theory of Probability (Spring 2018) Week 10 Brownian Motion and Conditional Probability 10.1 Standard Brownian Motion (SBM) Brownian motion is a stochastic process with both practical and theoretical

More information

Econ 204 Differential Equations. 1 Existence and Uniqueness of Solutions

Econ 204 Differential Equations. 1 Existence and Uniqueness of Solutions Econ 4 Differential Equations In this supplement, we use the methods we have developed so far to study differential equations. 1 Existence and Uniqueness of Solutions Definition 1 A differential equation

More information

SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS

SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties

More information

Introduction to Semidefinite Programming I: Basic properties a

Introduction to Semidefinite Programming I: Basic properties a Introduction to Semidefinite Programming I: Basic properties and variations on the Goemans-Williamson approximation algorithm for max-cut MFO seminar on Semidefinite Programming May 30, 2010 Semidefinite

More information

Chapter Two Elements of Linear Algebra

Chapter Two Elements of Linear Algebra Chapter Two Elements of Linear Algebra Previously, in chapter one, we have considered single first order differential equations involving a single unknown function. In the next chapter we will begin to

More information

7 Planar systems of linear ODE

7 Planar systems of linear ODE 7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,

More information

arxiv: v2 [math.pr] 4 Sep 2017

arxiv: v2 [math.pr] 4 Sep 2017 arxiv:1708.08576v2 [math.pr] 4 Sep 2017 On the Speed of an Excited Asymmetric Random Walk Mike Cinkoske, Joe Jackson, Claire Plunkett September 5, 2017 Abstract An excited random walk is a non-markovian

More information

Solutions Chapter 9. u. (c) u(t) = 1 e t + c 2 e 3 t! c 1 e t 3c 2 e 3 t. (v) (a) u(t) = c 1 e t cos 3t + c 2 e t sin 3t. (b) du

Solutions Chapter 9. u. (c) u(t) = 1 e t + c 2 e 3 t! c 1 e t 3c 2 e 3 t. (v) (a) u(t) = c 1 e t cos 3t + c 2 e t sin 3t. (b) du Solutions hapter 9 dode 9 asic Solution Techniques 9 hoose one or more of the following differential equations, and then: (a) Solve the equation directly (b) Write down its phase plane equivalent, and

More information

Strauss PDEs 2e: Section Exercise 4 Page 1 of 6

Strauss PDEs 2e: Section Exercise 4 Page 1 of 6 Strauss PDEs 2e: Section 5.3 - Exercise 4 Page of 6 Exercise 4 Consider the problem u t = ku xx for < x < l, with the boundary conditions u(, t) = U, u x (l, t) =, and the initial condition u(x, ) =, where

More information

Linear Algebra (wi1403lr) Lecture no.4

Linear Algebra (wi1403lr) Lecture no.4 Linear Algebra (wi1403lr) Lecture no.4 EWI / DIAM / Numerical Analysis group Matthias Möller 29/04/2014 M. Möller (EWI/NA group) LA (wi1403lr) 29/04/2014 1 / 28 Review of lecture no.3 1.5 Solution Sets

More information

4 Uniform convergence

4 Uniform convergence 4 Uniform convergence In the last few sections we have seen several functions which have been defined via series or integrals. We now want to develop tools that will allow us to show that these functions

More information

Econ Slides from Lecture 1

Econ Slides from Lecture 1 Econ 205 Sobel Econ 205 - Slides from Lecture 1 Joel Sobel August 23, 2010 Warning I can t start without assuming that something is common knowledge. You can find basic definitions of Sets and Set Operations

More information

Linear Algebra, 4th day, Thursday 7/1/04 REU Info:

Linear Algebra, 4th day, Thursday 7/1/04 REU Info: Linear Algebra, 4th day, Thursday 7/1/04 REU 004. Info http//people.cs.uchicago.edu/laci/reu04. Instructor Laszlo Babai Scribe Nick Gurski 1 Linear maps We shall study the notion of maps between vector

More information

= ϕ r cos θ. 0 cos ξ sin ξ and sin ξ cos ξ. sin ξ 0 cos ξ

= ϕ r cos θ. 0 cos ξ sin ξ and sin ξ cos ξ. sin ξ 0 cos ξ 8. The Banach-Tarski paradox May, 2012 The Banach-Tarski paradox is that a unit ball in Euclidean -space can be decomposed into finitely many parts which can then be reassembled to form two unit balls

More information

Topic 4 Notes Jeremy Orloff

Topic 4 Notes Jeremy Orloff Topic 4 Notes Jeremy Orloff 4 auchy s integral formula 4. Introduction auchy s theorem is a big theorem which we will use almost daily from here on out. Right away it will reveal a number of interesting

More information

Geometry and Motion, MA 134 Week 1

Geometry and Motion, MA 134 Week 1 Geometry and Motion, MA 134 Week 1 Mario J. Micallef Spring, 2007 Warning. These handouts are not intended to be complete lecture notes. They should be supplemented by your own notes and, importantly,

More information

Equations in Quadratic Form

Equations in Quadratic Form Equations in Quadratic Form MATH 101 College Algebra J. Robert Buchanan Department of Mathematics Summer 2012 Objectives In this lesson we will learn to: make substitutions that allow equations to be written

More information

Complex Trigonometric and Hyperbolic Functions (7A)

Complex Trigonometric and Hyperbolic Functions (7A) Complex Trigonometric and Hyperbolic Functions (7A) 07/08/015 Copyright (c) 011-015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation

More information

Optimization Theory. A Concise Introduction. Jiongmin Yong

Optimization Theory. A Concise Introduction. Jiongmin Yong October 11, 017 16:5 ws-book9x6 Book Title Optimization Theory 017-08-Lecture Notes page 1 1 Optimization Theory A Concise Introduction Jiongmin Yong Optimization Theory 017-08-Lecture Notes page Optimization

More information

Step 1. Step 2. Step 4. The corrected trial solution y with evaluated coefficients d 1, d 2,..., d k becomes the particular solution y p.

Step 1. Step 2. Step 4. The corrected trial solution y with evaluated coefficients d 1, d 2,..., d k becomes the particular solution y p. Definition Atoms A and B are related if and only if their successive derivatives share a common atom. Then x 3 is related to x and x 101, while x is unrelated to e x, xe x and x sin x. Atoms x sin x and

More information

Economics 204 Summer/Fall 2011 Lecture 5 Friday July 29, 2011

Economics 204 Summer/Fall 2011 Lecture 5 Friday July 29, 2011 Economics 204 Summer/Fall 2011 Lecture 5 Friday July 29, 2011 Section 2.6 (cont.) Properties of Real Functions Here we first study properties of functions from R to R, making use of the additional structure

More information

Algebra Workshops 10 and 11

Algebra Workshops 10 and 11 Algebra Workshops 1 and 11 Suggestion: For Workshop 1 please do questions 2,3 and 14. For the other questions, it s best to wait till the material is covered in lectures. Bilinear and Quadratic Forms on

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

UNDETERMINED COEFFICIENTS SUPERPOSITION APPROACH *

UNDETERMINED COEFFICIENTS SUPERPOSITION APPROACH * 4.4 UNDETERMINED COEFFICIENTS SUPERPOSITION APPROACH 19 Discussion Problems 59. Two roots of a cubic auxiliary equation with real coeffi cients are m 1 1 and m i. What is the corresponding homogeneous

More information

Partial differential equations (ACM30220)

Partial differential equations (ACM30220) (ACM3. A pot on a stove has a handle of length that can be modelled as a rod with diffusion constant D. The equation for the temperature in the rod is u t Du xx < x

More information

Lecture 3. Econ August 12

Lecture 3. Econ August 12 Lecture 3 Econ 2001 2015 August 12 Lecture 3 Outline 1 Metric and Metric Spaces 2 Norm and Normed Spaces 3 Sequences and Subsequences 4 Convergence 5 Monotone and Bounded Sequences Announcements: - Friday

More information

LECTURE 2 LINEAR REGRESSION MODEL AND OLS

LECTURE 2 LINEAR REGRESSION MODEL AND OLS SEPTEMBER 29, 2014 LECTURE 2 LINEAR REGRESSION MODEL AND OLS Definitions A common question in econometrics is to study the effect of one group of variables X i, usually called the regressors, on another

More information

Math 331 Homework Assignment Chapter 7 Page 1 of 9

Math 331 Homework Assignment Chapter 7 Page 1 of 9 Math Homework Assignment Chapter 7 Page of 9 Instructions: Please make sure to demonstrate every step in your calculations. Return your answers including this homework sheet back to the instructor as a

More information

Mathematics 104 Fall Term 2006 Solutions to Final Exam. sin(ln t) dt = e x sin(x) dx.

Mathematics 104 Fall Term 2006 Solutions to Final Exam. sin(ln t) dt = e x sin(x) dx. Mathematics 14 Fall Term 26 Solutions to Final Exam 1. Evaluate sin(ln t) dt. Solution. We first make the substitution t = e x, for which dt = e x. This gives sin(ln t) dt = e x sin(x). To evaluate the

More information

Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R

Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract We consider semilinear

More information

FIRST YEAR CALCULUS W W L CHEN

FIRST YEAR CALCULUS W W L CHEN FIRST YER CLCULUS W W L CHEN c W W L Chen, 994, 28. This chapter is available free to all individuals, on the understanding that it is not to be used for financial gain, and may be downloaded and/or photocopied,

More information

Iterative methods for positive definite linear systems with a complex shift

Iterative methods for positive definite linear systems with a complex shift Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution

More information

We denote the space of distributions on Ω by D ( Ω) 2.

We denote the space of distributions on Ω by D ( Ω) 2. Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study

More information

CAAM 454/554: Stationary Iterative Methods

CAAM 454/554: Stationary Iterative Methods CAAM 454/554: Stationary Iterative Methods Yin Zhang (draft) CAAM, Rice University, Houston, TX 77005 2007, Revised 2010 Abstract Stationary iterative methods for solving systems of linear equations are

More information

On Chern-Simons-Schrödinger equations including a vortex point

On Chern-Simons-Schrödinger equations including a vortex point On Chern-Simons-Schrödinger equations including a vortex point Alessio Pomponio Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Workshop in Nonlinear PDEs Brussels, September 7

More information

1 Review of last lecture and introduction

1 Review of last lecture and introduction Semidefinite Programming Lecture 10 OR 637 Spring 2008 April 16, 2008 (Wednesday) Instructor: Michael Jeremy Todd Scribe: Yogeshwer (Yogi) Sharma 1 Review of last lecture and introduction Let us first

More information

Functional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...

Functional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability... Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................

More information

Subsequences and the Bolzano-Weierstrass Theorem

Subsequences and the Bolzano-Weierstrass Theorem Subsequences and the Bolzano-Weierstrass Theorem A subsequence of a sequence x n ) n N is a particular sequence whose terms are selected among those of the mother sequence x n ) n N The study of subsequences

More information