dy = dx y dx Project 1 Consider the following differential equations: I. xy II. III.
|
|
- Caren Johnson
- 5 years ago
- Views:
Transcription
1 Project 1 Consider the following differential equations: dy = dy = y dy x + sin( x) y = x I. xy II. III. 1. Graph the direction field in the neighborhood of the origin. Graph approximate integral curves for 10 different initial conditions (you can use dfield8.m) 3. Solve the differential equation, using Euler's method (e.g. with euler.m), for three initial conditions and two values of the stepsize, h 4. If possible compute the error, either using the exact solution (if available) or a very accurate numerical solution (from matlab) Please submit a write-up of your project, in electronic form. It should include a statement of the problem, tabular results, graphical results, summary and conclusion. You should also include printouts of any m-files you have written, or modified. You should also give proper citations to any links you may have used. Do your own work.
2 I. Using matlab, we investigate the following differential equation dy = a) Direction Field. The direction field near the origin, , is shown below xy 5 x y b) Integral curves. 10 different initial conditions are shown in the figure below Note that y = 0 is a critical point for the differential equation, but is unstable.
3 c) Numerical solution. The exact solution for y ( 0) = 1 (found by integration factors or using seperability) is given by This can also be found using the Matlab command >> dsolve ( Dy=x*y, y(0)=1, x ) y = e The graph of the results of Euler s method, with h=0.1 and h=0.01, is compared to the exact solution below x Because of the convexity of the exact solution, the numerical solution given by Euler s method always lies below the graph of the solution. The exact solution for y(0)= (found by integration factors or using seperability) is given by y = e This can also be found using the Matlab command >> dsolve( Dy=x*y, y(0)=, x ) The graph of the results of Euler s method, with h=0.1 and h=0.01, is compared to the exact solution below x
4 The exact solution for y(0)=-1 (found by integration factors or using seperability) is given by This can also be found using the Matlab command >> dsolve( Dy=x*y, y(0)=-1, x ) y = e The graph of the results of Euler s method, with h=0.1 and h=0.01, is compared to the exact solution below x Because of the convexity of the exact solution, the numerical solution given by Euler s method always lies above the graph of the solution.
5 d) Error. Error, as a function of h. We fix the initial condition, and the interval [0,1], and then compute the error (defined as the maximum of the absolute value of the pointwise error over the interval) versus the step size h h=step size Error = max(abs(y-exact)) This results in an essentially linear graph of error vs. h. Summary. There is only one critical point for the differential equation dy = xy, given by y = 0 critical point is not stable. The exact solution can be found by the method of integrating factors; the method of seperability; or by using Matlab s built-in dsolve command. The exact solution is given y = Ce x. The by, where C is determined by the initial conditions. Euler s method converges linearly to the exact solution, in terms of step size, for a fixed interval.
6 II. Using matlab, we now investigate the following differential equation dy = a) Direction Field. The direction field near the origin, 5 x 5. 5 y 5, is shown below y b) Integral curves. 10 different initial conditions are shown in the figure below Note that y = 0 is a critical point for the differential equation, but is unstable.
7 c) Numerical solution. The exact solution for y(0)=1 (found by integration factors or using seperability) is given by 1 y = 1 x This can also be found using the Matlab command >> dsolve( Dy=y^, y(0)=1, x ) The graph of the results of Euler s method, with h=0.1 and h=0.01, is compared to the exact solution below Note, the exact solution has a singularity at x=1, so no matter how small the step size, Euler will never be close to the exact solution as x approaches 1! Also, because of convexity, the numerical solution given by Euler s method lies below the exact curve. The exact solution for y(0)= (found by integration factors or using seperability) is given by This can also be found using the Matlab command >> dsolve( Dy=y^, y(0)=, x ) 1 y = 1/ x The graph of the results of Euler s method, with h=0.1 and h=0.01, is compared to the exact solution below
8 The exact solution for y(0)=1/ (found by integration factors or using seperability) is given by This can also be found using the Matlab command >> dsolve( Dy=y^, y(0)=1/, x ) y 1 = The graph of the results of Euler s method, with h=0.1 and h=0.01, is compared to the exact solution below x Notice that since we are well away from the singularity (x=), Euler s method will converge to the exact solution.
9 d) Error. Error, as a function of h. We fix the initial condition (y(0)=1/), and the interval [0,1], and then compute the error (defined as the maximum of the absolute value of the pointwise error over the interval) versus the step size h h=step size Error = max(abs(y-exact)) This is essentially a linear graph Summary. There is only one critical point for the differential equation dy = y, given by y = 0 critical point is not stable. The exact solution can be found by the method of integrating factors; the method of seperability; or by using Matlab s built-in dsolve command. The exact solution is given by 1 y = C x. The, where C is determined by the initial conditions. Euler s method converges linearly to the exact solution, in terms of step size, for a fixed interval (as long as you are sufficiently far from the singularity). Without knowing the exact solution, it would be difficult to predict where the singularity would occur, even by an extremely accurate numerical solver. This is a fundamental problem of nonlinear differential equations.
10 III. The final differential equation we consider is or, in standard form, dy x + sin( x) y = x dy sin( x = x ) y x a) Direction Field. The direction field near the origin, 5 x 5. 5 y 5, is shown below b) Integral curves. 10 different initial conditions are shown in the figure below
11 c) Numerical solution. This equation has no closed form solution, or analytical solution (at least in terms of elementary functions). In order to get an accurate solution to compare to, we use the Matlab command ode45. This can be found using the Matlab command >> [t,y] = ode45(@rhs45,[1,10],[1]); where the listing of rhs.m is given in the Appendix. The graph of the results of Euler s method, with h=0.1, is compared to the numerical solution (ode45) is shown below below Although it may look like the differential equation has a problem at x=0, the term sin(x)/x is bounded (and infinitely smooth) at x=0. This is why Euler s method gives such good results. Summary. There are no critical points for the differential equation. There is no exact solution, in terms of elementary functions. Euler s method converges to the exact solution, in terms of step size, for a fixed interval. Because of the smooth behavior of sin(x)/x near the origin, there is no singularity near x=0. Conclusions. In this project we examine three different ordinary differential equations One linear (with an integrating factor), one linear (with no integrating factor), and one nonlinear. As long as the exact solution remains smooth and bounded, Euler s method has no trouble converging to the exact solution. Convergence is linear in terms of step size. When no exact solution is known, as in the third example, we can still compare Euler s method to a more sophisticated (and accurate) numerical method, such as ode45 which is built-in to Matlab. Matlab provides a power numerical and graphical tool for investigating the direction fields, integral (solution) curves, as well as approximate numerical solutions. In addition to allowing one to program custom numerical routines, Matlab comes with built-in numerical methods which work well.
12 Appendix I. Program listings Euler.m function [tout, yout] = euler(rhs, t0, y0, h, N) % euler('rhs',t0,y0,h,n) % % integrates a single ordinary differential equation % of the form % y'=f(t,y) % where % the initial value of t = t0 % and y(t0) = y0 % the stepsize = h % and the number of steps = N % Note: the right hand side must be defined in the file rhs.m % allocate space and initialize to vectors to zero tout = zeros(n+1,1); yout = zeros(n+1,1); % set initial values tout(1) = t0; yout(1) = y0; % apply Euler's method for i=1:n tout(i+1) = tout(i) + h; slope = feval(rhs,tout(i),yout(i)); yout(i+1) = yout(i) + slope*h; end rhs1.m function yprime = rhs1(x, y) % rhs returns a function which is the right hand side % of a first order differential equation % dy/ = y' = f(x,y) yprime = x*y; rhs.m function yprime = rhs1(x, y) % rhs returns a function which is the right hand side % of a first order differential equation % dy/ = y' = f(x,y) yprime = y^;
13 rhs3.m function yprime = rhs3(x, y) % rhs returns a function which is the right hand side % of a first order differential equation % dy/ = y' = f(x,y) yprime = x-sin(x)*y/x; rhs45.m function [ydot] = rhs45(t, y) % rhs returns a function which is the right hand side % of a first order differential equation % dy/ = y' = f(x,y) ydot = zeros(1,1); ydot(1) = t-sin(t)*y(1)/t;
MAT 275 Laboratory 3 Numerical Solutions by Euler and Improved Euler Methods (scalar equations)
MATLAB sessions: Laboratory 3 1 MAT 275 Laboratory 3 Numerical Solutions by Euler and Improved Euler Methods (scalar equations) In this session we look at basic numerical methods to help us understand
More informationEuler s Methods (a family of Runge- Ku9a methods)
Euler s Methods (a family of Runge- Ku9a methods) ODE IVP An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable: The equation is coupled with an
More informationODE IVP. An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable:
Euler s Methods ODE IVP An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable: The equation is coupled with an initial value/condition (i.e., value
More informationApplication 2.4 Implementing Euler's Method
Application 2.4 Implementing Euler's Method One's understanding of a numerical algorithm is sharpened by considering its implementation in the form of a calculator or computer program. Figure 2.4.13 in
More informationMAT 275 Laboratory 3 Numerical Solutions by Euler and Improved Euler Methods (scalar equations)
MAT 275 Laboratory 3 Numerical Solutions by Euler and Improved Euler Methods (scalar equations) In this session we look at basic numerical methods to help us understand the fundamentals of numerical approximations.
More informationMATH2071: LAB 2: Explicit ODE methods
MATH2071: LAB 2: Explicit ODE methods 1 Introduction Introduction Exercise 1 Euler s method review Exercise 2 The Euler Halfstep (RK2) Method Exercise 3 Runge-Kutta Methods Exercise 4 The Midpoint Method
More informationChE 400: Applied Chemical Engineering Calculations Tutorial 6: Numerical Solution of ODE Using Excel and Matlab
ChE 400: Applied Chemical Engineering Calculations Tutorial 6: Numerical Solution of ODE Using Excel and Matlab Tutorial 6: Numerical Solution of ODE Gerardine G. Botte This handout contains information
More informationIntroduction to Simulink
Introduction to Simulink There are several computer packages for finding solutions of differential equations, such as Maple, Mathematica, Maxima, MATLAB, etc. These systems provide both symbolic and numeric
More informationThe listing says y(1) = How did the computer know this?
18.03 Class 2, Feb 3, 2010 Numerical Methods [1] How do you know the value of e? [2] Euler's method [3] Sources of error [4] Higher order methods [1] The study of differential equations has three parts:.
More informationDirection Fields; Euler s Method
Direction Fields; Euler s Method It frequently happens that we cannot solve first order systems dy (, ) dx = f xy or corresponding initial value problems in terms of formulas. Remarkably, however, this
More informationWorkpackage 5 - Ordinary Differential Equations
Mathematics for I Workpackage 5 - Ordinary Differential Equations Introduction During this laboratory you will be introduced to some of Matlab s facilities for solving ordinary differential equations (ode).
More informationFind the specific function values. Complete parts (a) through (d) below. f (x,y,z) = x y y 2 + z = (Simplify your answer.) ID: 14.1.
. Find the specific function values. Complete parts (a) through (d) below. f (x,y,z) = x y y 2 + z 2 (a) f(2, 4,5) = (b) f 2,, 3 9 = (c) f 0,,0 2 (d) f(4,4,00) = = ID: 4..3 2. Given the function f(x,y)
More informationover The idea is to construct an algorithm to solve the IVP ODE (9.1)
Runge- Ku(a Methods Review of Heun s Method (Deriva:on from Integra:on) The idea is to construct an algorithm to solve the IVP ODE (9.1) over To obtain the solution point we can use the fundamental theorem
More information(Part - 1) P. Sam Johnson. April 14, Numerical Solution of. Ordinary Differential Equations. (Part - 1) Sam Johnson. NIT Karnataka.
P. April 14, 2015 1/51 Overview We discuss the following important methods of solving ordinary differential equations of first / second order. Picard s method of successive approximations (Method of successive
More informationApproximate Solutions to Differential Equations Slope Fields (graphical) and Euler s Method (numeric) by Dave Slomer
Approximate Solutions to Differential Equations Slope Fields (graphical) and Euler s Method (numeric) by Dave Slomer Leonhard Euler was a great Swiss mathematician. The base of the natural logarithm function,
More informationMathematical Methods and Modeling Laboratory class. Numerical Integration of Ordinary Differential Equations
Mathematical Methods and Modeling Laboratory class Numerical Integration of Ordinary Differential Equations Exact Solutions of ODEs Cauchy s Initial Value Problem in normal form: Recall: if f is locally
More informationOrdinary differential equations solving methods
Radim Hošek, A07237 radhost@students.zcu.cz Ordinary differential equations solving methods Problem: y = y2 (1) y = x y (2) y = sin ( + y 2 ) (3) Where it is possible we try to solve the equations analytically,
More informationSection 18-1: Graphical Representation of Linear Equations and Functions
Section 18-1: Graphical Representation of Linear Equations and Functions Prepare a table of solutions and locate the solutions on a coordinate system: f(x) = 2x 5 Learning Outcome 2 Write x + 3 = 5 as
More informationReview Initial Value Problems Euler s Method Summary
THE EULER METHOD P.V. Johnson School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 INITIAL VALUE PROBLEMS The Problem Posing a Problem 3 EULER S METHOD Method Errors 4 SUMMARY OUTLINE 1 REVIEW 2 INITIAL
More informationMA 114 Worksheet #17: Average value of a function
Spring 2019 MA 114 Worksheet 17 Thursday, 7 March 2019 MA 114 Worksheet #17: Average value of a function 1. Write down the equation for the average value of an integrable function f(x) on [a, b]. 2. Find
More informationDifferential Equations (92.236) Listing of Matlab Lab Exercises
Differential Equations (92.236) Listing of Matlab Lab Exercises This page contains a summary list of the Matlab lab exercises available for this course. We will have roughly 10 12 lab sessions that highlight
More informationAP Calculus BC Course Description
AP Calculus BC Course Description COURSE OUTLINE: The following topics define the AP Calculus BC course as it is taught over three trimesters, each consisting of twelve week grading periods. Limits and
More informationEuler and improved Euler using SAGE
Euler and improved Euler using SAGE The goal is to find an approximate solution to the problem y = f(x, y), y(a) = c, (1) where f(x, y) is some given function We shall try to approximate the value of the
More informationWhat you will learn today
What you will learn today Tangent Planes and Linear Approximation and the Gradient Vector Vector Functions 1/21 Recall in one-variable calculus, as we zoom in toward a point on a curve, the graph becomes
More informationApplication 7.6A The Runge-Kutta Method for 2-Dimensional Systems
Application 7.6A The Runge-Kutta Method for -Dimensional Systems Figure 7.6. in the text lists TI-85 and BASIC versions of the program RKDIM that implements the Runge-Kutta iteration k (,, ) = f tn xn
More informationODEs occur quite often in physics and astrophysics: Wave Equation in 1-D stellar structure equations hydrostatic equation in atmospheres orbits
Solving ODEs General Stuff ODEs occur quite often in physics and astrophysics: Wave Equation in 1-D stellar structure equations hydrostatic equation in atmospheres orbits need workhorse solvers to deal
More informationModule1: Numerical Solution of Ordinary Differential Equations. Lecture 6. Higher order Runge Kutta Methods
Module1: Numerical Solution of Ordinary Differential Equations Lecture 6 Higher order Runge Kutta Methods Keywords: higher order methods, functional evaluations, accuracy Higher order Runge Kutta Methods
More informationNew Mexico Tech Hyd 510
Numerics Motivation Modeling process (JLW) To construct a model we assemble and synthesize data and other information to formulate a conceptual model of the situation. The model is conditioned on the science
More informationMATH 19520/51 Class 6
MATH 19520/51 Class 6 Minh-Tam Trinh University of Chicago 2017-10-06 1 Review partial derivatives. 2 Review equations of planes. 3 Review tangent lines in single-variable calculus. 4 Tangent planes to
More informationc x y f() f (x) Determine the Determine the Approximate c : Replacin on the AP exam: under-approximation
Tangent Lines and Linear Approximations Students should be able to: Determine the slope of tangent line to a curve at a point Determine the equations of tangent lines and normal lines Approximate a value
More informationEast Penn School District Secondary Curriculum
East Penn School District Secondary Curriculum A Planned Course Statement for Analytic Geometry and Calculus (BC) AP Course # 360 Grade(s) 12 Department: Math ength of Period (mins.) 41 Total Clock Hours:
More informationUniversity of California, Berkeley
University of California, Berkeley FINAL EXAMINATION, Fall 2012 DURATION: 3 hours Department of Mathematics MATH 53 Multivariable Calculus Examiner: Sean Fitzpatrick Total: 100 points Family Name: Given
More informationover The idea is to construct an algorithm to solve the IVP ODE (8.1)
Runge- Ku(a Methods Review of Heun s Method (Deriva:on from Integra:on) The idea is to construct an algorithm to solve the IVP ODE (8.1) over To obtain the solution point we can use the fundamental theorem
More informationPolymath 6. Overview
Polymath 6 Overview Main Polymath Menu LEQ: Linear Equations Solver. Enter (in matrix form) and solve a new system of simultaneous linear equations. NLE: Nonlinear Equations Solver. Enter and solve a new
More informationUsing Modified Euler Method (MEM) for the Solution of some First Order Differential Equations with Initial Value Problems (IVPs).
Using Modified Euler Method (MEM) for the Solution of some First Order Differential Equations with Initial Value Problems (IVPs). D.I. Lanlege, Ph.D. * ; U.M. Garba, B.Sc.; and A. Aluebho, B.Sc. Department
More informationIntroduction to Computational Mathematics
Introduction to Computational Mathematics Introduction Computational Mathematics: Concerned with the design, analysis, and implementation of algorithms for the numerical solution of problems that have
More informationLAB #8 Numerical Methods
LAB #8 Numerical Methods Goal: The purpose of this lab is to explain how computers numerically approximate solutions to differential equations. Required tools: Matlab routine dfield ; numerical routines
More informationFirst of all, we need to know what it means for a parameterize curve to be differentiable. FACT:
CALCULUS WITH PARAMETERIZED CURVES In calculus I we learned how to differentiate and integrate functions. In the chapter covering the applications of the integral, we learned how to find the length of
More informationWe can conclude that if f is differentiable in an interval containing a, then. f(x) L(x) = f(a) + f (a)(x a).
= sin( x) = 8 Lecture :Linear Approximations and Differentials Consider a point on a smooth curve y = f(x), say P = (a, f(a)), If we draw a tangent line to the curve at the point P, we can see from the
More informationIntroduction to Motion
Date Partners Objectives: Introduction to Motion To investigate how motion appears on a position versus time graph To investigate how motion appears on a velocity versus time graph and the relationship
More informationTangent Lines and Linear Approximations Solutions
Solutions We have intentionally included more material than can be covered in most Student Study Sessions to account for groups that are able to answer the questions at a faster rate. Use your own judgment,
More informationENGI Parametric & Polar Curves Page 2-01
ENGI 3425 2. Parametric & Polar Curves Page 2-01 2. Parametric and Polar Curves Contents: 2.1 Parametric Vector Functions 2.2 Parametric Curve Sketching 2.3 Polar Coordinates r f 2.4 Polar Curve Sketching
More informationMathematics for chemical engineers
Mathematics for chemical engineers Drahoslava Janovská Department of mathematics Winter semester 2013-2014 Numerical solution of ordinary differential equations Initial value problem Outline 1 Introduction
More informationDouble Pendulum. Freddie Witherden. February 10, 2009
Double Pendulum Freddie Witherden February 10, 2009 Abstract We report on the numerical modelling of a double pendulum using C++. The system was found to be very sensitive to both the initial starting
More informationNAG Library Chapter Introduction. D04 Numerical Differentiation
D4 Numerical Differentiation Introduction D4 NAG Library Chapter Introduction D4 Numerical Differentiation Contents 1 Scope of the Chapter... 2 2 Background to the Problems... 2 2.1 Description of the
More informationFunctions of Two variables.
Functions of Two variables. Ferdinánd Filip filip.ferdinand@bgk.uni-obuda.hu siva.banki.hu/jegyzetek 27 February 217 Ferdinánd Filip 27 February 217 Functions of Two variables. 1 / 36 Table of contents
More informationObjectives. Materials
Activity 13 Objectives Understand what a slope field represents in terms of Create a slope field for a given differential equation Materials TI-84 Plus / TI-83 Plus Graph paper Introduction One of the
More informationApplications of MATLAB/Simulink for Process Dynamics and Control
Applications of MATLAB/Simulink for Process Dynamics and Control (This lecture was modified from slides provided by Professor Kirk Dolan and Wei Liao at MSU and Venkat Subramanian at WashU) Outline Introduction
More information8 Piecewise Polynomial Interpolation
Applied Math Notes by R. J. LeVeque 8 Piecewise Polynomial Interpolation 8. Pitfalls of high order interpolation Suppose we know the value of a function at several points on an interval and we wish to
More informationPolar (BC Only) They are necessary to find the derivative of a polar curve in x- and y-coordinates. The derivative
Polar (BC Only) Polar coordinates are another way of expressing points in a plane. Instead of being centered at an origin and moving horizontally or vertically, polar coordinates are centered at the pole
More informationAOE 5204: Homework Assignment 4 Due: Wednesday, 9/21 in class. Extended to Friday 9/23 Online Students:
AOE 5204: Homework Assignment 4 Due: Wednesday, 9/21 in class. Extended to Friday 9/23 Online Students: Email (cdhall@vt.edu) by 5 PM Suppose F b is initially aligned with F i and at t = 0 begins to rotate
More information4 Visualization and. Approximation
4 Visualization and Approximation b A slope field for the differential equation y tan(x + y) tan(x) tan(y). It is not always possible to write down an explicit formula for the solution to a differential
More informationA Study on Numerical Exact Solution of Euler, Improved Euler and Runge - Kutta Method
A Study on Numerical Exact Solution of Euler, Improved Euler and Runge - Kutta Method Mrs. P. Sampoornam P.K.R Arts College for Women, Gobichettipalayam Abstract: In this paper, I will discuss the Runge
More informationMath 128a - Programming Project 1 - ANSWERS. The main program is called ZeroCurve.m, and it calls three other functions:
Math 8a - Programming Project - ANSWERS The main program is called ZeroCurve.m, and it calls three other functions: Dist.m which just computes Euclidean distances between points, NextZero.m which, given
More informationQuasilinear First-Order PDEs
MODULE 2: FIRST-ORDER PARTIAL DIFFERENTIAL EQUATIONS 16 Lecture 3 Quasilinear First-Order PDEs A first order quasilinear PDE is of the form a(x, y, z) + b(x, y, z) x y = c(x, y, z). (1) Such equations
More informationObjectives. Materials
S Activity 4 Objectives Materials Understand what a slope field represents in terms of dy Create a slope field for a given differential equation T-84 Plus / T-83 Plus Graph paper ntroduction ntroduction
More informationParametric and Polar Curves
Chapter 2 Parametric and Polar Curves 2.1 Parametric Equations; Tangent Lines and Arc Length for Parametric Curves Parametric Equations So far we ve described a curve by giving an equation that the coordinates
More informationParametric and Polar Curves
Chapter 2 Parametric and Polar Curves 2.1 Parametric Equations; Tangent Lines and Arc Length for Parametric Curves Parametric Equations So far we ve described a curve by giving an equation that the coordinates
More informationDirectional Derivatives as Vectors
Directional Derivatives as Vectors John Ganci 1 Al Lehnen 2 1 Richland College Dallas, TX jganci@dcccd.edu 2 Madison Area Technical College Madison, WI alehnen@matcmadison.edu Statement of problem We are
More informationParameter Estimation in Differential Equations: A Numerical Study of Shooting Methods
Parameter Estimation in Differential Equations: A Numerical Study of Shooting Methods Franz Hamilton Faculty Advisor: Dr Timothy Sauer January 5, 2011 Abstract Differential equation modeling is central
More informationA First Course on Kinetics and Reaction Engineering Example S4.5
Example S4.5 Problem Purpose The purpose of this example is to illustrate how to use the MATLAB template file FitNumDifMR.m to perform a regression analysis for multiple response data with a model that
More informationOrdinary Differential Equations
Next: Partial Differential Equations Up: Numerical Analysis for Chemical Previous: Numerical Differentiation and Integration Subsections Runge-Kutta Methods Euler's Method Improvement of Euler's Method
More informationObjectives. Materials
Activity 6 Local Linearity, Differentiability, and Limits of Difference Quotients Objectives Connect the concept of local linearity to differentiability through numerical explorations of limits of difference
More informationCS 450 Numerical Analysis. Chapter 7: Interpolation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationtutorial i: Using matcont for numerical integration of ODEs
tutorial i: Using matcont for numerical integration of ODEs u.a. Kuznetsov Department of Mathematics Utrecht University Budapestlaan TA, Utrecht September, This session illustrates how to input a system
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationLecture VIII. Global Approximation Methods: I
Lecture VIII Global Approximation Methods: I Gianluca Violante New York University Quantitative Macroeconomics G. Violante, Global Methods p. 1 /29 Global function approximation Global methods: function
More informationDifferentiation. The Derivative and the Tangent Line Problem 10/9/2014. Copyright Cengage Learning. All rights reserved.
Differentiation Copyright Cengage Learning. All rights reserved. The Derivative and the Tangent Line Problem Copyright Cengage Learning. All rights reserved. 1 Objectives Find the slope of the tangent
More informationKeywords: TDC, Submarine, Torpedo, World War II
Self Position Keeper Modeling Abstract During World War II, all the submarines of all warring nations used analog computers to assist submarine commanders with aiming straight-running torpedoes (e.g. UK
More informationChapter 6. Curves and Surfaces. 6.1 Graphs as Surfaces
Chapter 6 Curves and Surfaces In Chapter 2 a plane is defined as the zero set of a linear function in R 3. It is expected a surface is the zero set of a differentiable function in R n. To motivate, graphs
More informationNumerical Aspects of Special Functions
Numerical Aspects of Special Functions Nico M. Temme In collaboration with Amparo Gil and Javier Segura, Santander, Spain. Nico.Temme@cwi.nl Centrum voor Wiskunde en Informatica (CWI), Amsterdam Numerics
More informationA Hybrid Magnetic Field Solver Using a Combined Finite Element/Boundary Element Field Solver
A Hybrid Magnetic Field Solver Using a Combined Finite Element/Boundary Element Field Solver Abstract - The dominant method to solve magnetic field problems is the finite element method. It has been used
More informationChapter 10 Homework: Parametric Equations and Polar Coordinates
Chapter 1 Homework: Parametric Equations and Polar Coordinates Name Homework 1.2 1. Consider the parametric equations x = t and y = 3 t. a. Construct a table of values for t =, 1, 2, 3, and 4 b. Plot the
More informationNAME: Section # SSN: X X X X
Math 155 FINAL EXAM A May 5, 2003 NAME: Section # SSN: X X X X Question Grade 1 5 (out of 25) 6 10 (out of 25) 11 (out of 20) 12 (out of 20) 13 (out of 10) 14 (out of 10) 15 (out of 16) 16 (out of 24)
More informationMath 126 Winter CHECK that your exam contains 8 problems.
Math 126 Winter 2016 Your Name Your Signature Student ID # Quiz Section Professor s Name TA s Name CHECK that your exam contains 8 problems. This exam is closed book. You may use one 8 1 11 sheet of hand-written
More informationf xx (x, y) = 6 + 6x f xy (x, y) = 0 f yy (x, y) = y In general, the quantity that we re interested in is
1. Let f(x, y) = 5 + 3x 2 + 3y 2 + 2y 3 + x 3. (a) Final all critical points of f. (b) Use the second derivatives test to classify the critical points you found in (a) as a local maximum, local minimum,
More information=.1( sin(2πx/24) y) y(0) = 70.
Differential Equations, Spring 2017 Computer Project 2, Due 11:30 am, Friday, March 10 The goal of this project is to implement the Euler Method and the Improved Euler Method, and to use these methods
More informationMath Lab 6: Powerful Fun with Power Series Representations of Functions Due noon Thu. Jan. 11 in class *note new due time, location for winter quarter
Matter & Motion Winter 2017 18 Name: Math Lab 6: Powerful Fun with Power Series Representations of Functions Due noon Thu. Jan. 11 in class *note new due time, location for winter quarter Goals: 1. Practice
More informationMath 133A, September 10: Numerical Methods (Euler & Runge-Kutta) Section 1: Numerical Methods
Math 133A, September 10: Numerical Methods (Euler & Runge-Kutta) Section 1: Numerical Methods We have seen examples of first-order differential equations which can, by one method of another, be solved
More informationCalculus II (Math 122) Final Exam, 11 December 2013
Name ID number Sections B Calculus II (Math 122) Final Exam, 11 December 2013 This is a closed book exam. Notes and calculators are not allowed. A table of trigonometric identities is attached. To receive
More informationEF432. Introduction to spagetti and meatballs
EF432 Introduction to spagetti and meatballs CSC 418/2504: Computer Graphics Course web site (includes course information sheet): http://www.dgp.toronto.edu/~karan/courses/418/fall2015 Instructor: Karan
More informationSystematic Calculation of Loop Points for Parametric Curves
Systematic Calculation of Loop Points for Parametric Curves Skip Thompson Department of Mathematics & Statistics Radford University Radford, VA 24142 thompson@radford.edu January 11, 2013 Abstract We consider
More informationHomework Set #2-3, Math 475B
Homework Set #2-3, Math 475B Part I: Matlab In the last semester you learned a number of essential features of MATLAB. 1. In this instance, you will learn to make 3D plots and contour plots of z = f(x,
More informationJournal of Engineering Research and Studies E-ISSN
Journal of Engineering Research and Studies E-ISS 0976-79 Research Article SPECTRAL SOLUTIO OF STEADY STATE CODUCTIO I ARBITRARY QUADRILATERAL DOMAIS Alavani Chitra R 1*, Joshi Pallavi A 1, S Pavitran
More informationMAE 384 Numerical Methods for Engineers
MAE 384 Numerical Methods for Engineers Instructor: Huei-Ping Huang office: ERC 359, email: hp.huang@asu.edu (Huei rhymes with way ) Tu/Th 9:00-10:15 PM WGHL 101 Textbook: Numerical Methods for Engineers
More informationComputer Graphics. - Rasterization - Philipp Slusallek
Computer Graphics - Rasterization - Philipp Slusallek Rasterization Definition Given some geometry (point, 2D line, circle, triangle, polygon, ), specify which pixels of a raster display each primitive
More informationTopic 2.3: Tangent Planes, Differentiability, and Linear Approximations. Textbook: Section 14.4
Topic 2.3: Tangent Planes, Differentiability, and Linear Approximations Textbook: Section 14.4 Warm-Up: Graph the Cone & the Paraboloid paraboloid f (x, y) = x 2 + y 2 cone g(x, y) = x 2 + y 2 Do you notice
More informationSYDE 312 UNIT 4: EXTRA QUADRATURE PROBLEMS
SYDE 3 UNIT 4: EXTRA QUADRATURE PROBLEMS Problem E x x + Use f(x) = /x x + You can t use the Matlab gaussquad file in the provided form, because it implements composite Gauss-Legendre quadrature with multiple
More information(1) Tangent Lines on Surfaces, (2) Partial Derivatives, (3) Notation and Higher Order Derivatives.
Section 11.3 Partial Derivatives (1) Tangent Lines on Surfaces, (2) Partial Derivatives, (3) Notation and Higher Order Derivatives. MATH 127 (Section 11.3) Partial Derivatives The University of Kansas
More information14.5 Directional Derivatives and the Gradient Vector
14.5 Directional Derivatives and the Gradient Vector 1. Directional Derivatives. Recall z = f (x, y) and the partial derivatives f x and f y are defined as f (x 0 + h, y 0 ) f (x 0, y 0 ) f x (x 0, y 0
More informationCode from Maple, MATLAB, and Fortran for Example 1 from "Extending Explicit and Linearly Implicit ODE Solvers for Index-1 DAEs."
Appendix A Code from Maple, MATLAB, and Fortran for Example 1 from "Extending Explicit and Linearly Implicit ODE Solvers for Index-1 DAEs." 1. Maple Code (using dsolve s rkf45) Use y1, y2, etc. for all
More informationGRAPHING WORKSHOP. A graph of an equation is an illustration of a set of points whose coordinates satisfy the equation.
GRAPHING WORKSHOP A graph of an equation is an illustration of a set of points whose coordinates satisfy the equation. The figure below shows a straight line drawn through the three points (2, 3), (-3,-2),
More informationMath 226A Homework 4 Due Monday, December 11th
Math 226A Homework 4 Due Monday, December 11th 1. (a) Show that the polynomial 2 n (T n+1 (x) T n 1 (x)), is the unique monic polynomial of degree n + 1 with roots at the Chebyshev points x k = cos ( )
More informationExact equations are first order DEs of the form M(x, y) + N(x, y) y' = 0 for which we can find a function f(x, φ(x)) so that
Section 2.6 Exact Equations (ONLY) Key Terms: Exact equations are first order DEs of the form M(x, y) + N(x, y) y' = 0 for which we can find a function f(x, φ(x)) so that The construction of f(x, φ(x))
More informationVolumes of Solids of Revolution Lecture #6 a
Volumes of Solids of Revolution Lecture #6 a Sphereoid Parabaloid Hyperboloid Whateveroid Volumes Calculating 3-D Space an Object Occupies Take a cross-sectional slice. Compute the area of the slice. Multiply
More informationPURE MATHEMATICS 212 Multivariable Calculus CONTENTS. Page. 1. Assignment Summary... i 2. Summary Assignments...2
PURE MATHEMATICS 212 Multivariable Calculus CONTENTS Page 1. Assignment Summary... i 2. Summary...1 3. Assignments...2 i PMTH212, Multivariable Calculus Assignment Summary 2010 Assignment Date to be Posted
More informationFinite difference methods
Finite difference methods Siltanen/Railo/Kaarnioja Spring 8 Applications of matrix computations Applications of matrix computations Finite difference methods Spring 8 / Introduction Finite difference methods
More informationParametric and Polar Curves
Chapter 2 Parametric and Polar Curves 2.1 Parametric Equations; Tangent Lines and Arc Length for Parametric Curves Parametric Equations So far we ve described a curve by giving an equation that the coordinates
More informationProceedings of the Third International DERIVE/TI-92 Conference
Using the TI-92 and TI-92 Plus to Explore Derivatives, Riemann Sums, and Differential Equations with Symbolic Manipulation, Interactive Geometry, Scripts, Regression, and Slope Fields Sally Thomas, Orange
More informationWinfried Auzinger, Dirk Praetorius (SS2012)
Computer Algebra using Maple Part V: Extensions, some special topics. Didactical aspects Winfried Auzinger, Dirk Praetorius (SS2012) restart: 1 Polynomial equations and systems For polynomial equations
More information