Calculus of Variations Summer Term 2014
|
|
- Grant Nash
- 6 years ago
- Views:
Transcription
1 Calculus of Variations Summer Term 2014 Lecture Juni 2014 c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
2 Purpose of Lesson Purpose of Lesson: To discuss numerical solutions of the variational problems To introduce Euler s Finite Difference Method and Ritz s Method. c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
3 9. Numerical Solutions c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
4 Numerical Solutions: The Euler-Lagrange equations may be hard to solve. Natural response is to find numerical methods. 1 Numerical solution of the Euler-Lagrange equations We won t consider these here (see other courses) 2 Euler s finite difference method 3 Ritz (Rayleigh-Ritz) In 2D: Kantorovich s method c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
5
6 Numerical approximation of integrals: Numerical Approximation use an arbitrary set of mesh points a = x 0 < x 1 < x 2 < < x n = b. approximate y (x i ) = y i+1 y i x i+1 x i = y i x i rectangle rule J[y] = b a n 1 F(x, y, y )dx F i=0 Ĵ[ ] is a function of the vector y = (y 1, y 2,..., y n ). ( x i, y i, y ) i x i = x Ĵ[y] i c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
7 Finite Difference Method (FDM) Treat this as a maximization of a function of n variables, so that we require J y i = 0 for all i = 1, 2,..., n. Typically use uniform grid so x i = x = b a n. c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
8 Simple Example Example 12.1 Find extremals for J[y] = with y(0) = 0 and y(1) = [ 1 2 y ] 2 y 2 y dx The Euler-Lagrange equation y y = 1. c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
9 Simple Example: direct solution Example 12.1 (direct solution) E-L equation: y y = 1 Solution to homogeneous equation y y = 0 is given by e λx giving characteristic equation so λ = ±1 Particular solution y = 1. λ 2 1 = 0, Final solution is y(x) = Ae x + Be x + 1. c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
10 Simple Example: direct solution Example 12.1 (direct solution) The boundary conditions y(0) = y(1) = 0 constrain A + B = 1 Ae + Be 1 = 1 so A = 1 e e 2 1 and B = e e2 e 2 1. Then the exact solution to the extremal problem is y(x) = 1 e e 2 1 ex + e e2 e 2 1 e x + 1. c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
11 Simple Example: Euler s FDM Example 12.1 (Euler s FDM) Find extremals for Euler s FDM: J[y] = 1 0 [ 1 2 y ] 2 y 2 y dx Take the grid x i = i/n, for i = 0, 1,..., n so So end points y 0 = 0 and y n = 0 x = 1/n and y i = y i+1 y i. y i = y i / x = n(y i+1 y i ) and y 2 i = n 2 ( y 2 i 2y i y i+1 + y 2 i+1). c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
12 Simple Example: Euler s FDM Example 12.1 (Euler s FDM) Find extremals for J[y] = Its FDM approximation is 1 n 1 J[y] = F (x i, y i, y i )dx = = i=0 n 1 i=0 n 1 i=0 0 [ 1 2 y ] 2 y 2 y dx 1 ( ) ( ) 2 n2 yi 2 2y i y i+1 + yi+1 2 x + yi 2 /2 y i x 1 ( ) 2 n yi 2 2y i y i+1 + yi y 2 i /2 y i n. c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
13 Simple Example: end-conditions Example 12.1 (end-conditions) We know the end conditions y(0) = y(1) = 0, which imply that y 0 = y n = 0. Include them into the objective using Lagrange multipliers H[y] = n 1 i=0 1 ( ) 2 n yi 2 2y i y i+1 + yi y i 2 /2 y i n + λ 0 y 0 + λ n y n. c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
14 Simple Example: Euler s FDM Example 12.1 (Euler s FDM) Taking derivatives, note that y i only appears in two terms of the FDM approximation H[y] = H[y] y i = n 1 i=0 1 ( ) 2 n yi 2 2y i y i+1 + yi y i 2 /2 y i n + λ 0 y 0 + λ n y n n(y 0 y 1 ) + y λ 0 for i = 0 n n(2y i y i+1 y i 1 ) + y i /n 1/n for i = 1,..., n 1 n(y n y n 1 ) + λ n for i = n We need to set the derivatives to all be zero, so we now have n = 3 linear equations, including y 0 = y n = 0, and n + 3 variables including the two Lagrange multipliers. We can solve this system numerically using, e.g., Maple. c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
15 Simple Example: Euler s FDM Example 12.1 (Euler s FDM) Example: n = 4, solve where A = Az = b and b = (0.00, 0.25, 0.25, 0.25, 0.25, 0.00, 0.00) T c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
16 Simple Example: Euler s FDM Example 12.1 (Euler s FDM) First n + 1 terms of z give y Last two terms of z give the Lagrange multipliers λ 0 and λ n. Solving the system we get for n = 4 y 1 = y 3 = , y 2 = c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
17 Simple Example: results Example 12.1 (results) c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
18 Simple Example: results Example 12.1 (results) c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
19 Convergence of Euler s FDM Convergence of Euler s FDM n 1 Ĵ[y] = F i=0 ( x i, y i, y i x i Only two terms in the sum involve y i, so Ĵ y i = ( F y i = 1 F x y + F y i x i 1, y i 1, y i 1 x ( x i 1, y i 1, y i 1 x i ( x i, y i, y i x = F F (x i, y i, y i y ) y i i ) 1 x ( ) x and y i = y i+1 y i ) + F y ) i ( x i, y i, y ) i x ( F y i x i, y i, y i x ( x i, y i, y i x ) F y i x ) x i 1, y i 1, y i 1 x ) c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
20 Convergence of Euler s FDM Convergence of Euler s FDM Ĵ = F (x i, y i, y y i y i ( F i ) y i x i, y i, y i x ) ( F y i x x i 1, y i 1, y i 1 x ) = 0. In limit n, then x 0, and so we get F y d ( ) F dx y = 0 which are the Euler-Lagrange equations. i.e., the finite difference solution converges to the solution of the Euler-Lagrange equations. c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
21 Comments Remarks There are lots of ways to improve Euler s FDM use a better method of numerical quadrature (integration) trapezoidal rule Simpson s rule Romberg s method use a non-uniform grid make it finer where there is more variation We can use a different approach that can be even better. c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
22 Ritz s Method Ritz s Method In Ritz s method (called Kantorovich s method where there is more than one independent variable), we approximate our functions (the extremal in particular) using a family of simple functions. Again we can reduce the problem into a standard multivariable maximization problem, but now we seek coefficients for our approximation. c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
23 Ritz s Method Assume we can approximate y(x) by y(x) = φ 0 (x) + c 1 φ 1 (x) + c 2 φ 2 (x) + + c n φ n (x) where we choose a convenient set of functions φ j (x) and find the values of c j which produce an extremal. For fixed end-points problem: Choose φ 0 (x) to satisfy the end conditions. Then φ j (x 0 ) = φ j (x 1 ) = 0 for j = 1, 2,..., n The φ can be chosen from standard sets of functions, e.g. power series, trigonometric functions, Bessel s functions, etc. (but must be linearly independent). c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
24 Ritz s Method Select { φ j } n j=0 Approximate y n (x) = φ 0 (x) + c 1 φ 1 (x) + c 2 φ 2 (x) + + c n φ n (x) Approximate J[y] J[y n ] = F(x, y n, y n)dx. x 0 x 1 Integrate to get J[y n ] = J n (c 1, c 2,..., c n ). J n is a known function of n variables, so we can maximize (or minimize) it as usual by J n c i = 0 for all i = 1, 2,..., n. c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
25 Upper Bounds Assume the extremal of interest is a minimum, then for the extremal J[y] < J[ŷ] for all ŷ within the neighborhood of y. Assume our approximating function y n is close enough to be in that neighborhood, then J[y] J[y n ] = J n [c] so the approximation provides an upper bound on the minimum J[y]. Another way to think about it is that we optimize on a smaller set of possible functions y, so we can t get quite as good a minimum. c Daria Apushkinskaya 2014 () Calculus of variations lecture Juni / 25
Calculus of Variations Summer Term 2016
Calculus of Variations Summer Term 2016 Lecture 12 Universität des Saarlandes 17. Juni 2016 c Daria Apushkinskaya (UdS) Calculus of variations lecture 12 17. Juni 2016 1 / 32 Purpose of Lesson Purpose
More informationVariational Methods & Optimal Control
Variational Methods & Optimal Control lecture 12 Matthew Roughan Discipline of Applied Mathematics School of Mathematical Sciences University of Adelaide April 14, 216
More informationCalculus of Variations Summer Term 2015
Calculus of Variations Summer Term 2015 Lecture 12 Universität des Saarlandes 17. Juni 2015 c Daria Apushkinskaya (UdS) Calculus of variations lecture 12 17. Juni 2015 1 / 31 Purpose of Lesson Purpose
More informationCalculus of Variations Summer Term 2014
Calculus of Variations Summer Term 2014 Lecture 5 7. Mai 2014 c Daria Apushkinskaya 2014 () Calculus of variations lecture 5 7. Mai 2014 1 / 25 Purpose of Lesson Purpose of Lesson: To discuss catenary
More informationCalculus of Variations Summer Term 2014
Calculus of Variations Summer Term 2014 Lecture 2 25. April 2014 c Daria Apushkinskaya 2014 () Calculus of variations lecture 2 25. April 2014 1 / 20 Purpose of Lesson Purpose of Lesson: To discuss the
More informationCalculus of Variations Summer Term 2014
Calculus of Variations Summer Term 2014 Lecture 9 23. Mai 2014 c Daria Apushkinskaya 2014 () Calculus of variations lecture 9 23. Mai 2014 1 / 23 Purpose of Lesson Purpose of Lesson: To consider several
More informationCalculus of Variations Summer Term 2014
Calculus of Variations Summer Term 2014 Lecture 1 28. April 2014 c Daria Apushkinskaya 2014 () Calculus of variations lecture 1 28. April 2014 1 / 22 Organisatorial Issues General Information Lecturer:
More informationCalculus of Variations Summer Term 2015
Calculus of Variations Summer Term 2015 Lecture 14 Universität des Saarlandes 24. Juni 2015 c Daria Apushkinskaya (UdS) Calculus of variations lecture 14 24. Juni 2015 1 / 20 Purpose of Lesson Purpose
More informationCalculus of Variations Summer Term 2016
Calculus of Variations Summer Term 2016 Lecture 14 Universität des Saarlandes 28. Juni 2016 c Daria Apushkinskaya (UdS) Calculus of variations lecture 14 28. Juni 2016 1 / 31 Purpose of Lesson Purpose
More informationVariational Methods and Optimal Control
Variational Methods and Optimal Control A/Prof. Matthew Roughan July 26, 2010 lecture01 Introduction What is the point of this course? Revision Example 1: The money pit. Example 2: Catenary: shape of a
More informationStat 451 Lecture Notes Numerical Integration
Stat 451 Lecture Notes 03 12 Numerical Integration Ryan Martin UIC www.math.uic.edu/~rgmartin 1 Based on Chapter 5 in Givens & Hoeting, and Chapters 4 & 18 of Lange 2 Updated: February 11, 2016 1 / 29
More informationMath 10C - Fall Final Exam
Math 1C - Fall 217 - Final Exam Problem 1. Consider the function f(x, y) = 1 x 2 (y 1) 2. (i) Draw the level curve through the point P (1, 2). Find the gradient of f at the point P and draw the gradient
More informationNumerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error. 2- Fixed point
Numerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error In this method we assume initial value of x, and substitute in the equation. Then modify x and continue till we
More informationDay 2 Notes: Riemann Sums In calculus, the result of f ( x)
AP Calculus Unit 6 Basic Integration & Applications Day 2 Notes: Riemann Sums In calculus, the result of f ( x) dx is a function that represents the anti-derivative of the function f(x). This is also sometimes
More informationSufficient Conditions for Finite-variable Constrained Minimization
Lecture 4 It is a small de tour but it is important to understand this before we move to calculus of variations. Sufficient Conditions for Finite-variable Constrained Minimization ME 256, Indian Institute
More informationPhysics 342 Lecture 2. Linear Algebra I. Lecture 2. Physics 342 Quantum Mechanics I
Physics 342 Lecture 2 Linear Algebra I Lecture 2 Physics 342 Quantum Mechanics I Wednesday, January 27th, 21 From separation of variables, we move to linear algebra Roughly speaking, this is the study
More informationFirst and Second Order ODEs
Civil Engineering 2 Mathematics Autumn 211 M. Ottobre First and Second Order ODEs Warning: all the handouts that I will provide during the course are in no way exhaustive, they are just short recaps. Notation
More informationMechanical Systems II. Method of Lagrange Multipliers
Mechanical Systems II. Method of Lagrange Multipliers Rafael Wisniewski Aalborg University Abstract. So far our approach to classical mechanics was limited to finding a critical point of a certain functional.
More informationPhysics 342 Lecture 2. Linear Algebra I. Lecture 2. Physics 342 Quantum Mechanics I
Physics 342 Lecture 2 Linear Algebra I Lecture 2 Physics 342 Quantum Mechanics I Wednesday, January 3th, 28 From separation of variables, we move to linear algebra Roughly speaking, this is the study of
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 information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationPractice Exam 2 (Solutions)
Math 5, Fall 7 Practice Exam (Solutions). Using the points f(x), f(x h) and f(x + h) we want to derive an approximation f (x) a f(x) + a f(x h) + a f(x + h). (a) Write the linear system that you must solve
More information12.0 Properties of orthogonal polynomials
12.0 Properties of orthogonal polynomials In this section we study orthogonal polynomials to use them for the construction of quadrature formulas investigate projections on polynomial spaces and their
More informationFinite-Elements Method 2
Finite-Elements Method 2 January 29, 2014 2 From Applied Numerical Analysis Gerald-Wheatley (2004), Chapter 9. Finite-Elements Method 3 Introduction Finite-element methods (FEM) are based on some mathematical
More informationLESSON 25: LAGRANGE MULTIPLIERS OCTOBER 30, 2017
LESSON 5: LAGRANGE MULTIPLIERS OCTOBER 30, 017 Lagrange multipliers is another method of finding minima and maxima of functions of more than one variable. In fact, many of the problems from the last homework
More information6 Lecture 6b: the Euler Maclaurin formula
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Fall 217 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 217 March 26, 218 6 Lecture 6b: the Euler Maclaurin formula
More informationMultivariable Calculus and Matrix Algebra-Summer 2017
Multivariable Calculus and Matrix Algebra-Summer 017 Homework 4 Solutions Note that the solutions below are for the latest version of the problems posted. For those of you who worked on an earlier version
More informationCalculus of Variations
ECE 68 Midterm Exam Solution April 1, 8 1 Calculus of Variations This exam is open book and open notes You may consult additional references You may even discuss the problems (with anyone), but you must
More informationCalculus of Variations Summer Term 2014
Calculus of Variations Summer Term 2014 Lecture 20 18. Juli 2014 c Daria Apushkinskaya 2014 () Calculus of variations lecture 20 18. Juli 2014 1 / 20 Purpose of Lesson Purpose of Lesson: To discuss several
More informationOptimization using Calculus. Optimization of Functions of Multiple Variables subject to Equality Constraints
Optimization using Calculus Optimization of Functions of Multiple Variables subject to Equality Constraints 1 Objectives Optimization of functions of multiple variables subjected to equality constraints
More informationIntroduction to the Calculus of Variations
236861 Numerical Geometry of Images Tutorial 1 Introduction to the Calculus of Variations Alex Bronstein c 2005 1 Calculus Calculus of variations 1. Function Functional f : R n R Example: f(x, y) =x 2
More informationMATH 280 Multivariate Calculus Fall Integrating a vector field over a surface
MATH 280 Multivariate Calculus Fall 2011 Definition Integrating a vector field over a surface We are given a vector field F in space and an oriented surface in the domain of F as shown in the figure below
More informationWeek 4: Calculus and Optimization (Jehle and Reny, Chapter A2)
Week 4: Calculus and Optimization (Jehle and Reny, Chapter A2) Tsun-Feng Chiang *School of Economics, Henan University, Kaifeng, China September 27, 2015 Microeconomic Theory Week 4: Calculus and Optimization
More informationHigher-order ordinary differential equations
Higher-order ordinary differential equations 1 A linear ODE of general order n has the form a n (x) dn y dx n +a n 1(x) dn 1 y dx n 1 + +a 1(x) dy dx +a 0(x)y = f(x). If f(x) = 0 then the equation is called
More information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More informationGeneralization to inequality constrained problem. Maximize
Lecture 11. 26 September 2006 Review of Lecture #10: Second order optimality conditions necessary condition, sufficient condition. If the necessary condition is violated the point cannot be a local minimum
More informationHow to Solve Linear Differential Equations
How to Solve Linear Differential Equations Definition: Euler Base Atom, Euler Solution Atom Independence of Atoms Construction of the General Solution from a List of Distinct Atoms Euler s Theorems Euler
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 1. Runge-Kutta Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Mathematics National Institute of Technology Durgapur Durgapur-71309 email: anita.buie@gmail.com 1 . Chapter 8 Numerical Solution of
More informationMA22S3 Summary Sheet: Ordinary Differential Equations
MA22S3 Summary Sheet: Ordinary Differential Equations December 14, 2017 Kreyszig s textbook is a suitable guide for this part of the module. Contents 1 Terminology 1 2 First order separable 2 2.1 Separable
More informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
More informationChapter 2: First Order DE 2.6 Exact DE and Integrating Fa
Chapter 2: First Order DE 2.6 Exact DE and Integrating Factor First Order DE Recall the general form of the First Order DEs (FODE): dy dx = f(x, y) (1) (In this section x is the independent variable; not
More informationNumerical Optimization Algorithms
Numerical Optimization Algorithms 1. Overview. Calculus of Variations 3. Linearized Supersonic Flow 4. Steepest Descent 5. Smoothed Steepest Descent Overview 1 Two Main Categories of Optimization Algorithms
More information1 Lesson 13: Methods of Integration
Lesson 3: Methods of Integration Chapter 6 Material: pages 273-294 in the textbook: Lesson 3 reviews integration by parts and presents integration via partial fraction decomposition as the third of the
More informationCalculus: Area. Mathematics 15: Lecture 22. Dan Sloughter. Furman University. November 12, 2006
Calculus: Area Mathematics 15: Lecture 22 Dan Sloughter Furman University November 12, 2006 Dan Sloughter (Furman University) Calculus: Area November 12, 2006 1 / 7 Area Note: formulas for the areas of
More information. (a) Express [ ] as a non-trivial linear combination of u = [ ], v = [ ] and w =[ ], if possible. Otherwise, give your comments. (b) Express +8x+9x a
TE Linear Algebra and Numerical Methods Tutorial Set : Two Hours. (a) Show that the product AA T is a symmetric matrix. (b) Show that any square matrix A can be written as the sum of a symmetric matrix
More informationSimple ODE Solvers - Derivation
Simple ODE Solvers - Derivation These notes provide derivations of some simple algorithms for generating, numerically, approximate solutions to the initial value problem y (t =f ( t, y(t y(t 0 =y 0 Here
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 informationLECTURE 14: REGULAR SINGULAR POINTS, EULER EQUATIONS
LECTURE 14: REGULAR SINGULAR POINTS, EULER EQUATIONS 1. Regular Singular Points During the past few lectures, we have been focusing on second order linear ODEs of the form y + p(x)y + q(x)y = g(x). Particularly,
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More information1.8. Integration using Tables and CAS
1.. INTEGRATION USING TABLES AND CAS 39 1.. Integration using Tables and CAS The use of tables of integrals and Computer Algebra Systems allow us to find integrals very quickly without having to perform
More informationHIGHER ORDER METHODS. There are two principal means to derive higher order methods. b j f(x n j,y n j )
HIGHER ORDER METHODS There are two principal means to derive higher order methods y n+1 = p j=0 a j y n j + h p j= 1 b j f(x n j,y n j ) (a) Method of Undetermined Coefficients (b) Numerical Integration
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationFirst Order Differential Equations
Chapter 2 First Order Differential Equations 2.1 9 10 CHAPTER 2. FIRST ORDER DIFFERENTIAL EQUATIONS 2.2 Separable Equations A first order differential equation = f(x, y) is called separable if f(x, y)
More information12d. Regular Singular Points
October 22, 2012 12d-1 12d. Regular Singular Points We have studied solutions to the linear second order differential equations of the form P (x)y + Q(x)y + R(x)y = 0 (1) in the cases with P, Q, R real
More informationMTH301 Calculus II Glossary For Final Term Exam Preparation
MTH301 Calculus II Glossary For Final Term Exam Preparation Glossary Absolute maximum : The output value of the highest point on a graph over a given input interval or over all possible input values. An
More informationMath6100 Day 8 Notes 6.1, 6.2 & 6.3, Area
Math6100 Day 8 Notes 6.1, 6.2 & 6.3, Area 6.1 Area of Polygonal Regions Let's first derive formulas for the area of these shapes. 1. Rectangle 2. Parallelogram 3. Triangle 4. Trapezoid 1 Ex 1: Find the
More informationMethod of Lagrange Multipliers
Method of Lagrange Multipliers A. Salih Department of Aerospace Engineering Indian Institute of Space Science and Technology, Thiruvananthapuram September 2013 Lagrange multiplier method is a technique
More informationJUST THE MATHS UNIT NUMBER NUMERICAL MATHEMATICS 6 (Numerical solution) of (ordinary differential equations (A)) A.J.Hobson
JUST THE MATHS UNIT NUMBER 17.6 NUMERICAL MATHEMATICS 6 (Numerical solution) of (ordinary differential equations (A)) by A.J.Hobson 17.6.1 Euler s unmodified method 17.6.2 Euler s modified method 17.6.3
More informationReview. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn
Review Taylor Series and Error Analysis Roots of Equations Linear Algebraic Equations Optimization Numerical Differentiation and Integration Ordinary Differential Equations Partial Differential Equations
More informationEssential Ordinary Differential Equations
MODULE 1: MATHEMATICAL PRELIMINARIES 10 Lecture 2 Essential Ordinary Differential Equations In this lecture, we recall some methods of solving first-order IVP in ODE (separable and linear) and homogeneous
More informationIntegration. Topic: Trapezoidal Rule. Major: General Engineering. Author: Autar Kaw, Charlie Barker.
Integration Topic: Trapezoidal Rule Major: General Engineering Author: Autar Kaw, Charlie Barker 1 What is Integration Integration: The process of measuring the area under a function plotted on a graph.
More informationSecond-Order Homogeneous Linear Equations with Constant Coefficients
15 Second-Order Homogeneous Linear Equations with Constant Coefficients A very important class of second-order homogeneous linear equations consists of those with constant coefficients; that is, those
More informationLESSON 21: DIFFERENTIALS OF MULTIVARIABLE FUNCTIONS OCTOBER 18, 2017
LESSON 21: DIFFERENTIALS OF MULTIVARIABLE FUNCTIONS OCTOBER 18, 2017 Today we do a quick review of differentials for functions of a single variable and then discuss how to extend this notion to functions
More informationLink lecture - Lagrange Multipliers
Link lecture - Lagrange Multipliers Lagrange multipliers provide a method for finding a stationary point of a function, say f(x, y) when the variables are subject to constraints, say of the form g(x, y)
More informationReading group: Calculus of Variations and Optimal Control Theory by Daniel Liberzon
: Calculus of Variations and Optimal Control Theory by Daniel Liberzon 16th March 2017 1 / 30 Content 1 2 Recall on finite-dimensional of a global minimum 3 Infinite-dimensional 4 2 / 30 Content 1 2 Recall
More informationMcGill University Math 325A: Differential Equations LECTURE 12: SOLUTIONS FOR EQUATIONS WITH CONSTANTS COEFFICIENTS (II)
McGill University Math 325A: Differential Equations LECTURE 12: SOLUTIONS FOR EQUATIONS WITH CONSTANTS COEFFICIENTS (II) HIGHER ORDER DIFFERENTIAL EQUATIONS (IV) 1 Introduction (Text: pp. 338-367, Chap.
More informationPhysics 200 Lecture 4. Integration. Lecture 4. Physics 200 Laboratory
Physics 2 Lecture 4 Integration Lecture 4 Physics 2 Laboratory Monday, February 21st, 211 Integration is the flip-side of differentiation in fact, it is often possible to write a differential equation
More informationLESSON 21: DIFFERENTIALS OF MULTIVARIABLE FUNCTIONS MATH FALL 2018
LESSON 21: DIFFERENTIALS OF MULTIVARIABLE FUNCTIONS MATH 16020 FALL 2018 ELLEN WELD 1. Quick Review of Differentials Ex 1. Consider the function f(x) x. We know that f(9) 9 3, but what is f(9.1) 9.1? Obviously,
More informationdt 2 The Order of a differential equation is the order of the highest derivative that occurs in the equation. Example The differential equation
Lecture 18 : Direction Fields and Euler s Method A Differential Equation is an equation relating an unknown function and one or more of its derivatives. Examples Population growth : dp dp = kp, or = kp
More informationMATH 1231 MATHEMATICS 1B Calculus Section 2: - ODEs.
MATH 1231 MATHEMATICS 1B 2007. For use in Dr Chris Tisdell s lectures: Tues 11 + Thur 10 in KBT Calculus Section 2: - ODEs. 1. Motivation 2. What you should already know 3. Types and orders of ODEs 4.
More informationMath Assignment 5
Math 2280 - Assignment 5 Dylan Zwick Fall 2013 Section 3.4-1, 5, 18, 21 Section 3.5-1, 11, 23, 28, 35, 47, 56 Section 3.6-1, 2, 9, 17, 24 1 Section 3.4 - Mechanical Vibrations 3.4.1 - Determine the period
More informationEcon Slides from Lecture 14
Econ 205 Sobel Econ 205 - Slides from Lecture 14 Joel Sobel September 10, 2010 Theorem ( Lagrange Multipliers ) Theorem If x solves max f (x) subject to G(x) = 0 then there exists λ such that Df (x ) =
More information17.2 Nonhomogeneous Linear Equations. 27 September 2007
17.2 Nonhomogeneous Linear Equations 27 September 2007 Nonhomogeneous Linear Equations The differential equation to be studied is of the form ay (x) + by (x) + cy(x) = G(x) (1) where a 0, b, c are given
More informationAPPM/MATH Problem Set 6 Solutions
APPM/MATH 460 Problem Set 6 Solutions This assignment is due by 4pm on Wednesday, November 6th You may either turn it in to me in class or in the box outside my office door (ECOT ) Minimal credit will
More informationMonday, 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 information01 Harmonic Oscillations
Utah State University DigitalCommons@USU Foundations of Wave Phenomena Library Digital Monographs 8-2014 01 Harmonic Oscillations Charles G. Torre Department of Physics, Utah State University, Charles.Torre@usu.edu
More informationNumerical differentiation
Numerical differentiation Paul Seidel 1801 Lecture Notes Fall 011 Suppose that we have a function f(x) which is not given by a formula but as a result of some measurement or simulation (computer experiment)
More informationReview for Exam 2. Review for Exam 2.
Review for Exam 2. 5 or 6 problems. No multiple choice questions. No notes, no books, no calculators. Problems similar to homeworks. Exam covers: Regular-singular points (5.5). Euler differential equation
More informationLagrange Multipliers
Lagrange Multipliers (Com S 477/577 Notes) Yan-Bin Jia Nov 9, 2017 1 Introduction We turn now to the study of minimization with constraints. More specifically, we will tackle the following problem: minimize
More informationCOURSE Numerical integration of functions
COURSE 6 3. Numerical integration of functions The need: for evaluating definite integrals of functions that has no explicit antiderivatives or whose antiderivatives are not easy to obtain. Let f : [a,
More informationHomogeneous Linear ODEs of Second Order Notes for Math 331
Homogeneous Linear ODEs of Second Order Notes for Math 331 Richard S. Ellis Department of Mathematics and Statistics University of Massachusetts Amherst, MA 01003 In this document we consider homogeneous
More informationCalculus of Variations
Structural Engineering Research Group (SERG) Summer Seminar Series #9 July, 4 Calculus of Variations Tzuyang Yu Associate Professor, Ph.D. Department of Civil and Environmental Engineering The University
More informationDynamic Programming. Macro 3: Lecture 2. Mark Huggett 2. 2 Georgetown. September, 2016
Macro 3: Lecture 2 Mark Huggett 2 2 Georgetown September, 2016 Three Maps: Review of Lecture 1 X = R 1 + and X grid = {x 1,..., x n } X where x i+1 > x i 1. T (v)(x) = max u(f (x) y) + βv(y) s.t. y Γ 1
More informationUnit #24 - Lagrange Multipliers Section 15.3
Unit #24 - Lagrange Multipliers Section 1.3 Some material from Calculus, Single and MultiVariable by Hughes-Hallett, Gleason, McCallum et. al. Copyright 200 by John Wiley & Sons, Inc. This material is
More informationHOMEWORK 7 SOLUTIONS
HOMEWORK 7 SOLUTIONS MA11: ADVANCED CALCULUS, HILARY 17 (1) Using the method of Lagrange multipliers, find the largest and smallest values of the function f(x, y) xy on the ellipse x + y 1. Solution: The
More informationEC487 Advanced Microeconomics, Part I: Lecture 2
EC487 Advanced Microeconomics, Part I: Lecture 2 Leonardo Felli 32L.LG.04 6 October, 2017 Properties of the Profit Function Recall the following property of the profit function π(p, w) = max x p f (x)
More informationMath RE - Calculus II Antiderivatives and the Indefinite Integral Page 1 of 5
Math 201-203-RE - Calculus II Antiderivatives and the Indefinite Integral Page 1 of 5 What is the Antiderivative? In a derivative problem, a function f(x) is given and you find the derivative f (x) using
More informationOrdinary 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 informationOrdinary differential equations Notes for FYS3140
Ordinary differential equations Notes for FYS3140 Susanne Viefers, Dept of Physics, University of Oslo April 4, 2018 Abstract Ordinary differential equations show up in many places in physics, and these
More informationPreface. 2 Linear Equations and Eigenvalue Problem 22
Contents Preface xv 1 Errors in Computation 1 1.1 Introduction 1 1.2 Floating Point Representation of Number 1 1.3 Binary Numbers 2 1.3.1 Binary number representation in computer 3 1.4 Significant Digits
More informationCh 5.4: Regular Singular Points
Ch 5.4: Regular Singular Points! Recall that the point x 0 is an ordinary point of the equation if p(x) = Q(x)/P(x) and q(x)= R(x)/P(x) are analytic at at x 0. Otherwise x 0 is a singular point.! Thus,
More informationMathematical Economics. Lecture Notes (in extracts)
Prof. Dr. Frank Werner Faculty of Mathematics Institute of Mathematical Optimization (IMO) http://math.uni-magdeburg.de/ werner/math-ec-new.html Mathematical Economics Lecture Notes (in extracts) Winter
More informationPDE and Boundary-Value Problems Winter Term 2014/2015
PDE and Boundary-Value Problems Winter Term 214/215 Lecture 17 Saarland University 22. Januar 215 c Daria Apushkinskaya (UdS) PDE and BVP lecture 17 22. Januar 215 1 / 28 Purpose of Lesson To show how
More informationIntroduction to First Order Equations Sections
A B I L E N E C H R I S T I A N U N I V E R S I T Y Department of Mathematics Introduction to First Order Equations Sections 2.1-2.3 Dr. John Ehrke Department of Mathematics Fall 2012 Course Goals The
More informationExercises * on Linear Algebra
Exercises * on Linear Algebra Laurenz Wiskott Institut für Neuroinformatik Ruhr-Universität Bochum, Germany, EU 4 February 7 Contents Vector spaces 4. Definition...............................................
More informationLeast Squares Regression
Least Squares Regression Chemical Engineering 2450 - Numerical Methods Given N data points x i, y i, i 1 N, and a function that we wish to fit to these data points, fx, we define S as the sum of the squared
More informationPhysics 129a Calculus of Variations Frank Porter Revision
Physics 129a Calculus of Variations 71113 Frank Porter Revision 171116 1 Introduction Many problems in physics have to do with extrema. When the problem involves finding a function that satisfies some
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationIntegration Using Tables and Summary of Techniques
Integration Using Tables and Summary of Techniques Philippe B. Laval KSU Today Philippe B. Laval (KSU) Summary Today 1 / 13 Introduction We wrap up integration techniques by discussing the following topics:
More information