Symbolic Solution of higher order equations
|
|
- Oscar Antony Bennett
- 5 years ago
- Views:
Transcription
1 Math Assignment 4 - Higher Order Equations and Systems of Equations Due: Monday, April 16. Nothing accepted after Tuesday, April 17. This is worth 15 points. 10% points off for being late. You may work by yourself or in pairs. If you work in pairs it is only necessary to turn in one set of answers for each pair. Put both your names on it. Do Problems 1-6 in the discussion below. Copy and paste the input and output for each problem from the Command Window (or plot window) to a Microsoft Word file. Do this after you are satisfied with your input and output. Label the work for each problem in the Word file: Problem 1, Problem 2, etc. Print a copy of this file to hand in. Or use a Live Script. Symbolic Solution of higher order equations MATLAB can find explicit solutions to some higher order equations using dsolve, but unless the equation has constant coefficients or a very particular form, solutions cannot be written explicitly using elementary functions. Consider the equation y'' + 3y' + 2y = 0. Here is the general solution. >> syms y(t) >> dsolve(diff(y,2) + 3*diff(y, 1) + 2*y == 0) C1*exp(-t) + C2*exp(-2*t) We can give partial initial conditions >> dsolve(diff(y,2) + 3*diff(y, 1) + 2*y == 0, y(0) = 2) C5*exp(-2*t) - exp(-t)*(c5-2) or complete initial conditions 1
2 >> Dy = diff(y, t); >> dsolve(diff(y,2) + 3*diff(y, 1) + 2*y == 0, y(0) = 2, Dy(0) == 1) 5*exp(-t) - 3*exp(-2*t) Problem 1: Find the general solution of y'' +2y' +17y = 0. Find the solution that satisfies the initial conditions y(0) = 1 and y'(0) = -1. Plot this solution for 0 t 10. Linear equations with variable coefficients are typically unsolvable, or have solutions expressed using special functions defined for this purpose. For example, the solution to the equation y'' + ty = 0 involves the Airy functions. >> dsolve(diff(y,2) + t*y == 0) C8*airy(0, -t) + C9*airy(2, -t) Problem 2: Find the general solution of t 2 y'' - 12y = 0. Find the solution that satisfies the initial conditions y(1) = 3 and y'(1) = 5. Plot this solution for 0.1 < t 4. Here is the nonlinear equation y'' + y' y = 0. >> dsolve(diff(y,2) + diff(y,1)/y == 0) Warning: Explicit solution could not be found; implicit solution returned. solve(int(1/(c10 - log(y)), y, 'IgnoreSpecialCases', true, 'IgnoreAnalyticConstraints', true) - t - C12 == 0, y) Problem 3: Find the general solution of 4t 2 y'' + 4ty' + (4t 2-1)y = 0. 2
3 Symbolic solution of systems One can solve simple systems of linear differential equations with constant coefficients by hand either using substitution, elimination with differential operators or by writing the system in matrix form and solving using eigenvalues and eigenvectors. In MATLAB one can use dsolve to solve a constant coefficient linear system. We write the system as a vector of equations. For example, to solve x = y 2 with the initial conditions x(0) = 1 y(0) = 2 we do the following. y = 3 2x >> syms x(t) y(t) >> [x(t), y(t)] = dsolve([diff(x,t) == y-2, diff(y,t) == 3-2*x], x(0) == 1, y(0) == 0) x(t) = 3/2 - (3*cos(atan(2*2^(1/2)) - 2^(1/2)*t))/2 y(t) = 2 - (3*8^(1/2)*cos(atan(2^(1/2)/4) + 2^(1/2)*t))/4 There are three different plots that reveal different properties of the solutions. The first two are just plots of x and y vs t as t varies over some range. Here are these two for the solution in the example above. >> tspan = linspace(0, 10, 100); >> plot(tspan, x(tspan)) >> title('x = x(t)') >> ylabel('x') >> plot(tspan, y(tspan)) >> title('y = y(t)') 3
4 As you can see both from the formulas and the graphs, both x and y are periodic functions of t. The third graph plots the ordered pairs (x(t), y(t)) as t varies over some range. This graph is a curve in the xy plane called the trajectory. Here is the trajectory for this solution. >> plot(x(tspan), y(tspan)) >> title('trajectory: (x(t), y(t))') >> xlabel('x') >> axis equal >> axis([-0.5, 3.5, -0.5, 4.5]) The trajectory is the ellipse at the right. The command axis equal made the scales on the two axes the same. The trajectories of the solutions of systems of differential equations can have a variety of interesting shapes. x = 2x + y e2t Problem 4:Find the general solution of y = x + 2y. Find the solution that satisfies the initial conditions x(0) = 1 y(0) = -1. Make the following three plots. First, plot x versus t for 0 t 10. Then, plot y versus t for 0 t 10. Finally, plot y versus x for 0 t 10 to create the trajectory. Numerical approximations We can approximate solutions to systems of equations with ode45 just like we did with a single equation. Recall ode45 is a modified Runge-Kutta method which is good for most equations. 4
5 x = xy + 2 x(0) = - 2 Consider the system y = y + 2x + 1 with the initial conditions y(0) = 1.5. Suppose we want an approximate solution on the interval 0 t 3. Before using ode45 we need to write the system using vector notation, i.e. we write the two equations in terms of x1 and x2 instead of x and y. Letting x1 = x and x2 = y, the equations become x 1 = x 1x2 + 2 x2 = x2 + 2x solution as follows. x1(0) = - 2 and the initial conditions become x 2(0) = 1.5. Then we generate a numerical >> [tsol, xsol] = ode45(@(t,x) [x(1)*x(2)+2; x(2)+2*x(1)+1], [0 3], [-2 1.5]); The equations are listed as entries in a column vector, and the initial conditions x1(0) = 2, x2(0) = 1.5 are given as a vector as well. Again, there are three plots: x vs t, y vs t and finally the trajectory where we plot (x(t), y(t)) for t in some interval >> plot(tsol, xsol(:,1)) >> title('x = x(t)') >> ylabel('x') >> plot(tsol, xsol(:,2)) >> title('y = y(t)') >> plot(xsol(:,1), xsol(:,2)) >> title('trajectory: (x(t), y(t))') >> xlabel('x') >> axis equal This produces the graphs below. 5
6 x = 2x - 2xy Problem 5: Find a numerical approximation to the solution of y = xy - y that satisfies the initial conditions x(0) = 1 y(0) = 3. Plot x and y versus t for 0 t 10 and also the corresponding trajectory in the xy plane. To approximate the solution of a higher order equation we need to rewrite it a first order system by introducing the new variables y1 = y and y2 = y'. For example, suppose we want to approximate the solution to y'' + 3y' + 2y = 0 that satisfies the initial conditions y(0) = 2 and y'(0) = 1. We let y1 = y and y2 = y'. Then the equation becomes y2' + 3y2 + 2y1 = 0. Combining with the definition y2 = y' we get the system y1 = y2 y2 = - 2y 1-3y 2 together with the initial conditions y 1(0) = 2 y 2(0) = 1. Suppose we want an approximation on the interval 0 t 20. Then we can input it into ode45 as follows. >> [tsol, ysol] = ode45(@(t,y) [y(2); - 2*y(1) - 3*y(2)], [0 20], [2 1]); Here is a plot of the approximate solution. >> plot(tsol, ysol(:,1)) Problem 6: Find a numerical approximation to the solution of the van der Pol equation y'' - (1 y 2 )y' + y = 0 that satisfies the initial conditions y(0) = 1 and y'(0) = 0. Plot y versus t for 0 t 30. 6
WILEY. Differential Equations with MATLAB (Third Edition) Brian R. Hunt Ronald L. Lipsman John E. Osborn Jonathan M. Rosenberg
Differential Equations with MATLAB (Third Edition) Updated for MATLAB 2011b (7.13), Simulink 7.8, and Symbolic Math Toolbox 5.7 Brian R. Hunt Ronald L. Lipsman John E. Osborn Jonathan M. Rosenberg All
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 informationMath Assignment 3 - Linear Algebra
Math 216 - Assignment 3 - Linear Algebra Due: Tuesday, March 27. Nothing accepted after Thursday, March 29. This is worth 15 points. 10% points off for being late. You may work by yourself or in pairs.
More informationLesson 14: Van der Pol Circuit and ode23s
Lesson 4: Van der Pol Circuit and ode3s 4. Applied Problem. A series LRC circuit when coupled via mutual inductance with a triode circuit can generate a sequence of pulsing currents that have very rapid
More informationAPPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015.
APPM 23: Final Exam :3am :pm, May, 25. ON THE FRONT OF YOUR BLUEBOOK write: ) your name, 2) your student ID number, 3) lecture section, 4) your instructor s name, and 5) a grading table for eight questions.
More informationFall 2003 Math 308/ Numerical Methods 3.6, 3.7, 5.3 Runge-Kutta Methods Mon, 27/Oct c 2003, Art Belmonte
Fall Math 8/ Numerical Methods.6,.7,. Runge-Kutta Methods Mon, 7/Oct c, Art Belmonte Summary Geometrical idea Runge-Kutta methods numerically approximate the solution of y = f (t, y), y(a) = y by using
More informationLinear Systems of ODE: Nullclines, Eigenvector lines and trajectories
Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 2013 Outline
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 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 informationMath 216 First Midterm 19 October, 2017
Math 6 First Midterm 9 October, 7 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 informationMath 2250 Maple Project 2, March Tacoma Narrows. CLASSTIME VERSION A-E, F-K, L-R, S-Z
Math 2250 Maple Project 2, March 2004. Tacoma Narrows. NAME CLASSTIME VERSION A-E, F-K, L-R, S-Z Circle the version - see problem 2.1. There are six (6) problems in this project. Please answer the questions
More informationODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class
ODE Homework Due Wed. 9 August 2009; At the beginning of the class. (a) Solve Lẏ + Ry = E sin(ωt) with y(0) = k () L, R, E, ω are positive constants. (b) What is the limit of the solution as ω 0? (c) Is
More informationAMS 27L LAB #6 Winter 2009
AMS 27L LAB #6 Winter 2009 Symbolically Solving Differential Equations Objectives: 1. To learn about the MATLAB Symbolic Solver 2. To expand knowledge of solutions to Diff-EQs 1 Symbolically Solving Differential
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 information2015 School Year. Graduate School Entrance Examination Problem Booklet. Mathematics
2015 School Year Graduate School Entrance Examination Problem Booklet Mathematics Examination Time: 10:00 to 12:30 Instructions 1. Do not open this problem booklet until the start of the examination is
More informationMath 216 First Midterm 18 October, 2018
Math 16 First Midterm 18 October, 018 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 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 informationApplied Calculus. Review Problems for the Final Exam
Math135 Study Guide 1 Math 131/135/194, Fall 2004 Applied Calculus Daniel Kaplan, Macalester College Review Problems for the Final Exam Problem 1../DE/102b.tex Problem 3../DE/107.tex Consider the pair
More information(Linear equations) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Linear equations) Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots
More information2 Solving Ordinary Differential Equations Using MATLAB
Penn State Erie, The Behrend College School of Engineering E E 383 Signals and Control Lab Spring 2008 Lab 3 System Responses January 31, 2008 Due: February 7, 2008 Number of Lab Periods: 1 1 Objective
More informationSYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES
SYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES by Mark H. Holmes Yuklun Au J. W. Stayman Department of Mathematical Sciences Rensselaer Polytechnic Institute, Troy, NY, 12180 Abstract
More informationHOMEWORK 3: Phase portraits for Mathmatical Models, Analysis and Simulation, Fall 2010 Report due Mon Oct 4, Maximum score 7.0 pts.
HOMEWORK 3: Phase portraits for Mathmatical Models, Analysis and Simulation, Fall 2010 Report due Mon Oct 4, 2010. Maximum score 7.0 pts. Three problems are to be solved in this homework assignment. The
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 information1 The relation between a second order linear ode and a system of two rst order linear odes
Math 1280 Spring, 2010 1 The relation between a second order linear ode and a system of two rst order linear odes In Chapter 3 of the text you learn to solve some second order linear ode's, such as x 00
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 informationChapter 1: Introduction
Chapter 1: Introduction Definition: A differential equation is an equation involving the derivative of a function. If the function depends on a single variable, then only ordinary derivatives appear 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 informationLS.5 Theory of Linear Systems
LS.5 Theory of Linear Systems 1. General linear ODE systems and independent solutions. We have studied the homogeneous system of ODE s with constant coefficients, (1) x = Ax, where A is an n n matrix of
More informationMatlab Course. Anna Kristine Wåhlin. Department of Geophysics, University of Oslo. January Matlab Course p.1/??
Matlab Course Anna Kristine Wåhlin Department of Geophysics, University of Oslo January 2003 Matlab Course p.1/?? Numerical estimate of the derivative An estimate of the time derivative of dataset at time
More informationFinal Exam Sample Problems, Math 246, Spring 2018
Final Exam Sample Problems, Math 246, Spring 2018 1) Consider the differential equation dy dt = 9 y2 )y 2. a) Find all of its stationary points and classify their stability. b) Sketch its phase-line portrait
More informationLinear Systems of ODE: Nullclines, Eigenvector lines and trajectories
Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 203 Outline
More informationSolutions to Math 53 First Exam April 20, 2010
Solutions to Math 53 First Exam April 0, 00. (5 points) Match the direction fields below with their differential equations. Also indicate which two equations do not have matches. No justification is necessary.
More informationLecture 8: Calculus and Differential Equations
Lecture 8: Calculus and Differential Equations Dr. Mohammed Hawa Electrical Engineering Department University of Jordan EE201: Computer Applications. See Textbook Chapter 9. Numerical Methods MATLAB provides
More informationLecture 8: Calculus and Differential Equations
Lecture 8: Calculus and Differential Equations Dr. Mohammed Hawa Electrical Engineering Department University of Jordan EE21: Computer Applications. See Textbook Chapter 9. Numerical Methods MATLAB provides
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 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 informationECE2111 Signals and Systems UMD, Spring 2013 Experiment 1: Representation and manipulation of basic signals in MATLAB
ECE2111 Signals and Systems UMD, Spring 2013 Experiment 1: Representation and manipulation of basic signals in MATLAB MATLAB is a tool for doing numerical computations with matrices and vectors. It can
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 informationMath 266: Phase Plane Portrait
Math 266: Phase Plane Portrait Long Jin Purdue, Spring 2018 Review: Phase line for an autonomous equation For a single autonomous equation y = f (y) we used a phase line to illustrate the equilibrium solutions
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 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 information(1 + 2y)y = x. ( x. The right-hand side is a standard integral, so in the end we have the implicit solution. y(x) + y 2 (x) = x2 2 +C.
Midterm 1 33B-1 015 October 1 Find the exact solution of the initial value problem. Indicate the interval of existence. y = x, y( 1) = 0. 1 + y Solution. We observe that the equation is separable, and
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 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 informationMATH 251 Examination II July 28, Name: Student Number: Section:
MATH 251 Examination II July 28, 2008 Name: Student Number: Section: This exam has 9 questions for a total of 100 points. In order to obtain full credit for partial credit problems, all work must be shown.
More informationLecture 7 Symbolic Computations
Lecture 7 Symbolic Computations The focus of this course is on numerical computations, i.e. calculations, usually approximations, with floating point numbers. However, Matlab can also do symbolic computations,
More informationAPPM 2460 CHAOTIC DYNAMICS
APPM 2460 CHAOTIC DYNAMICS 1. Introduction Today we ll continue our exploration of dynamical systems, focusing in particular upon systems who exhibit a type of strange behavior known as chaos. We will
More informationSIGNALS AND LINEAR SYSTEMS LABORATORY EELE Experiment (2) Introduction to MATLAB - Part (2) Prepared by:
The Islamic University of Gaza Faculty of Engineering Electrical Engineering Department SIGNALS AND LINEAR SYSTEMS LABORATORY EELE 110 Experiment () Introduction to MATLAB - Part () Prepared by: Eng. Mohammed
More informationNonconstant Coefficients
Chapter 7 Nonconstant Coefficients We return to second-order linear ODEs, but with nonconstant coefficients. That is, we consider (7.1) y + p(t)y + q(t)y = 0, with not both p(t) and q(t) constant. The
More informationMath 331 Homework Assignment Chapter 7 Page 1 of 9
Math Homework Assignment Chapter 7 Page of 9 Instructions: Please make sure to demonstrate every step in your calculations. Return your answers including this homework sheet back to the instructor as a
More informationExtreme Values and Positive/ Negative Definite Matrix Conditions
Extreme Values and Positive/ Negative Definite Matrix Conditions James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 8, 016 Outline 1
More information3.3. SYSTEMS OF ODES 1. y 0 " 2y" y 0 + 2y = x1. x2 x3. x = y(t) = c 1 e t + c 2 e t + c 3 e 2t. _x = A x + f; x(0) = x 0.
.. SYSTEMS OF ODES. Systems of ODEs MATH 94 FALL 98 PRELIM # 94FA8PQ.tex.. a) Convert the third order dierential equation into a rst oder system _x = A x, with y " y" y + y = x = @ x x x b) The equation
More informationMath 310 Introduction to Ordinary Differential Equations Final Examination August 9, Instructor: John Stockie
Make sure this exam has 15 pages. Math 310 Introduction to Ordinary Differential Equations inal Examination August 9, 2006 Instructor: John Stockie Name: (Please Print) Student Number: Special Instructions
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 informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationExample 1 (Characteristic Equation, Eigenvalue, and Eigenvector)
Matlab Lab 3 Example 1 (Characteristic Equation, Eigenvalue, and Eigenvector) A polynomial equation is uniquely determined by the coefficients of the monomial terms. For example, the quadratic equation
More informationPolytechnic Institute of NYU MA 2132 Final Practice Answers Fall 2012
Polytechnic Institute of NYU MA Final Practice Answers Fall Studying from past or sample exams is NOT recommended. If you do, it should be only AFTER you know how to do all of the homework and worksheet
More information(Mathematical Operations with Arrays) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Mathematical Operations with Arrays) Contents Getting Started Matrices Creating Arrays Linear equations Mathematical Operations with Arrays Using Script
More informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
More informationHigher Order Taylor Methods
Higher Order Taylor Methods Marcelo Julio Alvisio & Lisa Marie Danz May 6, 2007 Introduction Differential equations are one of the building blocks in science or engineering. Scientists aim to obtain numerical
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 informationSolutions of the Sample Problems for the First In-Class Exam Math 246, Fall 2017, Professor David Levermore
Solutions of the Sample Problems for the First In-Class Exam Math 246, Fall 207, Professor David Levermore () (a) Give the integral being evaluated by the following Matlab command. int( x/(+xˆ4), x,0,inf)
More informationSolving systems of first order equations with ode Systems of first order differential equations.
A M S 20 MA TLA B NO T E S U C S C Solving systems of first order equations with ode45 c 2015, Yonatan Katznelson The M A T L A B numerical solver, ode45 is designed to work with first order differential
More informationFINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.
FINAL EXAM MATH0 Theory of Ordinary Differential Equations There are 5 problems on 2 pages. Spring 2009. 25 points Consider the linear plane autonomous system x + y x y. Find a fundamental matrix of the
More informationmultistep methods Last modified: November 28, 2017 Recall that we are interested in the numerical solution of the initial value problem (IVP):
MATH 351 Fall 217 multistep methods http://www.phys.uconn.edu/ rozman/courses/m351_17f/ Last modified: November 28, 217 Recall that we are interested in the numerical solution of the initial value problem
More informationSystems of Linear ODEs
P a g e 1 Systems of Linear ODEs Systems of ordinary differential equations can be solved in much the same way as discrete dynamical systems if the differential equations are linear. We will focus here
More informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
More informationSOLUTIONS to Exercises from Optimization
SOLUTIONS to Exercises from Optimization. Use the bisection method to find the root correct to 6 decimal places: 3x 3 + x 2 = x + 5 SOLUTION: For the root finding algorithm, we need to rewrite the equation
More informationFIRST-ORDER ORDINARY DIFFERENTIAL EQUATIONS II: Graphical and Numerical Methods David Levermore Department of Mathematics University of Maryland
FIRST-ORDER ORDINARY DIFFERENTIAL EQUATIONS II: Graphical and Numerical Methods David Levermore Department of Mathematics University of Maryland 9 January 0 Because the presentation of this material in
More informationMath 5BI: Problem Set 6 Gradient dynamical systems
Math 5BI: Problem Set 6 Gradient dynamical systems April 25, 2007 Recall that if f(x) = f(x 1, x 2,..., x n ) is a smooth function of n variables, the gradient of f is the vector field f(x) = ( f)(x 1,
More informationu n 2 4 u n 36 u n 1, n 1.
Exercise 1 Let (u n ) be the sequence defined by Set v n = u n 1 x+ u n and f (x) = 4 x. 1. Solve the equations f (x) = 1 and f (x) =. u 0 = 0, n Z +, u n+1 = u n + 4 u n.. Prove that if u n < 1, then
More informationSection 5. Graphing Systems
Section 5. Graphing Systems 5A. The Phase Plane 5A-1. Find the critical points of each of the following non-linear autonomous systems. x = x 2 y 2 x = 1 x + y a) b) y = x xy y = y + 2x 2 5A-2. Write each
More informationMath 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv
Math 1270 Honors ODE I Fall, 2008 Class notes # 1 We have learned how to study nonlinear systems x 0 = F (x; y) y 0 = G (x; y) (1) by linearizing around equilibrium points. If (x 0 ; y 0 ) is an equilibrium
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 informationMathematical Computing
IMT2b2β Department of Mathematics University of Ruhuna A.W.L. Pubudu Thilan Differential Equations Types of Differential Equations Differential equations can basically be classified as ordinary differential
More informationYou may use a calculator, but you must show all your work in order to receive credit.
Math 2410-010/015 Exam II April 7 th, 2017 Name: Instructions: Key Answer each question to the best of your ability. All answers must be written clearly. Be sure to erase or cross out any work that you
More informationMECH : a Primer for Matlab s ode suite of functions
Objectives MECH 4-563: a Primer for Matlab s ode suite of functions. Review the fundamentals of initial value problems and why numerical integration methods are needed.. Introduce the ode suite of numerical
More informationDON T PANIC! If you get stuck, take a deep breath and go on to the next question. Come back to the question you left if you have time at the end.
Math 307A, Midterm 1 Spring 2013 Name: Instructions. DON T PANIC! If you get stuck, take a deep breath and go on to the next question. Come back to the question you left if you have time at the end. There
More informationThe University of British Columbia Final Examination - April 20, 2009 Mathematics 152 All Sections. Closed book examination. No calculators.
The University of British Columbia Final Examination - April 20, 2009 Mathematics 152 All Sections Closed book examination. No calculators. Time: 2.5 hours Last Name First Signature Student Number Section
More informationA A x i x j i j (i, j) (j, i) Let. Compute the value of for and
7.2 - Quadratic Forms quadratic form on is a function defined on whose value at a vector in can be computed by an expression of the form, where is an symmetric matrix. The matrix R n Q R n x R n Q(x) =
More informationSolution: (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 informationMATH 100 Introduction to the Profession
MATH 100 Introduction to the Profession Differential Equations in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2012 fasshauer@iit.edu MATH 100 ITP 1 What
More informationFirst-Order ODEs. Chapter Separable Equations. We consider in this chapter differential equations of the form dy (1.1)
Chapter 1 First-Order ODEs We consider in this chapter differential equations of the form dy (1.1) = F (t, y), where F (t, y) is a known smooth function. We wish to solve for y(t). Equation (1.1) is called
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 informationLaboratory 1. Solving differential equations with nonzero initial conditions
Laboratory 1 Solving differential equations with nonzero initial conditions 1. Purpose of the exercise: - learning symbolic and numerical methods of differential equations solving with MATLAB - using Simulink
More informationMidterm Exam 2. October 27, Time: 1 hour, 30 minutes. San Francisco State University ECON 715 Fall Name:
San Francisco State University Michael Bar ECON 75 Fall 206 Midterm Exam 2 October 27, 206 Time: hour, 30 minutes Name: Instructions:. One double-sided sheet with any content is allowed. 2. Calculators
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationFORTRAN 77 Lesson 7. Reese Haywood Ayan Paul
FORTRAN 77 Lesson 7 Reese Haywood Ayan Paul 1 Numerical Integration A general all purpose Numerical Integration technique is the Trapezoidal Method. For most functions, this method is fast and accurate.
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 informationProperties of Linear Transformations from R n to R m
Properties of Linear Transformations from R n to R m MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Topic Overview Relationship between the properties of a matrix transformation
More informationSection 2.4: Add and Subtract Rational Expressions
CHAPTER Section.: Add and Subtract Rational Expressions Section.: Add and Subtract Rational Expressions Objective: Add and subtract rational expressions with like and different denominators. You will recall
More information1. Type your solutions. This homework is mainly a programming assignment.
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5243 INTRODUCTION TO CYBER-PHYSICAL SYSTEMS H O M E W O R K S # 6 + 7 Ahmad F. Taha October 22, 2015 READ Homework Instructions: 1. Type your solutions. This homework
More informationStability of critical points in Linear Systems of Ordinary Differential Equations (SODE)
Stability of critical points in Linear Systems of Ordinary Differential Equations (SODE) In this chapter Mathematica will be used to study the stability of the critical points of twin equation linear SODE.
More informationMATH 44041/64041 Applied Dynamical Systems, 2018
MATH 4441/6441 Applied Dynamical Systems, 218 Answers to Assessed Coursework 1 1 Method based on series expansion of the matrix exponential Let the coefficient matrix of α the linear system be A, then
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 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 informationComputer lab for MAN460
Computer lab for MAN460 (version 20th April 2006, corrected 20 May) Prerequisites Matlab is rather user friendly, and to do the first exercises, it is enough to write an m-file consisting of two lines,
More informationChapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.
More informationSolution. It is evaluating the definite integral r 1 + r 4 dr. where you can replace r by any other variable.
Solutions of Sample Problems for First In-Class Exam Math 246, Fall 202, Professor David Levermore () (a) Give the integral being evaluated by the following MATLAB command. int( x/(+xˆ4), x,0,inf) Solution.
More informationIntroduction to Ordinary Differential Equations with Mathematica
ALFRED GRAY MICHAEL MEZZINO MARKA. PINSKY Introduction to Ordinary Differential Equations with Mathematica An Integrated Multimedia Approach %JmT} Web-Enhanced Includes CD-ROM TABLE OF CONTENTS Preface
More information