2D1240 Numerical Methods II / André Jaun, NADA, KTH NADA
|
|
- Patience Ferguson
- 5 years ago
- Views:
Transcription
1 ND 2D24 Numerical Methods II / ndré Jaun, ND, KTH 2. INITIL- / OUNDRY VLUE PROLEMS 2.. Remember: what we saw during the last lesson Euler and Runge-Kutta methods for solving ODEs Predator-prey model using a system of ODEs nharmonic pendulum writing 2 a order ODE as a system with two ODEs nd 2..2 Overview: what you will learn today oundary value problem Solution by shooting and finite differences 5 theory questions for the exam
2 2.2 Initial- / oundary value problems (NM 8.7, H.) Take Newton s second law of motion for example! ψ! ψ x v v, x = f(t; x; x ) = This second order problem is specified with 2 additional constraints. Initial value problems satisfy 2 initial conditions (I) at the beginning of the interval m F x(t ) = x ; v(t ) = v Knowing the initial position and velocity of a bullet, where is it going to land? oundary value problems satisfy 2 boundary conditions () at the boundaries = x a ; x(b) = x b x(a) What is the trajectory of a bullet shot a from to b? oundary conditions can be specified as the x(a) value (Dirichlet), the derivative (a) (Neuman) or a linear combination of the two x(a) + ffx (a) = fi (Robin). x
3 t=[ ]; v=2.5; options=odeset('events','on','outputsel',[2],'refine',5) >> a=pi/6; X=[. v*cos(a) v*sin(a)]; [T,Y]=ode23(shoot,t,X,options); >> a=pi/3; X=[. v*cos(a) v*sin(a)]; [T2,Y2]=ode23(shoot,t,X,options); >> x=:.:.5; plot(x,exp(x.^2), Y(:,),Y(:,2),'--',Y2(:,),Y2(:,2),':') >> 2.3 Shooting method for boundary value problems (NM 8.7., H.3) The method solves initial value problems, adjusting the I to approach the. Shoot from a with a guess velocity; if the bullet falls short, increase the velocity, etc Example determine the canon angle to shoot from the hill (; ) to the target (; =e). varargout = shoot(t,x,flag); function.4 Shooting from the top of a Gaussian hill flag switch '', case.2 varargout } = dydt(t,x); 'events', [varargout :3}]= events(t,x); case error(['unknown flag.']); otherwise, end f=dydt(t,x) function X(4); ; -9.8]; f=[x(3); f=(vx,vy,ax,ay) [value,isterminal,direction]=events(t,x) function vertical coordinate x(2) shapemountain=exp(-x()^2); ground elevation = X(2)-shapeMountain; value = ; isterminal present height stop if ground.5.5 horizontal coordinate x() direction = -; going down Here events(t,x) is used to stop the integration when the bullet hits the ground; experimentally, the canon angle is found to be somewhere ff 2 [ß=6; ß=3].
4 b=pi/3; tol=pi/8; t=[ ]; v=2.5; xtarget=; a=pi/6; options=odeset('events','on','outputsel',[2],'refine',5). v*cos(a) v*sin(a)];[t,y]=ode23(shoot,t,x,options);fa=y(end,)-xtarget; X=[. v*cos(b) v*sin(b)];[t,y]=ode23(shoot,t,x,options);fb=y(end,)-xtarget; X=[ abs(b-a)>tol while x=b-fb*(b-a)/(fb-fa); The solution is computed by solving a non-linear equation f(ff) = x end (ff) x b = : the secant method is most effective when f (ff) is not known explicitly. In Matlab use the same driver shoot() to search with the secant method: fa=fb; b=x; fb=fx; a=b; x b] sol=8/pi*[a.4 Shooting from the top of a Gaussian hill X=[. v*cos(x) v*sin(x)];[t,y]=ode23(shoot,t,x,options);fx=y(end,)-xtarget; end.2 x b a vertical coordinate x(2) ==> alpha = 48. +/-.3 degree.5.5 horizontal coordinate x() Linear boundary value problems can be solved with a system using only two trials.
5 i+ 2y + i i y p y i+ y y i i x a 2 + h ha p xi ) ( e i d i y i + ( + p xi ) z } i f b i 2.4 Finite differences for boundary value problems (NM 8.7.2, H.4) The method divides the interval where the solution is sought into equidistant pieces using finite difference (FD) to approximate derivatives with unknowns y i = y(x i ) x i = x a + i h; h = (x b x a )=N i ß y i+ y i y ; y i ß y i+ 2y i + y i h 2 2h here for example with a O(h precision ). Linear example: solve the boundary value problem in the interval 2 2 [; 4] x y + a p xy (6 y)x = ; y() = ;y(4) = for parameter values a 2 f 3; 2; ; ; ; 2g. Discretize using a uniform mesh x i = x a + i h, x a =, x b = 4, h = 4=N where N is the number of intervals. Insert the FD approximation and evaluate for every mesh point x i inside 2h This yields a tridiagonal linear y (E + D + F)y = b system (6 y i )x i = 8i = 2 :::;N i (2 h 2 x i ) z } y i+ = 6h 2 x i z } y ha 8i 2 z } 2
6 y 2 y n y n y y 2 y n y n y y(a) 2 b oundary conditions are usually imposed by replacing the corresponding equation ::: 2 d 2 f 2 ::: e e 3 d 3 f 3 ::: y 3. = b 3. ::: e n d n f n ::: n b y(b) and substituting the columns to keep the matrix symmetric if it was in the beginning ::: d 2 f 2 ::: y 3. e 3 d 3 f 3 ::: = y(a) 2 e 2 y(a) b b 3. ::: e n d n ::: n f n y(b) b y(b) The system is can be solved directly using tridia or solvers for sparse systems.
7 plot(x,y), subplot(3,2,k); ylabel(l) l=sprintf('a=i',a), In Matlab, this yields X=[; 4]; Y=[; -]; k=:6 for a=k-4; n=; h=diff(x)/(n-); x=linspace(x(),x(2),n)'; e=-.5*h*a*sqrt(x); boundaries parameter a n- intervals assemble tridiagonal d=-(2-h^2*x); matrix f=+.5*h*a*sqrt(x); d [; f(:n-)]],-:,n,n); =spdiags([[e(2:n);] assemble RHS b=6*h^2*x; (,2)=; b()=y(); (,)=.; (n,n-)=; b(n)=y(2); (n,n)=.; y= b; solve a= 3 a= a= y +ay x /2 (6 y)x=, y()=, y(4)= a= 2 a= a= end Peer Teaching (2 minutes to think, explain to your neighbour and vote) omputational cost. Examine the number of operations required to solve a linear boundary value problem in dimension. Which method is more efficient? ψ shooting " about equal # cannot say! finite differences
8 i+ 2y i + y i y 2 + cy 2 i (6 y i)x i = 8i = 2 :::;N h ( G(yi ) = (2 h 2 x i )y i + ch 2 y 2 i Non-linear example: solve the boundary value problem in the interval x 2 [; 4] + cy 2 (6 y)x = with y() = ;y(4) = ; for c = : y Discretize using finite differences on a uniform mesh b i = 6h 2 x i This non-linear system of the form F(y) = is solved after inserting the boundary conditions using the Newton method y new = y old + ffiy, withjffiy = f(y old ). The Jacobian matrix J i;j = F i =y j is tridiagonal y i + G(y i ) + y i+ b i = with (y ) ::: G G (y 2 ) :::..... J = ::: G (y n ) ::: G (y n ) with G (y i ) = (2 h 2 x i ) + 2ch 2 y i. What should we take for the initial guess?
9 randolinj, y"+c*y^2=(6-y)x; y()=;y(4)=- bandmatrix, non-linear system with Newton's method clf clear, xend= 4; x=[ xend]; x=; yend=-; y=[ ]; plot(x,y); hold on y=; + cy^2 =(6-y)x, y()=, y(4)=- ') title('y" ntot=4; h=xend/ntot; Solve linear problem with c= for an initial guess c= one=ones(n-,); b=6*h^2*x; dia=-(2-h^2*x); b(n)=b(n)-yend; b()=b()-y; y=tridia(dia,one,one,b); x; xend]; Y=[y; y; yend]; plot(x,y,'--') X=[; Method: gradually increase the non-linearity With strongly non-linear equations, first switch-off the non-linear term and calculate a linear solution for the inital guess. Then gradually switch-on the nonlinearity, using increasingly non-linear solutions to start the Newton iterations. In Matlab this can be implemented as n=ntot-; x=h*(:n)'; [ymax,ix]=max(y); text(x(ix),ymax,'c='), drawnow
10 ymax [ymax,ix]=max(y); num2str(c)]), drawnow text(x(ix),ymax,['c=' end plot(x,y,'r') Gradually increase c using Newton's method with previous solutions as initial guess 4 y" + cy 2 =(6 y)x, y()=, y(4)= cvec=[.25.5.]; nr=:3 for err=; c=cvec(nr), experimentally iter=; 3 2 c=. c=.25 c=.5 err>e-4 & iter<5 while Gprim=dia+2*c*h^2*y; G=dia.*y+c*h^2*y.^2; y(:n-)]+g+[y(2:n);yend]-6*h^2*x; F=[y; J dy=-f(y) dy=-tridia(gprim,one,one,f); err=norm(dy,inf)/norm(y,inf) iter=iter+; y=y+dy; end iter==5, disp('no convergence'),break,end if 2 c= Y=[y; y; yend]; plot(x,y,':') The non-linearity is very strong as can be seen from the plot.
2D1240 Numerical Methods II / André Jaun, NADA, KTH NADA
A NADA D4 Numerical Methods II / André Jaun, NADA, KTH 8. PRECISION INTEGRATION, NON-LINEAR EQUATIONS 8.. Remember: what we saw during the last lesson Numerical error propagation, convergence and Richardson
More information(6.2) -The Shooting Method for Nonlinear Problems
(6.2) -The Shooting Method for Nonlinear Problems Consider the boundary value problems (BVPs) for the second order differential equation of the form (*) y f x,y,y, a x b, y a and y b. When f x,y,y is linear
More informationNumerical Methods - Boundary Value Problems for ODEs
Numerical Methods - Boundary Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Boundary Value Problems for ODEs 2013 1 / 14 Outline 1 Boundary Value
More informationAPPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified
More informationChap. 20: Initial-Value Problems
Chap. 20: Initial-Value Problems Ordinary Differential Equations Goal: to solve differential equations of the form: dy dt f t, y The methods in this chapter are all one-step methods and have the general
More informationSYSTEMS OF ODES. mg sin ( (x)) dx 2 =
SYSTEMS OF ODES Consider the pendulum shown below. Assume the rod is of neglible mass, that the pendulum is of mass m, and that the rod is of length `. Assume the pendulum moves in the plane shown, and
More informationMecE 390 Final examination, Winter 2014
MecE 390 Final examination, Winter 2014 Directions: (i) a double-sided 8.5 11 formula sheet is permitted, (ii) no calculators are permitted, (iii) the exam is 80 minutes in duration; please turn your paper
More informationComputation Fluid Dynamics
Computation Fluid Dynamics CFD I Jitesh Gajjar Maths Dept Manchester University Computation Fluid Dynamics p.1/189 Garbage In, Garbage Out We will begin with a discussion of errors. Useful to understand
More informationOptimization and Calculus
Optimization and Calculus To begin, there is a close relationship between finding the roots to a function and optimizing a function. In the former case, we solve for x. In the latter, we solve: g(x) =
More informationInitial Value Problems
Numerical Analysis, lecture 13: Initial Value Problems (textbook sections 10.1-4, 10.7) differential equations standard form existence & uniqueness y 0 y 2 solution methods x 0 x 1 h h x 2 y1 Euler, Heun,
More informationPractice Problems for Final Exam
Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationOrdinary Differential equations
Chapter 1 Ordinary Differential equations 1.1 Introduction to Ordinary Differential equation In mathematics, an ordinary differential equation (ODE) is an equation which contains functions of only one
More informationMultiplication of Polynomials
Summary 391 Chapter 5 SUMMARY Section 5.1 A polynomial in x is defined by a finite sum of terms of the form ax n, where a is a real number and n is a whole number. a is the coefficient of the term. n is
More informationComputational Methods CMSC/AMSC/MAPL 460. Solving nonlinear equations and zero finding. Finding zeroes of functions
Computational Methods CMSC/AMSC/MAPL 460 Solving nonlinear equations and zero finding Ramani Duraiswami, Dept. of Computer Science Where does it arise? Finding zeroes of functions Solving functional equations
More informationNumerical Solutions of Boundary Value Problems
Numerical Solutions of Boundary Value Problems In these note we will consider the solution to boundary value problems of the form y = f(x,y,y ), a < x < b, () y(a) = A, (2) y(b) = B. (3) We will consider
More informationOrdinary Differential Equations (ode)
Ordinary Differential Equations (ode) Numerical Methods for Solving Initial condition (ic) problems and Boundary value problems (bvp) What is an ODE? =,,...,, yx, dx dx dx dx n n 1 n d y d y d y In general,
More informationExercises, module A (ODEs, numerical integration etc)
FYTN HT18 Dept. of Astronomy and Theoretical Physics Lund University, Sweden Exercises, module A (ODEs, numerical integration etc) 1. Use Euler s method to solve y (x) = y(x), y() = 1. (a) Determine y
More informationDifferential Equations FMNN10 Graded Project #1 c G Söderlind 2017
Differential Equations FMNN10 Graded Project #1 c G Söderlind 2017 Published 2017-10-30. Instruction in computer lab 2017-11-02/08/09. Project report due date: Monday 2017-11-13 at 10:00. Goals. The goal
More informationNumerical Analysis II. Problem Sheet 12
P. Grohs S. Hosseini Ž. Kereta Spring Term 2015 Numerical Analysis II ETH Zürich D-MATH Problem Sheet 12 Problem 12.1 The Exponential Runge-Kutta Method Consider the exponential Runge Kutta single-step
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationOrdinary Differential Equations (ODEs)
Ordinary Differential Equations (ODEs) 1 Computer Simulations Why is computation becoming so important in physics? One reason is that most of our analytical tools such as differential calculus are best
More informationAn Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations
An Overly Simplified and Brief Review of Differential Equation Solution Methods We will be dealing with initial or boundary value problems. A typical initial value problem has the form y y 0 y(0) 1 A typical
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More informationThe Fundamental Theorem of Calculus
The Fundamental Theorem of Calculus Objectives Evaluate a definite integral using the Fundamental Theorem of Calculus. Understand and use the Mean Value Theorem for Integrals. Find the average value of
More informationNumerical Methods for Embedded Optimization and Optimal Control. Exercises
Summer Course Numerical Methods for Embedded Optimization and Optimal Control Exercises Moritz Diehl, Daniel Axehill and Lars Eriksson June 2011 Introduction This collection of exercises is intended to
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationGraded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo
2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical
More informationCoupled differential equations
Coupled differential equations Example: dy 1 dy 11 1 1 1 1 1 a y a y b x a y a y b x Consider the case with b1 b d y1 a11 a1 y1 y a1 ay dy y y e y 4/1/18 One way to address this sort of problem, is to
More informationExam 2 Study Guide: MATH 2080: Summer I 2016
Exam Study Guide: MATH 080: Summer I 016 Dr. Peterson June 7 016 First Order Problems Solve the following IVP s by inspection (i.e. guessing). Sketch a careful graph of each solution. (a) u u; u(0) 0.
More informationChapter 11 ORDINARY DIFFERENTIAL EQUATIONS
Chapter 11 ORDINARY DIFFERENTIAL EQUATIONS The general form of a first order differential equations is = f(x, y) with initial condition y(a) = y a We seek the solution y = y(x) for x > a This is shown
More informationHomework 2 - Solutions MA/CS 375, Fall 2005
Homework 2 - Solutions MA/CS 375, Fall 2005 1. Use the bisection method, Newton s method, and the Matlab R function fzero to compute a positive real number x satisfying: sinh x = cos x. For each of the
More informationBOUNDARY VALUE PROBLEMS
BOUNDARY VALUE PROBLEMS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 BOUNDARY VALUE PROBLEMS 3 NEWTONS SHOOTING METHOD 4 SUMMARY OUTLINE 1 REVIEW 2 BOUNDARY VALUE PROBLEMS 3 NEWTONS SHOOTING
More informationModeling & Simulation 2018 Lecture 12. Simulations
Modeling & Simulation 2018 Lecture 12. Simulations Claudio Altafini Automatic Control, ISY Linköping University, Sweden Summary of lecture 7-11 1 / 32 Models of complex systems physical interconnections,
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 information1. Consider the initial value problem: find y(t) such that. y = y 2 t, y(0) = 1.
Engineering Mathematics CHEN30101 solutions to sheet 3 1. Consider the initial value problem: find y(t) such that y = y 2 t, y(0) = 1. Take a step size h = 0.1 and verify that the forward Euler approximation
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationThe family of Runge Kutta methods with two intermediate evaluations is defined by
AM 205: lecture 13 Last time: Numerical solution of ordinary differential equations Today: Additional ODE methods, boundary value problems Thursday s lecture will be given by Thomas Fai Assignment 3 will
More informationNumerical Methods. Root Finding
Numerical Methods Solving Non Linear 1-Dimensional Equations Root Finding Given a real valued function f of one variable (say ), the idea is to find an such that: f() 0 1 Root Finding Eamples Find real
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 75 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP. Use MATLAB solvers for solving higher order ODEs and systems
More informationOutline. Additional Nonlinear Systems. Abstract. Finding Equilibrium Points Numerically. Newton s Method
Outline Finding Equilibrium Points Numerically Additional Nonlinear Systems James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University June 13, 2017
More informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
More informationButcher tableau Can summarize an s + 1 stage Runge Kutta method using a triangular grid of coefficients
AM 205: lecture 13 Last time: ODE convergence and stability, Runge Kutta methods Today: the Butcher tableau, multi-step methods, boundary value problems Butcher tableau Can summarize an s + 1 stage Runge
More informationODE Runge-Kutta methods
ODE Runge-Kutta methods The theory (very short excerpts from lectures) First-order initial value problem We want to approximate the solution Y(x) of a system of first-order ordinary differential equations
More informationThe Shooting Method for Boundary Value Problems
1 The Shooting Method for Boundary Value Problems Consider a boundary value problem of the form y = f(x, y, y ), a x b, y(a) = α, y(b) = β. (1.1) One natural way to approach this problem is to study the
More informationSolutions to Section 1.1
Solutions to Section True-False Review: FALSE A derivative must involve some derivative of the function y f(x), not necessarily the first derivative TRUE The initial conditions accompanying a differential
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationOrdinary Differential Equations: Initial Value problems (IVP)
Chapter Ordinary Differential Equations: Initial Value problems (IVP) Many engineering applications can be modeled as differential equations (DE) In this book, our emphasis is about how to use computer
More informationMATH 2053 Calculus I Review for the Final Exam
MATH 05 Calculus I Review for the Final Exam (x+ x) 9 x 9 1. Find the limit: lim x 0. x. Find the limit: lim x + x x (x ).. Find lim x (x 5) = L, find such that f(x) L < 0.01 whenever 0 < x
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MATLAB sessions: Laboratory 4 MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for
More informationNumerical Methods for the Solution of Differential Equations
Numerical Methods for the Solution of Differential Equations Markus Grasmair Vienna, winter term 2011 2012 Analytical Solutions of Ordinary Differential Equations 1. Find the general solution of the differential
More informationCS520: numerical ODEs (Ch.2)
.. CS520: numerical ODEs (Ch.2) Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca people.cs.ubc.ca/ ascher/520.html Uri Ascher (UBC) CPSC 520: ODEs (Ch. 2) Fall
More informationLesson 9: Predator-Prey and ode45
Lesson 9: Predator-Prey and ode45 9.1 Applied Problem. In this lesson we will allow for more than one population where they depend on each other. One population could be the predator such as a fox, and
More informationGetting Started With The Predator - Prey Model: Nullclines
Getting Started With The Predator - Prey Model: Nullclines James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 28, 2013 Outline The Predator
More informationProblem set 7 Math 207A, Fall 2011 Solutions
Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for solving higher order ODEs and
More informationINTRODUCTION TO NUMERICAL ANALYSIS
INTRODUCTION TO NUMERICAL ANALYSIS Cho, Hyoung Kyu Department of Nuclear Engineering Seoul National University 3. SOLVING NONLINEAR EQUATIONS 3.1 Background 3.2 Estimation of errors in numerical solutions
More informationComputers, Lies and the Fishing Season
1/47 Computers, Lies and the Fishing Season Liz Arnold May 21, 23 Introduction Computers, lies and the fishing season takes a look at computer software programs. As mathematicians, we depend on computers
More informationAlgebra II Through Competitions Chapter 7 Function Composition and Operations
. FUNCTIONS. Definition A function is a relationship between the independent variable x and dependent variable y. Each value of x corresponds exactly one value of y. Note two different values of x can
More informationConsider the following example of a linear system:
LINEAR SYSTEMS Consider the following example of a linear system: Its unique solution is x + 2x 2 + 3x 3 = 5 x + x 3 = 3 3x + x 2 + 3x 3 = 3 x =, x 2 = 0, x 3 = 2 In general we want to solve n equations
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationNumerical methods. Examples with solution
Numerical methods Examples with solution CONTENTS Contents. Nonlinear Equations 3 The bisection method............................ 4 Newton s method.............................. 8. Linear Systems LU-factorization..............................
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationDRAFT - Math 101 Lecture Note - Dr. Said Algarni
2 Limits 2.1 The Tangent Problems The word tangent is derived from the Latin word tangens, which means touching. A tangent line to a curve is a line that touches the curve and a secant line is a line that
More informationThe Derivative. Appendix B. B.1 The Derivative of f. Mappings from IR to IR
Appendix B The Derivative B.1 The Derivative of f In this chapter, we give a short summary of the derivative. Specifically, we want to compare/contrast how the derivative appears for functions whose domain
More informationInitial value problems for ordinary differential equations
AMSC/CMSC 660 Scientific Computing I Fall 2008 UNIT 5: Numerical Solution of Ordinary Differential Equations Part 1 Dianne P. O Leary c 2008 The Plan Initial value problems (ivps) for ordinary differential
More informationSolving Nonlinear Equations
Solving Nonlinear Equations Jijian Fan Department of Economics University of California, Santa Cruz Oct 13 2014 Overview NUMERICALLY solving nonlinear equation Four methods Bisection Function iteration
More informationMotivation: We have already seen an example of a system of nonlinear equations when we studied Gaussian integration (p.8 of integration notes)
AMSC/CMSC 460 Computational Methods, Fall 2007 UNIT 5: Nonlinear Equations Dianne P. O Leary c 2001, 2002, 2007 Solving Nonlinear Equations and Optimization Problems Read Chapter 8. Skip Section 8.1.1.
More informationMon Jan Improved acceleration models: linear and quadratic drag forces. Announcements: Warm-up Exercise:
Math 2250-004 Week 4 notes We will not necessarily finish the material from a given day's notes on that day. We may also add or subtract some material as the week progresses, but these notes represent
More information3.2 A2 - Just Like Derivatives but Backwards
3. A - Just Like Derivatives but Backwards The Definite Integral In the previous lesson, you saw that as the number of rectangles got larger and larger, the values of Ln, Mn, and Rn all grew closer and
More informationMaterials Genome Assessment
Materials Genome Assessment Lecture 12 : Quantum Mechanics Theory: Differential equations with boundary conditions Schrodinger s cat Particle in a Box Programming: Modification of the Pendulum program
More informationMATH 3795 Lecture 13. Numerical Solution of Nonlinear Equations in R N.
MATH 3795 Lecture 13. Numerical Solution of Nonlinear Equations in R N. Dmitriy Leykekhman Fall 2008 Goals Learn about different methods for the solution of F (x) = 0, their advantages and disadvantages.
More informationBasic Theory of Differential Equations
page 104 104 CHAPTER 1 First-Order Differential Equations 16. The following initial-value problem arises in the analysis of a cable suspended between two fixed points y = 1 a 1 + (y ) 2, y(0) = a, y (0)
More informationSolutions to Homework # 1 Math 381, Rice University, Fall (x y) y 2 = 0. Part (b). We make a convenient change of variables:
Hildebrand, Ch. 8, # : Part (a). We compute Subtracting, we eliminate f... Solutions to Homework # Math 38, Rice University, Fall 2003 x = f(x + y) + (x y)f (x + y) y = f(x + y) + (x y)f (x + y). x = 2f(x
More informationMath 551 Homework Assignment 3 Page 1 of 6
Math 551 Homework Assignment 3 Page 1 of 6 Name and section: ID number: E-mail: 1. Consider Newton s method for finding + α with α > 0 by finding the positive root of f(x) = x 2 α = 0. Assuming that x
More informationThe Corrected Trial Solution in the Method of Undetermined Coefficients
Definition of Related Atoms The Basic Trial Solution Method Symbols Superposition Annihilator Polynomial for f(x) Annihilator Equation for f(x) The Corrected Trial Solution in the Method of Undetermined
More informationModule Two: Differential Calculus(continued) synopsis of results and problems (student copy)
Module Two: Differential Calculus(continued) synopsis of results and problems (student copy) Srikanth K S 1 Syllabus Taylor s and Maclaurin s theorems for function of one variable(statement only)- problems.
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationAstronomy 8824: Numerical Methods Notes 2 Ordinary Differential Equations
Astronomy 8824: Numerical Methods Notes 2 Ordinary Differential Equations Reading: Numerical Recipes, chapter on Integration of Ordinary Differential Equations (which is ch. 15, 16, or 17 depending on
More informationScientific Computing with Case Studies SIAM Press, Lecture Notes for Unit V Solution of
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 1 Dianne P. O Leary c 2008 1 The
More informationLecture 7. Root finding I. 1 Introduction. 2 Graphical solution
1 Introduction Lecture 7 Root finding I For our present purposes, root finding is the process of finding a real value of x which solves the equation f (x)=0. Since the equation g x =h x can be rewritten
More informationOptimal Control - Homework Exercise 3
Optimal Control - Homework Exercise 3 December 17, 2010 In this exercise two different problems will be considered, first the so called Zermelo problem the problem is to steer a boat in streaming water,
More informationCLASS NOTES Models, Algorithms and Data: Introduction to computing 2018
CLASS NOTES Models, Algorithms and Data: Introduction to computing 2018 Petros Koumoutsakos, Jens Honore Walther (Last update: April 16, 2018) IMPORTANT DISCLAIMERS 1. REFERENCES: Much of the material
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More information. For each initial condition y(0) = y 0, there exists a. unique solution. In fact, given any point (x, y), there is a unique curve through this point,
1.2. Direction Fields: Graphical Representation of the ODE and its Solution Section Objective(s): Constructing Direction Fields. Interpreting Direction Fields. Definition 1.2.1. A first order ODE of the
More information2D1240 Numerical Methods II / André Jaun, NADA, KTH NADA
A NADA 2D24 Numericl Methods II / André Jun, NADA, KTH. NON-LINEAR LEAST SQUARE FITTING, OPTIMIZATION.. Remember: wht we sw during the lst lesson Itertive solution for sprse sstems (Jcobi, Guss Seidel,
More informationSHOOTING METHOD IN SOLVING BOUNDARY VALUE PROBLEM
www.arpapress.com/volumes/vol21issue1/ijrras_21_1_02.pdf SHOOTING METHOD IN SOLVING BOUNDARY VALUE PROBLEM Badradeen Adam 1 & Mohsin H. A. Hashim 2 1 Department of Mathematics, Faculty of Education, University
More informationLecture 2: Review of Prerequisites. Table of contents
Math 348 Fall 217 Lecture 2: Review of Prerequisites Disclaimer. As we have a textbook, this lecture note is for guidance and supplement only. It should not be relied on when preparing for exams. In this
More informationElementary ODE Review
Elementary ODE Review First Order ODEs First Order Equations Ordinary differential equations of the fm y F(x, y) () are called first der dinary differential equations. There are a variety of techniques
More informationPhysics 110 Homework Solutions Week #6 - Wednesday
Physics 110 Homework Solutions Week #6 - Wednesday Friday, May3, 2013 Chapter 6 Questions - none Multiple-Choice 66 C 67 D 68 B 69 C Problems 612 It s velocity as the ball hits the ground is found from
More informationSolutions to Exam 2, Math 10560
Solutions to Exam, Math 6. Which of the following expressions gives the partial fraction decomposition of the function x + x + f(x = (x (x (x +? Solution: Notice that (x is not an irreducile factor. If
More informationNumerical Methods for ODEs. Lectures for PSU Summer Programs Xiantao Li
Numerical Methods for ODEs Lectures for PSU Summer Programs Xiantao Li Outline Introduction Some Challenges Numerical methods for ODEs Stiff ODEs Accuracy Constrained dynamics Stability Coarse-graining
More information1. Method 1: bisection. The bisection methods starts from two points a 0 and b 0 such that
Chapter 4 Nonlinear equations 4.1 Root finding Consider the problem of solving any nonlinear relation g(x) = h(x) in the real variable x. We rephrase this problem as one of finding the zero (root) of a
More informationResearch Computing with Python, Lecture 7, Numerical Integration and Solving Ordinary Differential Equations
Research Computing with Python, Lecture 7, Numerical Integration and Solving Ordinary Differential Equations Ramses van Zon SciNet HPC Consortium November 25, 2014 Ramses van Zon (SciNet HPC Consortium)Research
More informationReview Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester
HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE
More informationChapter 8. Exploring Polynomial Functions. Jennifer Huss
Chapter 8 Exploring Polynomial Functions Jennifer Huss 8-1 Polynomial Functions The degree of a polynomial is determined by the greatest exponent when there is only one variable (x) in the polynomial Polynomial
More informationy 2y = 4 x, Name Form Solution method
An Introduction to Higher-Order Differential Equations Up to this point in the class, we have only specifically studied solution techniques for first-order differential equations, i.e. equations whose
More informationIntroductory Numerical Analysis
Introductory Numerical Analysis Lecture Notes December 16, 017 Contents 1 Introduction to 1 11 Floating Point Numbers 1 1 Computational Errors 13 Algorithm 3 14 Calculus Review 3 Root Finding 5 1 Bisection
More informationMath Numerical Analysis Mid-Term Test Solutions
Math 400 - Numerical Analysis Mid-Term Test Solutions. Short Answers (a) A sufficient and necessary condition for the bisection method to find a root of f(x) on the interval [a,b] is f(a)f(b) < 0 or f(a)
More informationMath 2a Prac Lectures on Differential Equations
Math 2a Prac Lectures on Differential Equations Prof. Dinakar Ramakrishnan 272 Sloan, 253-37 Caltech Office Hours: Fridays 4 5 PM Based on notes taken in class by Stephanie Laga, with a few added comments
More information