1. Consider the initial value problem: find y(t) such that. y = y 2 t, y(0) = 1.
|
|
- Amice Norris
- 5 years ago
- Views:
Transcription
1 Engineering Mathematics CHEN30101 solutions to sheet 3 1. Consider the initial value problem: find y(t) such that y = y 2 t, y(0) = 1. Take a step size h = 0.1 and verify that the forward Euler approximation to y(0.1) is y 1 = 1.0 and that the approximations to subsequent points are y 2 = 1.01, y 3 = and y 4 = Compare these estimates with the exact values y(0.1), y(0.2), y(0.3) and y(0.4) obtained by solving the ODE algebraically. Running the forward Euler (fe) process the estimates y j y(t j ) are tabulated below. The results are rounded to four decimal places. h t j y j f(y j,t j ) y(t j ) Note that the exact solution is the function y(t) = 2 2 t Bytakingastepsizeofh = 0.5, computetheforward Eulerapproximation of z(5), where z(t) satisfies the initial value problem: z = 1 z, z(0) = 1. 50(1+t) Running (fe) the approximate solution z j is tabulated below. The results are rounded to four decimal places. h t j z j f(z j,t j ) z(t j ) Note that the exact solution is the function z(t) = (1+t) 1/50. 1
2 3. Consider the initial value problem: find y(t) such that y = π 1 y 2, y(0) = 0. Take a step size h = 0.1 and show that computing the forward Euler approximation to y(1) is problematic. Running the forward Euler process gives... h t j y j f(y j,t j ) If y = , y = π 1 y 2 cannot be calculated and the method stops. 4. Consider the initial value problem: find y(x) such that y = 1 2 ex y, y(0) = 0. Compute the forward Euler approximation of y(1) using (a) a single step of size h = 1, (b) two steps each of size h = 0.5, (c) four steps each of size h = 0.25, and (d) ten steps each of size h = 0.1. Compare your estimates with the exact value. (a) Running (fe) for h = 1 gives so y(1) 0.5. (b) Running (fe) for h = 0.5 gives so y(1) (c) Running (fe) for h = 0.25 gives so y(1) h x j y j f(y j,x j ) h x j y j f(y j,x j ) h x j y j f(y j,x j )
3 (d) Running (fe) for h = 0.1 gives so y(1) h x j y j f(y j,x j ) The solution is y(x) = ln ( ) 1+e x 2 so y(1) = to four decimal places. The accuracy of the forward Euler approximation improves when the step size h is reduced. 5. By taking a step size of h = 0.5, compute the modified Euler approximation of z(2), where z(t) satisfies the initial value problem: z = 1 z, z(0) = 1. 50(1+t) Running the modified Euler (me) process the approximate solution y j is tabulated below. Note that t m = t j + h 2 and z m = z j + h 2 f(z j,t j ). h t j z j f(z j,t j ) t m z m f(z m,t m ) z(t j ) These results are clearly much more accurate than those obtained for the forward Euler approximation scheme (cf. question 2). 6. Consider the initial value problem: find y(x) such that y = 17 12x 3y, y(0) = 9. Using four steps of size h = 0.25, estimate y(1) using (i) the forward Euler method, (ii) the modified Euler method. Compare your estimates with the exact value. 3
4 The exact solution is the function y(x) = 2e 3x + 7 4x. Running the forward Euler process gives... h x j y j f(y j,x j ) Running the modified Euler process gives... h x j y j f(y j,x j ) x m y m f(y m,x m ) y(x j ) Neither solution is particularly accurate but the modified Euler method gives more accurate results than the forward Euler method. 7. Heun s method for approximating the solution of an initial value problem is given by K 1 = f(x n,y n ), K 2 = f(x n +h,y n +hk 1 ) y n+1 = y n + h 2 (K 1 +K 2 ). Carry out two steps of the method (with step size h = 0.5) to compute an estimate of the solution y(x) of the problem y = y x, y(0) = 20. Running the process with x = x n +h and y = y n +hk 1 gives... h x n y n K 1 x y K 2 y(t n ) Consider the initial value problem: find y(t) such that y = 1+lny, y(0) = 1. Carry out a step of the 4th-order Runge Kutta method (with step size h = 0.1) and verify that K 1 = 1, K 2 = , K 3 = and K 4 = so that the approximation to y(0.1) is given by y 1 =
5 For the next step of the method, show that K 1 = , K 2 = , K 3 = and K 4 = so that y 2 = Continue for one more step and verify that y 3 = Running the Runge Kutta (rk4) process gives... h = 0.1 t = 0.00 y 0 = K 1 = first temporary point t = 0.05 y m = K 2 = second temporary point t = 0.05 y m = K 3 = third temporary point t = 0.10 y = K 4 = new point t = 0.10 y 1 = h = 0.1 t = 0.10 y 1 = K 1 = first temporary point t = 0.15 y m = K 2 = second temporary point t = 0.15 y m = K 3 = third temporary point t = 0.20 y = K 4 = new point t = 0.20 y 2 = h = 0.1 t = 0.20 y 2 = K 1 = first temporary point t = 0.25 y m = K 2 = second temporary point t = 0.25 y m = K 3 = third temporary point t = 0.30 y = K 4 = new point t = 0.30 y 3 = Considerfluiddrainingoutof atank accordingtotheequation y = k y with k = 0.1 and y = 20 at t = 0. Use the 4th-order Runge Kutta method with step size h = 10 to estimate the values of y(t) (the height, in metres, of fluid in the tank) for t = 10, 20 and 30 measured in seconds. Running the (rk4) process gives... first step h = 10 t = 0 y 0 = K 1 = first temporary point t = 5 y m = K 2 = second temporary point t = 5 y m = K 3 = third temporary point t = 10 y = K 4 = new point t = 10 y 1 = h = 10 t = 10 y 1 = K 1 = first temporary point t = 15 y m = K 2 = second temporary point t = 15 y m = K 3 = third temporary point t = 20 y = K 4 = new point t = 20 y 2 = h = 10 t = 20 y 2 = K 1 = first temporary point t = 25 y m = K 2 = second temporary point t = 25 y m = K 3 = third temporary point t = 30 y = K 4 = new point t = 30 y 3 =
6 10. Consider the initial value problem: find y(x) such that y = 2xy 2 y(0) = 1. Carry out two steps of each of the following methods for solving the equation. (a) forward Euler with h = 1 (b) forward Euler with h = 0.5 (c) modified Euler with h = 1 (d) modified Euler with h = 0.5 (e) Runge Kutta with h = 1 (f) Runge Kutta with h = 0.5. Compare your results with the exact solution. (a) Running (fe) with h = 1... Unrealistic results. (b) Running (fe) with h = h x j y j f(y j,x j ) y(x j ) h x j y j f(y j,x j ) y(x j ) The final value is correct (by coincidence!). (c) Running (me) with h = 1... h x j y j f(y j,x j ) x m y m f(y m,x m ) y(x j ) All subsequent values will be zero. Not a realistic solution. (d) Running (me) with h = h x j y j f(y j,x j ) x m y m f(y m,x m ) y(x j ) These values are not particularly accurate but are pointing in the right direction. 6
7 (e) Running (rk4) with h = 1... h = 1.0 x = 0.00 y 0 = K 1 = first temporary point x = 0.50 y m = K 2 = second temporary point x = 0.50 y m = K 3 = third temporary point x = 1.00 y = K 4 = new point x = 1.00 y 1 = h = 1.0 x = 1.00 y 1 = K 1 = first temporary point x = 1.50 y m = K 2 = second temporary point x = 1.50 y m = K 3 = third temporary point x = 2.00 y = K 4 = new point x = 2.00 y = These values are not particularly accurate but are again pointing in the correct direction. (f) Running (rk4) with h = h = 0.5 x = 0.00 y 0 = K 1 = first temporary point x = 0.25 y m = K 2 = second temporary point x = 0.25 y m = K 3 = third temporary point x = 0.50 y = K 4 = new point x = 0.50 y 1 = h = 0.5 x = 0.50 y 1 = K 1 = first temporary point x = 0.75 y m = K 2 = second temporary point x = 0.75 y m = K 3 = third temporary point x = 1.00 y = K 4 = new point x = 1.00 y 2 = These values are much more accurate. 11. Express the following high-order equations as systems of first-order equations. (a) d2 y 2 +12y = 0; y(0) = 1, dy (0) = 0. (b) d2 y 2 +2dy 4y = x; y(0) = 0, dy (0) = 0. (c) d3 y 3 dy +y = dy x3, y(0) = 1, (0) = 0, d 2 y 2(0) = 0. (d) d2 p dt 2 = 1 dp (p+t) dt ; p(0) = 1, dp (0) = 2. dt (e) the coupled system d 2 y dt 2 2dz dt = y +z; y(t 0) = y 0, d 2 z dt 2 +2dy dt = y z; z(t 0) = z 0, dy dt (t 0) = v A, dz dt (t 0) = v B. 7
8 (a) Introducing z = dy dz so dy = z, y(0) = 1, dz = 12y, z(0) = 0. (b) Introducing z = dy dz so dy = z, y(0) = 0, dz = x+4y 2z, z(0) = 0. dz and w = (c) Introducing z = dy dy = z, y(0) = 1, dz = w, z(0) = 0, dw = x3 +z y, w(0) = 0. (d) Introducing q = dp dt dp dt = q, p(0) = 1, = d2 y gives the following coupled system. 2 = d2 y gives the following coupled system. 2 = d2 y gives the following system. 2 dq so dt = d2 p gives the following system. dt2 dq dt = q, q(0) = 2. p+t (e) Finally, introducing x = dy dz and w = gives the following system dt dt of four equations with associated initial conditions. dy dt = x, y(t 0) = y 0, dt = y +z +2w, x(t 0) = v A, dz dt = w, z(t 0) = z 0, dw dt = y z 2x, w(t 0) = v B. 8
9 12. Consider the coupled initial value problem dt = 3y, dy dt = 3x, x(0) = 1, y(0) = 0. Approximate the solution of this system using the Runge Kutta method with step size h = 0.1. Show that at the first step, the K values are (0,3), ( 0.45,3), ( 0.45,2.9325), ( ,2.865) giving new values x 1 and y 1 of and Show that the second step gives values x 2 and y 2 of and Compare these results with the exact values given by x = cos3t, y = sin3t. Running (rk4) with h = h = 0.1 t = 0.00 x 0 = y 0 = K f1 = K g1 = first temporary point t = 0.05 x m = y m = K f2 = K g2 = second temporary point t = 0.05 x m = y m = K f3 = K g3 = third temporary point t = 0.10 x = y = K f4 = K g4 = new point t = 0.10 x 1 = y 1 = h = 0.1 t = 0.10 x 1 = y 1 = K f1 = K g1 = first temporary point t = 0.15 x m = y m = K f2 = K g2 = second temporary point t = 0.15 x m = y m = K f3 = K g3 = third temporary point t = 0.20 x = y = K f4 = K g4 = new point t = 0.20 x 2 = y 2 = The exact solutions are the functions x(t) = cos3t, y(t) = sin3t. Thus when when t = 0.1, we have x = and y = and when t = 0.2 we have x = and y = We see that the numerical approximations are correct to 4 but not 5 decimal places. 13. Consider the coupled initial value problem dt = 1 3 y, dy dt = 6x3, x(0) = 1; y(0) = 3. Compute an approximation to the solution vector (x,y) at the final time t = 0.5 using the Runge Kutta method with step size h = 0.1. Then repeat the calculation for h = Compare the numerical results with the exact values given by x = 1/(t+1), y = 3/(t+1) 2. Editing rhs and running rk4 gives the following output: >> [time, zzrk4, nfeval] = rk4([1;-3],0.5,5); RK4 iteration in 2 dimensions... t x y all done 9
10 >> [time, zzrk4, nfeval] = rk4([1;-3],0.5,10); RK4 iteration in 2 dimensions... t x y all done Theexactvaluesarex(0.5) = andy(0.5) = The approximation errors for h = 0.1 are x(0.5) x 5 = and y(0.5) y 5 = For the smaller value h = 0.05, we have x(0.5) x 10 = and y(0.5) y 10 = The numerical results are much more accurate for the smaller value of h. 14. Consider the second-order initial value problem d 2 y dt 2 +y +ǫy3 = 0; y(0) = 1, dy dt (0) = 0 with the value ǫ = 0.1. (This is an example of a Duffing oscillator.) Compute three steps of Runge Kutta with a step size h = 0.1. The first step is to express the second-order problem as a coupled firstorder system. dy dt = v; y(0) = 1, dv dt = y 0.1y3 ; v(0) = 0. Editing rhs and running rk4 gives the following output: >> [time, zzrk4, nfeval] = rk4([1;0],0.3,3); RK4 iteration in 2 dimensions... t x y all done 10
11 15. A sphere floating in water (in the absence of friction and with its centre y above the water level) experiences a force due to gravity of 4 3 πr3 ρg and a force due to buoyancy of 1 3 ρ waterπ ( 2r 3 +y(y 2 3r 2 ) ) g. It s motion is determined by the ODE d 2 y dt 2 = g+gρ water 4ρ ] [2+ y3 r 3 3y r A 0.4 metre radius sphere with density 0.7 of that of water is initially at rest with its centre 25 cm below the level of the water. You may assume that g = 10 ms 2. (a) Express the equation of motion as a system of two first-order equations and state the appropriate initial conditions. (b) Take a step size of h = 0.1 and compute two steps of Runge Kutta to estimate the position of the centre of the sphere when t = 0.1 and when t = 0.2. The first step is to express the second-order problem as a coupled firstorder system: dy dt = v with y(0) = 0.25 dv dt = g +gρ water 4ρ [2+ y3 r 3 3y r ] with v(0) = 0 Next, we substitute the given values for the parameters to give dy dt dv 25 = 10+ dt 7 with v(0) = 0 = v with y(0) = 0.25 [ y 3 7.5y ] = y y Finally, editing rhs and running rk4 gives the following output: >> [time, zzrk4, nfeval] = rk4([-0.25;0],0.2,2); RK4 iteration in 2 dimensions... t x y all done The sphere is quickly floating back to the surface! Running rk4 for longer times suggests that the sphere tends to bob up and down in the water. This is illustrated by the following output. 11
12 >> [time, zzrk4, nfeval] = rk4([-0.25;0],10,100); RK4 iteration in 2 dimensions... t x y all done >> plot(time,zzrk4(1,:), -kx ), axis square 16. Compute a centered difference approximation to the BVP x 2 d2 y 2 4xdy +4y = x3, x (1,2); y(1) = 1, y(2) = 6, using a grid of size h = Compare your result with the exact solution values y(1.25) = 4.771, y(1.5) = , y(1.75) = We start by dividing the interval into 4 pieces so that x 0 = 1, x 1 = 1.25, x 2 = 1.5, x 3 = 1.75, x 4 = 2, and then we define the grid values y 0 = y(x 0 ) = 1; y 1 y(x 1 ), y 2 y(x 2 ), y 3 y(x 3 ), y 4 = y(x 4 ) = 6. 12
13 We also define the two finite difference approximations y (x j ) = 1 h 2(y j 1 2y j +y j+1 ), y (x j ) = 1 2h (y j+1 y j 1 ). To simplify notation define y 1 = A, y 2 = B and y 3 = C. The difference approximations at the point x 1 = 1.25 are y (x j ) B = 2B 2 and y (x j ) B 2A+1 = 32A + 16B Substituting these expressions into the ODE gives ( 32A+16B +16) (2B 2)+4A = that is, 46A+35B = The difference approximations at the point x 2 = 1.5 are y (x j ) C A 0.5 = 2C 2A and y (x j ) C 2B+A 0.25 = 16A 32B +16C. Substituting these 2 expressions into the ODE gives (16A 32B +16C)+4 1.5(2C 2A)+4B = that is, 24A 68B +48C = The difference approximations at the point x 3 = 1.75 are y (x j ) 6 B 0.5 = 12 2B and y (x j ) 6 2C+B = 96 32C + 16B. Substituting these expressions into the ODE gives (96 32C +16B) (12 2B)+4C = that is 35B 94C = Solving the 3 linear equations gives the results A = (= y 1 ), B = (= y 2 ) and C = (= y 3 ). Comparing with the exact values we see that the errors are of order of Compute a centered difference approximation to the BVP d 2 y 2 + π2 4 y = ǫy3, x (0,1); y(0) = 1, y(1) = 0, for thevalueǫ = 0.1, usingagrid ofsize h = 0.5. Compareyour computed result with the exact solution at the interior grid point. Dividing the interval into 2 pieces we introduce the unknown grid point value y 1 y(0.5). To simplify notation we let A = y 1. The difference approximations at the point x 1 = 0.5 are y (x j ) = 1 and y (x j ) 0 2A+1 = 4 8A. Substituting these expressions into the ODE gives 4 8A+(π 2 /4)A = 0.1A 3. This is a nonlinear equation. To solve it we set up a Newton iteration with f(a) = 4+(π 2 /4 8)A 0.1A 3 so that f (A) = π 2 / A 2. The iteration converges to the solution A = (= y 1 ). A f(a) f (A)
14 The exact solution solution to the BVP when ǫ = 0 is the function y(x) = cos( πx 2 ) so that y(0.5) = The difference approximation is surprisingly accurate. 18. An iron bar is kept at temperatures T 0 and T 1 at its two ends. The radiation of heat at point x is proportional its temperature AT(x) and this effect is balanced by thermal diffusion, that is, κ d2 T. Thus, when 2 the temperature T(x) has reached equilibrium, it satisfies the BVP κ d2 T 2 = AT, x (0,1); T(0) = T 0, T(1) = T 1. Write down a centered difference approximation to T(x) for values of x of 0.25, 0.5 and 0.75 and solve the resulting linear system in the case where A = 1, κ = 1, T 0 = 1 and T 1 = 2. Multiplying by minus one and substituting parameter values gives d 2 T T = 0 witht(0) = 1 andt(1) = 2. 2 Let s definegridvaluessothatt(0.25) P, T(0.5) QandT(0.75) R. Thedifferenceapproximationatthepointx 1 = 0.25isT (x j ) Q 2P = 16Q 32P +16. Substituting this expression into the ODE gives 16Q 33P +16 = 0. (a) Thedifferenceapproximationatthepointx 2 = 0.5isT (x j ) R 2Q+P = 16R 32Q+16P. Substituting this expression into the ODE gives 16P 33Q+16R = 0. (b) Thedifferenceapproximationatthepointx 3 = 0.75isT (x j ) 2 2R+Q = 32 32R+16Q. Substituting this expression into the ODE gives 16Q 33R+32 = 0. (c) Solving the 3 linear equations (a) (c) gives the results P = , Q = and R = Note that the exact solution to the BVP is T = 2e 1 e 2 1 ex + e2 2e e 2 1 e x. Thustheexact values aret(0.25) = , T(0.5) = , T(0.75) = so the numerical values are in good agreement. 14
15 19. Imagine two metal bars, one with ends maintained at temperature T = 600 and the other with a temperature maintained at T = 400. As well as thermal diffusion and general radiation, there is heat transferred from one bar to the other. The bars thus satisfy a coupled system of ODEs 0.01 d2 T T (T 1 T 2 ) = 0, 0.01 d2 T T (T 2 T 1 ) = 0, subject to T 1 (0) = T 1 (1) = 600 and T 2 (0) = T 2 (1) = 400. Take a grid with points x 0 = 0, x 1 = 0.25, x 2 = 0.5, x 3 = 0.75 and x 4 = 1 and compute a centered difference approximation to T 1 and T 2. First we introduce new functions y = 0.01T 1 and z = 0.01T 2 and substitute into the ODEs. This gives the rescaled equations: y 1.5y +0.5z = 0 z +0.5y 1.5z = 0 subject to y(0) = y(1) = 6 and z(0) = z(1) = 4. Let s define grid values so that A y(0.25), B y(0.5), C y(0.75), D z(0.25), E z(0.5) and F z(0.75). Substituting the second derivative difference approximation into the ODEs evaluated at the grid points gives a system of 6 algebraic equations: 6 2A+B A+0.5D = 0 67A+32B +D +192 = 0 C 2B +A B +0.5E = 0 32A 67B +32C +E = 0 B 2C C +0.5F = 0 32B 67C +F +192 = 0 4 2D +E D +0.5A = 0 A 67D +32E +128 = 0 D 2E +F E +0.5B = 0 B +32D 67E +32F = 0 E 2F F +0.5C = 0 C +32E 67F +128 = 0 Solving these equations gives A = 5.42, B = 5.23, C = 5.42, D = 3.73, E = 3.64, F = The corresponding temperature estimates are x = 0 x = 0.25 x = 0.5 x = 0.75 x = 1 Bar Bar
APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified
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 informationPhysics 115/242 Comparison of methods for integrating the simple harmonic oscillator.
Physics 115/4 Comparison of methods for integrating the simple harmonic oscillator. Peter Young I. THE SIMPLE HARMONIC OSCILLATOR The energy (sometimes called the Hamiltonian ) of the simple harmonic oscillator
More informationChap. 20: Initial-Value Problems
Chap. 20: Initial-Value Problems Ordinary Differential Equations Goal: to solve differential equations of the form: dy dt f t, y The methods in this chapter are all one-step methods and have the general
More informationIntegration is the reverse of the process of differentiation. In the usual notation. k dx = kx + c. kx dx = 1 2 kx2 + c.
PHYS122 - Electricity and Magnetism Integration Reminder Integration is the reverse of the process of differentiation. In the usual notation f (x)dx = f(x) + constant The derivative of the RHS gives you
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 information. For each initial condition y(0) = y 0, there exists a. unique solution. In fact, given any point (x, y), there is a unique curve through this point,
1.2. Direction Fields: Graphical Representation of the ODE and its Solution Section Objective(s): Constructing Direction Fields. Interpreting Direction Fields. Definition 1.2.1. A first order ODE of the
More informationPartial Derivatives. w = f(x, y, z).
Partial Derivatives 1 Functions of Several Variables So far we have focused our attention of functions of one variable. These functions model situations in which a variable depends on another independent
More information4 Partial Differentiation
4 Partial Differentiation Many equations in engineering, physics and mathematics tie together more than two variables. For example Ohm s Law (V = IR) and the equation for an ideal gas, PV = nrt, which
More informationFirst order differential equations
First order differential equations Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra by Goode and Annin Samy T. First
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 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 informationProblem set 7 Math 207A, Fall 2011 Solutions
Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase
More informationSolution: (a) Before opening the parachute, the differential equation is given by: dv dt. = v. v(0) = 0
Math 2250 Lab 4 Name/Unid: 1. (25 points) A man bails out of an airplane at the altitute of 12,000 ft, falls freely for 20 s, then opens his parachute. Assuming linear air resistance ρv ft/s 2, taking
More informationAP Calculus Testbank (Chapter 6) (Mr. Surowski)
AP Calculus Testbank (Chapter 6) (Mr. Surowski) Part I. Multiple-Choice Questions 1. Suppose that f is an odd differentiable function. Then (A) f(1); (B) f (1) (C) f(1) f( 1) (D) 0 (E). 1 1 xf (x) =. The
More informationProblem Set. Assignment #1. Math 3350, Spring Feb. 6, 2004 ANSWERS
Problem Set Assignment #1 Math 3350, Spring 2004 Feb. 6, 2004 ANSWERS i Problem 1. [Section 1.4, Problem 4] A rocket is shot straight up. During the initial stages of flight is has acceleration 7t m /s
More informationDefinition of differential equations and their classification. Methods of solution of first-order differential equations
Introduction to differential equations: overview Definition of differential equations and their classification Solutions of differential equations Initial value problems Existence and uniqueness Mathematical
More informationCompound Damped Pendulum: An Example
Compound Damped Pendulum: An Example Temple H. Fay Department of Mathematical Technology 1 Tshwane University of Technology Pretoria, South Africa thf ay@hotmail:com Abstract: In this article, we use an
More informationENGI Partial Differentiation Page y f x
ENGI 344 4 Partial Differentiation Page 4-0 4. Partial Differentiation For functions of one variable, be found unambiguously by differentiation: y f x, the rate of change of the dependent variable can
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationM343 Homework 3 Enrique Areyan May 17, 2013
M343 Homework 3 Enrique Areyan May 17, 013 Section.6 3. Consider the equation: (3x xy + )dx + (6y x + 3)dy = 0. Let M(x, y) = 3x xy + and N(x, y) = 6y x + 3. Since: y = x = N We can conclude that this
More informationGeometry and Motion Selected answers to Sections A and C Dwight Barkley 2016
MA34 Geometry and Motion Selected answers to Sections A and C Dwight Barkley 26 Example Sheet d n+ = d n cot θ n r θ n r = Θθ n i. 2. 3. 4. Possible answers include: and with opposite orientation: 5..
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationDifferential equations
Differential equations Math 27 Spring 2008 In-term exam February 5th. Solutions This exam contains fourteen problems numbered through 4. Problems 3 are multiple choice problems, which each count 6% of
More informationOrdinary Differential Equations (ODEs)
Ordinary Differential Equations (ODEs) NRiC Chapter 16. ODEs involve derivatives wrt one independent variable, e.g. time t. ODEs can always be reduced to a set of firstorder equations (involving only first
More informationSixth Term Examination Papers 9475 MATHEMATICS 3
Sixth Term Examination Papers 9475 MATHEMATICS 3 Morning WEDNESDAY 26 JUNE 2013 Time: 3 hours Additional Materials: Answer Booklet Formulae Booklet INSTRUCTIONS TO CANDIDATES Please read this page carefully,
More informationTP 3:Runge-Kutta Methods-Solar System-The Method of Least Squares
TP :Runge-Kutta Methods-Solar System-The Method of Least Squares December 8, 2009 1 Runge-Kutta Method The problem is still trying to solve the first order differential equation dy = f(y, x). (1) dx In
More informationPolytechnic Institute of NYU MA 2132 Final Practice Answers Fall 2012
Polytechnic Institute of NYU MA Final Practice Answers Fall Studying from past or sample exams is NOT recommended. If you do, it should be only AFTER you know how to do all of the homework and worksheet
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More informationMath 308, Sections 301, 302, Summer 2008 Review before Test I 06/09/2008
Math 308, Sections 301, 302, Summer 2008 Review before Test I 06/09/2008 Chapter 1. Introduction Section 1.1 Background Definition Equation that contains some derivatives of an unknown function is called
More informationdx n a 1(x) dy
HIGHER ORDER DIFFERENTIAL EQUATIONS Theory of linear equations Initial-value and boundary-value problem nth-order initial value problem is Solve: a n (x) dn y dx n + a n 1(x) dn 1 y dx n 1 +... + a 1(x)
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 informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationSMA 208: Ordinary differential equations I
SMA 208: Ordinary differential equations I First Order differential equations Lecturer: Dr. Philip Ngare (Contacts: pngare@uonbi.ac.ke, Tue 12-2 PM) School of Mathematics, University of Nairobi Feb 26,
More informationBasic Theory of Differential Equations
page 104 104 CHAPTER 1 First-Order Differential Equations 16. The following initial-value problem arises in the analysis of a cable suspended between two fixed points y = 1 a 1 + (y ) 2, y(0) = a, y (0)
More information1 Introduction IPICSE-2016
(06) DOI: 0.05/ matecconf/06860006 IPICSE-06 Numerical algorithm for solving of nonlinear problems of structural mechanics based on the continuation method in combination with the dynamic relaxation method
More information1. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal
. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal 3 9 matrix D such that A = P DP, for A =. 3 4 3 (a) P = 4, D =. 3 (b) P = 4, D =. (c) P = 4 8 4, D =. 3 (d) P
More informationCoupled differential equations
Coupled differential equations Example: dy 1 dy 11 1 1 1 1 1 a y a y b x a y a y b x Consider the case with b1 b d y1 a11 a1 y1 y a1 ay dy y y e y 4/1/18 One way to address this sort of problem, is to
More informationChapter1. Ordinary Differential Equations
Chapter1. Ordinary Differential Equations In the sciences and engineering, mathematical models are developed to aid in the understanding of physical phenomena. These models often yield an equation that
More informationFirst-Order Differential Equations
CHAPTER 1 First-Order Differential Equations 1. Diff Eqns and Math Models Know what it means for a function to be a solution to a differential equation. In order to figure out if y = y(x) is a solution
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
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 informationCalculus Math 21B, Winter 2009 Final Exam: Solutions
Calculus Math B, Winter 9 Final Exam: Solutions. (a) Express the area of the region enclosed between the x-axis and the curve y = x 4 x for x as a definite integral. (b) Find the area by evaluating the
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationAdvection, Conservation, Conserved Physical Quantities, Wave Equations
EP711 Supplementary Material Thursday, September 4, 2014 Advection, Conservation, Conserved Physical Quantities, Wave Equations Jonathan B. Snively!Embry-Riddle Aeronautical University Contents EP711 Supplementary
More informationNumerical Methods for Initial Value Problems; Harmonic Oscillators
Lab 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 information1.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 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 informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More informationMB4018 Differential equations
MB4018 Differential equations Part II http://www.staff.ul.ie/natalia/mb4018.html Prof. Natalia Kopteva Spring 2015 MB4018 (Spring 2015) Differential equations Part II 0 / 69 Section 1 Second-Order Linear
More informationMAT 211 Final Exam. Spring Jennings. Show your work!
MAT 211 Final Exam. pring 215. Jennings. how your work! Hessian D = f xx f yy (f xy ) 2 (for optimization). Polar coordinates x = r cos(θ), y = r sin(θ), da = r dr dθ. ylindrical coordinates x = r cos(θ),
More informationSimple ODE Solvers - Derivation
Simple ODE Solvers - Derivation These notes provide derivations of some simple algorithms for generating, numerically, approximate solutions to the initial value problem y (t =f ( t, y(t y(t 0 =y 0 Here
More informationMath 34B. Practice Exam, 3 hrs. March 15, 2012
Math 34B Practice Exam, 3 hrs March 15, 2012 9.3.4c Compute the indefinite integral: 10 x+9 dx = 9.3.4c Compute the indefinite integral: 10 x+9 dx = = 10 x 10 9 dx = 10 9 10 x dx = 10 9 e ln 10x dx 9.3.4c
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationMultistep Methods for IVPs. t 0 < t < T
Multistep Methods for IVPs We are still considering the IVP dy dt = f(t,y) t 0 < t < T y(t 0 ) = y 0 So far we have looked at Euler s method, which was a first order method and Runge Kutta (RK) methods
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationModels of ocean circulation are all based on the equations of motion.
Equations of motion Models of ocean circulation are all based on the equations of motion. Only in simple cases the equations of motion can be solved analytically, usually they must be solved numerically.
More informationSME 3023 Applied Numerical Methods
UNIVERSITI TEKNOLOGI MALAYSIA SME 3023 Applied Numerical Methods Ordinary Differential Equations Abu Hasan Abdullah Faculty of Mechanical Engineering Sept 2012 Abu Hasan Abdullah (FME) SME 3023 Applied
More information(1) Rate of change: A swimming pool is emptying at a constant rate of 90 gal/min.
CHAPTER 1 Introduction 1. Bacground Models of physical situations from Calculus (1) Rate of change: A swimming pool is emptying at a constant rate of 90 gal/min. With V = volume in gallons and t = time
More informationMATH 32A: MIDTERM 1 REVIEW. 1. Vectors. v v = 1 22
MATH 3A: MIDTERM 1 REVIEW JOE HUGHES 1. Let v = 3,, 3. a. Find e v. Solution: v = 9 + 4 + 9 =, so 1. Vectors e v = 1 v v = 1 3,, 3 b. Find the vectors parallel to v which lie on the sphere of radius two
More informationBoiling Heat Transfer and Two-Phase Flow Fall 2012 Rayleigh Bubble Dynamics. Derivation of Rayleigh and Rayleigh-Plesset Equations:
Boiling Heat Transfer and Two-Phase Flow Fall 2012 Rayleigh Bubble Dynamics Derivation of Rayleigh and Rayleigh-Plesset Equations: Let us derive the Rayleigh-Plesset Equation as the Rayleigh equation can
More informationThe Fundamental Theorem of Calculus: Suppose f continuous on [a, b]. 1.) If G(x) = x. f(t)dt = F (b) F (a) where F is any antiderivative
1 Calulus pre-requisites you must know. Derivative = slope of tangent line = rate. Integral = area between curve and x-axis (where area can be negative). The Fundamental Theorem of Calculus: Suppose f
More informationIntermediate Tier - Algebra revision
Intermediate Tier - Algebra revision Contents : Collecting like terms Multiplying terms together Indices Expanding single brackets Expanding double brackets Substitution Solving equations Finding nth term
More informationIntroduction to First Order Equations Sections
A B I L E N E C H R I S T I A N U N I V E R S I T Y Department of Mathematics Introduction to First Order Equations Sections 2.1-2.3 Dr. John Ehrke Department of Mathematics Fall 2012 Course Goals The
More informationWorksheet 8 Sample Solutions
Technische Universität München WS 2016/17 Lehrstuhl für Informatik V Scientific Computing Univ.-Prof. Dr. M. Bader 19.12.2016/21.12.2016 M.Sc. S. Seckler, M.Sc. D. Jarema Worksheet 8 Sample Solutions Ordinary
More informationFirst order differential Equations
Chapter 10 First order differential Equations 10.1 What is a Differential Equation? A differential equation is an equation involving an unknown function and its derivatives. A general differential equation
More informationOrdinary Differential Equations (ODEs)
c01.tex 8/10/2010 22: 55 Page 1 PART A Ordinary Differential Equations (ODEs) Chap. 1 First-Order ODEs Sec. 1.1 Basic Concepts. Modeling To get a good start into this chapter and this section, quickly
More informationMath 266: Ordinary Differential Equations
Math 266: Ordinary Differential Equations Long Jin Purdue University, Spring 2018 Basic information Lectures: MWF 8:30-9:20(111)/9:30-10:20(121), UNIV 103 Instructor: Long Jin (long249@purdue.edu) Office
More informationCalculus III. Math 233 Spring Final exam May 3rd. Suggested solutions
alculus III Math 33 pring 7 Final exam May 3rd. uggested solutions This exam contains twenty problems numbered 1 through. All problems are multiple choice problems, and each counts 5% of your total score.
More informationSOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)
SOLUTIONS TO THE 18.02 FINAL EXAM BJORN POONEN December 14, 2010, 9:00am-12:00 (3 hours) 1) For each of (a)-(e) below: If the statement is true, write TRUE. If the statement is false, write FALSE. (Please
More informationSKMM 3023 Applied Numerical Methods
UNIVERSITI TEKNOLOGI MALAYSIA SKMM 3023 Applied Numerical Methods Ordinary Differential Equations ibn Abdullah Faculty of Mechanical Engineering Òº ÙÐÐ ÚºÒÙÐÐ ¾¼½ SKMM 3023 Applied Numerical Methods Ordinary
More informationOrdinary Differential Equations
Chapter 13 Ordinary Differential Equations We motivated the problem of interpolation in Chapter 11 by transitioning from analzying to finding functions. That is, in problems like interpolation and regression,
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationYear 8. Semester 2 Revisions
Semester Revisions Year 8 Semester Two Revisions 01 Students are advised to review previously completed tests, assignments and homework followed by this semester one revision booklet. Be sure to seek extra
More informationExercises for Multivariable Differential Calculus XM521
This document lists all the exercises for XM521. The Type I (True/False) exercises will be given, and should be answered, online immediately following each lecture. The Type III exercises are to be done
More informationODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class
ODE Homework Due Wed. 9 August 2009; At the beginning of the class. (a) Solve Lẏ + Ry = E sin(ωt) with y(0) = k () L, R, E, ω are positive constants. (b) What is the limit of the solution as ω 0? (c) Is
More informationENGI Partial Differentiation Page y f x
ENGI 3424 4 Partial Differentiation Page 4-01 4. Partial Differentiation For functions of one variable, be found unambiguously by differentiation: y f x, the rate of change of the dependent variable can
More informationSect2.1. Any linear equation:
Sect2.1. Any linear equation: Divide a 0 (t) on both sides a 0 (t) dt +a 1(t)y = g(t) dt + a 1(t) a 0 (t) y = g(t) a 0 (t) Choose p(t) = a 1(t) a 0 (t) Rewrite it in standard form and ḡ(t) = g(t) a 0 (t)
More informationME3.6 Sheet 2 Answers
ME36 Sheet 2 Answers (i) We have dx/dt Ý6x 2xy Ý 8 fx,y,dy/dt y 2 Ý x 2 gx,y The critical points satisfy g y x and f Ý6x 2xy Ý 8 Substitute y x to get Ý6x 2x 2 Ý 8 x Ý 4x x 4 or Ý If instead we substitute
More informationMcGill University April 16, Advanced Calculus for Engineers
McGill University April 16, 2014 Faculty of cience Final examination Advanced Calculus for Engineers Math 264 April 16, 2014 Time: 6PM-9PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer
More informationPractice Exam 1 Solutions
Practice Exam 1 Solutions 1a. Let S be the region bounded by y = x 3, y = 1, and x. Find the area of S. What is the volume of the solid obtained by rotating S about the line y = 1? Area A = Volume 1 1
More informationSample Final Questions: Solutions Math 21B, Winter y ( y 1)(1 + y)) = A y + B
Sample Final Questions: Solutions Math 2B, Winter 23. Evaluate the following integrals: tan a) y y dy; b) x dx; c) 3 x 2 + x dx. a) We use partial fractions: y y 3 = y y ) + y)) = A y + B y + C y +. Putting
More informationMAT 275 Test 1 SOLUTIONS, FORM A
MAT 75 Test SOLUTIONS, FORM A The differential equation xy e x y + y 3 = x 7 is D neither linear nor homogeneous Solution: Linearity is ruinied by the y 3 term; homogeneity is ruined by the x 7 on the
More informationSolutions to the Review Questions
Solutions to the Review Questions Short Answer/True or False. True or False, and explain: (a) If y = y + 2t, then 0 = y + 2t is an equilibrium solution. False: (a) Equilibrium solutions are only defined
More informationSection 7.4 Runge-Kutta Methods
Section 7.4 Runge-Kutta Methods Key terms: Taylor methods Taylor series Runge-Kutta; methods linear combinations of function values at intermediate points Alternatives to second order Taylor methods Fourth
More informationFirst Year Physics: Prelims CP1. Classical Mechanics: Prof. Neville Harnew. Problem Set III : Projectiles, rocket motion and motion in E & B fields
HT017 First Year Physics: Prelims CP1 Classical Mechanics: Prof Neville Harnew Problem Set III : Projectiles, rocket motion and motion in E & B fields Questions 1-10 are standard examples Questions 11-1
More informationMTHE 227 Problem Set 2 Solutions
MTHE 7 Problem Set Solutions 1 (Great Circles). The intersection of a sphere with a plane passing through its center is called a great circle. Let Γ be the great circle that is the intersection of the
More informationMathematics for chemical engineers. Numerical solution of ordinary differential equations
Mathematics for chemical engineers Drahoslava Janovská Numerical solution of ordinary differential equations Initial value problem Winter Semester 2015-2016 Outline 1 Introduction 2 One step methods Euler
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
More informationIntroduction to Differential Equations
Chapter 1 Introduction to Differential Equations 1.1 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 informationDifferential Equations
Differential Equations Definitions Finite Differences Taylor Series based Methods: Euler Method Runge-Kutta Methods Improved Euler, Midpoint methods Runge Kutta (2nd, 4th order) methods Predictor-Corrector
More informationRunga-Kutta Schemes. Exact evolution over a small time step: Expand both sides in a small time increment: d(δt) F (x(t + δt),t+ δt) Ft + FF x ( t)
Runga-Kutta Schemes Exact evolution over a small time step: x(t + t) =x(t)+ t 0 d(δt) F (x(t + δt),t+ δt) Expand both sides in a small time increment: x(t + t) =x(t)+x (t) t + 1 2 x (t)( t) 2 + 1 6 x (t)+
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Linear Multistep methods Gustaf Söderlind and Carmen Arévalo Numerical Analysis, Lund University Textbooks: A First Course in the
More informationLAURENTIAN UNIVERSITY UNIVERSITÉ LAURENTIENNE
Page 1 of 15 LAURENTIAN UNIVERSITY UNIVERSITÉ LAURENTIENNE Wednesday, December 13 th Course and No. 2006 MATH 2066 EL Date................................... Cours et no........................... Total
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 informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 1. Runge-Kutta Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Mathematics National Institute of Technology Durgapur Durgapur-71309 email: anita.buie@gmail.com 1 . Chapter 8 Numerical Solution of
More informationIntroductory Differential Equations
Introductory Differential Equations Lecture Notes June 3, 208 Contents Introduction Terminology and Examples 2 Classification of Differential Equations 4 2 First Order ODEs 5 2 Separable ODEs 5 22 First
More informationOrdinary Differential Equations
Ordinary Differential Equations Part 4 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Ordinary Differential Equations p. 1/23 Two-point Boundary Value Problems NRiC 17. BCs specified at two
More information