Introduction to ODEs in Scilab
|
|
- Eugene Gilbert
- 6 years ago
- Views:
Transcription
1 Introduction to Aditya Sengupta Indian Institute of Technology Bombay April 15th, 2010, Smt. Indira Gandhi College of Engineering
2 Scilab can be used to model and simulate a variety of systems, such as: 1 Ordinary Differential Equations 2 Boundary Value Problems 3 Difference Equations 4 Differential Algebraic Equations We will deal with Ordinary Differential Equations in this talk.
3 Scilab can be used to model and simulate a variety of systems, such as: 1 Ordinary Differential Equations 2 Boundary Value Problems 3 Difference Equations 4 Differential Algebraic Equations We will deal with Ordinary Differential Equations in this talk.
4 We will do two things: 1 Model the ODE in a way Scilab can understand. 2 Solve the system for a given set of initial values.
5 Modeling the system We will model the system as a first-order equation: ẏ = f (t, y) Note: Scilab tools assume the differential equation to have been written as first order system. Some models are initially written in terms of higher-order derivatives, but they can always be rewritten as first-order systems by the introduction of additional variables
6 Modeling the system We will model the system as a first-order equation: ẏ = f (t, y) Note: Scilab tools assume the differential equation to have been written as first order system. Some models are initially written in terms of higher-order derivatives, but they can always be rewritten as first-order systems by the introduction of additional variables
7 Modeling the system We will model the system as a first-order equation: ẏ = f (t, y) Note: Scilab tools assume the differential equation to have been written as first order system. Some models are initially written in terms of higher-order derivatives, but they can always be rewritten as first-order systems by the introduction of additional variables
8 Example Let us consider the simple system: dx dt = sin(2t) We can model this system using this code: 1 f u n c t i o n dx = f ( t, x ) 2 dx = s i n (2 t ) ; 3 e n d f u n c t i o n
9 Example Let us consider the simple system: dx dt = sin(2t) We can model this system using this code: 1 f u n c t i o n dx = f ( t, x ) 2 dx = s i n (2 t ) ; 3 e n d f u n c t i o n
10 Solution We know that the solution is supposed to be x = 1 2 cos(2t) + c where c is a constant that depends on the initial value of the problem.
11 Depending on the initial value, the plot will look like this:
12 Solving an ODE in Scilab The simulation tool we will use for solving is the ode function The simplest calling sequence for ode is: y=ode(y0, t0, t, f) where y0 is the initial value at t0 and t contains the points in time at which the solution is to be determined. f is the function corresponding to ẏ = f (t, y)
13 For our example, we will take the initial value to be y0 = -0.5 at t0 = 0. Let us evaluate the ODE from t = 0:0.1:5. The code is: 1 t0 = 0 2 x0 = t = 0 : 0. 1 : 5 ; 4 x = ode ( x0, t0, t, f ) ; 5 p l o t 2 d ( t, x )
14 For our example, we will take the initial value to be y0 = -0.5 at t0 = 0. Let us evaluate the ODE from t = 0:0.1:5. The code is: 1 t0 = 0 2 x0 = t = 0 : 0. 1 : 5 ; 4 x = ode ( x0, t0, t, f ) ; 5 p l o t 2 d ( t, x )
15 You should get a graph that looks like this:
16 Higher Order Derivatives When we have ODEs formulated in terms of higher order derivatives, we need to rewrite them as first-order systems. We do this by using variables to fill in the intermediate order derivaties. For example, let us consider the system: d 2 y dt 2 = sin(2t) whose one solution we can easily guess to be y = (1/4)sin(2t)
17 Higher Order Derivatives When we have ODEs formulated in terms of higher order derivatives, we need to rewrite them as first-order systems. We do this by using variables to fill in the intermediate order derivaties. For example, let us consider the system: d 2 y dt 2 = sin(2t) whose one solution we can easily guess to be y = (1/4)sin(2t)
18 We convert the second order equation into two first order equations: dy/dt = z dz/dt = sin(2t) Therefore, we have the ode in the form: where dx and x are vectors: dx/dt = f (t, x) x = [z; sin(2t)] dx = [dy/dt; dz/dt] We then proceed to replace z, dy/dt, and dz/dt with vector components x(2), dx(1), and dx(2)
19 We convert the second order equation into two first order equations: dy/dt = z dz/dt = sin(2t) Therefore, we have the ode in the form: where dx and x are vectors: dx/dt = f (t, x) x = [z; sin(2t)] dx = [dy/dt; dz/dt] We then proceed to replace z, dy/dt, and dz/dt with vector components x(2), dx(1), and dx(2)
20 We convert the second order equation into two first order equations: dy/dt = z dz/dt = sin(2t) Therefore, we have the ode in the form: where dx and x are vectors: dx/dt = f (t, x) x = [z; sin(2t)] dx = [dy/dt; dz/dt] We then proceed to replace z, dy/dt, and dz/dt with vector components x(2), dx(1), and dx(2)
21 We model the system thus: 1 f u n c t i o n dx = f ( t, x ) 2 dx ( 1 ) = x ( 2 ) 3 dx ( 2 ) = s i n (2 t ) 4 e n d f u n c t i o n and simulate the ODE thus: 1 t = 0 : : 4 %pi ; 2 3 y=ode ( [ 0 ; 1/2], 0, t, f ) ; 4 // Note t h e i m p o r t a n c e o f g i v i n g c o r r e c t s t a r t i n g v a l u e s. Try to put a l t e r n a t e s t a r t i n g v a l u e s and s e e t h e d i f f e r e n c e. 5 6 p l o t 2 d ( t, [ y ( 1, : ) y ( 2, : ) ] ) 7 // The c u r v e i n b l a c k i s t h e f i n a l s o l u t i o n. The o t h e r c u r v e i s f o r i l l u s t r a t i o n to show t h e i n t e r m e d i a t e s t e p.
22 We model the system thus: 1 f u n c t i o n dx = f ( t, x ) 2 dx ( 1 ) = x ( 2 ) 3 dx ( 2 ) = s i n (2 t ) 4 e n d f u n c t i o n and simulate the ODE thus: 1 t = 0 : : 4 %pi ; 2 3 y=ode ( [ 0 ; 1/2], 0, t, f ) ; 4 // Note t h e i m p o r t a n c e o f g i v i n g c o r r e c t s t a r t i n g v a l u e s. Try to put a l t e r n a t e s t a r t i n g v a l u e s and s e e t h e d i f f e r e n c e. 5 6 p l o t 2 d ( t, [ y ( 1, : ) y ( 2, : ) ] ) 7 // The c u r v e i n b l a c k i s t h e f i n a l s o l u t i o n. The o t h e r c u r v e i s f o r i l l u s t r a t i o n to show t h e i n t e r m e d i a t e s t e p.
23 Root Finding Sometimes- we just want to simulate a differential equation up to the time that a specific event occurs. For example, an engine being revved until it reaches a particular speed- after which the gear is to be changed. For such circumstances, we need to define a quantity that signals the occurance of the event.
24 Root Finding Sometimes- we just want to simulate a differential equation up to the time that a specific event occurs. For example, an engine being revved until it reaches a particular speed- after which the gear is to be changed. For such circumstances, we need to define a quantity that signals the occurance of the event.
25 Root Finding Sometimes- we just want to simulate a differential equation up to the time that a specific event occurs. For example, an engine being revved until it reaches a particular speed- after which the gear is to be changed. For such circumstances, we need to define a quantity that signals the occurance of the event.
26 In Scilab we use the ode root function, which is called thus: [y, rd] = ode("root", y0, t0, t, f, ng, g) where g is a function that becomes zero valued when the constraining event occurs and ng is the size of g. rd is a vector that contains the stopping time as its first element.
27 In Scilab we use the ode root function, which is called thus: [y, rd] = ode("root", y0, t0, t, f, ng, g) where g is a function that becomes zero valued when the constraining event occurs and ng is the size of g. rd is a vector that contains the stopping time as its first element.
28 Example Let us consider the example of the engine that is revved. We wish to constrain the revving of the engine till it reaches a certain point. We build a first order approximation of an engine using the following code (call it engine.sci): 1 f u n c t i o n d_revs = engine ( t, revs ) 2 d_revs = %eˆ( revs ) 3 e n d f u n c t i o n We can simulate the behaviour of the engine when it is unconstrained using the following code: 1 e x e c engine. sci 2 revs = ode ( 0, 0, 0 : 0. 1 : 1 0, engine ) ; 3 p l o t 2 d ( 0 : 0. 1 : 1 0, revs )
29 Example Let us consider the example of the engine that is revved. We wish to constrain the revving of the engine till it reaches a certain point. We build a first order approximation of an engine using the following code (call it engine.sci): 1 f u n c t i o n d_revs = engine ( t, revs ) 2 d_revs = %eˆ( revs ) 3 e n d f u n c t i o n We can simulate the behaviour of the engine when it is unconstrained using the following code: 1 e x e c engine. sci 2 revs = ode ( 0, 0, 0 : 0. 1 : 1 0, engine ) ; 3 p l o t 2 d ( 0 : 0. 1 : 1 0, revs )
30 Example Let us consider the example of the engine that is revved. We wish to constrain the revving of the engine till it reaches a certain point. We build a first order approximation of an engine using the following code (call it engine.sci): 1 f u n c t i o n d_revs = engine ( t, revs ) 2 d_revs = %eˆ( revs ) 3 e n d f u n c t i o n We can simulate the behaviour of the engine when it is unconstrained using the following code: 1 e x e c engine. sci 2 revs = ode ( 0, 0, 0 : 0. 1 : 1 0, engine ) ; 3 p l o t 2 d ( 0 : 0. 1 : 1 0, revs )
31 We then write the constraining function (call it gearbox.sci): 1 f u n c t i o n stop = gearbox ( t, revs ) 2 stop = 1. 5 revs //We choose to s t o p t h e e n g i n e when t h e r e v s r e a c h t h e v a l u e 1. 5 ( You can choose any o t h e r v a l u e ) 3 e n d f u n c t i o n We then simulate the behaviour of the engine when it is constrained as above. 1 e x e c engine. sci 2 e x e c gearbox. sci 3 [ revs, stop_time ]= ode ( r o o t, 0, 0, 0 : 0. 1 : 1 0, engine, 1, gearbox ) ; 4 p l o t 2 d ( [ 0 : 0. 1 : stop_time ( 1 ), stop_time ( 1 ) ], revs )
32 We then write the constraining function (call it gearbox.sci): 1 f u n c t i o n stop = gearbox ( t, revs ) 2 stop = 1. 5 revs //We choose to s t o p t h e e n g i n e when t h e r e v s r e a c h t h e v a l u e 1. 5 ( You can choose any o t h e r v a l u e ) 3 e n d f u n c t i o n We then simulate the behaviour of the engine when it is constrained as above. 1 e x e c engine. sci 2 e x e c gearbox. sci 3 [ revs, stop_time ]= ode ( r o o t, 0, 0, 0 : 0. 1 : 1 0, engine, 1, gearbox ) ; 4 p l o t 2 d ( [ 0 : 0. 1 : stop_time ( 1 ), stop_time ( 1 ) ], revs )
33 Compare the two graphs- can you see where the simulation was halted in the second case?
34 Linear Systems Since they appear so often, there are special functions for modeling and simulating linear systems. For instance, you can create a linear system thus: 1 s = p o l y ( 0, s ) 2 sys = s y s l i n ( c, 1/( s+1) ) and simulate it thus: 1 t = 0 : 0. 1 : 1 0 ; 2 y = csim ( s t e p, t, sys ) ; 3 p l o t 2 d ( t, y ) ; 4 5 z = csim ( s i n (5 t ), t, sys ) ; 6 p l o t 2 d ( t, z ) ; 7 8 bode ( sys, , 100) ;
35 Linear Systems Since they appear so often, there are special functions for modeling and simulating linear systems. For instance, you can create a linear system thus: 1 s = p o l y ( 0, s ) 2 sys = s y s l i n ( c, 1/( s+1) ) and simulate it thus: 1 t = 0 : 0. 1 : 1 0 ; 2 y = csim ( s t e p, t, sys ) ; 3 p l o t 2 d ( t, y ) ; 4 5 z = csim ( s i n (5 t ), t, sys ) ; 6 p l o t 2 d ( t, z ) ; 7 8 bode ( sys, , 100) ;
36 Now try to: Model and simulate your own systems Use the help command to find more options
37 Thanks
Speed and Velocity: Recall from Calc 1: If f (t) gives the position of an object at time t, then. velocity at time t = f (t) speed at time t = f (t)
Speed and Velocity: Recall from Calc 1: If f (t) gives the position of an object at time t, then velocity at time t = f (t) speed at time t = f (t) Math 36-Multi (Sklensky) In-Class Work January 8, 013
More informationLinear Algebra and Optimization Using Scilab
U. Patil deepakp@ee.iitb.ac.in Indian Institute of Technology, Bombay November 2, 2009 System of Linear Equations Consider x 1 + x 2 + x 3 = 10 3x 1 + x 2 + 2x 3 = 5 x 1 + x 2 x 3 = 1. System of Linear
More informationComputer Problems for Methods of Solving Ordinary Differential Equations
Computer Problems for Methods of Solving Ordinary Differential Equations 1. Have a computer make a phase portrait for the system dx/dt = x + y, dy/dt = 2y. Clearly indicate critical points and separatrices.
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 informationHands-on Session: Linear Algebra and Ordinary Differential Equations with SCILAB
Hands-on Session: Linear Algebra and Ordinary Differential Equations with SCILAB Scilab and Its Applications to Global Optimization and Fractional Differential Equations SGGS IE & T, Nanded, April 23-25,
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 informationLinear Algebra, Optimization and Solving Ordinary Differential Equations Using Scilab
, and Solving Ordinary Differential Equations Using Scilab U. Patil deepakp@ee.iitb.ac.in Indian Institute of Technology, Bombay November 8, 2009, and Solving Ordinary Differential , and Solving Ordinary
More informationMathematics Engineering Calculus III Fall 13 Test #1
Mathematics 2153-02 Engineering Calculus III Fall 13 Test #1 Instructor: Dr. Alexandra Shlapentokh (1) Which of the following statements is always true? (a) If x = f(t), y = g(t) and f (1) = 0, then dy/dx(1)
More informationLinear algebra and differential equations (Math 54): Lecture 20
Linear algebra and differential equations (Math 54): Lecture 20 Vivek Shende April 7, 2016 Hello and welcome to class! Last time We started discussing differential equations. We found a complete set of
More information3.2 Systems of Two First Order Linear DE s
Agenda Section 3.2 Reminders Lab 1 write-up due 9/26 or 9/28 Lab 2 prelab due 9/26 or 9/28 WebHW due 9/29 Office hours Tues, Thurs 1-2 pm (5852 East Hall) MathLab office hour Sun 7-8 pm (MathLab) 3.2 Systems
More informationA quadratic expression is a mathematical expression that can be written in the form 2
118 CHAPTER Algebra.6 FACTORING AND THE QUADRATIC EQUATION Textbook Reference Section 5. CLAST OBJECTIVES Factor a quadratic expression Find the roots of a quadratic equation A quadratic expression is
More informationFirst Order Linear Ordinary Differential Equations
First Order Linear Ordinary Differential Equations The most general first order linear ODE is an equation of the form p t dy dt q t y t f t. 1 Herepqarecalledcoefficients f is referred to as the forcing
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 information5.2 LENGTHS OF CURVES & AREAS OF SURFACES OF REVOLUTION
5.2 Arc Length & Surface Area Contemporary Calculus 1 5.2 LENGTHS OF CURVES & AREAS OF SURFACES OF REVOLUTION This section introduces two additional geometric applications of integration: finding the length
More informationLecture 9. Systems of Two First Order Linear ODEs
Math 245 - Mathematics of Physics and Engineering I Lecture 9. Systems of Two First Order Linear ODEs January 30, 2012 Konstantin Zuev (USC) Math 245, Lecture 9 January 30, 2012 1 / 15 Agenda General Form
More informationSolving Differential Equations Using MATLAB
Solving Differential Equations Using MATLAB Abraham Asfaw aasfaw.student@manhattan.edu November 28, 2011 1 Introduction In this lecture, we will follow up on lecture 2 with a discussion of solutions to
More informationLecture 1. Scott Pauls 1 3/28/07. Dartmouth College. Math 23, Spring Scott Pauls. Administrivia. Today s material.
Lecture 1 1 1 Department of Mathematics Dartmouth College 3/28/07 Outline Course Overview http://www.math.dartmouth.edu/~m23s07 Matlab Ordinary differential equations Definition An ordinary differential
More information37. f(t) sin 2t cos 2t 38. f(t) cos 2 t. 39. f(t) sin(4t 5) 40.
28 CHAPTER 7 THE LAPLACE TRANSFORM EXERCISES 7 In Problems 8 use Definition 7 to find {f(t)} 2 3 4 5 6 7 8 9 f (t),, f (t) 4,, f (t) t,, f (t) 2t,, f (t) sin t,, f (t), cos t, t t t 2 t 2 t t t t t t t
More informationExam 1 Review SOLUTIONS
1. True or False (and give a short reason): Exam 1 Review SOLUTIONS (a) If the parametric curve x = f(t), y = g(t) satisfies g (1) = 0, then it has a horizontal tangent line when t = 1. FALSE: To make
More informationSubstitutions and by Parts, Area Between Curves. Goals: The Method of Substitution Areas Integration by Parts
Week #7: Substitutions and by Parts, Area Between Curves Goals: The Method of Substitution Areas Integration by Parts 1 Week 7 The Indefinite Integral The Fundamental Theorem of Calculus, b a f(x) dx =
More informationMA 102 Mathematics II Lecture Feb, 2015
MA 102 Mathematics II Lecture 1 20 Feb, 2015 Differential Equations An equation containing derivatives is called a differential equation. The origin of differential equations Many of the laws of nature
More informationMATH 2413 TEST ON CHAPTER 4 ANSWER ALL QUESTIONS. TIME 1.5 HRS.
MATH 1 TEST ON CHAPTER ANSWER ALL QUESTIONS. TIME 1. HRS. M1c Multiple Choice Identify the choice that best completes the statement or answers the question. 1. Use the summation formulas to rewrite the
More informationTable of contents. d 2 y dx 2, As the equation is linear, these quantities can only be involved in the following manner:
M ath 0 1 E S 1 W inter 0 1 0 Last Updated: January, 01 0 Solving Second Order Linear ODEs Disclaimer: This lecture note tries to provide an alternative approach to the material in Sections 4. 4. 7 and
More informationSOME EXAMPLES OF THE USE OF COMPUTER ALGEBRA SYSTEMS IN INTRODUCTORY COURSES
SOME EXAMPLES OF THE USE OF COMPUTER ALGEBRA SYSTEMS IN INTRODUCTORY COURSES Abi Fattahi Whittier College Department of Mathematics Whittier, CA 90605 afattahi@whittier.edu In this presentation, we will
More informationSection 1.6 Inverse Functions
0 Chapter 1 Section 1.6 Inverse Functions A fashion designer is travelling to Milan for a fashion show. He asks his assistant, Betty, what 7 degrees Fahrenheit is in Celsius, and after a quick search on
More informationMATH 307: Problem Set #3 Solutions
: Problem Set #3 Solutions Due on: May 3, 2015 Problem 1 Autonomous Equations Recall that an equilibrium solution of an autonomous equation is called stable if solutions lying on both sides of it tend
More informationParametric Equations, Function Composition and the Chain Rule: A Worksheet
Parametric Equations, Function Composition and the Chain Rule: A Worksheet Prof.Rebecca Goldin Oct. 8, 003 1 Parametric Equations We have seen that the graph of a function f(x) of one variable consists
More informationSection 1.4. Meaning of Slope for Equations, Graphs, and Tables
Section 1.4 Meaning of Slope for Equations, Graphs, and Tables Finding Slope from a Linear Equation Finding Slope from a Linear Equation Example Find the slope of the line Solution Create a table using
More informationEEL2216 Control Theory CT1: PID Controller Design
EEL6 Control Theory CT: PID Controller Design. Objectives (i) To design proportional-integral-derivative (PID) controller for closed loop control. (ii) To evaluate the performance of different controllers
More informationMAT187H1F Lec0101 Burbulla
Spring 2017 Second Order Linear Homogeneous Differential Equation DE: A(x) d 2 y dx 2 + B(x)dy dx + C(x)y = 0 This equation is called second order because it includes the second derivative of y; it is
More informationChemical Engineering Applications in Scilab
Chemical Engineering Applications in Scilab Prashant Dave Indian Institute of Technology Bombay (IIT Bombay) Introduction In Chemical Engineering the type of problems that occur are Modeling and Simulation
More informationRiemann Sums. Outline. James K. Peterson. September 15, Riemann Sums. Riemann Sums In MatLab
Riemann Sums James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 15, 2013 Outline Riemann Sums Riemann Sums In MatLab Abstract This
More informationOrdinary Differential Equations
Chapter 10 Ordinary Differential Equations 10.1 Introduction Relationship between rate of change of variables rather than variables themselves gives rise to differential equations. Mathematical formulation
More informationModeling Rarefaction and Shock waves
April 30, 2013 Inroduction Inroduction Rarefaction and shock waves are combinations of two wave fronts created from the initial disturbance of the medium. Inroduction Rarefaction and shock waves are combinations
More information(Refer Slide Time: 00:32)
Nonlinear Dynamical Systems Prof. Madhu. N. Belur and Prof. Harish. K. Pillai Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture - 12 Scilab simulation of Lotka Volterra
More informationSection 2.1 (First Order) Linear DEs; Method of Integrating Factors. General first order linear DEs Standard Form; y'(t) + p(t) y = g(t)
Section 2.1 (First Order) Linear DEs; Method of Integrating Factors Key Terms/Ideas: General first order linear DEs Standard Form; y'(t) + p(t) y = g(t) Integrating factor; a function μ(t) that transforms
More informationDIFFERENTIATION RULES
3 DIFFERENTIATION RULES DIFFERENTIATION RULES We have: Seen how to interpret derivatives as slopes and rates of change Seen how to estimate derivatives of functions given by tables of values Learned how
More informationOrdinary Differential Equations
Ordinary Differential Equations for Engineers and Scientists Gregg Waterman Oregon Institute of Technology c 2017 Gregg Waterman This work is licensed under the Creative Commons Attribution 4.0 International
More informationA. Incorrect! Apply the rational root test to determine if any rational roots exist.
College Algebra - Problem Drill 13: Zeros of Polynomial Functions No. 1 of 10 1. Determine which statement is true given f() = 3 + 4. A. f() is irreducible. B. f() has no real roots. C. There is a root
More information1.3 Subtraction of Rational Numbers
www.ck2.org Chapter. Real Numbers Operations.3 Subtraction of Rational Numbers Learning objectives Find additive inverses. Subtract rational numbers. Evaluate change using a variable expression. Solve
More informationO.K. But what if the chicken didn t have access to a teleporter.
The intermediate value theorem, and performing algebra on its. This is a dual topic lecture. : The Intermediate value theorem First we should remember what it means to be a continuous function: A function
More informationSelected Solutions: 3.5 (Undetermined Coefficients)
Selected Solutions: 3.5 (Undetermined Coefficients) In Exercises 1-10, we want to apply the ideas from the table to specific DEs, and solve for the coefficients as well. If you prefer, you might start
More informationLimits: How to approach them?
Limits: How to approach them? The purpose of this guide is to show you the many ways to solve it problems. These depend on many factors. The best way to do this is by working out a few eamples. In particular,
More informationSection 2.1 Differential Equation and Solutions
Section 2.1 Differential Equation and Solutions Key Terms: Ordinary Differential Equation (ODE) Independent Variable Order of a DE Partial Differential Equation (PDE) Normal Form Solution General Solution
More informationTime Response of Systems
Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =
More informationReview for Ma 221 Final Exam
Review for Ma 22 Final Exam The Ma 22 Final Exam from December 995.a) Solve the initial value problem 2xcosy 3x2 y dx x 3 x 2 sin y y dy 0 y 0 2 The equation is first order, for which we have techniques
More informationAdvanced Calculus Questions
Advanced Calculus Questions What is here? This is a(n evolving) collection of challenging calculus problems. Be warned - some of these questions will go beyond the scope of this course. Particularly difficult
More informationOrdinary Differential Equations. Monday, October 10, 11
Ordinary Differential Equations Monday, October 10, 11 Problems involving ODEs can always be reduced to a set of first order differential equations. For example, By introducing a new variable z, this can
More informationa k 0, then k + 1 = 2 lim 1 + 1
Math 7 - Midterm - Form A - Page From the desk of C. Davis Buenger. https://people.math.osu.edu/buenger.8/ Problem a) [3 pts] If lim a k = then a k converges. False: The divergence test states that if
More informationAPPLICATIONS OF INTEGRATION
6 APPLICATIONS OF INTEGRATION APPLICATIONS OF INTEGRATION 6.5 Average Value of a Function In this section, we will learn about: Applying integration to find out the average value of a function. AVERAGE
More informationRegression I - the least squares line
Regression I - the least squares line The difference between correlation and regression. Correlation describes the relationship between two variables, where neither variable is independent or used to predict.
More informationSection 7.2 Addition and Subtraction Identities. In this section, we begin expanding our repertoire of trigonometric identities.
Section 7. Addition and Subtraction Identities 47 Section 7. Addition and Subtraction Identities In this section, we begin expanding our repertoire of trigonometric identities. Identities The sum and difference
More informationPhysics 202 Laboratory 3. Root-Finding 1. Laboratory 3. Physics 202 Laboratory
Physics 202 Laboratory 3 Root-Finding 1 Laboratory 3 Physics 202 Laboratory The fundamental question answered by this week s lab work will be: Given a function F (x), find some/all of the values {x i }
More informationx(t) = t[u(t 1) u(t 2)] + 1[u(t 2) u(t 3)]
ECE30 Summer II, 2006 Exam, Blue Version July 2, 2006 Name: Solution Score: 00/00 You must show all of your work for full credit. Calculators may NOT be used.. (5 points) x(t) = tu(t ) + ( t)u(t 2) u(t
More informationUse of Differential Equations In Modeling and Simulation of CSTR
Use of Differential Equations In Modeling and Simulation of CSTR JIRI VOJTESEK, PETR DOSTAL Department of Process Control, Faculty of Applied Informatics Tomas Bata University in Zlin nám. T. G. Masaryka
More informationA Glimpse at Scipy FOSSEE. June Abstract This document shows a glimpse of the features of Scipy that will be explored during this course.
A Glimpse at Scipy FOSSEE June 010 Abstract This document shows a glimpse of the features of Scipy that will be explored during this course. 1 Introduction SciPy is open-source software for mathematics,
More informationOften, in this class, we will analyze a closed-loop feedback control system, and end up with an equation of the form
ME 32, Spring 25, UC Berkeley, A. Packard 55 7 Review of SLODEs Throughout this section, if y denotes a function (of time, say), then y [k or y (k) denotes the k th derivative of the function y, y [k =
More informationWeek 7: Binary Outcomes (Scott Long Chapter 3 Part 2)
Week 7: (Scott Long Chapter 3 Part 2) Tsun-Feng Chiang* *School of Economics, Henan University, Kaifeng, China April 29, 2014 1 / 38 ML Estimation for Probit and Logit ML Estimation for Probit and Logit
More informationInfinite series, improper integrals, and Taylor series
Chapter 2 Infinite series, improper integrals, and Taylor series 2. Introduction to series In studying calculus, we have explored a variety of functions. Among the most basic are polynomials, i.e. functions
More information9.4 CALCULUS AND PARAMETRIC EQUATIONS
9.4 Calculus with Parametric Equations Contemporary Calculus 1 9.4 CALCULUS AND PARAMETRIC EQUATIONS The previous section discussed parametric equations, their graphs, and some of their uses for visualizing
More informationWeek 2 Techniques of Integration
Week Techniques of Integration Richard Earl Mathematical Institute, Oxford, OX LB, October Abstract Integration by Parts. Substitution. Rational Functions. Partial Fractions. Trigonometric Substitutions.
More information(Refer Slide Time: 02:11 to 04:19)
Digital Signal Processing Prof. S. C. Dutta Roy Department of Electrical Engineering Indian Institute of Technology, Delhi Lecture - 24 Analog Chebyshev LPF Design This is the 24 th lecture on DSP and
More information1 Differential. Equations. A differential equation is any equation that involves a derivative. For example, Newton s second law F ma
1 Differential Equations b The strange attractor for a Sprott system consisting of three quadratic differential equations. 1 A differential equation is any equation that involves a derivative. For example,
More information3.5: Euler s Method (cont d) and 3.7: Population Modeling
3.5: Euler s Method (cont d) and 3.7: Population Modeling Mathematics 3 Lecture 19 Dartmouth College February 15, 2010 Typeset by FoilTEX Example 1 Let s consider a very simple first-order IVP: dy dx =
More informationComputational Fluid Dynamics Prof. Dr. SumanChakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computational Fluid Dynamics Prof. Dr. SumanChakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. #11 Fundamentals of Discretization: Finite Difference
More informationMTH4100 Calculus I. Week 6 (Thomas Calculus Sections 3.5 to 4.2) Rainer Klages. School of Mathematical Sciences Queen Mary, University of London
MTH4100 Calculus I Week 6 (Thomas Calculus Sections 3.5 to 4.2) Rainer Klages School of Mathematical Sciences Queen Mary, University of London Autumn 2008 R. Klages (QMUL) MTH4100 Calculus 1 Week 6 1 /
More informationNumerical techniques to solve equations
Programming for Applications in Geomatics, Physical Geography and Ecosystem Science (NGEN13) Numerical techniques to solve equations vaughan.phillips@nateko.lu.se Vaughan Phillips Associate Professor,
More informationLecture IX. Definition 1 A non-singular Sturm 1 -Liouville 2 problem consists of a second order linear differential equation of the form.
Lecture IX Abstract When solving PDEs it is often necessary to represent the solution in terms of a series of orthogonal functions. One way to obtain an orthogonal family of functions is by solving a particular
More informationAdvanced Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras
Advanced Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras Lecture - 3 Simplex Method for Bounded Variables We discuss the simplex algorithm
More informationDIFFERENTIATION RULES
3 DIFFERENTIATION RULES DIFFERENTIATION RULES Before starting this section, you might need to review the trigonometric functions. DIFFERENTIATION RULES In particular, it is important to remember that,
More informationStandard Form of Conics
When we teach linear equations in Algebra1, we teach the simplest linear function (the mother function) as y = x. We then usually lead to the understanding of the effects of the slope and the y-intercept
More informationTAYLOR POLYNOMIALS DARYL DEFORD
TAYLOR POLYNOMIALS DARYL DEFORD 1. Introduction We have seen in class that Taylor polynomials provide us with a valuable tool for approximating many different types of functions. However, in order to really
More informationModelling and Mathematical Methods in Process and Chemical Engineering
FS 07 February 0, 07 Modelling and Mathematical Methods in Process and Chemical Engineering Solution Series. Systems of linear algebraic equations: adjoint, determinant, inverse The adjoint of a (square)
More informationMATH 307 Introduction to Differential Equations Autumn 2017 Midterm Exam Monday November
MATH 307 Introduction to Differential Equations Autumn 2017 Midterm Exam Monday November 6 2017 Name: Student ID Number: I understand it is against the rules to cheat or engage in other academic misconduct
More informationCURVILINEAR MOTION: GENERAL & RECTANGULAR COMPONENTS
CURVILINEAR MOTION: GENERAL & RECTANGULAR COMPONENTS Today s Objectives: Students will be able to: 1. Describe the motion of a particle traveling along a curved path. 2. Relate kinematic quantities in
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 informationParameter Estimation of Mathematical Models Described by Differential Equations
Parameter Estimation p.1/12 Parameter Estimation of Mathematical Models Described by Differential Equations Hossein Zivaripiran Department of Computer Science University of Toronto Outline Parameter Estimation
More informationComputational Techniques Prof. Dr. Niket Kaisare Department of Chemical Engineering Indian Institute of Technology, Madras
Computational Techniques Prof. Dr. Niket Kaisare Department of Chemical Engineering Indian Institute of Technology, Madras Module No. # 07 Lecture No. # 05 Ordinary Differential Equations (Refer Slide
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 informationUnit VIII Matrices, Eigenvalues and Eigenvectors
Unit VIII Matrices, Eigenvalues and Eigenvectors These materials were developed for use at Grinnell College and neither Grinnell College nor the author, Mark Schneider, assume any responsibility for their
More informationCURVILINEAR MOTION: RECTANGULAR COMPONENTS (Sections )
CURVILINEAR MOTION: RECTANGULAR COMPONENTS (Sections 12.4-12.5) Today s Objectives: Students will be able to: a) Describe the motion of a particle traveling along a curved path. b) Relate kinematic quantities
More informationOrdinary Differential Equations (ODEs)
Chapter 13 Ordinary Differential Equations (ODEs) We briefly review how to solve some of the most standard ODEs. 13.1 First Order Equations 13.1.1 Separable Equations A first-order ordinary differential
More informationPredator - Prey Model Trajectories and the nonlinear conservation law
Predator - Prey Model Trajectories and the nonlinear conservation law James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 28, 2013 Outline
More informationMath 111: Calculus. David Perkinson
Math : Calculus David Perkinson Fall 207 Contents Week, Monday: Introduction: derivatives, integrals, and the fundamental theorem. 5 Week, Wednesday: Average speed, instantaneous speed. Definition of the
More informationES 111 Mathematical Methods in the Earth Sciences Lecture Outline 15 - Tues 20th Nov 2018 First and Higher Order Differential Equations
ES 111 Mathematical Methods in the Earth Sciences Lecture Outline 15 - Tues 20th Nov 2018 First and Higher Order Differential Equations Integrating Factor Here is a powerful technique which will work (only!)
More informationMath 3B: Lecture 14. Noah White. February 13, 2016
Math 3B: Lecture 14 Noah White February 13, 2016 Last time Accumulated change problems Last time Accumulated change problems Adding up a value if it is changing over time Last time Accumulated change problems
More informationNormal Modes, Wave Motion and the Wave Equation Hilary Term 2011 Lecturer: F Hautmann
Normal Modes, Wave Motion and the Wave Equation Hilary Term 2011 Lecturer: F Hautmann Part A: Normal modes ( 4 lectures) Part B: Waves ( 8 lectures) printed lecture notes slides will be posted on lecture
More informationThis file is /conf/snippets/setheader.pg you can use it as a model for creating files which introduce each problem set.
Victoria Howle WeBWorK assignment number WW06 is due : 03/30/2012 at 08:07pm CDT. The (* replace with url for the course home page *) for the course contains the syllabus, grading policy and other information.
More informationDynamics of Ocean Structures Prof. Dr. Srinivasan Chandrasekaran Department of Ocean Engineering Indian Institute of Technology, Madras
Dynamics of Ocean Structures Prof. Dr. Srinivasan Chandrasekaran Department of Ocean Engineering Indian Institute of Technology, Madras Lecture 25 Continuous System In the last class, in this, we will
More informationthen the substitution z = ax + by + c reduces this equation to the separable one.
7 Substitutions II Some of the topics in this lecture are optional and will not be tested at the exams. However, for a curious student it should be useful to learn a few extra things about ordinary differential
More informationSecond-Order Homogeneous Linear Equations with Constant Coefficients
15 Second-Order Homogeneous Linear Equations with Constant Coefficients A very important class of second-order homogeneous linear equations consists of those with constant coefficients; that is, those
More informationProperties of Real Numbers
Properties of Real Numbers Essential Understanding. Relationships that are always true for real numbers are called properties, which are rules used to rewrite and compare expressions. Two algebraic expressions
More informationExercise 4) Newton's law of cooling is a model for how objects are heated or cooled by the temperature of an ambient medium surrounding them.
Exercise 4) Newton's law of cooling is a model for how objects are heated or cooled by the temperature of an ambient medium surrounding them. In this model, the bo temperature T = T t changes at a rate
More information1.4 Techniques of Integration
.4 Techniques of Integration Recall the following strategy for evaluating definite integrals, which arose from the Fundamental Theorem of Calculus (see Section.3). To calculate b a f(x) dx. Find a function
More informationPuzzle 1 Puzzle 2 Puzzle 3 Puzzle 4 Puzzle 5 /10 /10 /10 /10 /10
MATH-65 Puzzle Collection 6 Nov 8 :pm-:pm Name:... 3 :pm Wumaier :pm Njus 5 :pm Wumaier 6 :pm Njus 7 :pm Wumaier 8 :pm Njus This puzzle collection is closed book and closed notes. NO calculators are allowed
More informationSpring 2015, MA 252, Calculus II, Final Exam Preview Solutions
Spring 5, MA 5, Calculus II, Final Exam Preview Solutions I will put the following formulas on the front of the final exam, to speed up certain problems. You do not need to put them on your index card,
More informationCLEP College Algebra - Problem Drill 21: Solving and Graphing Linear Inequalities
CLEP College Algebra - Problem Drill 21: Solving and Graphing Linear Inequalities No. 1 of 10 1. Which inequality represents the statement three more than seven times a real number is greater than or equal
More informationMath 23: Differential Equations (Winter 2017) Midterm Exam Solutions
Math 3: Differential Equations (Winter 017) Midterm Exam Solutions 1. [0 points] or FALSE? You do not need to justify your answer. (a) [3 points] Critical points or equilibrium points for a first order
More information(Refer Slide Time: 01:17)
Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay Lecture No. 7 Heat Conduction 4 Today we are going to look at some one dimensional
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 information