Project 2: Using linear systems for numerical solution of boundary value problems

Size: px
Start display at page:

Download "Project 2: Using linear systems for numerical solution of boundary value problems"

Transcription

1 LINEAR ALGEBRA, MATH 124 Instructor: Dr. T.I. Lakoba Project 2: Using linear systems for numerical solution of boundary value problems Goal Introduce one of the most important applications of Linear Algebra to Engineering. General requirements You may work alone or with one other person. If you work with someone else, hand in one answer sheet with both of your names on it. No groups bigger than two. No collaboration between groups. Write your answers on the answer sheet provided in the last few pages of this document. Staple 1 all the paper showing your neatly presented 2 work to the answer sheet. Introduction Most phenomena in Physics and Engineering are modeled by differential equations. You were introduced to simple differential equations of the form dy/dx = f(x, y) in Calculus II; perhaps, you also took a course (MATH 230 or MATH 271) focused on methods for solving such equations. A boundary value problem (BVP), whose approximate solution is illustrated in this Project, is just a differential equation supplied with two boundary values for the unknown variable y. BVPs occur in a great variety of applications, such as (i) oscillation of membranes and strings (think of music!), (ii) sagging of beams and other structures under mechanical load, (iii) radiation by atoms and molecules, and (iv) hurricane modeling, to name just a few. Most BVPs arising in modern applications of science and engineering are so complex that they cannot be solved exactly with pen and paper. Then one has to resort to solving them approximately on a computer. There is a whole area of applied mathematics which deals solely with this kind of solution of differential equations. It is called Numerical Analysis, and Linear Algebra lies at its foundation. Numerical analysis is very much driven by the need to solve practical problems: There are quite a few hi-tech companies which develop commercial software for solving differential equations on a computer. In this Project, you will learn how to apply your knowledge of linear systems to solving one simple BVP. Linear systems enter this problem from two different angles. First, you will set up and solve a linear system to find an approximate formula for the second derivative of a function. Second, you will use that formula to set up another linear system, which describes the numerical solution of the BVP. Then you will use MATLAB to solve the latter linear system on a computer. 1 Your grade will be reduced by 5% if you hand in a pile of non-stapled sheets. 2 I will reduce your grade by an amount left to my discretion in each particular case if your work is presented in a messy way and I have to waste time deciphering it. 1

2 T y i ( i=8 ) cable yh (x,y) 0 T W(x) x i ( i=8 ) L 0 x Problem description and method of solution i = A flexible power cable has one end staked to the ground and the other end fastened H ft up a pole L ft away, as shown in the left figure above. Let y be the cable s elevation at a point located x horizontal units from the stake. Our goal will be to find y(x). We will do so by (approximately) solving the BVP satisfied by y(x). This BVP consists of a second-order differential equation and two boundary conditions. The differential equation is: d 2 y dx = W (x), (1) 2 T where the notations W and T are as follows (see also the left figure above): W (x) is the cable s weight per unit length at a point (x, y(x)); W may depend on x if, for example, part of the cable is iced; T is the (constant) tension of the cable. We assume that both W (x) and T are known; their values will be given later on. Since the ends of the cable are pinned at specific locations, the following boundary conditions must supplement the differential equation: y(0) = 0 and y(l) = H. (2) The above equations (1) and (2) form a BVP for the cable s elevation, y = y(x), as a function of the distance, x, from the stake. The BVP (1) & (2) is sufficiently simple so that one can solve it exactly using just pen and paper. However, as you read in the Introduction, most differential equations of practical interest cannot be solved exactly. To solve them approximately, one has to use an approach that we illustrate below for the simple BVP (1) & (2). This approach ultimately relies on Linear Algebra. The key idea behind finding an approximate solution to BVP (1) & (2) is this. Instead of looking for the cable s elevation, y(x), for all values of x within the interval [0, L], limit your goal to finding y only at discrete values of x on that interval. To this end, let h = L/n and consider discrete points This is illustrated in the right figure above for n = Note that x i = i h, where i = 0, 1,..., n. (3a) x 0 = 0 and x n = L. (3b) 3 Incidentally, notice that the points x i subdivide the interval [0, L] into n subintervals of length h = L/n: [0, h], [h, 2h], [2h, 3h],..., [(n 1)h, L]. This may look familiar: You saw exactly the same approach when you first learned definite integration and Riemann sums in Calculus I and later learned various applications of it, such as finding areas and volumes, in Calculus II. 2

3 Further, denote y(x i ) = y i. So now, finding an approximate solution to the BVP (1) & (2) will amount to finding only the finitely many values y i, rather than the continuous function y(x) for all values of x on [0, L]. Clearly, if h is taken to be very small, there is, practically, almost no difference between the continous variables x and y, on one hand, and the discrete sets {x 0, x 1, x 2,..., x n } and {y 0, y 1, y 2,..., y n }, on the other. Yet, this discretization trick will allow you to replace the BVP given by Equations (1) and (2), which you do not know how to solve, by a linear system, which you know how to solve. The steps of transition from the BVP to a linear system are described below. (3c) Step-by-step instructions 1. Your goal in this part is to derive an approximate formula for the second derivative, y (x), that involves only the discrete values of y(x) at certain points. It is shown in higher-level courses like Numerical Analysis that for any sufficiently smooth function y(x), the second derivative y (x) at a point x = a (for some arbitrary a) can be approximated using the values of y(x) at three nearby points: x = a h, x = a, and x = a + h: y (a) A 1 y(a h) + A 0 y(a) + A 1 y(a + h), (4) where A 1, A 0, A 1 are some constant coefficients. Find A 1, A 0, A 1 in Equation (4), using as an example either pp of the textbook or/and Appendix A found after these Instructions. Write your formula (i.e., Equation (4) with your numerical values for the coefficients A 1, A 0, A 1 ) on the answer sheet and attach to it the paper with your derivation. Check your answer by comparing it with the solution of Exercise 21 of Section Let now a be one of the points x i defined in Equation (3a) above: a = x i. For a nearby point a h, we then have: a h = x i h = ih h = (i 1)h = x i 1. According to Equation (3c), we can write y(a) = y(x i ) = y i, and then similarly: y(a h) = y(x i 1 ) = y i 1. Use these observations to rewrite the approximate formula you have obtained for y (a) (that is, for y (x i )) in terms of y i 1, y i, etc. Next, at each interior point x i of the interval [0, L] (i.e., when 0 < i < n), approximate the differential equation (1) using the approximate formula for y (x i ) that you have just written down. To this end, simply substitute that formula into the left-hand side of Equation (1). On the right-hand side of (1), use similar notations, i.e., W (x i ) = W i. Write the resulting equation at the point x i in the answer sheet. Thus, the outcome of this step should have the form ( )y i 1 + ( )y i + ( )y i+1 = (something known), (5) and your job is to supply the coefficients or expressions inside the parentheses. 3

