Solution: (a) Before opening the parachute, the differential equation is given by: dv dt. = v. v(0) = 0
|
|
- Rosaline Gray
- 6 years ago
- Views:
Transcription
1 Math 2250 Lab 4 Name/Unid: 1. (35 points) Leslie Leroy Irvin bails out of an airplane at the altitude of 16,000 ft, falls freely for 20 s, then opens his parachute. Assuming linear air resistance ρv ft/s 2, taking ρ=0.18 without the parachute and ρ= 1.8 with the parachute. (a) Write down the ODE that describes the velocity v before he opens the parachute and solve the IVP. (b) Find the velocity v at time 1 s, the velocity when he just opens the parachute and the terminal velocity. (c) Find the height y(t) when he opens the parachute. (d) Write down the ODE describing the velocity after he opens the parachute, write the new initial condition and solve the IVP. (e) Write down the equation for the height depending of time, use the new initial condition for the height. (f) How long will it take him to reach the ground (from the moment he bails out of the plane)? You may want to solve for t using some software. Solution: (a) Before opening the parachute, the differential equation is given by: dv dt = v v(0) = 0 Solving the ODE by integrating factor and using the initial condition, we get: P (t) = 0.18 (1) Q(t) = 32 (2) ρ(t) = e 0.18dt = e 0.18t (3) e 0.18t (v (t) + 0.2v) = 32e 0.18t (4) (v(t)e 0.18t ) = 32e 0.18t (5) v(t)e 0.18t = e 0.18t + C (6) v(0) = 0 = C = (7) v(t) = (e 0.18t 1) (8)
2 (b) Then, we can evaluate v(t) at time t=1 and t=20 (20s is when the man opens the parachute) and t (for the terminal velocity) v(1) = ft/s, v(20) = ft/s, v termi = ft/s (c) For finding the height, integrate v(t) so we get: y(t) = e 0.18t t + c Using the initial condition y(0)=16,000 ft. c= Thus He opens the parachute at t=20s y(t) = e 0.18t t + 16, y(20) = f t (d) After opening the parachute, the differential equation change to: dv = v dt v(0) = (e) Solving again by integrating factor and using the initial condition, we get: Integrating v(t) for getting the height, v(t) = 17.7e 1.8t y(t) = 9.83e 1.8t t + c Using the initial condition y(0) = 13, , the equation for the height is: y(t) = 9.83e 1.8t t (f) To reach the ground means y = 0 Thus solving the last equation for t, y = 0 when t = 86.83s The total of time descending is t = = s. Page 2
3 2. (25 points) Consider the same IVP as in problem 1 (part a) Apply the following methods to approximate the solution function on the time interval 0 t 1. Do these computations by hand (using a calculator or other technology to do each step, but not using programmed code). You may wish to use programmed code (numericswn.m) to check your work before you hand it in. (a) Euler Method: use step size h = 0.5, i.e. 2 time steps. (b) Improved Euler Method: use step size h = 0.5 again. (c) Runge Kutta Method: use step size h = 1.0, i.e. just 1 time step. (d) Compare these approximations for v(1) with the exact solution at t = 1 and comment about the relative accuracy of the three techniques, by computing the relative errors v approx v exact v exact in each case.(the step sizes are so large that you don t expect any of the techniques to be extremely accurate.) Solution: h=1, time interval [0,3], v(0) = 0 (a) Euler Method, v n+1 = v n + h f(t n, v n ) v 0 = 0 = k = = 32 = v 1 = v h = 16 v 1 = 16 = k = ( 16) = = v 2 = h = (b) Improved Euler. step 1: v 0 = 0 = k 1 = 32, k 2 = = v 1 = ( ) = step 2: v 1 = = k 1 = k 2 = = v 2 = (c) R K Method k 1 = f(x n, y n )k 2 = f(x n + h 2, y n + h 2 k 1)k 3 = f(x n + h 2, y n + h 2 k 2)k 4 = f(x n+1, y n + hk 3 )v n+1 1 step only: v 0 = 0 = k1 = 32 k2 = k3 = k4 = = v 1 = v (k 1 + 2k 2 + 2k 3 + k 4 ) = Page 3
4 (d) Results: Using part b) on problem 1, the exact velocity at time t = 1 is v(1) = The relative accuracies for the three methods are Eulers: = = Improved Euler s: = = Runge-Kutta: = = Conclusion: Runge-Kutta is much more accurate. 3. (40 points) If we drop a package from a helicopter with a parachute attached, the wind resistance provided by the parachute is not really linear, assume for a certain type of parachute the acceleration equation is modeled by (Assume velocity in ft/s) v (t) = v 0.16v 1.6 v(0) = 0 (a) What is the terminal velocity in this model? Notice that the slope function v 0.16v 1.6 is decreasing function of v for v > 0 that its value is 32 when v=0, and that its limit as v is. Thus the slope function has exactly one root and you can find it with the Matlab (or other software) solve command. (b) This is a differential equation which does NOT have an elementary solution. Use numericswn.m (You may modify the Matlab code according your need)to estimate the solution values at t=3 and t=6 seconds with the methods Euler, Improved Euler and Runge Kutta with step sizes h=0.3 and h=0.03. Write the approximate solution for each method with each h and each time t=3 and t=6. Plot your results and comment about apparent accuracy of the three methods with these different choices of time step size. (c) Since we don t know the exact solution function in this problem, there is not a quick direct way to see how accurate our numerical approximate solutions are. There are error estimates that some of you will learn in later numerical analysis classes (analogous to error estimates you may have learned for Taylor series), but for the purposes of this lab undertake the following crude test: Compare the numerical approximations for velocity at t = 3 and t = 6, for Runge-Kutta with h = 0.3 and h = If these approximations are close (and the approximate solution graphs appear to line up exactly), then it s likely that they are also close to the exact solution. Page 4
5 (d) Based on your work above, approximately what fraction of the terminal velocity is attained at t = 3 and t = 6 seconds? (e) What does the asymptote on the plot mean? What happen if you change the initial condition to v(0)= 5 and v(0)= 40? Solution: (a) The terminal velocity is found by calculating the root of the right hand side (acceleration equals to zero), using a software, with Maple solve code ( solve(32.2 v.15 v 1.5 = 0, v); ) you get ft/s. Alternatively, computing the Runge Kutta at h = 0.01 for 20 seconds predicts the stable velocity. (b) Plot of the methods: Euler, Improved Euler and RK for h= Numerical Solution. h= v Euler method Improved Euler Runge Kutta method t For h=0.3 the values of v on t=3 are the following: Euler v(3)= Improved Euler v(3)= RK v(3)= For h=0.3 the values of v on t=6 are the following: Euler v(6)= Improved Euler v(6)= RK v(6)= Plot of the methods: Euler, Improved Euler and RK for h=0.03 Page 5
6 30 Numerical Solution. h= v Euler method Improved Euler Runge Kutta method t For h=0.03 the values of v on t=3 are the following: Euler v(3)= Improved Euler v(3)= RK v(3)= For h=0.03 the values of v on t=6 are the following: Euler v(6)= Improved Euler v(6)= RK v(6)= We can see that on our second graph the three methods give us similar numerical solutions. (c) At t = 3 for Runge-Kutta with h = 0.3 and h = 0.03 the values of the velocity are v(3) = and v(3) = respectively. At t = 6 for Runge-Kutta with h = 0.3 and h = 0.03 the values of the velocity are v(6) = and v(6) = respectively. Notice these values are very closed to each other, so it is likely that they are also close to the exact solution. (d) At t = 3, 99.58% of the terminal velocity is attained. At t = 6, 99.99% of the terminal velocity is attained. (e) The asymptote on the plot represents the terminal velocity, if we change the initial condition to v(0)= 5 the graph of the velocity is increasing until attain the terminal velocity and if v(0)= 40 then the graph of the velocity is decreasing until attain the terminal velocity. Page 6
Solution: (a) Before opening the parachute, the differential equation is given by: dv dt. = v. v(0) = 0
Math 2250 Lab 4 Name/Unid: 1. (25 points) A man bails out of an airplane at the altitute of 12,000 ft, falls freely for 20 s, then opens his parachute. Assuming linear air resistance ρv ft/s 2, taking
More informationChapter 9b: Numerical Methods for Calculus and Differential Equations. Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers
Chapter 9b: Numerical Methods for Calculus and Differential Equations Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers Acceleration Initial-Value Problems Consider a skydiver
More informationSection 7.4 Runge-Kutta Methods
Section 7.4 Runge-Kutta Methods Key terms: Taylor methods Taylor series Runge-Kutta; methods linear combinations of function values at intermediate points Alternatives to second order Taylor methods Fourth
More informationNumerical method for approximating the solution of an IVP. Euler Algorithm (the simplest approximation method)
Section 2.7 Euler s Method (Computer Approximation) Key Terms/ Ideas: Numerical method for approximating the solution of an IVP Linear Approximation; Tangent Line Euler Algorithm (the simplest approximation
More informationCHAPTER 1 INTRODUCTION TO NUMERICAL METHOD
CHAPTER 1 INTRODUCTION TO NUMERICAL METHOD Presenter: Dr. Zalilah Sharer 2018 School of Chemical and Energy Engineering Universiti Teknologi Malaysia 16 September 2018 Chemical Engineering, Computer &
More informationMath Lecture 9
Math 2280 - Lecture 9 Dylan Zwick Fall 2013 In the last two lectures we ve talked about differential equations for modeling populations Today we ll return to the theme we touched upon in our second lecture,
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 informationMathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Mathematical Models MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Ordinary Differential Equations The topic of ordinary differential equations (ODEs) is
More informationMathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Spring Department of Mathematics
Mathematical Models MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 Ordinary Differential Equations The topic of ordinary differential equations (ODEs)
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 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 informationLecture 17: Ordinary Differential Equation II. First Order (continued)
Lecture 17: Ordinary Differential Equation II. First Order (continued) 1. Key points Maple commands dsolve dsolve[interactive] dsolve(numeric) 2. Linear first order ODE: y' = q(x) - p(x) y In general,
More informationAPPLICATIONS OF DIFFERENTIATION
4 APPLICATIONS OF DIFFERENTIATION APPLICATIONS OF DIFFERENTIATION 4.9 Antiderivatives In this section, we will learn about: Antiderivatives and how they are useful in solving certain scientific problems.
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationSolving systems of ODEs with Matlab
Solving systems of ODEs with Matlab James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 20, 2013 Outline 1 Systems of ODEs 2 Setting Up
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 information2.1 Differential Equations and Solutions. Blerina Xhabli
2.1 Math 3331 Differential Equations 2.1 Differential Equations and Solutions Blerina Xhabli Department of Mathematics, University of Houston blerina@math.uh.edu math.uh.edu/ blerina/teaching.html Blerina
More informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationOrdinary Differential Equations (ODEs)
Ordinary Differential Equations (ODEs) NRiC Chapter 16. ODEs involve derivatives wrt one independent variable, e.g. time t. ODEs can always be reduced to a set of firstorder equations (involving only first
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 informationNumerical methods for solving ODEs
Chapter 2 Numerical methods for solving ODEs We will study two methods for finding approximate solutions of ODEs. Such methods may be used for (at least) two reasons the ODE does not have an exact solution
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 informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No
More informationFirst Order ODEs (cont). Modeling with First Order ODEs
First Order ODEs (cont). Modeling with First Order ODEs September 11 15, 2017 Bernoulli s ODEs Yuliya Gorb Definition A first order ODE is called a Bernoulli s equation iff it is written in the form y
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 informationMath 2250 Lab 3 Due Date : 2/2/2017
Math 2250 Lab Due Date : 2/2/2017 Name: UID: Unless stated otherwise, show all your work and explain your reasoning. You are allowed to use any results from lecture or the text as long as they are referenced
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 informationEuler s Method, cont d
Jim Lambers MAT 461/561 Spring Semester 009-10 Lecture 3 Notes These notes correspond to Sections 5. and 5.4 in the text. Euler s Method, cont d We conclude our discussion of Euler s method with an example
More informationYou may not use your books, notes; calculators are highly recommended.
Math 301 Winter 2013-14 Midterm 1 02/06/2014 Time Limit: 60 Minutes Name (Print): Instructor This exam contains 8 pages (including this cover page) and 6 problems. Check to see if any pages are missing.
More informationINTRODUCTION TO COMPUTER METHODS FOR O.D.E.
INTRODUCTION TO COMPUTER METHODS FOR O.D.E. 0. Introduction. The goal of this handout is to introduce some of the ideas behind the basic computer algorithms to approximate solutions to differential equations.
More informationMath 116 Second Midterm March 19, 2012
Math 6 Second Midterm March 9, 202 Name: Instructor: Section:. Do not open this exam until you are told to do so. 2. This exam has pages including this cover. There are 9 problems. Note that the problems
More informationMath 308 Week 8 Solutions
Math 38 Week 8 Solutions There is a solution manual to Chapter 4 online: www.pearsoncustom.com/tamu math/. This online solutions manual contains solutions to some of the suggested problems. Here are solutions
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 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 informationMath 307 A - Spring 2015 Final Exam June 10, 2015
Name: Math 307 A - Spring 2015 Final Exam June 10, 2015 Student ID Number: There are 8 pages of questions. In addition, the last page is the basic Laplace transform table. Make sure your exam contains
More informationNumerical Methods for Initial Value Problems; Harmonic Oscillators
Lab 1 Numerical Methods for Initial Value Problems; Harmonic Oscillators Lab Objective: Implement several basic numerical methods for initial value problems (IVPs), and use them to study harmonic oscillators.
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationNumerical method for approximating the solution of an IVP. Euler Algorithm (the simplest approximation method)
Section 2.7 Euler s Method (Computer Approximation) Key Terms/ Ideas: Numerical method for approximating the solution of an IVP Linear Approximation; Tangent Line Euler Algorithm (the simplest approximation
More informationThese will be no tutorials for Math on Tuesday April 26.
Worksheet The purpose of this worksheet is 1. To understand how the differential equation describing simple harmonic motion is derived. 2. To explore how to predict what the solution to this differential
More informationNumerical Solution of Differential Equations
1 Numerical Solution of Differential Equations A differential equation (or "DE") contains derivatives or differentials. In a differential equation the unknown is a function, and the differential equation
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 informationENGI 3424 First Order ODEs Page 1-01
ENGI 344 First Order ODEs Page 1-01 1. Ordinary Differential Equations Equations involving only one independent variable and one or more dependent variables, together with their derivatives with respect
More informationMath Lecture 9
Math 2280 - Lecture 9 Dylan Zwick Spring 2013 In the last two lectures we ve talked about differential equations for modeling populations. Today we ll return to the theme we touched upon in our second
More informationKINEMATICS IN ONE DIMENSION p. 1
KINEMATICS IN ONE DIMENSION p. 1 Motion involves a change in position. Position can be indicated by an x-coordinate on a number line. ex/ A bumblebee flies along a number line... x = 2 when t = 1 sec 2
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 informationMath 216 Final Exam 24 April, 2017
Math 216 Final Exam 24 April, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationNumerical Methods for Initial Value Problems; Harmonic Oscillators
1 Numerical Methods for Initial Value Problems; Harmonic Oscillators Lab Objective: Implement several basic numerical methods for initial value problems (IVPs), and use them to study harmonic oscillators.
More informationDifferential Equation (DE): An equation relating an unknown function and one or more of its derivatives.
Lexicon Differential Equation (DE): An equation relating an unknown function and one or more of its derivatives. Ordinary Differential Equation (ODE): A differential equation that contains only ordinary
More informationMultistep Methods for IVPs. t 0 < t < T
Multistep Methods for IVPs We are still considering the IVP dy dt = f(t,y) t 0 < t < T y(t 0 ) = y 0 So far we have looked at Euler s method, which was a first order method and Runge Kutta (RK) methods
More informationCHAPTER 80 NUMERICAL METHODS FOR FIRST-ORDER DIFFERENTIAL EQUATIONS
CHAPTER 8 NUMERICAL METHODS FOR FIRST-ORDER DIFFERENTIAL EQUATIONS EXERCISE 33 Page 834. Use Euler s method to obtain a numerical solution of the differential equation d d 3, with the initial conditions
More informationChapter 6 - Ordinary Differential Equations
Chapter 6 - Ordinary Differential Equations 7.1 Solving Initial-Value Problems In this chapter, we will be interested in the solution of ordinary differential equations. Ordinary differential equations
More informationSolving Differential Equations: First Steps
30 ORDINARY DIFFERENTIAL EQUATIONS 3 Solving Differential Equations Solving Differential Equations: First Steps Now we start answering the question which is the theme of this book given a differential
More informationMAE 305 Engineering Mathematics I
MAE 305 Engineering Mathematics I Princeton University Assignment # 3 September 26, 1997 Due on Friday, 2PM, October 3, 1997 1. Where we are. We have learned how to deal with a single first order ODE for
More information1 Antiderivatives graphically and numerically
Math B - Calculus by Hughes-Hallett, et al. Chapter 6 - Constructing antiderivatives Prepared by Jason Gaddis Antiderivatives graphically and numerically Definition.. The antiderivative of a function f
More informationMiSP Force and Gravity Worksheet #3, L3
MiSP Force and Gravity Worksheet #3, L3 Name Date FORCE AND ACCELERATION Introduction Today you will view a video of a typical skydive. The skydiver had a skydiving altimeter mounted in a special box with
More informationSolution. Your sketch should show both x and y axes marked from 5 to 5 and circles of radius 1, 3, and 5 centered at the origin.
Solutions of the Sample Problems for the First In-Class Exam Math 246, Fall 208, Professor David Levermore () (a) Sketch the graph that would be produced by the following Matlab command. fplot(@(t) 2/t,
More informationRevision History Date Version Description Author 03/06/ Initial Version Sergio Miguel Martin. Table of Contents
Icarus Physics Engine Manual 1 - Basic Engine Operation Version 1.0 Revision History Date Version Description Author 03/06/2010 1.0 Initial Version Sergio Miguel Martin Table of Contents 1. Introduction
More informationInitial-Value Problems for ODEs. Introduction to Linear Multistep Methods
Initial-Value Problems for ODEs Introduction to Linear Multistep Methods Numerical Analysis (9th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More informationPhysics 584 Computational Methods
Physics 584 Computational Methods Introduction to Matlab and Numerical Solutions to Ordinary Differential Equations Ryan Ogliore April 18 th, 2016 Lecture Outline Introduction to Matlab Numerical Solutions
More informationVariable Step Size Differential Equation Solvers
Math55: Differential Equations 1/30 Variable Step Size Differential Equation Solvers Jason Brewer and George Little Introduction The purpose of developing numerical methods is to approximate the solution
More informationMath Homework 3 Solutions. (1 y sin x) dx + (cos x) dy = 0. = sin x =
2.6 #10: Determine if the equation is exact. If so, solve it. Math 315-01 Homework 3 Solutions (1 y sin x) dx + (cos x) dy = 0 Solution: Let P (x, y) = 1 y sin x and Q(x, y) = cos x. Note P = sin x = Q
More information(1) Rate of change: A swimming pool is emptying at a constant rate of 90 gal/min.
CHAPTER 1 Introduction 1. Bacground Models of physical situations from Calculus (1) Rate of change: A swimming pool is emptying at a constant rate of 90 gal/min. With V = volume in gallons and t = time
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 informationExperiment 2. F r e e F a l l
Suggested Reading for this Lab Experiment F r e e F a l l Taylor, Section.6, and standard deviation rule in Taylor handout. Review Chapters 3 & 4, Read Sections 8.1-8.6. You will also need some procedures
More informationMath 252 Fall 2002 Supplement on Euler s Method
Math 5 Fall 00 Supplement on Euler s Method Introduction. The textbook seems overly enthusiastic about Euler s method. These notes aim to present a more realistic treatment of the value of the method and
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 informationMath 2250 Lab 3 Due September 18, Name:
Math 2250 Lab 3 Due September 18, 2014 Name: 1. Consider the differential equation dx dt = 4x4 x 2 (a) Find the equilbria and draw the phase diagram. (b) Classify the equilibria as stable, unstable, or
More informationCalculus for the Life Sciences
Improved Calculus for the Life Sciences ntial Equations Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center
More informationDifferential Equations and the Parachute Problem
Differential Equations and the Parachute Problem Ronald Phoebus and Cole Reilly May 10, 2004 Abstract The parachute problem is a classical first semester differential equations problem often introduced
More information1 Systems of First Order IVP
cs412: introduction to numerical analysis 12/09/10 Lecture 24: Systems of First Order Differential Equations Instructor: Professor Amos Ron Scribes: Yunpeng Li, Mark Cowlishaw, Nathanael Fillmore 1 Systems
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 informationHomogeneous Equations with Constant Coefficients
Homogeneous Equations with Constant Coefficients MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 General Second Order ODE Second order ODEs have the form
More informationChapter 2 Notes, Kohler & Johnson 2e
Contents 2 First Order Differential Equations 2 2.1 First Order Equations - Existence and Uniqueness Theorems......... 2 2.2 Linear First Order Differential Equations.................... 5 2.2.1 First
More informationBoyce/DiPrima/Meade 11 th ed, Ch 1.1: Basic Mathematical Models; Direction Fields
Boyce/DiPrima/Meade 11 th ed, Ch 1.1: Basic Mathematical Models; Direction Fields Elementary Differential Equations and Boundary Value Problems, 11 th edition, by William E. Boyce, Richard C. DiPrima,
More informationThese will be no tutorials for Math on Tuesday April 26.
Worksheet The purpose of this worksheet is 1. To understand how the differential equation describing simple harmonic motion is derived. 2. To explore how to predict what the solution to this differential
More informationNotation Nodes are data points at which functional values are available or at which you wish to compute functional values At the nodes fx i
LECTURE 6 NUMERICAL DIFFERENTIATION To find discrete approximations to differentiation (since computers can only deal with functional values at discrete points) Uses of numerical differentiation To represent
More informationMath 116 Second Midterm March 19, 2012
Math 6 Second Midterm March 9, 22 Name: EXAM SOLUTIONS Instructor: Section:. Do not open this exam until you are told to do so. 2. This exam has pages including this cover. There are 9 problems. Note that
More informationWe begin exploring Euler s method by looking at direction fields. Consider the direction field below.
Emma Reid- MA 66, Lesson 9 (SU 17) Euler s Method (.7) So far in this course, we have seen some very special types of first order ODEs. We ve seen methods to solve linear, separable, homogeneous, Bernoulli,
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 informationSolving Linear Systems of ODEs with Matlab
Solving Linear Systems of ODEs with Matlab James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 27, 2013 Outline Linear Systems Numerically
More informationInitial value problems for ordinary differential equations
Initial value problems for ordinary differential equations Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 IVP
More informationPhysics 101 Discussion Week 3 Explanation (2011)
Physics 101 Discussion Week 3 Explanation (2011) D3-1. Velocity and Acceleration A. 1: average velocity. Q1. What is the definition of the average velocity v? Let r(t) be the total displacement vector
More informationLab 4: The Simple Pendulum
Page 1 Technical Math II Lab 4: Simple Pendulum Lab 4: The Simple Pendulum Purpose: To investigate the relationship between the length of a simple pendulum and the time it takes to complete a full swing.
More informationBozeman Public Schools Mathematics Curriculum Calculus
Bozeman Public Schools Mathematics Curriculum Calculus Process Standards: Throughout all content standards described below, students use appropriate technology and engage in the mathematical processes
More informationMTH 452/552 Homework 3
MTH 452/552 Homework 3 Do either 1 or 2. 1. (40 points) [Use of ode113 and ode45] This problem can be solved by a modifying the m-files odesample.m and odesampletest.m available from the author s webpage.
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 informationSection , #5. Let Q be the amount of salt in oz in the tank. The scenario can be modeled by a differential equation.
Section.3.5.3, #5. Let Q be the amount of salt in oz in the tank. The scenario can be modeled by a differential equation dq = 1 4 (1 + sin(t) ) + Q, Q(0) = 50. (1) 100 (a) The differential equation given
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 informationOrdinary Differential Equations
Ordinary Differential Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro September 19, 2014 1 / 55 Motivation
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
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 informationSolutions to Homework 3
Solutions to Homework 3 Section 3.4, Repeated Roots; Reduction of Order Q 1). Find the general solution to 2y + y = 0. Answer: The charactertic equation : r 2 2r + 1 = 0, solving it we get r = 1 as a repeated
More informationAntiderivatives. Definition A function, F, is said to be an antiderivative of a function, f, on an interval, I, if. F x f x for all x I.
Antiderivatives Definition A function, F, is said to be an antiderivative of a function, f, on an interval, I, if F x f x for all x I. Theorem If F is an antiderivative of f on I, then every function of
More information(2) Let f(x) = a 2 x if x<2, 4 2x 2 ifx 2. (b) Find the lim f(x). (c) Find all values of a that make f continuous at 2. Justify your answer.
(1) Let f(x) = x x 2 9. (a) Find the domain of f. (b) Write an equation for each vertical asymptote of the graph of f. (c) Write an equation for each horizontal asymptote of the graph of f. (d) Is f odd,
More informationMath 266: Ordinary Differential Equations
Math 266: Ordinary Differential Equations Long Jin Purdue University, Spring 2018 Basic information Lectures: MWF 8:30-9:20(111)/9:30-10:20(121), UNIV 103 Instructor: Long Jin (long249@purdue.edu) Office
More informationIntroduction to Differential Equations
> 22M:16 Fall 05 J. Simon ##################### Introduction to Differential Equations NOTE: This handout is a supplement to your text. There are several homework problems, part of the HW that is due on
More informationfor the initial position and velocity. Find formulas for v t and x t. b) Assuming x 0 = 0 and v 0 O 0, show that the maximum value of x t is x max
Math 2250-4 Wed Aug 28 Finish 1.2. Differential equations of the form y# x = f x. Then begin section 1.3, slope fields., We still have two exercises to finish from Tuesday's notes: the formula for the
More informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Transient Response of a Chemical Reactor Concentration of a substance in a chemical reactor
More information