Introduction to computational modelling
|
|
- Shannon Scott
- 5 years ago
- Views:
Transcription
1 Introduction to computational modelling Lecture 6 : Differential equations Physics: Pendulum Algorithm: Taylor s method Programming: Code explained Instructor : Cedric Weber Course : 4CCP1000
2 Schedule Class/Week Chapter Topic Milestones 1 Monte Carlo UNIX system / Fortran 2 Monte Carlo Fibonacci sequence 3 Monte Carlo Random variables 4 Monte Carlo Central Limit Theorem 5 Monte Carlo Monte Carlo integration Milestone 1 6 Differential equations The Pendulum 7 Differential equations A Quantum Particle in a box 8 Differential equations The Tacoma bridge Milestone 2 9 Linear Algebra System of equations 10 Linear Algebra Matrix operations Milestone 3 2
3 What did you learn last time? 1. Performing an integral with Monte Carlo ý 2. Using a subroutine with the call statement ý 3. Using a random generator to produce random values ý 4. Convert integers to real numbers ý 5. Generate and plot two columns files with xmgrace ý 3
4 Where are we going? 100% * Lecture 1-5 : * Statistics and integrals 90% 80% 70% 60% * Lecture 6 : * The pendulum : we will combine functions and arrays 50% 40% 30% 20% 10% Milestone 2 : Solving differential equations 0% week 1 week 2 week 3 week 4 week 5 week 6 week 7 week 8 week 9 week 10 4
5 5
6 Today s experiment Pendulum! Step 1 : Let s write down Forces R ² At time t=0, pendulum is dropped from angle q* with initial zero velocity ² Length R, gravity constant g=9.81 ² Forces: ² e n : ² e t : normal direction tangential direction 6
7 Equation of motion for a circular motion Step 2 : Let s define the acceleration d ² Newton : F = M a = M v(t) [a=acceleration, v=velocity, M=mass] dt ² point moving on a circle of radius R can be described by polar coordinates (R,q) ² In the Cartesian basis, a point on the circle : ² Carthesian Position: OP(t) = ( R cos( q(t) ), R sin( q(t) ) ) ² velocity: d v(t) = dt OP (t) d d v (t) = OP ( ) d dt ² v(t) = ( - R sin(q) dq/dt, R cos(q) dq/dt ) ² The velocity is along the tangential direction!
8 Newton s law Ø Acceleration : a = d dt v Ø we obtained: v(t) = ( - R sin(q) dq/dt, R cos(q) dq/dt ) Ø a = d/dt ( - R sin(q) dq/dt, R cos(q) dq/dt ) = [ Problem 1 ] Ø You will find that : Ø a = R (d 2 q(t)/dt 2 ) e t - R (dq/dt) 2 e n Ø e t =tangential direction, e n =normal direction Ø Acceleration along the tangential direction: at = R (d 2 q/dt 2 ) 8
9 Differential equation F = M a Newton law along tangential direction : * Force = - Mg sin(q(t)) * Acceleration = R (d 2 q(t) / dt 2 ) d 2 (t) dt 2 = g sin( (t)) R * We are looking for the function q(t) which satisfies this equation * With initial condition : (t = 0) = d d (t = 0) = 0 * To simplify the equation we use R=g : d 2 (t) dt 2 = sin( (t)) 9
10 How can we solve this? d 2 dx 2 = sin( ) * Let s start with small oscillations * In this limit, we can use : 0 sin( ) * Ha! And the differential equation simplifies to : * Guess for the solution? [Fill in] d 2 dx 2 = (x) =... * For the general case, we need a computer 10
11 11
12 A Differential equation: we know the derivative of a function at each point, but we do not know the function 12
13 ² Brook Taylor ( ) born in Edmonton ² Mathematician (Taylor series, foundation of derivation calculus) ² Differential equation: d y(x) =f(x, y) dx ² Question: Taylors s method ² If I start with the initial condition : y(x=a)=a, ² à what is the value b of y(x=b)? ² All we know is: where we are starting from (y(a)), and the derivative of y(x) at each point along the way up to x=b 13
14 * Taylor Algorithm : Taylor s method 1. Start at x=a, we know y(a) 2. Estimate y(a+h), where h is a small number 3. y(a+h) is obtained by using the slope (derivative) of the function y(x) at the point x=a 4. Do the same for x=a+h, starting from the obtained y(a+h), and do the estimation for y(a+2h), repeat until x=b. y i+1 = y i + hf(x i,y i ) * Iterative procedure (do loop) * The smaller the interval h, the smaller the error h 14
15 Discretization and error * For the program, we will use a fixed number of interval N * The larger N, the smaller the obtained error * a and b are given, this is the interval in which we want to solve the differential equation * The value y(a) is given by the initial condition h =(b a)/n x i = a + ih y i = y(a + ih) yi 0 = y 0 (a + ih) =f(x i,y i ) 15
16 Taylor s method: second order * The pendulum equation is a second order differential equation (double derivative), we need to modify the algorithm d 2 y dx 2 = sin(y) y(0) = y y 0 (0) = 0 * Simple idea: let s define a new variable: * We get the equation : Initial conditions dy dx = z dz dx = sin(y) * We obtain a system of coupled first order differential equations z(0) = 0 y(0) = y 16
17 Algorithm : F1(y,z) and F2(y,z) are given * More generally : dy dx = F 1(y, z) dz dx = F 2(y, z) y(0) = A z(0) = B 1. Start with : y(0)=a, z(0)=b 2. Compute : F1(y(0),z(0)) F2(y(0),z(0)) 3. Estimate y(h),z(h) : y(h) = y(0) + h*f1(y(0),z(0)) z(h) = z(0) + h*f2(y(0),z(0)) 4. Compute : F1(y(h),z(h)) F2(y(h),z(h)) 5. Estimate y(2h),z(2h) : y(2h) = y(h) + h*f1(y(h),z(h)) z(2h) = z(h) + h*f2(y(h),z(h)) 6... Repeat until y( N*h ) and z( N*h ) are obtained 17
18 Let s see what we get * pendulum angle as a function of time, with an initial dropping angle q i =p/4 Period When y(x) is zero the pendulum is in the down right position Periodic oscillation (no dissipation, perpetual movement). 18
19 19
20 Main structure of the code MODULE module library! contains! function F1(a,b) implicit none real(8) :: F1 real(8) :: a,b F1= [ FILL IN ] end function function F2(a,b) implicit none real(8) :: F2 real(8) :: a,b F2= [ FILL IN ] end function end module! PROGRAM program pendulum use library! implicit none integer,parameter :: N= Integer :: j real(8) :: h,a,b,alpha Real(8) :: x(0:n),y(0:n),z(0:n) a=0.0 ; b=30.0 ; h = [ FILL IN ] write(*,*) 'please enter initial angle' [FILL IN] x(0)=a y(0)=alpha z(0)=0.0 do j = 1, N y(j) = y(j-1) + h * F1( y(j-1),z(j-1) ) z(j) = z(j-1) + h * F2( y(j-1),z(j-1) ) x(j)= j*h End do end program 20
21 Reminder: Function Function name is F1 a and b are the arguments (inputs) The function is also defined as a variable which will contain the return value of the function module library! contains! function F1(a,b) Implicit none is an implicit none optional statement: if present, all real(8) :: F1 variables used in the real(8) :: a,b scope of the function F1= [ FILL IN ] have to be declared/ end function defined end module! 21
22 Part 1 Part 2 Part 3 22 program pendulum use library! implicit none integer, parameter :: N= integer :: j real(8) :: h,a,b,alpha real(8) :: x(0:n),y(0:n),z(0:n) a=0.0 ; b=30.0 ; h = [ FILL IN ] write(*,*) 'please enter initial angle' [FILL IN] x(0)=a y(0)=alpha z(0)=0.0 do j = 1, N y(j) = y(j-1) + h * F1( y(j-1), z(j-1) ) z(j) = z(j-1) + h * F2( y(j-1), z(j-1) ) x(j)= j*h End do end program Program structure Part 1 : use module and implicit none Part 2 : declaration of variables, the list of all objects you are going to work with (real(8) for double precision rational numbers and integer(4) for whole numbers) Part 3 : the calculations, we use variables to perform additions/ multiplications
23 program pendulum use library implicit none Integer(4),parameter :: N= Integer(4) :: j real(8) :: h,a,b,alpha real(8) :: x(0:n),y(0:n),z(0:n) Part 2 : variable declaration 1) Parameter statement (follows integer(4) or real(8) statements) : it will assign once and for all a value to this variable, which cannot be further changed. Basically a parameter variable is identical to replacing everywhere in your code N with in this example. 2) x(0:n) : the bracket after x denotes that x is an array (a vector). The general notation to declare an array is : x(i1:i2), where i1 is the first index of the vector, and i2 the last index of the vector. The dimension or size of this vector is hence i2-i1+1 The notation x(i2) is also valid, in this case the compiler assumes that i1=1 23
24 program pendulum use library implicit none integer,parameter :: N= integer :: j real(8) :: h,a,b,x(0:n),y(0:n),z(0:n),alpha a=0.0 ; b=30.0 ; h = [ FILL IN ] write(*,*) 'please enter initial angle' [FILL IN] x(0)=a y(0)=alpha z(0)=0.0 do j = 1, N y(j) = y(j-1) + h * F1( y(j-1),z(j-1) ) z(j) = z(j-1) + h * F2( y(j-1),z(j-1) ) x(j)= j*h End do end program We apply here the Taylor s method. The first thing to do is to find the small parameter h. The curve y(x) will be obtained for x=0,h,2h,3h If there are N points along y(x), and x goes from a to b, what is h? We then ask the initial angle (where we drop the pendulum from). How can you read the initial angle from the keyboard? 24
25 program pendulum use library implicit none integer,parameter :: N= integer :: j real(8) ::h,a,b,x(0:n),y(0:n),z(0:n),alpha x(0)=a y(0)=alpha z(0)=0.0 do j = 1, N y(j) = y(j-1) + h * F1( y(j-1),z(j-1) ) z(j) = z(j-1) + h * F2( y(j-1),z(j-1) ) x(j)= j*h End do end program x(0) : x was declared as an array, the indices of the array x range from 0 to N. x : array for the time, x(0)=first time t=0, x(1) second time t=h, x(2) third time t=2h y : array for the solution of the differential equation z: array for the function z(x)=dy/dx y(0)=alpha and z(0)=0 are the initial conditions of the problem 25
26 program pendulum use library implicit none integer,parameter :: N= integer :: j real(8) ::h,a,b,x(0:n),y(0:n),z(0:n),alpha x(0)=a y(0)=alpha z(0)=0.0 do j = 1, N y(j) = y(j-1) + h * F1( y(j-1),z(j-1) ) z(j) = z(j-1) + h * F2( y(j-1),z(j-1) ) x(j)= j*h End do end program We apply the algorithm (see page 15), and we compute the function y(h) and z(h) from the functions known at t=0 (y(0) and z(0)). y(h) = y(0) + h F1(y(0),z(0)) z(h) = z(0) + h F2(y(0),z(0)) This is a call to the function defined in the module. The argument a of the function F1(a,b) will take the value y(j-1), where y(j- 1) is the element j- 1 of the array y. First iteration of the do loop, j=1, so y(j- 1)=y(0), and z(j- 1)=z(0), these are the initial conditions. 26
27 program pendulum use library implicit none integer,parameter :: N= integer :: j real(8) ::h,a,b,x(0:n),y(0:n),z(0:n),alpha x(0)=a y(0)=alpha z(0)=0.0 do j = 1, N y(j) = y(j-1) + h * F1( y(j-1),z(j-1) ) z(j) = z(j-1) + h * F2( y(j-1),z(j-1) ) x(j) = j*h [FILL IN] End do end program For j=1: from the initial condition y(0) and z(0), we just found y(1) and z(1). Y(1) is the solution of the differential equation at time x=h. We also store the time positions at which we compute y(x) in the array x. Final step: we need to write the solution of the differential equation y(x) into a file, such that we can plot it with gnuplot or xmgrace. How can we do that? Fill the blank during the problem session. 27
28 Friday, problem session 1. Problem 1 : derive analytically the differential equations of the pendulum. 2. Problem 2 : calculate the solution y(x) of the differential equation for various initial condition, and plot the result. 28
Computational modeling
Computational modeling Lecture 1 : Linear algebra - Matrix operations Examination next week: How to get prepared Theory and programming: Matrix operations Instructor : Cedric Weber Course : 4CCP1 Schedule
More informationComputational modeling
Computational modeling Lecture 4 : Central Limit Theorem Theory: Normal distribution Programming: Arrays Instructor : Cedric Weber Course : 4CCP1000 Schedule Class/Week Chapter Topic Milestones 1 Monte
More informationMaterials Genome Assessment
Materials Genome Assessment Lecture 12 : Quantum Mechanics Theory: Differential equations with boundary conditions Schrodinger s cat Particle in a Box Programming: Modification of the Pendulum program
More informationdt 2 The Order of a differential equation is the order of the highest derivative that occurs in the equation. Example The differential equation
Lecture 18 : Direction Fields and Euler s Method A Differential Equation is an equation relating an unknown function and one or more of its derivatives. Examples Population growth : dp dp = kp, or = kp
More informationPractice / Lecture 8 : Differential equations Quantum Particle in a box
Problem 1 : Quantum Particle in a box In this problem we would like to modify the code of the pendulum, written in the last session, to solve the Schrodinger equation. Problem 1.a : d 2 y(x) dx 2 = 2m(V
More informationOrdinary Differential Equations (ODEs)
Ordinary Differential Equations (ODEs) 1 Computer Simulations Why is computation becoming so important in physics? One reason is that most of our analytical tools such as differential calculus are best
More informationEngineering Mechanics Prof. U. S. Dixit Department of Mechanical Engineering Indian Institute of Technology, Guwahati Kinematics
Engineering Mechanics Prof. U. S. Dixit Department of Mechanical Engineering Indian Institute of Technology, Guwahati Kinematics Module 10 - Lecture 24 Kinematics of a particle moving on a curve Today,
More informationSept , 17, 23, 29, 37, 41, 45, 47, , 5, 13, 17, 19, 29, 33. Exam Sept 26. Covers Sept 30-Oct 4.
MATH 23, FALL 2013 Text: Calculus, Early Transcendentals or Multivariable Calculus, 7th edition, Stewart, Brooks/Cole. We will cover chapters 12 through 16, so the multivariable volume will be fine. WebAssign
More informationTutorial-1, MA 108 (Linear Algebra)
Tutorial-1, MA 108 (Linear Algebra) 1. Verify that the function is a solution of the differential equation on some interval, for any choice of the arbitrary constants appearing in the function. (a) y =
More informationLimits and Continuity. 2 lim. x x x 3. lim x. lim. sinq. 5. Find the horizontal asymptote (s) of. Summer Packet AP Calculus BC Page 4
Limits and Continuity t+ 1. lim t - t + 4. lim x x x x + - 9-18 x-. lim x 0 4-x- x 4. sinq lim - q q 5. Find the horizontal asymptote (s) of 7x-18 f ( x) = x+ 8 Summer Packet AP Calculus BC Page 4 6. x
More informationab = c a If the coefficients a,b and c are real then either α and β are real or α and β are complex conjugates
Further Pure Summary Notes. Roots of Quadratic Equations For a quadratic equation ax + bx + c = 0 with roots α and β Sum of the roots Product of roots a + b = b a ab = c a If the coefficients a,b and c
More informationMath 113 Winter 2005 Key
Name Student Number Section Number Instructor Math Winter 005 Key Departmental Final Exam Instructions: The time limit is hours. Problem consists of short answer questions. Problems through are multiple
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 informationParametric Curves. Calculus 2 Lia Vas
Calculus Lia Vas Parametric Curves In the past, we mostly worked with curves in the form y = f(x). However, this format does not encompass all the curves one encounters in applications. For example, consider
More informationMATH 2554 (Calculus I)
MATH 2554 (Calculus I) Dr. Ashley K. University of Arkansas February 21, 2015 Table of Contents Week 6 1 Week 6: 16-20 February 3.5 Derivatives as Rates of Change 3.6 The Chain Rule 3.7 Implicit Differentiation
More informationMTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.
MTH4101 CALCULUS II REVISION NOTES 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) 1.1 Introduction Types of numbers (natural, integers, rationals, reals) The need to solve quadratic equations:
More informationPHYS 1114, Lecture 33, April 10 Contents:
PHYS 1114, Lecture 33, April 10 Contents: 1 This class is o cially cancelled, and has been replaced by the common exam Tuesday, April 11, 5:30 PM. A review and Q&A session is scheduled instead during class
More informationPurdue University Study Guide for MA Credit Exam
Purdue University Study Guide for MA 60 Credit Exam Students who pass the credit exam will gain credit in MA60. The credit exam is a twohour long exam with 5 multiple choice questions. No books or notes
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 informationIntegration, Separation of Variables
Week #1 : Integration, Separation of Variables Goals: Introduce differential equations. Review integration techniques. Solve first-order DEs using separation of variables. 1 Sources of Differential Equations
More informationChapter 11 ORDINARY DIFFERENTIAL EQUATIONS
Chapter 11 ORDINARY DIFFERENTIAL EQUATIONS The general form of a first order differential equations is = f(x, y) with initial condition y(a) = y a We seek the solution y = y(x) for x > a This is shown
More informationMotion in a Plane Uniform Circular Motion
Lecture 11 Chapter 8 Physics I Motion in a Plane Uniform Circular Motion Course website: http://faculty.uml.edu/andriy_danylov/teaching/physicsi IN THIS CHAPTER, you will learn to solve problems about
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 informationExam 4 SCORE. MA 114 Exam 4 Spring Section and/or TA:
Exam 4 Name: Section and/or TA: Last Four Digits of Student ID: Do not remove this answer page you will return the whole exam. You will be allowed two hours to complete this test. No books or notes may
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationNormal Force. W = mg cos(θ) Normal force F N = mg cos(θ) F N
Normal Force W = mg cos(θ) Normal force F N = mg cos(θ) Note there is no weight force parallel/down the include. The car is not pressing on anything causing a force in that direction. If there were a person
More informationFinal exam practice 4 (solutions) UCLA: Math 3B, Winter 2019
Final exam practice 4 (solutions) Instructor: Noah White Date: UCLA: Math 3B, Winter 2019 This exam has 7 questions, for a total of 80 points. Please print your working and answers neatly. Write your solutions
More informationSection 4.3 Vector Fields
Section 4.3 Vector Fields DEFINITION: A vector field in R n is a map F : A R n R n that assigns to each point x in its domain A a vector F(x). If n = 2, F is called a vector field in the plane, and if
More informationMATH 1241 Common Final Exam Fall 2010
MATH 1241 Common Final Exam Fall 2010 Please print the following information: Name: Instructor: Student ID: Section/Time: The MATH 1241 Final Exam consists of three parts. You have three hours for the
More informationPage Problem Score Max Score a 8 12b a b 10 14c 6 6
Fall 14 MTH 34 FINAL EXAM December 8, 14 Name: PID: Section: Instructor: DO NOT WRITE BELOW THIS LINE. Go to the next page. Page Problem Score Max Score 1 5 5 1 3 5 4 5 5 5 6 5 7 5 8 5 9 5 1 5 11 1 3 1a
More informationMultiple Choice Answers. MA 113 Calculus I Spring 2018 Exam 2 Tuesday, 6 March Question
MA 113 Calculus I Spring 2018 Exam 2 Tuesday, 6 March 2018 Name: Section: Last 4 digits of student ID #: This exam has 12 multiple choice questions (five points each) and 4 free response questions (ten
More information18.02 Multivariable Calculus Fall 2007
MIT OpenourseWare http://ocw.mit.edu 18.02 Multivariable alculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 18.02 Lecture 21. Test for
More informationBrushing up on Basic skills. Calculus AB (*problems for BC)
Brushing up on Basic skills To get you ready for Calculus AB (*problems for BC) Name: Directions: Use a pencil and the space provided next to each question to show all work. The purpose of this packet
More information4.9 Anti-derivatives. Definition. An anti-derivative of a function f is a function F such that F (x) = f (x) for all x.
4.9 Anti-derivatives Anti-differentiation is exactly what it sounds like: the opposite of differentiation. That is, given a function f, can we find a function F whose derivative is f. Definition. An anti-derivative
More informationMatthew S. Norton PHY 112 Lab Report CL-4 Damped Oscillator: Feynman-Newton Method for Solving First-Order Differential Equations March 2, 2006
Matthew S. Norton PHY 112 Lab Report CL-4 Damped Oscillator: Feynman-Newton Method for Solving First-Order Differential Equations March 2, 2006 Norton 0 Norton 1 Abstract In this lab, I used the Feynman-Newton
More informationReview session Midterm 1
AS.110.109: Calculus II (Eng) Review session Midterm 1 Yi Wang, Johns Hopkins University Fall 2018 7.1: Integration by parts Basic integration method: u-sub, integration table Integration By Parts formula
More information1 Motion of a single particle - Linear momentum, work and energy principle
1 Motion of a single particle - Linear momentum, work and energy principle 1.1 In-class problem A block of mass m slides down a frictionless incline (see Fig.). The block is released at height h above
More informationMOTION IN TWO OR THREE DIMENSIONS
MOTION IN TWO OR THREE DIMENSIONS 3 Sections Covered 3.1 : Position & velocity vectors 3.2 : The acceleration vector 3.3 : Projectile motion 3.4 : Motion in a circle 3.5 : Relative velocity 3.1 Position
More informationThe Principle of Least Action
The Principle of Least Action In their never-ending search for general principles, from which various laws of Physics could be derived, physicists, and most notably theoretical physicists, have often made
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationMATH 228: Calculus III (FALL 2016) Sample Problems for FINAL EXAM SOLUTIONS
MATH 228: Calculus III (FALL 216) Sample Problems for FINAL EXAM SOLUTIONS MATH 228 Page 2 Problem 1. (2pts) Evaluate the line integral C xy dx + (x + y) dy along the parabola y x2 from ( 1, 1) to (2,
More information18.02 Multivariable Calculus Fall 2007
MIT OpenourseWare http://ocw.mit.edu 8.02 Multivariable alculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 8.02 Lecture 8. hange of variables.
More informationENGI 4430 Parametric Vector Functions Page dt dt dt
ENGI 4430 Parametric Vector Functions Page 2-01 2. Parametric Vector Functions (continued) Any non-zero vector r can be decomposed into its magnitude r and its direction: r rrˆ, where r r 0 Tangent Vector:
More informationMathematics Specialist Units 3 & 4 Program 2018
Mathematics Specialist Units 3 & 4 Program 018 Week Content Assessments Complex numbers Cartesian Forms Term 1 3.1.1 review real and imaginary parts Re(z) and Im(z) of a complex number z Week 1 3.1. review
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 informationWave Phenomena Physics 15c
Wave Phenomena Physics 15c Lecture Harmonic Oscillators (H&L Sections 1.4 1.6, Chapter 3) Administravia! Problem Set #1! Due on Thursday next week! Lab schedule has been set! See Course Web " Laboratory
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 informationMath 308 Exam I Practice Problems
Math 308 Exam I Practice Problems This review should not be used as your sole source of preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..
More informationPHYSICS I. Lecture 1. Charudatt Kadolkar. Jul-Nov IIT Guwahati
PHYSICS I Lecture 1 Charudatt Kadolkar IIT Guwahati Jul-Nov 2014 Section 1 Introduction to the Course Syllabus Topics Classical Mechanics: Kinetic Energy rest mass energy Syllabus Topics Classical Mechanics:
More informationSolutions to old Exam 3 problems
Solutions to old Exam 3 problems Hi students! I am putting this version of my review for the Final exam review here on the web site, place and time to be announced. Enjoy!! Best, Bill Meeks PS. There are
More informationMidterm Exam 2. Wednesday, May 7 Temple, Winter 2018
Name: Student ID#: Section: Midterm Exam 2 Wednesday, May 7 Temple, Winter 2018 Show your work on every problem. orrect answers with no supporting work will not receive full credit. Be organized and use
More informationTIME SCHEDULE. Module Topic Periods I. 1.1 Functions, Limits & Continuity Differentiation. 20 II 2.1 Applications of Differentiation 27 III
SUBJECT TITLE : TECHNICAL MATHEMATICS I (Common for All Engineering / Technology programmes) SUBJECT CODE : PERIODS/WEEK : 6 PERIODS/SEMESTER : 8 CREDITS : 6 TIME SCHEDULE Module Topic Periods I. Functions,
More informationNumerical Optimization Algorithms
Numerical Optimization Algorithms 1. Overview. Calculus of Variations 3. Linearized Supersonic Flow 4. Steepest Descent 5. Smoothed Steepest Descent Overview 1 Two Main Categories of Optimization Algorithms
More informationExam 3 MATH Calculus I
Trinity College December 03, 2015 MATH 131-01 Calculus I By signing below, you attest that you have neither given nor received help of any kind on this exam. Signature: Printed Name: Instructions: Show
More informationLAURENTIAN UNIVERSITY UNIVERSITÉ LAURENTIENNE
Page 1 of 15 LAURENTIAN UNIVERSITY UNIVERSITÉ LAURENTIENNE Friday, December 14 th Course and No. 2007 MATH 2066 EL Date................................... Cours et no........................... Total no.
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 informationENGI Parametric Vector Functions Page 5-01
ENGI 3425 5. Parametric Vector Functions Page 5-01 5. Parametric Vector Functions Contents: 5.1 Arc Length (Cartesian parametric and plane polar) 5.2 Surfaces of Revolution 5.3 Area under a Parametric
More informationINTRODUCTION, FOUNDATIONS
1 INTRODUCTION, FOUNDATIONS ELM1222 Numerical Analysis Some of the contents are adopted from Laurene V. Fausett, Applied Numerical Analysis using MATLAB. Prentice Hall Inc., 1999 2 Today s lecture Information
More informationNumerical Methods for Initial Value Problems; Harmonic Oscillators
1 Numerical Methods for Initial Value Problems; Harmonic Oscillators Lab Objective: Implement several basic numerical methods for initial value problems (IVPs), and use them to study harmonic oscillators.
More informationSlope Fields, Differential Equations... and Maybe Euler s Method, Too
Slope Fields, Differential Equations... and Maybe Euler s Method, Too Chuck Garner, Ph.D. Department of Mathematics Rockdale Magnet School for Science and Technology October 10, 2014 / Georgia DOE AP Content
More informationPaper Reference. Further Pure Mathematics FP1 Advanced/Advanced Subsidiary. Friday 26 January 2007 Afternoon Time: 1 hour 30 minutes
Centre No. Candidate No. Paper Reference(s) 6674/01 Edexcel GCE Further Pure Mathematics FP1 Advanced/Advanced Subsidiary Friday 6 January 007 Afternoon Time: 1 hour 30 minutes Materials required for examination
More informationSYSTEMS OF ODES. mg sin ( (x)) dx 2 =
SYSTEMS OF ODES Consider the pendulum shown below. Assume the rod is of neglible mass, that the pendulum is of mass m, and that the rod is of length `. Assume the pendulum moves in the plane shown, and
More informationWhat will you learn?
Section 2.2 Basic Differentiation Rules & Rates of Change Calc What will you learn? Find the derivative using the Constant Rule Find the derivative using the Power Rule Find the derivative using the Constant
More informationLecture 9: Taylor Series
Math 8 Instructor: Padraic Bartlett Lecture 9: Taylor Series Week 9 Caltech 212 1 Taylor Polynomials and Series When we first introduced the idea of the derivative, one of the motivations we offered was
More informationPaper Reference. Paper Reference(s) 6668/01 Edexcel GCE Further Pure Mathematics FP2 Advanced/Advanced Subsidiary
Centre No. Candidate No. Paper Reference 6 6 6 8 0 1 Surname Signature Paper Reference(s) 6668/01 Edexcel GCE Further Pure Mathematics FP2 Advanced/Advanced Subsidiary Thursday 24 June 2010 Morning Time:
More informationFinal Exam. Monday March 19, 3:30-5:30pm MAT 21D, Temple, Winter 2018
Name: Student ID#: Section: Final Exam Monday March 19, 3:30-5:30pm MAT 21D, Temple, Winter 2018 Show your work on every problem. orrect answers with no supporting work will not receive full credit. Be
More informationResearch Computing with Python, Lecture 7, Numerical Integration and Solving Ordinary Differential Equations
Research Computing with Python, Lecture 7, Numerical Integration and Solving Ordinary Differential Equations Ramses van Zon SciNet HPC Consortium November 25, 2014 Ramses van Zon (SciNet HPC Consortium)Research
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Chapter 1 Introduction and Basic Terminology Most of the phenomena studied in the sciences and engineering involve processes that change with time. For example, it is well known
More informationLecture 7 - Separable Equations
Lecture 7 - Separable Equations Separable equations is a very special type of differential equations where you can separate the terms involving only y on one side of the equation and terms involving only
More information- - - - - - - - - - - - - - - - - - DISCLAIMER - - - - - - - - - - - - - - - - - - General Information: This midterm is a sample midterm. This means: The sample midterm contains problems that are of similar,
More informationHomework for Math , Spring 2013
Homework for Math 3220 2, Spring 2013 A. Treibergs, Instructor April 10, 2013 Our text is by Joseph L. Taylor, Foundations of Analysis, American Mathematical Society, Pure and Applied Undergraduate Texts
More informationd` = 1+( dy , which is part of the cone.
7.5 Surface area When we did areas, the basic slices were rectangles, with A = h x or h y. When we did volumes of revolution, the basic slices came from revolving rectangles around an axis. Depending on
More informationHere are brief notes about topics covered in class on complex numbers, focusing on what is not covered in the textbook.
Phys374, Spring 2008, Prof. Ted Jacobson Department of Physics, University of Maryland Complex numbers version 5/21/08 Here are brief notes about topics covered in class on complex numbers, focusing on
More informationPhysics 141. Lecture 18. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 18, Page 1
Physics 141. Lecture 18. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 18, Page 1 Physics 141. Lecture 18. Course Information. Topics to be discussed today: A
More informationWave Phenomena Physics 15c
Wave Phenomena Physics 15c Lecture Harmonic Oscillators (H&L Sections 1.4 1.6, Chapter 3) Administravia! Problem Set #1! Due on Thursday next week! Lab schedule has been set! See Course Web " Laboratory
More informationInfinite Series. 1 Introduction. 2 General discussion on convergence
Infinite Series 1 Introduction I will only cover a few topics in this lecture, choosing to discuss those which I have used over the years. The text covers substantially more material and is available for
More informationSudoku Puzzle A.P. Exam (Part B) Questions are from the 1997 and 1998 A.P. Exams A Puzzle by David Pleacher
Sudoku Puzzle A.P. Exam (Part B) Questions are from the 1997 and 1998 A.P. Exams A Puzzle by David Pleacher Solve the 4 multiple-choice problems below. A graphing calculator is required for some questions
More informationFinal exam practice 4 UCLA: Math 3B, Winter 2019
Instructor: Noah White Date: Final exam practice 4 UCLA: Math 3B, Winter 2019 This exam has 7 questions, for a total of 80 points. Please print your working and answers neatly. Write your solutions in
More informationODE Background: Differential (1A) Young Won Lim 12/29/15
ODE Background: Differential (1A Copyright (c 2011-2015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version
More informationPhysics 201, Lecture 28
Physics 01, Lecture 8 Today s Topics n Oscillations (Ch 15) n n n More Simple Harmonic Oscillation n Review: Mathematical Representation n Eamples: Simple Pendulum, Physical pendulum Damped Oscillation
More informationAnalysis Methods in Atmospheric and Oceanic Science
Analysis Methods in Atmospheric and Oceanic Science AOSC 652 Week 7, Day 1 13 Oct 2014 1 Student projects: 20% of the final grade: you will receive a numerical score for the project and final grade will
More informationIntroduction Derivation General formula List of series Convergence Applications Test SERIES 4 INU0114/514 (MATHS 1)
MACLAURIN SERIES SERIES 4 INU0114/514 (MATHS 1) Dr Adrian Jannetta MIMA CMath FRAS Maclaurin Series 1/ 21 Adrian Jannetta Recap: Binomial Series Recall that some functions can be rewritten as a power series
More informationThis practice exam is intended to help you prepare for the final exam for MTH 142 Calculus II.
MTH 142 Practice Exam Chapters 9-11 Calculus II With Analytic Geometry Fall 2011 - University of Rhode Island This practice exam is intended to help you prepare for the final exam for MTH 142 Calculus
More informationTake-Home Exam 1: pick up on Thursday, June 8, return Monday,
SYLLABUS FOR 18.089 1. Overview This course is a review of calculus. We will start with a week-long review of single variable calculus, and move on for the remaining five weeks to multivariable calculus.
More informationModule M4.3 Further differentiation
F L E X I B L E L E A R N I N G A P P R O A C H T O P H Y S I C S Module M4. Opening items 2 Derivatives of composite functions 2. A function of a function 2.2 The chain rule 2. Differentiating implicit
More informationNATIONAL BOARD FOR HIGHER MATHEMATICS. Research Scholarships Screening Test. Saturday, February 2, Time Allowed: Two Hours Maximum Marks: 40
NATIONAL BOARD FOR HIGHER MATHEMATICS Research Scholarships Screening Test Saturday, February 2, 2008 Time Allowed: Two Hours Maximum Marks: 40 Please read, carefully, the instructions on the following
More informationChapter 10 Conics, Parametric Equations, and Polar Coordinates Conics and Calculus
Chapter 10 Conics, Parametric Equations, and Polar Coordinates 10.1 Conics and Calculus 1. Parabola A parabola is the set of all points x, y ( ) that are equidistant from a fixed line and a fixed point
More information2.2 The derivative as a Function
2.2 The derivative as a Function Recall: The derivative of a function f at a fixed number a: f a f a+h f(a) = lim h 0 h Definition (Derivative of f) For any number x, the derivative of f is f x f x+h f(x)
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationMath 116 Final Exam April 24, 2017
On my honor, as a student, I have neither given nor received unauthorized aid on this academic work. Initials: Do not write in this area Your Initials Only: Math 6 Final Exam April 24, 207 Your U-M ID
More informationPursuit Curves. Molly Severdia May 16, 2008
Pursuit Curves Molly Severdia May 16, 2008 Abstract This paper will discuss the differential equations which describe curves of pure pursuit, in which the pursuer s velocity vector is always pointing directly
More informationEnergy in a Simple Harmonic Oscillator. Class 30. Simple Harmonic Motion
Simple Harmonic Motion Class 30 Here is a simulation of a mass hanging from a spring. This is a case of stable equilibrium in which there is a large extension in which the restoring force is linear in
More informationPhysics 201, Lecture 8
Physics 01, Lecture 8 Today s Topics q Physics 01, Review 1 q Important Notes: v v v v This review is not designed to be complete on its own. It is not meant to replace your own preparation efforts Exercises
More informationExtended Introduction to Computer Science CS1001.py. Lecture 8 part A: Finding Zeroes of Real Functions: Newton Raphson Iteration
Extended Introduction to Computer Science CS1001.py Lecture 8 part A: Finding Zeroes of Real Functions: Newton Raphson Iteration Instructors: Benny Chor, Amir Rubinstein Teaching Assistants: Yael Baran,
More informationSECTION A. f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes.
SECTION A 1. State the maximal domain and range of the function f(x) = ln(x). Sketch the graph of y = f(x), indicating the coordinates of any points where the graph crosses the axes. 2. By evaluating f(0),
More informationVector Calculus, Maths II
Section A Vector Calculus, Maths II REVISION (VECTORS) 1. Position vector of a point P(x, y, z) is given as + y and its magnitude by 2. The scalar components of a vector are its direction ratios, and represent
More informationMath 233. Practice Problems Chapter 15. i j k
Math 233. Practice Problems hapter 15 1. ompute the curl and divergence of the vector field F given by F (4 cos(x 2 ) 2y)i + (4 sin(y 2 ) + 6x)j + (6x 2 y 6x + 4e 3z )k olution: The curl of F is computed
More informationPhysics Exam 1 Formulas
INSTRUCTIONS: Write your NAME on the front of the blue exam booklet. The exam is closed book, and you may have only pens/pencils and a calculator (no stored equations or programs and no graphing). Show
More informationWeek 12: Optimisation and Course Review.
Week 12: Optimisation and Course Review. MA161/MA1161: Semester 1 Calculus. Prof. Götz Pfeiffer School of Mathematics, Statistics and Applied Mathematics NUI Galway November 21-22, 2016 Assignments. Problem
More informationFinal Exam Review Exercise Set A, Math 1551, Fall 2017
Final Exam Review Exercise Set A, Math 1551, Fall 2017 This review set gives a list of topics that we explored throughout this course, as well as a few practice problems at the end of the document. A complete
More information