4 3. Here you will practice setting up a linear system to approximate Equation (1) by using a small number n of discretization points. Later you will need to use a large n. Draw the interval [0, L], subdivide it into 5 subintervals, and label the resulting points x 0, x 1,..., x 5 (see Equations (3a) and (3b)). Write your equation (5) for i = 1. Then repeat this for each interior point x i of [0, L]. Note that each of your equations couples the value y i to the values y i 1 and y i+1. Whenever you require the values y 0 and y n, these are found from the boundary conditions (2). Thus, the equations at the interior points x i, i = 1,..., 4 form a linear system for the values y i. Write this linear system for the unknown values y 1, y 2, y 3, y 4 in a matrix form. (Note that W i and T are known in this problem.) Use the notation h for L/5, etc Set up and solve a linear system (5) for the case n = 20 using MATLAB (see Appendix B after these Instructions). Use the values L = 100 ft, H = 12 ft, T = 3000 lb, and { 6 lb/ft (iced cable), 0 x 50 ft W (x) = 3 lb/ft (clean cable), 50 < x 100 ft. Write the numerical values of your solution at several locations along the cable, as requested in the answer sheet. Round these values to 4 decimal places. Also, use MATLAB to plot the graph y = y(x) which you obtained on the interval 0 x L (i.e., include the end points). Label the axes. For example, to label the x-axis, type: xlabel( whatever you want to use for the x-label, fontsize,12). Print your computer output, including the plot of y(x), and staple it to the answer sheet. Extra credit (partial credit is given only if the solution is mostly correct) Consider a broken power line, where one end of the cable is still fastened to the pole but the other (free) end lies on the ground. At the point of contact with the ground, the cable s slope is zero, and therefore the condition y (0) = 0 (6) must be added to the boundary conditions (2). (Note that in this case, the tension T at the point of contact of the cable with the ground is provided solely by the static friction between the cable and the ground.) The resulting BVP, consisting of Equations (1), (2), and (6), has a solution only for one particular value of L. This L is the distance from the pole to the point where the broken cable first touches the ground. Find L by setting up and solving a linear system for this case. Use the values: T = 500 lb, W (x) = 6 lb/ft (assume the cable is uniformly iced), and the other parameters as before. Approximate the first derivative in the boundary condition (6) by the formula For definiteness, use h = 1 ft. y (a) y(a + h) y(a) h. (7) Hint: One (but not the only) way to solve this problem is as follows. Use only two (which ones?) of the three boundary conditions (2) and (6) to view this as an initial value problem (recall Calculus II). Then use the third boundary condition to find L. 4 That is, if you are to write L 2 /25, replace it with h 2 in your formula; similarly, replace 5/L with 1/h, and so on. This will become very handy in the next Step, where the small number n = 5 will be increased to a more realistic value. 4

5 Appendix A Here I present an example of finding the coefficients in a formula similar to Equation (4) of this Project. An analogous derivation is found on pp of the text book. Recall from Calculus I that the first derivative of a function is defined as y 1 (a) = lim (y(a + h) y(a)), h 0 h which means that y (a) 1 (y(a + h) y(a)) for sufficiently small h. (A1) h Note that Equation (A1) is satisfied exactly if y(x) = const and y(x) = x (verify this by substituting these two y(x) into (A1)!). If, however, y = x 2, the right-hand side of Equation (A1) yields 2a + h, which equals y (a) = 2a only approximately when h is small. Then we ask: Can one modify the expression on the right-hand side of Equation (A1) so that the new expression would exactly equal y (x) not only for y = const and y = x, but also for y = x 2? The short answer is yes. The details are provided below. We look for the required expression in the form given by the right-hand side of Equation (4): 5 y (a) A 1 y(a h) + A 0 y(a) + A 1 y(a + h). (A2) Recall that we want this equation to be satisfied exactly rather than approximately whenever y = 1, y = x, and y = x 2. Then we simply substitute each of these three functions into Equation (A2): when y(x) = 1 : 1 = 0 = A A A 1 1 (A3a) Note that the three 1 s on the r.h.s. occurred because y(a h) = y(a) = y(a + h) = 1 when y(x) 1. when y(x) = x : x = 1 = A 1 (a h) + A 0 (a) + A 1 (a + h) (A3b) Similarly to the previous note, the coefficients of A 1, A 0, A 1 came from the fact that y(a h) = a h etc. for y(x) = x. when y(x) = x 2 : (x 2 ) at x=a = 2a = A 1 (a h) 2 + A 0 (a) 2 + A 1 (a + h) 2 (A3c) The linear system (A3) for the unknown coefficients A 1, A 0, A 1 can be simplified once we notice that we have required it to be true for any x = a. In particular, it must be true for a = 0. Substituting a = 0 into Equations (A3) and transforming the resulting linear system to the reduced echelon form, we find (verify): A 1 = 1/(2h), A 0 = 0, A 1 = 1/(2h). (A4) The values for A 1, A 0, A 1 in Equation (4) are obtained similarly, also by requiring that that equation be exact for y = 1, y = x, and y = x 2. Appendix B Here I include some hints about MATLAB that may help you in successfully completing this Project. 5 You may wonder why the expressions on the right-hand sides in Equations (A2) and (4) have the same form, even though the derivatives on the left-hand sides are different. The answer is given in courses like Numerical Analysis, which is offered at UVM as Math 237 and

6 Read Appendix A in the text book (you do not need the material from subsections A.4, A.5, A.10, and A.11) and MATLAB Primer by K. Sigmon, posted on the class web site (you do not need the material from pp and 18 21). You may, but do not have to, read other MATLAB resources posted there. To open a script file for the first time, start MATLAB, then go to File in the main menu. Select New and then M-file. This opens a window where you can write your script. Save it in some directory on your \M: drive 6. When you want to run your script at any later time, make sure you change the directory (by typing cd m:\directoryname in the MATLAB command window) to the directory where your script file is saved. The coefficient matrix for solving the linear system for y 1, y 2,..., y n 1 has the dimension (n 1) (n 1). Thus, for n = 20, it contains 361 entries. However, only the elements on the main diagonal and the two diagonals closest to it are nonzero (you may easily see this after you complete Step 3 of this assignment). Typing all the 361 entries one by one into the matrix is not a good option. A good option is to use the for-loop of MATLAB (see the Primer, Sec. 6). For example, the loop that defines matrix is: for i=1:2 A(i,i)=5; A(i,i+1)=6; end A(3,3)=7; A few notes are in order. A = First, note that the entries which you do not specify are set to zero by default in MATLAB. Second, note that in the above loop, we defined nonzero entries only of the first two rows of matrix A. This is because the indices of such entries in these rows follow a certain pattern. The indices of the nonzero entry in the third row do not follow exactly the same pattern, and hence we had to define that entry outside the loop. Likewise, in the matrix that you have to set up for this Project, indices of nonzero entries follow the same pattern in most, but not in all, rows. Therefore, first identify those rows with such a pattern and define nonzero entries in these rows within a loop. Then, define nonzero entries whose indices do not follow a patern, outside the loop. Third, in the vector on the right-hand side of your linear system, you should also determine which entries follow a pattern and define them within a loop (or loops). Then, define the entries which do not follow a pattern outside the loop. Finally, make sure that the matrix you have constructed is indeed the one you intended to construct. You may do so by displaying that matrix for some n less than 20. This, on one hand, will not change the logic of your script and, on the other, will produce a matrix of a smaller dimension, which you should be able to inspect on the screen. 6 The same letter, M, in the words M-file and M: drive is a coincidence. There is no relation between MATLAB s M-files and the CEMS s M: drive. 6

7 To solve a linear system, one has many options in MATLAB. One is by using the already familiar command rref(b), where B is the augmented matrix. Another way is to use the inverse matrix A 1, for which MATLAB s notation is inv(a). This, however, may be slower than other options if A has large dimensions. A much better way is to use the backslash notation which MATLAB s founders invented: see Sec. 3 in MATLAB Primer or simply google backslash matlab. This comment is about plotting y(x) on the interval [0, L]. Note that your solution y 1,..., y n 1 does not include the boundary values y(0) and y(l) because they were known, not solved for. Therefore, you need to append these values to your solution vector. See Sec. A.4 in Appendix A of the textbook or Sec. 5 in the MATLAB Primer. Here is an example of how to append values to a row vector: a = [1 2]; a = [0 a 3] Note that since your solution {y 1, y 2,... } is a column (rather than row) vector, you will need to slightly modify the above example. Finally, you will need to define a vector x in order to plot your solution versus it. This can be done with a for-loop, but a much better way is described in Sec. 10 of the Primer; a related example is also found in Sec. 11 there. If you would like to learn more about numerical methods, Take the course MATH 337. Acknowledgement This project is based on an idea by Professor Scott C. Fulton (Clarkson University) and some material from Chapter 8 of the book Numerical Analysis: A Practical Approach by M. Maron and R. Lopez, 3rd Ed., Wadsworth (Belmont, 1991). 7

8 LINEAR ALGEBRA, MATH 124 Answer Sheet for Project 2 Reminder: You MUST staple pages with your neatly presented work to get full credit. Name(s): 1. (27 points) Approximate formula for y (a): 2. (15 points) Approximation to the differential equation (1) at the interior point x i : 3. (27 points) Linear system (in matrix form) for the case n = 5: 4. (31 points) Solution for the elevation y i computed for n = 20 and rounded to four decimal places: i x i, ft y i, ft Attach your code and a Matlab-produced plot of your solution; do not forget to label the axes. Extra credit (27 points) Attach your work. L = ft 8

Using Microsoft Excel

Using Microsoft Excel Using Microsoft Excel Objective: Students will gain familiarity with using Excel to record data, display data properly, use built-in formulae to do calculations, and plot and fit data with linear functions.

More information

Lab 2 Worksheet. Problems. Problem 1: Geometry and Linear Equations

Lab 2 Worksheet. Problems. Problem 1: Geometry and Linear Equations Lab 2 Worksheet Problems Problem : Geometry and Linear Equations Linear algebra is, first and foremost, the study of systems of linear equations. You are going to encounter linear systems frequently in

More information

Calculus at Rutgers. Course descriptions

Calculus at Rutgers. Course descriptions Calculus at Rutgers This edition of Jon Rogawski s text, Calculus Early Transcendentals, is intended for students to use in the three-semester calculus sequence Math 151/152/251 beginning with Math 151

More information

Experimental Uncertainty (Error) and Data Analysis

Experimental Uncertainty (Error) and Data Analysis E X P E R I M E N T 1 Experimental Uncertainty (Error) and Data Analysis INTRODUCTION AND OBJECTIVES Laboratory investigations involve taking measurements of physical quantities, and the process of taking

More information

New Mexico Tech Hyd 510

New Mexico Tech Hyd 510 Vectors vector - has magnitude and direction (e.g. velocity, specific discharge, hydraulic gradient) scalar - has magnitude only (e.g. porosity, specific yield, storage coefficient) unit vector - a unit

More information

Basic Linear Algebra in MATLAB

Basic Linear Algebra in MATLAB Basic Linear Algebra in MATLAB 9.29 Optional Lecture 2 In the last optional lecture we learned the the basic type in MATLAB is a matrix of double precision floating point numbers. You learned a number

More information

Experiment: Oscillations of a Mass on a Spring

Experiment: Oscillations of a Mass on a Spring Physics NYC F17 Objective: Theory: Experiment: Oscillations of a Mass on a Spring A: to verify Hooke s law for a spring and measure its elasticity constant. B: to check the relationship between the period

More information

Experimental Uncertainty (Error) and Data Analysis

Experimental Uncertainty (Error) and Data Analysis Experimental Uncertainty (Error) and Data Analysis Advance Study Assignment Please contact Dr. Reuven at yreuven@mhrd.org if you have any questions Read the Theory part of the experiment (pages 2-14) and

More information

Newton's 2 nd Law. . Your end results should only be interms of m

Newton's 2 nd Law. . Your end results should only be interms of m Newton's nd Law Introduction: In today's lab you will demonstrate the validity of Newton's Laws in predicting the motion of a simple mechanical system. The system that you will investigate consists of

More information

Mon Jan Improved acceleration models: linear and quadratic drag forces. Announcements: Warm-up Exercise:

Mon 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 information

Temperature measurement

Temperature measurement Luleå University of Technology Johan Carlson Last revision: July 22, 2009 Measurement Technology and Uncertainty Analysis - E7021E Lab 3 Temperature measurement Introduction In this lab you are given a

More information

Mathematica Project 3

Mathematica Project 3 Mathematica Project 3 Name: Section: Date: On your class s Sakai site, your instructor has placed 5 Mathematica notebooks. Please use the following table to determine which file you should select based

More information

LAB 2 - ONE DIMENSIONAL MOTION

LAB 2 - ONE DIMENSIONAL MOTION Name Date Partners L02-1 LAB 2 - ONE DIMENSIONAL MOTION OBJECTIVES Slow and steady wins the race. Aesop s fable: The Hare and the Tortoise To learn how to use a motion detector and gain more familiarity

More information

1.2. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy

1.2. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy .. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy = f(x, y). In this section we aim to understand the solution

More information

Computer projects for Mathematical Statistics, MA 486. Some practical hints for doing computer projects with MATLAB:

Computer projects for Mathematical Statistics, MA 486. Some practical hints for doing computer projects with MATLAB: Computer projects for Mathematical Statistics, MA 486. Some practical hints for doing computer projects with MATLAB: You can save your project to a text file (on a floppy disk or CD or on your web page),

More information

Mathematics 1104B. Systems of Equations and Inequalities, and Matrices. Study Guide. Text: Mathematics 11. Alexander and Kelly; Addison-Wesley, 1998.

Mathematics 1104B. Systems of Equations and Inequalities, and Matrices. Study Guide. Text: Mathematics 11. Alexander and Kelly; Addison-Wesley, 1998. Adult Basic Education Mathematics Systems of Equations and Inequalities, and Matrices Prerequisites: Mathematics 1104A, Mathematics 1104B Credit Value: 1 Text: Mathematics 11. Alexander and Kelly; Addison-Wesley,

More information

Student s Printed Name: _KEY Grading Guidelines CUID:

Student s Printed Name: _KEY Grading Guidelines CUID: Student s Printed Name: _KEY Grading Guidelines CUID: Instructor: Section # : You are not permitted to use a calculator on any portion of this test. You are not allowed to use any textbook, notes, cell

More information

Static and Kinetic Friction

Static and Kinetic Friction Experiment 12 If you try to slide a heavy box resting on the floor, you may find it difficult to get the box moving. Static friction is the force that is acting against the box. If you apply a light horizontal

More information

Linear Motion with Constant Acceleration

Linear Motion with Constant Acceleration Linear Motion 1 Linear Motion with Constant Acceleration Overview: First you will attempt to walk backward with a constant acceleration, monitoring your motion with the ultrasonic motion detector. Then

More information

INTEGRALS. In Chapter 2, we used the tangent and velocity problems to introduce the derivative the central idea in differential calculus.

INTEGRALS. In Chapter 2, we used the tangent and velocity problems to introduce the derivative the central idea in differential calculus. INTEGRALS 5 INTEGRALS In Chapter 2, we used the tangent and velocity problems to introduce the derivative the central idea in differential calculus. INTEGRALS In much the same way, this chapter starts

More information

The Revolution of the Moons of Jupiter

The Revolution of the Moons of Jupiter The Revolution of the Moons of Jupiter Overview: During this lab session you will make use of a CLEA (Contemporary Laboratory Experiences in Astronomy) computer program generously developed and supplied

More information

Introduction to Computer Tools and Uncertainties

Introduction to Computer Tools and Uncertainties Experiment 1 Introduction to Computer Tools and Uncertainties 1.1 Objectives To become familiar with the computer programs and utilities that will be used throughout the semester. To become familiar with

More information

Instructor Notes for Chapters 3 & 4

Instructor Notes for Chapters 3 & 4 Algebra for Calculus Fall 0 Section 3. Complex Numbers Goal for students: Instructor Notes for Chapters 3 & 4 perform computations involving complex numbers You might want to review the quadratic formula

More information

Newton s Second Law. Sample

Newton s Second Law. Sample Newton s Second Law Experiment 4 INTRODUCTION In your discussion of Newton s first law, you learned that when the sum of the forces acting on an object is zero, its velocity does not change. However, when

More information

Student s Printed Name:

Student s Printed Name: Student s Printed Name: Instructor: CUID: Section # : You are not permitted to use a calculator on any part of this test. You are not allowed to use any textbook, notes, cell phone, laptop, PDA, or any

More information

Math 308 Midterm Answers and Comments July 18, Part A. Short answer questions

Math 308 Midterm Answers and Comments July 18, Part A. Short answer questions Math 308 Midterm Answers and Comments July 18, 2011 Part A. Short answer questions (1) Compute the determinant of the matrix a 3 3 1 1 2. 1 a 3 The determinant is 2a 2 12. Comments: Everyone seemed to

More information

LABORATORY 1: KINEMATICS written by Melissa J. Wafer '95 June 1993

LABORATORY 1: KINEMATICS written by Melissa J. Wafer '95 June 1993 LABORATORY 1: KINEMATICS written by Melissa J. Wafer '95 June 1993 The purpose of this exercise is to re-enforce what you have learned about kinematics in class and to familiarize you with computer resources

More information

Chapter 8 Linear Algebraic Equations

Chapter 8 Linear Algebraic Equations PowerPoint to accompany Introduction to MATLAB for Engineers, Third Edition William J. Palm III Chapter 8 Linear Algebraic Equations Copyright 2010. The McGraw-Hill Companies, Inc. This work is only for

More information

MATH 333: Partial Differential Equations

MATH 333: Partial Differential Equations MATH 333: Partial Differential Equations Problem Set 9, Final version Due Date: Tues., Nov. 29, 2011 Relevant sources: Farlow s book: Lessons 9, 37 39 MacCluer s book: Chapter 3 44 Show that the Poisson

More information

Daily Update. Dondi Ellis. January 27, 2015

Daily Update. Dondi Ellis. January 27, 2015 Daily Update Dondi Ellis January 27, 2015 CLASS 1: Introduction and Section 1.1 REMINDERS: Assignment: Read sections 1.1 and 1.2, and Student Guide (online). Buy a TI-84 or other graphing calculator satisfying

More information

Solving Equations by Adding and Subtracting

Solving Equations by Adding and Subtracting SECTION 2.1 Solving Equations by Adding and Subtracting 2.1 OBJECTIVES 1. Determine whether a given number is a solution for an equation 2. Use the addition property to solve equations 3. Determine whether

More information

Question: Total. Points:

Question: Total. Points: MATH 308 May 23, 2011 Final Exam Name: ID: Question: 1 2 3 4 5 6 7 8 9 Total Points: 0 20 20 20 20 20 20 20 20 160 Score: There are 9 problems on 9 pages in this exam (not counting the cover sheet). Make

More information

Lesson 5b Solving Quadratic Equations

Lesson 5b Solving Quadratic Equations Lesson 5b Solving Quadratic Equations In this lesson, we will continue our work with Quadratics in this lesson and will learn several methods for solving quadratic equations. The first section will introduce

More information

Student s Printed Name: KEY_&_Grading Guidelines_CUID:

Student s Printed Name: KEY_&_Grading Guidelines_CUID: Student s Printed Name: KEY_&_Grading Guidelines_CUID: Instructor: Section # : You are not permitted to use a calculator on any portion of this test. You are not allowed to use any textbook, notes, cell

More information

MATH 1242 FINAL EXAM Spring,

MATH 1242 FINAL EXAM Spring, MATH 242 FINAL EXAM Spring, 200 Part I (MULTIPLE CHOICE, NO CALCULATORS).. Find 2 4x3 dx. (a) 28 (b) 5 (c) 0 (d) 36 (e) 7 2. Find 2 cos t dt. (a) 2 sin t + C (b) 2 sin t + C (c) 2 cos t + C (d) 2 cos t

More information

Introduction. Pre-Lab Questions: Physics 1CL PERIODIC MOTION - PART II Fall 2009

Introduction. Pre-Lab Questions: Physics 1CL PERIODIC MOTION - PART II Fall 2009 Introduction This is the second of two labs on simple harmonic motion (SHM). In the first lab you studied elastic forces and elastic energy, and you measured the net force on a pendulum bob held at an

More information

Sample Questions PREPARING FOR THE AP (AB) CALCULUS EXAMINATION. tangent line, a+h. a+h

Sample Questions PREPARING FOR THE AP (AB) CALCULUS EXAMINATION. tangent line, a+h. a+h Sample Questions PREPARING FOR THE AP (AB) CALCULUS EXAMINATION B B A B tangent line,, a f '(a) = lim h 0 f(a + h) f(a) h a+h a b b f(x) dx = lim [f(x ) x + f(x ) x + f(x ) x +...+ f(x ) x ] n a n B B

More information

9.1 - Systems of Linear Equations: Two Variables

9.1 - Systems of Linear Equations: Two Variables 9.1 - Systems of Linear Equations: Two Variables Recall that a system of equations consists of two or more equations each with two or more variables. A solution to a system in two variables is an ordered

More information

Making Piecewise Functions Continuous and Differentiable by Dave Slomer

Making Piecewise Functions Continuous and Differentiable by Dave Slomer Making Piecewise Functions Continuous and Differentiable by Dave Slomer Piecewise-defined functions are applied in areas such as Computer Assisted Drawing (CAD). Many piecewise functions in textbooks are

More information

Lesson 4 Linear Functions and Applications

Lesson 4 Linear Functions and Applications In this lesson, we take a close look at Linear Functions and how real world situations can be modeled using Linear Functions. We study the relationship between Average Rate of Change and Slope and how

More information

Lecture 5b: Starting Matlab

Lecture 5b: Starting Matlab Lecture 5b: Starting Matlab James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University August 7, 2013 Outline 1 Resources 2 Starting Matlab 3 Homework

More information

Motion II. Goals and Introduction

Motion II. Goals and Introduction Motion II Goals and Introduction As you have probably already seen in lecture or homework, and if you ve performed the experiment Motion I, it is important to develop a strong understanding of how to model

More information

MATLAB Project 2: MATH240, Spring 2013

MATLAB Project 2: MATH240, Spring 2013 1. Method MATLAB Project 2: MATH240, Spring 2013 This page is more information which can be helpful for your MATLAB work, including some new commands. You re responsible for knowing what s been done before.

More information

Solving Equations by Factoring. Solve the quadratic equation x 2 16 by factoring. We write the equation in standard form: x

Solving Equations by Factoring. Solve the quadratic equation x 2 16 by factoring. We write the equation in standard form: x 11.1 E x a m p l e 1 714SECTION 11.1 OBJECTIVES 1. Solve quadratic equations by using the square root method 2. Solve quadratic equations by completing the square Here, we factor the quadratic member of

More information

Math 309 Notes and Homework for Days 4-6

Math 309 Notes and Homework for Days 4-6 Math 309 Notes and Homework for Days 4-6 Day 4 Read Section 1.2 and the notes below. The following is the main definition of the course. Definition. A vector space is a set V (whose elements are called

More information

The Theory of Second Order Linear Differential Equations 1 Michael C. Sullivan Math Department Southern Illinois University

The Theory of Second Order Linear Differential Equations 1 Michael C. Sullivan Math Department Southern Illinois University The Theory of Second Order Linear Differential Equations 1 Michael C. Sullivan Math Department Southern Illinois University These notes are intended as a supplement to section 3.2 of the textbook Elementary

More information

MATH 18.01, FALL PROBLEM SET # 2

MATH 18.01, FALL PROBLEM SET # 2 MATH 18.01, FALL 2012 - PROBLEM SET # 2 Professor: Jared Speck Due: by Thursday 4:00pm on 9-20-12 (in the boxes outside of Room 2-255 during the day; stick it under the door if the room is locked; write

More information

CSCE 155N Fall Homework Assignment 2: Stress-Strain Curve. Assigned: September 11, 2012 Due: October 02, 2012

CSCE 155N Fall Homework Assignment 2: Stress-Strain Curve. Assigned: September 11, 2012 Due: October 02, 2012 CSCE 155N Fall 2012 Homework Assignment 2: Stress-Strain Curve Assigned: September 11, 2012 Due: October 02, 2012 Note: This assignment is to be completed individually - collaboration is strictly prohibited.

More information

Problem Worth Score Total 14

Problem Worth Score Total 14 MATH 241, Fall 14 Extra Credit Preparation for Final Name: INSTRUCTIONS: Write legibly. Indicate your answer clearly. Revise and clean up solutions. Do not cross anything out. Rewrite the page, I will

More information

Objectives. Materials

Objectives. Materials Activity 8 Exploring Infinite Series Objectives Identify a geometric series Determine convergence and sum of geometric series Identify a series that satisfies the alternating series test Use a graphing

More information

SECTION 1-5 Integer Exponents

SECTION 1-5 Integer Exponents 42 Basic Algebraic Operations In Problems 45 52, imagine that the indicated solutions were given to you by a student whom you were tutoring in this class. (A) Is the solution correct? If the solution is

More information

5.6 Logarithmic and Exponential Equations

5.6 Logarithmic and Exponential Equations SECTION 5.6 Logarithmic and Exponential Equations 305 5.6 Logarithmic and Exponential Equations PREPARING FOR THIS SECTION Before getting started, review the following: Solving Equations Using a Graphing

More information

MA 125 CALCULUS I SPRING 2007 April 27, 2007 FINAL EXAM. Name (Print last name first):... Student ID Number (last four digits):...

MA 125 CALCULUS I SPRING 2007 April 27, 2007 FINAL EXAM. Name (Print last name first):... Student ID Number (last four digits):... CALCULUS I, FINAL EXAM 1 MA 125 CALCULUS I SPRING 2007 April 27, 2007 FINAL EXAM Name (Print last name first):............................................. Student ID Number (last four digits):........................

More information

Math Practice Final - solutions

Math Practice Final - solutions Math 151 - Practice Final - solutions 2 1-2 -1 0 1 2 3 Problem 1 Indicate the following from looking at the graph of f(x) above. All answers are small integers, ±, or DNE for does not exist. a) lim x 1

More information

Forces and Newton s Second Law

Forces and Newton s Second Law Forces and Newton s Second Law Goals and Introduction Newton s laws of motion describe several possible effects of forces acting upon objects. In particular, Newton s second law of motion says that when

More information

2: SIMPLE HARMONIC MOTION

2: SIMPLE HARMONIC MOTION 2: SIMPLE HARMONIC MOTION Motion of a Mass Hanging from a Spring If you hang a mass from a spring, stretch it slightly, and let go, the mass will go up and down over and over again. That is, you will get

More information

AP Physics B SUMMER ASSIGNMENT Summer, Instructor: Rhett Butler

AP Physics B SUMMER ASSIGNMENT Summer, Instructor: Rhett Butler AP Physics B SUMMER ASSIGNMENT Summer, 2013 Instructor: Rhett Butler Welcome to AP Physics! I am delighted to be in the position to help further your understanding of the physical world around us. I also

More information

LAB 3 - VELOCITY AND ACCELERATION

LAB 3 - VELOCITY AND ACCELERATION Name Date Partners L03-1 LAB 3 - VELOCITY AND ACCELERATION OBJECTIVES A cheetah can accelerate from 0 to 50 miles per hour in 6.4 seconds. Encyclopedia of the Animal World A Jaguar can accelerate from

More information

College Physics 201 Graphing Review and Analysis

College Physics 201 Graphing Review and Analysis Northern Michigan University Graphing Review and Analysis Page-1 College Physics 201 Graphing Review and Analysis Materials: Graph paper, ruler, pencil and calculator. 1 Purpose The goal of these laboratory

More information

Integrals. D. DeTurck. January 1, University of Pennsylvania. D. DeTurck Math A: Integrals 1 / 61

Integrals. D. DeTurck. January 1, University of Pennsylvania. D. DeTurck Math A: Integrals 1 / 61 Integrals D. DeTurck University of Pennsylvania January 1, 2018 D. DeTurck Math 104 002 2018A: Integrals 1 / 61 Integrals Start with dx this means a little bit of x or a little change in x If we add up

More information

Ingredients of Multivariable Change: Models, Graphs, Rates & 9.1 Cross-Sectional Models and Multivariable Functions

Ingredients of Multivariable Change: Models, Graphs, Rates & 9.1 Cross-Sectional Models and Multivariable Functions Chapter 9 Ingredients of Multivariable Change: Models, Graphs, Rates & 9.1 Cross-Sectional Models and Multivariable Functions For a multivariable function with two input variables described by data given

More information

Student s Printed Name:

Student s Printed Name: Student s Printed Name: Instructor: CUID: Section # : You are not permitted to use a calculator on any portion of this test. You are not allowed to use any textbook, notes, cell phone, laptop, PDA, smart

More information

Contents 16. Higher Degree Equations

Contents 16. Higher Degree Equations Contents 16. Higher Degree Equations 2 16.3 Finding Roots of Higher Degree Equations................. 2 Example 16.15............................... 2 Example 16.16............................... 2 Example

More information

If you need more room, use the backs of the pages and indicate that you have done so.

If you need more room, use the backs of the pages and indicate that you have done so. Math 125 Final Exam Winter 2018 Your Name Your Signature Student ID # Quiz Section Professor s Name TA s Name Turn off and stow away all cell phones, watches, pagers, music players, and other similar devices.

More information

Math 115 Calculus Exam I October 6, 1999

Math 115 Calculus Exam I October 6, 1999 Math 5 Calculus Exam I October 6, 999 Department of Mathematics University of Michigan 999 Name: Instructor: Signature: Section: General Instructions: Do not open this exam until you are told to begin.

More information

PHY 123 Lab 9 Simple Harmonic Motion

PHY 123 Lab 9 Simple Harmonic Motion PHY 123 Lab 9 Simple Harmonic Motion (updated 11/17/16) The purpose of this lab is to study simple harmonic motion of a system consisting of a mass attached to a spring. You will establish the relationship

More information

AP Calculus BC Summer Assignment Mrs. Comeau

AP Calculus BC Summer Assignment Mrs. Comeau AP Calculus BC Summer Assignment 2015-2016 Mrs. Comeau Please complete this assignment DUE: the first day of class, SEPTEMBER 2nd. Email me if you have questions, or need help over the summer. I would

More information

Some Notes on Linear Algebra

Some Notes on Linear Algebra Some Notes on Linear Algebra prepared for a first course in differential equations Thomas L Scofield Department of Mathematics and Statistics Calvin College 1998 1 The purpose of these notes is to present

More information

PHYSICS LAB FREE FALL. Date: GRADE: PHYSICS DEPARTMENT JAMES MADISON UNIVERSITY

PHYSICS LAB FREE FALL. Date: GRADE: PHYSICS DEPARTMENT JAMES MADISON UNIVERSITY PHYSICS LAB FREE FALL Printed Names: Signatures: Date: Lab Section: Instructor: GRADE: PHYSICS DEPARTMENT JAMES MADISON UNIVERSITY Revision August 2003 Free Fall FREE FALL Part A Error Analysis of Reaction

More information

ARE YOU READY FOR CALCULUS?

ARE YOU READY FOR CALCULUS? ARE YOU READY FOR CALCULUS? Congratulations! You made it to Calculus AB! Instructions 1. Please complete the packet (see below), which will be due the day of registration. This packet will help you review

More information

Lesson 4b More on LOGARITHMS

Lesson 4b More on LOGARITHMS Lesson 4b More on LOGARITHMS Learning Packet Student Name Due Date Class Time/Day Submission Date THIS BOX FOR INSTRUCTOR GRADING USE ONLY Mini-Lesson is complete and information presented is as found

More information

7.1 Indefinite Integrals Calculus

7.1 Indefinite Integrals Calculus 7.1 Indefinite Integrals Calculus Learning Objectives A student will be able to: Find antiderivatives of functions. Represent antiderivatives. Interpret the constant of integration graphically. Solve differential

More information

Lab 10 Logarithmic Functions

Lab 10 Logarithmic Functions Name Student ID # Instructor Lab Period Date Due Lab 10 Logarithmic Functions Objectives 1. To develop an ability to visualize the graphs of logarithmic functions. 2. To become familiar with applications

More information

Lab 2: Static Response, Cantilevered Beam

Lab 2: Static Response, Cantilevered Beam Contents 1 Lab 2: Static Response, Cantilevered Beam 3 1.1 Objectives.......................................... 3 1.2 Scalars, Vectors and Matrices (Allen Downey)...................... 3 1.2.1 Attribution.....................................

More information

Purdue University Study Guide for MA Credit Exam

Purdue 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 information

Introduction. Pre-Lab Questions: Physics 1CL PERIODIC MOTION - PART II Spring 2009

Introduction. Pre-Lab Questions: Physics 1CL PERIODIC MOTION - PART II Spring 2009 Introduction This is the second of two labs on simple harmonic motion (SHM). In the first lab you studied elastic forces and elastic energy, and you measured the net force on a pendulum bob held at an

More information

Intermediate Algebra. Gregg Waterman Oregon Institute of Technology

Intermediate Algebra. Gregg Waterman Oregon Institute of Technology Intermediate Algebra Gregg Waterman Oregon Institute of Technology c 2017 Gregg Waterman This work is licensed under the Creative Commons Attribution 4.0 International license. The essence of the license

More information

Math Week 1 notes

Math Week 1 notes Math 2250-004 Week 1 notes We will not necessarily finish the material from a given day's notes on that day. Or on an amazing day we may get farther than I've predicted. We may also add or subtract some

More information

12.2 Simplifying Radical Expressions

12.2 Simplifying Radical Expressions x n a a m 1 1 1 1 Locker LESSON 1. Simplifying Radical Expressions Texas Math Standards The student is expected to: A.7.G Rewrite radical expressions that contain variables to equivalent forms. Mathematical

More information

How to Write a Laboratory Report

How to Write a Laboratory Report How to Write a Laboratory Report For each experiment you will submit a laboratory report. Laboratory reports are to be turned in at the beginning of the lab period, one week following the completion of

More information

Math 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 # 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 information

Page Points Score Total: 100

Page Points Score Total: 100 Math 1130 Spring 2019 Sample Midterm 2a 2/28/19 Name (Print): Username.#: Lecturer: Rec. Instructor: Rec. Time: This exam contains 10 pages (including this cover page) and 9 problems. Check to see if any

More information

Calculus is Cool. Math 160 Calculus for Physical Scientists I Exam 1 September 18, 2014, 5:00-6:50 pm. NAME: Instructor: Time your class meets:

Calculus is Cool. Math 160 Calculus for Physical Scientists I Exam 1 September 18, 2014, 5:00-6:50 pm. NAME: Instructor: Time your class meets: NAME: Instructor: Time your class meets: Math 160 Calculus for Physical Scientists I Exam 1 September 18, 2014, 5:00-6:50 pm How can it be that mathematics, being after all a product of human thought independent

More information

MATH10212 Linear Algebra B Homework Week 5

MATH10212 Linear Algebra B Homework Week 5 MATH Linear Algebra B Homework Week 5 Students are strongly advised to acquire a copy of the Textbook: D C Lay Linear Algebra its Applications Pearson 6 (or other editions) Normally homework assignments

More information

Math 320, spring 2011 before the first midterm

Math 320, spring 2011 before the first midterm Math 320, spring 2011 before the first midterm Typical Exam Problems 1 Consider the linear system of equations 2x 1 + 3x 2 2x 3 + x 4 = y 1 x 1 + 3x 2 2x 3 + 2x 4 = y 2 x 1 + 2x 3 x 4 = y 3 where x 1,,

More information

FINAL - PART 1 MATH 150 SPRING 2017 KUNIYUKI PART 1: 135 POINTS, PART 2: 115 POINTS, TOTAL: 250 POINTS No notes, books, or calculators allowed.

FINAL - PART 1 MATH 150 SPRING 2017 KUNIYUKI PART 1: 135 POINTS, PART 2: 115 POINTS, TOTAL: 250 POINTS No notes, books, or calculators allowed. Math 150 Name: FINAL - PART 1 MATH 150 SPRING 2017 KUNIYUKI PART 1: 135 POINTS, PART 2: 115 POINTS, TOTAL: 250 POINTS No notes, books, or calculators allowed. 135 points: 45 problems, 3 pts. each. You

More information

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam. Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam. Signature: Instructor s Name and Section: (Circle Your Section)

More information

Learning Packet. Lesson 5b Solving Quadratic Equations THIS BOX FOR INSTRUCTOR GRADING USE ONLY

Learning Packet. Lesson 5b Solving Quadratic Equations THIS BOX FOR INSTRUCTOR GRADING USE ONLY Learning Packet Student Name Due Date Class Time/Day Submission Date THIS BOX FOR INSTRUCTOR GRADING USE ONLY Mini-Lesson is complete and information presented is as found on media links (0 5 pts) Comments:

More information

Math 1526 Excel Lab 2 Summer 2012

Math 1526 Excel Lab 2 Summer 2012 Math 1526 Excel Lab 2 Summer 2012 Riemann Sums, Trapezoidal Rule and Simpson's Rule: In this lab you will learn how to recover information from rate of change data. For instance, if you have data for marginal

More information

36-309/749 Math Review 2014

36-309/749 Math Review 2014 36-309/749 Math Review 2014 The math content of 36-309 is not high. We will use algebra, including logs. We will not use calculus or matrix algebra. This optional handout is intended to help those students

More information

PHY 123 Lab 10-Simple Harmonic Motion

PHY 123 Lab 10-Simple Harmonic Motion 1 To print higher-resolution math symbols, click the Hi-Res Fonts for Printing button on the jsmath control panel. PHY 123 Lab 10-Simple Harmonic Motion The purpose of this lab is to study simple harmonic

More information

Midterm 1 revision source for MATH 227, Introduction to Linear Algebra

Midterm 1 revision source for MATH 227, Introduction to Linear Algebra Midterm revision source for MATH 227, Introduction to Linear Algebra 5 March 29, LJB page 2: Some notes on the Pearson correlation coefficient page 3: Practice Midterm Exam page 4: Spring 27 Midterm page

More information

MATH 60 Course Notebook Chapter #1

MATH 60 Course Notebook Chapter #1 MATH 60 Course Notebook Chapter #1 Integers and Real Numbers Before we start the journey into Algebra, we need to understand more about the numbers and number concepts, which form the foundation of Algebra.

More information

Comprehensive Exam in Real Analysis Fall 2006 Thursday September 14, :00-11:30am INSTRUCTIONS

Comprehensive Exam in Real Analysis Fall 2006 Thursday September 14, :00-11:30am INSTRUCTIONS Exam Packet # Comprehensive Exam in Real Analysis Fall 2006 Thursday September 14, 2006 9:00-11:30am Name (please print): Student ID: INSTRUCTIONS (1) The examination is divided into three sections to

More information

LAB 2: INTRODUCTION TO MOTION

LAB 2: INTRODUCTION TO MOTION Lab 2 - Introduction to Motion 3 Name Date Partners LAB 2: INTRODUCTION TO MOTION Slow and steady wins the race. Aesop s fable: The Hare and the Tortoise Objectives To explore how various motions are represented

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Experiment 03: Work and Energy

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Experiment 03: Work and Energy MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.01 Fall Term 2010 Experiment 03: Work and Energy Purpose of the Experiment: In this experiment you allow a cart to roll down an inclined

More information

REVIEW PROBLEMS FOR MIDTERM I MATH 2373, SPRING 2019 UNIVERSITY OF MINNESOTA ANSWER KEY

REVIEW PROBLEMS FOR MIDTERM I MATH 2373, SPRING 2019 UNIVERSITY OF MINNESOTA ANSWER KEY REVIEW PROBLEMS FOR MIDTERM I MATH 2373, SPRING 209 UNIVERSITY OF MINNESOTA ANSWER KEY This list of problems is not guaranteed to be a complete review. For a complete review make sure that you know how

More information

Math Refresher Answer Sheet (NOTE: Only this answer sheet and the following graph will be evaluated)

Math Refresher Answer Sheet (NOTE: Only this answer sheet and the following graph will be evaluated) Name: Score: / 50 Math Refresher Answer Sheet (NOTE: Only this answer sheet and the following graph will be evaluated) 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. MAKE SURE CALCULATOR

More information

AMS 27L LAB #6 Winter 2009

AMS 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 information

Solution: By inspection, the standard matrix of T is: A = Where, Ae 1 = 3. , and Ae 3 = 4. , Ae 2 =

Solution: By inspection, the standard matrix of T is: A = Where, Ae 1 = 3. , and Ae 3 = 4. , Ae 2 = This is a typical assignment, but you may not be familiar with the material. You should also be aware that many schools only give two exams, but also collect homework which is usually worth a small part

More information