# Ordinary Differential Equations

Save this PDF as:

Size: px
Start display at page:

## Transcription

1 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, the unknown is a function f, and the job of the algorithm is to fill in missing data. We continue this discussion by considering similar problems involving filling in function values. Here, our unknown continues to be a function f, but rather than simply guessing missing values we would like to solve more complex design problems. For example, consider the following problems: Find f approximating some other function f 0 but satisfying additional criteria (smoothness, continuity, low-frequency, etc.). Simulate some dynamical or physical relationship as f (t) where t is time. Find f with similar values to f 0 but certain properties in common with a different function g 0. In each of these cases, our unknown is a function f, but our criteria for success is more involved than matches a given set of data points. The theories of ordinary differential equations (ODEs) and partial differential equations (PDEs) study the case where we wish to find a function f ( x) based on information about or relationships between its derivatives. Notice we already have solved a simple version of this problem in our discussion of quadrature: Given f (t), methods for quadrature provide ways of approximating f (t) using integration. In this chapter, we will consider the case of an ordinary differential equations and in particular initial value problems. Here, the unknown is a function f (t) : R R n, and we given an equation satisfied by f and its derivatives as well as f (0); our goal is to predict f (t) for t > 0. We will provide several examples of ODEs appearing in the computer science literature and then will proceed to describe common solution techniques. As a note, we will use the notation f to denote the derivative d f /dt of f : [0, ) R n. Our goal will be to find f (t) given relationships between t, f (t), f (t), f (t), and so on. 1

2 13.1 Motivation ODEs appear in nearly any part of scientific example, and it is not difficult to encounter practical situations requiring their solution. For instance, the basic laws of physical motion are given by an ODE: Example 13.1 (Newton s Second Law). Continuing from 5.1.2, recall that Newton s Second Law of Motion states that F = ma, that is, the total force on an object is equal to its mass times its acceleration. If we simulate n particles simultaneously, then we can think of combining all their positions into a vector x R 3n. Similarly, we can write a function F(t, x, x ) R 3n taking the time, positions of the particles, and their velocities and returning the total force on each particle. This function can take into account interrelationships between particles (e.g. gravitational forces or springs), external effects like wind resistance (which depends on x ), external forces varying with time t, and so on. Then, to find the positions of all the particles as functions of time, we wish to solve the equation x = F(t, x, x )/m. We usually are given the positions and velocities of all the particles at time t = 0 as a starting condition. Example 13.2 (Protein folding). On a smaller scale, the equations governing motions of molecules also are ordinary differential equations. One particularly challenging case is that of protein folding, in which the geometry structure of a protein is predicted by simulating intermolecular forces over time. These forces take many often nonlinear forms that continue to challenge researchers in computational biology. Example 13.3 (Gradient descent). Suppose we are wishing to minimize an energy function E( x) over all x. We learned in Chapter 8 that E( x) points in the direction E decreases the most at a given x, so we did line search along that direction from x to minimize E locally. An alternative option popular in certain theories is to solve an ODE of the form x = E( x); in other words, think of x as a function of time x(t) that is attempting to decrease E by walking downhill. For example, suppose we wish to solve A x = b for symmetric positive definite A. We know from that this is equivalent to minimizing E( x) 1 2 x A x b x + c. Thus, we could attempt to solve the ODE x = f ( x) = b A x. As t, we expect x(t) to better and better satisfy the linear system. Example 13.4 (Crowd simulation). Suppose we are writing video game software requiring realistic simulation of virtual crowds of humans, animals, spaceships, and the like. One strategy for generating plausible motion, illustrated in Figure NUMBER, is to use differential equations. Here, the velocity of a member of the crowd is determined as a function of its environment; for example, in human crowds the proximity of other humans, distance to obstacles, and so on can affect the direction a given agent is moving. These rules can be simple, but in the aggregate their interaction is complex. Stable integrators for differential equations underlie this machinery, since we do not wish to have noticeably unrealistic or unphysical behavior Theory of ODEs A full treatment of the theory of ordinary differential equations is outside the scope of our discussion, and we refer the reader to CITE for more details. This aside, we mention some highlights here that will be relevant to our development in future sections. 2

3 Basic Notions The most general ODE initial value problem takes the following form: Find f (t) : R R n Satisfying F[t, f (t), f (t), f (t),..., f (k) (t)] = 0 Given f (0), f (0), f (0),..., f (k 1) (0) Here, F is some relationship between f and all its derivatives; we use f (l) to denote the l-th derivative of f. We can think of ODEs as determining evolution of f over time t; we know f and its derivatives at time zero and wish to predict it moving forward. ODEs take many forms even in a single variable. For instance, denote y = f (t) and suppose y R 1. Then, examples of ODEs include: y = 1 + cos t: This ODE can be solved by integrating both sides e.g. using quadrature methods y = ay: This ODE is linear in y y = ay + e t : This ODE is time and position-dependent y + 3y y = t: This ODE involves multiple derivatives of y y sin y = e ty : This ODE is nonlinear in y and t. Obviously the most general ODEs can be challenging to solve. We will restrict most of our discussion to the case of explicit ODEs, in which the highest-order derivative can be isolated: Definition 13.1 (Explicit ODE). An ODE is explicit if can be written in the form f (k) (t) = F[t, f (t), f (t), f (t),..., f (k 1) (t)]. For example, an explicit form of Newton s second law is x (t) = 1 m a(t, x(t), x (t)). Surprisingly, generalizing the trick in 5.1.2, in fact any explicit ODE can be converted to a first-order equation f (t) = F[t, f (t)], where f has multidimensional output. This observation implies that we need no more than one derivative in our treatment of ODE algorithms. To see this relationship, we simply recall that d2 y/dt 2 = d /dt( dy /dt). Thus, we can define an intermediate variable z dy /dt, and understand d2 y/dt 2 as dz /dt with the constraint z = dy /dt. More generally, if we wish to solve the explicit problem f (k) (t) = F[t, f (t), f (t), f (t),..., f (k 1) (t)], where f : R R n, then we define g(t) : R R kn using the first-order ODE: d dt g 1 (t) g 2 (t). g k 1 (t) g k (t) = g 2 (t) g 3 (t). g k (t) F[t, g 1 (t), g 2 (t),..., g k 1 (t)] 3

4 Here, we denote g i (t) : R R n to contain n components of g. Then, g 1 (t) satisfies the original ODE. To see so, we just check that our equation above implies g 2 (t) = g 1 (t), g 3(t) = g 2 (t) = g 1 (t), and so on. Thus, making these substitutions shows that the final row encodes the original ODE. The trick above will simplify our notation, but some care should be taken to understand that this approach does not trivialize computations. In particular, in may cases our function f (t) will only have a single output, but the ODE will be in several derivatives. We replace this case with one derivative and several outputs. Example 13.5 (ODE expansion). Suppose we wish to solve y = 3y 2y + y where y(t) : R R. This equation is equivalent to: y y d z = z dt w w Just as our trick above allows us to consider only first-order ODEs, we can restrict our notation even more to autonomous ODEs. These equations are of the form f (t) = F[ f (t)], that is, F no longer depends on t. To do so, we could define g(t) ( f (t) ḡ(t) Then, we can solve the following ODE for g instead: ( g f (t) = ) ( (t) F[ f (t), ḡ(t)] ḡ = (t) 1 In particular, ḡ(t) = t assuming we take ḡ(0) = 0. It is possible to visualize the behavior of ODEs in many ways, illustrated in Figure NUMBER. For instance, if the unknown f (t) is a function of a single variable, then we can think of F[ f (t)] as providing the slope of f (t), as shown in Figure NUMBER. Alternatively, if f (t) has output in R 2, we no longer can visualize the dependence on time t, but we can draw phase space, which shows the tangent of f (t) at each (x, y) R Existence and Uniqueness Before we can proceed to discretizations of the initial value problem, we should briefly acknowledge that not all differential equation problems are solvable. Furthermore, some differential equations admit multiple solutions. Example 13.6 (Unsolvable ODE). Consider the equation y = 2y/t, with y(0) = 0 given; notice we are not dividing by zero because y(0) is prescribed. Rewriting as 1 dy y dt = 2 t and integrating with respect to t on both sides shows: ). ln y = 2 ln t + c Or equivalently, y = Ct 2 for some C R. Notice that y(0) = 0 in this expression, contradicting our initial conditions. Thus, this ODE has no solution with the given initial conditions. 4 ).

5 Example 13.7 (Nonunique solutions). Now, consider the same ODE with y(0) = 0. Consider y(t) given by y(t) = Ct 2 for any C R. Then, y (t) = 2Ct. Thus, 2y t = 2Ct2 t = 2Ct = y (t), showing that the ODE is solved by this function regardless of C. nonunique. Thus, solutions of this problem are Thankfully, there is a rich theory characterizing behavior and stability of solutions to differential equations. Our development in the next chapter will have a stronger set of conditions needed for existence of a solution, but in fact under weak conditions on f it is possible to show that an ODE f (t) = F[ f (t)] has a solution. For instance, one such theorem guarantees local existence of a solution: Theorem 13.1 (Local existence and uniqueness). Suppose F is continuous and Lipschitz, that is, F[ y] F[ x] 2 L y x 2 for some L. Then, the ODE f (t) = F[ f (t)] admits exactly one solution for all t 0 regardless of initial conditions. In our subsequent development, we will assume that the ODE we are attempting to solve satisfies the conditions of such a theorem; this assumption is fairly realistic in that at least locally there would have to be fairly degenerate behavior to break such weak assumptions Model Equations One way to gain intuition for the behavior of ODEs is to examine behavior of solutions to some simple model equations that can be solved in closed form. These equations represent linearizations of more practical equations, and thus locally they model the type of behavior we can expect. We start with ODEs in a single variable. Given our simplifications in , the simplest equation we might expect to work with would be y = F[y], where y(t) : R R. Taking a linear approximation would yield equations of type y = ay + b. Substituting ȳ y + b /a shows: ȳ = y = ay + b = a(ȳ b /a) + b = aȳ. Thus, in our model equations the constant b simply induces a shift, and for our phenomenological study in this section we can assume b = 0. By the argument above, we locally can understand behavior of y = F[y] by studying the linear equation y = ay. In fact, applying standard arguments from calculus shows that y(t) = Ce at. Obviously, there are three cases, illustrated in Figure NUMBER: 1. a > 0: In this case solutions get larger and larger; in fact, if y(t) and ŷ(t) both satisfy the ODE with slightly different starting conditions, as t they diverge. 2. a = 0: The system in this case is solved by constants; solutions with different starting points stay the same distance apart. 3. a < 0: Then, all solutions of the ODE approach 0 as t. 5

6 We say cases 2 and 3 are stable, in the sense that perturbing y(0) slightly yields solutions that get close and close over time; case 1 is unstable, since a small mistake in specifying the input parameter y(0) will be amplified as time t advances. Unstable ODEs generate ill-posed computational problems; without careful consideration we cannot expect numerical methods to generate usable solutions in this case, since even theoretical outputs are so sensitive to perturbations of the input. On the other hand, stable problems are well-posed since small mistakes in y(0) get diminished over time. Advancing to multiple dimensions, we could study the linearized equation y = A y. As explained in 5.1.2, if y 1,, y k are eigenvectors of A with eigenvalues λ 1,..., λ k and y(0) = c 1 y c k y k, then y(t) = c 1 e λ 1t y c k e λ kt y k. In other words, the eigenvalues of A take the place of a in our one-dimensional example. From this result, it is not hard to intuit that a multivariable system is stable exactly when its spectral radius is less than one. In reality we wish to solve y = F[ y] for general functions F. Assuming F is differentiable, we can write F[ y] F[ y 0 ] + J F ( y 0 )( y y 0 ), yielding the model equation above after a shift. Thus, for short periods of time we expect behavior similar to the model equation. Additionally, the conditions in Theorem 13.1 can be viewed as a bound on the behavior of J F, providing a connection to less localized theories of ODE Time-Stepping Schemes We now proceed to describe several methods for solving the nonlinear ODE y = F[ y] for potentially nonlinear functions F. In general, given a time step h, our methods will be used to generate estimates of y(t + h) given y(t). Applying these methods iteratively generates estimates of y 0 y(t), y 1 y(t + h), y 2 y(t + 2h), y 3 y(t + 3h), and so on. Notice that since F has no t dependence the mechanism for generating each additional step is the same as the first, so for the most part we will only need to describe a single step of these methods. We call methods for generating approximations of y(t) integrators, reflecting the fact that they are integrating out the derivatives in the input equation. Of key importance to our consideration is the idea of stability. Just as ODEs can be stable or unstable, so can discretizations. For example, if h is too large, some schemes will accumulate error at an exponential rate; contrastingly, other methods are stable in that even if h is large the solutions will remain bounded. Stability, however, can compete with accuracy; often time stable schemes are bad approximations of y(t), even if they are guaranteed not to have wild behavior Forward Euler Our first ODE strategy comes from our construction of the forward differencing scheme in : F[ y k ] = y (t) = y k+1 y k h + O(h) 6

7 Solving this relationship for y k+1 shows y k+1 = y k + hf[ y k ] + O(h 2 ) y k + hf[ y k ]. Thus, the forward Euler scheme applies the formula on the right to estimate y k+1. It is one of the most efficient strategies for time-stepping, since it simply evaluates F and adds a multiple of the result to y k. For this reason, we call it an explicit method, that is, there is an explicit formula for y k+1 in terms of y k and F. Analyzing the accuracy of this method is fairly straightforward. Notice that our approximation of y k+1 is O(h 2 ), so each step induces quadratic error. We call this error localized truncation error because it is the error induced by a single step; the word truncation refers to the fact that we truncated a Taylor series to obtain this formula. Of course, our iterate y k already may be inaccurate thanks to accumulated truncation errors from previous iterations. If we integrate from t 0 to t with O( 1 /h) steps, then our total error looks like O(h); this estimate represents global truncation error, and thus we usually write that the forward Euler scheme is first-order accurate. The stability of this method requires somewhat more consideration. In our discussion, we will work out the stability of methods in the one-variable case y = ay, with the intuition that similar statements carry over to multidimensional equations by replacing a with the spectral radius. In this case, we know y k+1 = y k + ahy k = (1 + ah)y k. In other words, y k = (1 + ah) k y 0. Thus, the integrator is stable when 1 + ah 1, since otherwise y k exponentially. Assuming a < 0 (otherwise the problem is ill-posed), we can simplify: 1 + ah ah 1 2 ah 0 0 h 2 a, since a < 0 Thus, forward Euler admits a time step restriction for stability given by our final condition on h. In other words, the output of forward Euler can explode even when y = ay is stable if h is not small enough. Figure NUMBER illustrates what happens when this condition is obeyed or violated. In multiple dimensions, we can replace this restriction with an analogous one using the spectral radius of A. For nonlinear ODEs this formula gives a guide for stability at least locally in time; globally h may have to be adjusted if the Jacobian of F becomes worse conditioned Backward Euler Similarly, we could have applied the backward differencing scheme at y k+1 to design an ODE integrator: F[ y k+1 ] = y (t) = y k+1 y k h + O(h) Thus, we solve the following potentially nonlinear system of equations for y k+1 : y k = y k+1 hf[ y k+1 ]. Because we have to solve this equation for y k+1, backward Euler is an implicit integrator. 7

8 This method is first-order accurate like forward Euler by an identical proof. The stability of this method, however, contrasts considerably with forward Euler. Once again considering the model equation y = ay, we write: y k = y k+1 hay k+1 = y k+1 = y k 1 ha. Paralleling our previous argument, backward Euler is stable under the following condition: ha 1 1 ha 1 ha 1 or 1 ha 1 h 2 a or h 0, for a < 0 Obviously we always take h 0, so backward Euler is unconditionally stable. Of course, even if backward Euler is stable it is not necessarily accurate. If h is too large, y k will approach zero far too quickly. When simulating cloth and other physical materials that require lots of high-frequency detail to be realistic, backward Euler may not be an effective choice. Furthermore, we have to invert F[ ] to solve for y k+1. Example 13.8 (Backward Euler). Suppose we wish to solve y = A y for A R n n. Then, to find y k+1 we solve the following system: Trapezoidal Method y k = y k+1 ha y k+1 = y k+1 = (I n n ha) 1 y k. Suppose y k is known at time and y k+1 represents the value at time +1 = + h. Suppose we also know y k+ 1/2 halfway in between these two steps. Then, by our derivation of the centered differencing we know: y k+1 = y k + hf[ y k+ 1/2] + O(h 3 ) From our derivation of the trapezoidal rule: F[ y k+1 ] + F[ y k ] 2 = F[ y k+ 1/2] + O(h 2 ) Substituting this relationship yields our first second-order integration scheme, the trapezoid method for integrating ODEs: y k+1 = y k + h F[ y k+1] + F[ y k ] 2 Like backward Euler, this method is implicit since we must solve this equation for y k+1. Once again carrying out stability analysis on y = ay, we find in this case time steps of the trapezoidal method solve y k+1 = y k ha(y k+1 + y k ) In other words, ( y k = ha ) k ha y 0. 8

9 The method is thus stable when ha ha < 1. It is easy to see that this inequality holds whenever a < 0 and h > 0, showing that the trapezoid method is unconditionally stable. Despite its higher order of accuracy with maintained stability, the trapezoid method, however, has some drawbacks that make it less popular than backward Euler. In particular, consider the ratio R y k+1 y k = ha ha When a < 0, for large enough h this ratio eventually becomes negative; in fact, as h, we have R 1. Thus, as illustrated in Figure NUMBER, if time steps are too large, the trapezoidal method of integration tends to exhibit undesirable oscillatory behavior that is not at all like what we might expect for solutions of y = ay Runge-Kutta Methods A class of integrators can be derived by making the following observation: tk + t y k+1 = y k + y (t) dt by the Fundamental Theorem of Calculus = y k + tk + t F[ y(t)] dt Of course, using this formula outright does not work for formulating a method for time-stepping since we do nonow y(t), but careful application of our quadrature formulae from the previous chapter can generate feasible strategies. For example, suppose we apply the trapezoidal method for integration. Then, we find: y k+1 = y k + h 2 (F[ y k] + F[ y k+1 ]) + O(h 3 ) This is the formula we wrote for the trapezoidal method in If we do not wish to solve for y k+1 implicitly, however, we must find an expression to approximate F[ y k+1 ]. Using Euler s method, however, we know that y k+1 = y k + hf[ y k ] + O(h 2 ). Making this substitution for y k+1 does not affect the order of approximation of the trapezoidal time step above, so we can write: y k+1 = y k + h 2 (F[ y k] + F[ y k + hf[ y k ]]) + O(h 3 ) Ignoring the O(h 3 ) terms yields a new integration strategy known as Heun s method, which is second-order accurate and explicit. 9

10 If we study stability behavior of Heun s method for y = ay for a < 0, we know: Thus, the method is stable when y k+1 = y k + h 2 (ay k + a(y k + hay k )) = (1 + h2 ) a(2 + ha) y k = (1 + ha + 12 ) h2 a 2 y k ha h2 a ha + h 2 a 2 0 The inequality on the right shows h 2, and the one on the left is always true for h > 0 and a a < 0, so the stability condition is h 2 a. Heun s method is an example of a Runge-Kutta method derived by applying quadrature methods to the integral above and substituting Euler steps into F[ ]. Forward Euler is a first-order accurate Runge-Kutta method, and Heun s method is second-order. A popular fourth-order Runge- Kutta method (abbreviated RK4 ) is given by: y k+1 = y k + h 6 ( k1 + 2 k2 + 2 k3 + k4 ) where k1 = F [ y k ] k2 = F [ y k + 12 ] h k1 k3 = F [ y k + 12 ] h k2 ] k4 = F [ y k + h k3 This method can be derived by applying Simpson s rule for quadrature. Runge-Kutta methods are popular because they are explicit and thus easy to evaluate while providing high degrees of accuracy. The cost of this accuracy, however, is that F[ ] must be evaluated more times. Furthermore, Runge-Kutta strategies can be extended to implicit methods that can solve stiff equations Exponential Integrators One class of integrators that achieves strong accuracy when F[ ] is approximately linear is to use our solution to the model equation explicitly. In particular, if we were solving the ODE y = A y, using eigenvectors of A (or any other method) we could find an explicit solution y(t) as explained in We usually write y k+1 = e Ah y k, where e Ah encodes our exponentiation of the eigenvalues (in fact we can find a matrix e Ah from this expression that solves the ODE to time h). Now, if we write y = A y + G[ y], 10

11 where G is a nonlinear but small function, we can achieve fairly high accuracy by integrating the A part explicitly and then approximating the nonlinear G part separately. For example, the first-order exponential integrator applies forward Euler to the nonlinear G term: y k+1 = e Ah y k A 1 (1 e Ah )G[ y k ] Analysis revealing the advantages of this method is more complex than what we have written, but intuitively it is clear that these methods will behave particularly well when G is small Multivalue Methods The transformations in enabled us to simplify notation in the previous section considerably by reducing all explicit ODEs to the form y = F[ y]. In fact, while all explicit ODEs can be written this way, it is not clear that they always should. In particular, when we reduced k-th order ODEs to first-order ODEs, we introduced a number of variables representing the first through k 1-st derivatives of the desired output. In fact, in our final solution we only care about the zeroth derivative, that is, the function itself, so orders of accuracy on the temporary variables are less important. From this perspective, consider the Taylor series y( + h) = y( ) + h y ( ) + h2 2 y ( ) + O(h 3 ) If we only know y up to O(h 2 ), this does not affect our approximation, since y gets multiplied by h. Similarly, if we only know y up to O(h), this approximation will not affect the Taylor series terms above because it will get multiplied by h2 /2. Thus, we now consider multivalue methods, designed to integrate y (k) (t) = F[t, y (t), y (t),..., y (k 1) (t)] with different-order accuracy for different derivatives of the function y. Given the importance of Newton s second law F = ma, we will restrict to the case y = F[t, y, y ]; many extensions exist for the less common k-th order case. We introduce a velocity vector v(t) = y (t) and an acceleration vector a. By our previous reduction, we wish to solve the following first-order system: y (t) = v(t) v (t) = a(t) a(t) = F[t, y(t), v(t)] Our goal is to derive an integrator specifically tailored to this system Newmark Schemes We begin by deriving the famous class of Newmark integrators. 1 Denote y k, v k, and a k as the position, velocity, and acceleration vectors at time ; our goal is to advance to time +1 + h. 1 We follow the development in pdf. 11

12 Use y(t), v(t), and a(t) to denote the functions of time assuming we start at. Then, obviously we can write v k+1 = v k + a(t) dt We also can write y k+1 as an integral involving a(t), by following a few steps: y k+1 = y k + v(t) dt = y k + [t v(t)] +1 = y k + +1 v k+1 v k t a(t) dt after integration by parts = y k + h v k + +1 v k+1 +1 v k t a(t) dt by expanding the difference term = y k + h v k + +1 ( v k+1 v k ) = y k + h v k + +1 a(t) dt = y k + h v k + (+1 t) a(t) dt t a(t) dt by adding and subtracting h v k t a(t) dt after factoring t a(t) dt since v (t) = a(t) Suppose we choose τ [, +1 ]. Then, we can write expressions for a k and a k+1 using the Taylor series about τ: a k = a(τ) + a (τ)( τ) + O(h 2 ) a k+1 = a(τ) + a (τ)(+1 τ) + O(h 2 ) For any constant γ R, if we scale the first equation by 1 γ and the second by γ and sum the results, we find: a(τ) = (1 γ) a k + γ a k+1 + a (τ)((γ 1)( τ) γ(+1 τ)) + O(h 2 ) = (1 γ) a k + γ a k+1 + a (τ)(τ hγ ) + O(h 2 ) after substituting +1 = + h In the end, we wish to integrate a from to +1 to get the change in velocity. Thus, we compute: a(τ) dτ = (1 γ)h a k + γh a k+1 + = (1 γ)h a k + γh a k+1 + O(h 2 ), a (τ)(τ hγ ) dτ + O(h 3 ) where the second step holds because the integrand is O(h) and the interval of integration is of width h. In other words, we now know: v k+1 = v k + (1 γ)h a k + γh a k+1 + O(h 2 ) To make a similar approximation for y k+1, we can write (+1 t) a(t) dt = (+1 τ)((1 γ) a k + γ a k+1 + a (τ)(τ hγ )) dτ + O(h 3 ) = 1 2 (1 γ)h2 a k γh2 a k+1 + O(h 2 ) by a similar argument. 12

13 Thus, we can use our earlier relationship to show: y k+1 = y k + h v k + (+1 t) a(t) dt from before t ( k ) 1 = y k + h v k + 2 β h 2 a k + βh 2 a k+1 + O(h 2 ) Here, we use β instead of γ (and absorb a factor of two in the process) because the γ we choose for approximating y k+1 does not have to be the same as the one we choose for approximating v k+1. After all this integration, we have derived the class of Newmark schemes, with two input parameters γ and β, which has first-order accuracy by the proof above: ( ) 1 y k+1 = y k + h v k + 2 β h 2 a k + βh 2 a k+1 v k+1 = v k + (1 γ)h a k + γh a k+1 a k = F[, y k, v k ] Different choices of β and γ lead to different schemes. For instance, consider the following examples: β = γ = 0 gives the constant acceleration integrator: y k+1 = y k + h v k h2 a k v k+1 = v k + h a k This integrator is explicit and holds exactly when the acceleration is a constant function. β = 1 /2, γ = 1 gives the constant implicit acceleration integrator: y k+1 = y k + h v k h2 a k+1 v k+1 = v k + h a k+1 Here, the velocity is stepped implicitly using backward Euler, giving first-order accuracy. The update of y, however, can be written y k+1 = y k h( v k + v k+1 ), showing that it locally is updated by the midpoint rule; this is our first example of a scheme where the velocity and position updates have different orders of accuracy. Even so, it is possible to show that this technique, however, is globally first-order accurate in y. β = 1 /4, γ = 1 /2 gives the following second-order trapezoidal scheme after some algebra: x k+1 = x k h( v k + v k+1 ) v k+1 = v k h( a k + a k+1 ) 13

14 β = 0, γ = 1 /2 gives a second-order accurate central differencing scheme. In the canonical form, we have x k+1 = x k + h v k h2 a k v k+1 = v k h( a k + a k+1 ) The method earns its name because simplifying the equations above leads to the alternative form: v k+1 = y k+2 y k 2h a k+1 = y k+1 2 y k+1 + y k h 2 It is possible to show that Newmark s methods are unconditionally stable when 4β > 2γ > 1 and that second-order accuracy occurs exactly when γ = 1 /2 (CHECK) Staggered Grid A different way to achieve second-order accuracy in y is to use centered differences about the time + 1/2 + h /2: y k+1 = y k + h v k+ 1/2 Rather than use Taylor arguments to try to move v k+ 1/2, we can simply store velocities v at half points on the grid of time steps. Then, we can use a similar update to step forward the velocities: v k+ 3/2 = v k+ 1/2 + h a k+1. Notice that this update actually is second-order accurate for x as well, since if we substitute our expressions for v k+ 1/2 and v k+ 3/2 we can write: a k+1 = 1 h 2 ( y k+2 2 y k+1 + y k ) Finally, a simple approximation suffices for the acceleration term since it is a higher-order term: [ a k+1 = F +1, x k+1, 1 ] 2 ( v k+1/2 + v k+ 3/2) This expression can be substituted into the expression for v k+ 3/2. When F[ ] has no dependence on v, the method is fully explicit: y k+1 = y k + h v k+ 1/2 a k+1 = F[+1, y k+1 ] v k+ 3/2 = v k+ 1/2 + h a k+1 This is known as the leapfrog method of integration, thanks to the staggered grid of times and the fact that each midpoint is used to update the next velocity or position. 14

15 Otherwise, if the velocity update has dependence on v then the method becomes implicit. Often times, dependence on velocity is symmetric; for instance, wind resistance simply changes sign if you reverse the direction you are moving. This property can lead to symmetric matrices in the implicit step for updating velocities, making it possible to use conjugate gradients and related fast iterative methods to solve To Do Define stiff ODE Give table of time stepping methods for F[t; y] Use y notation more consistently 13.6 Problems TVD RK Multistep/multivalue methods a la Heath Verlet integration Symplectic integrators 15

### Ordinary 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

### Ordinary Differential Equations I

Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 32 Theme of Last Few

### Lecture V: The game-engine loop & Time Integration

Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change

### Numerical solution of ODEs

Numerical solution of ODEs Arne Morten Kvarving Department of Mathematical Sciences Norwegian University of Science and Technology November 5 2007 Problem and solution strategy We want to find an approximation

### Introduction to Initial Value Problems

Chapter 2 Introduction to Initial Value Problems The purpose of this chapter is to study the simplest numerical methods for approximating the solution to a first order initial value problem (IVP). Because

### Review for Exam 2 Ben Wang and Mark Styczynski

Review for Exam Ben Wang and Mark Styczynski This is a rough approximation of what we went over in the review session. This is actually more detailed in portions than what we went over. Also, please note

### Section 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

### Numerical Differential Equations: IVP

Chapter 11 Numerical Differential Equations: IVP **** 4/16/13 EC (Incomplete) 11.1 Initial Value Problem for Ordinary Differential Equations We consider the problem of numerically solving a differential

### ORDINARY DIFFERENTIAL EQUATIONS

ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 4884 NOVEMBER 9, 7 Summary This is an introduction to ordinary differential equations We

### PowerPoints 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

### Ordinary differential equation II

Ordinary Differential Equations ISC-5315 1 Ordinary differential equation II 1 Some Basic Methods 1.1 Backward Euler method (implicit method) The algorithm looks like this: y n = y n 1 + hf n (1) In contrast

### Review Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester

HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE

### Multistage Methods I: Runge-Kutta Methods

Multistage Methods I: Runge-Kutta Methods Varun Shankar January, 0 Introduction Previously, we saw that explicit multistep methods (AB methods) have shrinking stability regions as their orders are increased.

### CS205b/CME306. Lecture 4. x v. + t

CS05b/CME306 Lecture 4 Time Integration We now consider seeral popular approaches for integrating an ODE Forward Euler Forward Euler eolution takes on the form n+ = n + Because forward Euler is unstable

### 4 Stability analysis of finite-difference methods for ODEs

MATH 337, by T. Lakoba, University of Vermont 36 4 Stability analysis of finite-difference methods for ODEs 4.1 Consistency, stability, and convergence of a numerical method; Main Theorem In this Lecture

### CHAPTER 5: Linear Multistep Methods

CHAPTER 5: Linear Multistep Methods Multistep: use information from many steps Higher order possible with fewer function evaluations than with RK. Convenient error estimates. Changing stepsize or order

### DIFFERENTIAL EQUATIONS

DIFFERENTIAL EQUATIONS Basic Concepts Paul Dawkins Table of Contents Preface... Basic Concepts... 1 Introduction... 1 Definitions... Direction Fields... 8 Final Thoughts...19 007 Paul Dawkins i http://tutorial.math.lamar.edu/terms.aspx

### Introductory Numerical Analysis

Introductory Numerical Analysis Lecture Notes December 16, 017 Contents 1 Introduction to 1 11 Floating Point Numbers 1 1 Computational Errors 13 Algorithm 3 14 Calculus Review 3 Root Finding 5 1 Bisection

### Exact and Approximate Numbers:

Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.

### An Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations

An Overly Simplified and Brief Review of Differential Equation Solution Methods We will be dealing with initial or boundary value problems. A typical initial value problem has the form y y 0 y(0) 1 A typical

### Ph 22.1 Return of the ODEs: higher-order methods

Ph 22.1 Return of the ODEs: higher-order methods -v20130111- Introduction This week we are going to build on the experience that you gathered in the Ph20, and program more advanced (and accurate!) solvers

### Physically Based Modeling: Principles and Practice Differential Equation Basics

Physically Based Modeling: Principles and Practice Differential Equation Basics Andrew Witkin and David Baraff Robotics Institute Carnegie Mellon University Please note: This document is 1997 by Andrew

### Fourth 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)

### Chapter 9 Implicit Methods for Linear and Nonlinear Systems of ODEs

Chapter 9 Implicit Methods for Linear and Nonlinear Systems of ODEs In the previous chapter, we investigated stiffness in ODEs. Recall that an ODE is stiff if it exhibits behavior on widelyvarying timescales.

### Existence and uniqueness of solutions for nonlinear ODEs

Chapter 4 Existence and uniqueness of solutions for nonlinear ODEs In this chapter we consider the existence and uniqueness of solutions for the initial value problem for general nonlinear ODEs. Recall

### Jim Lambers MAT 772 Fall Semester Lecture 21 Notes

Jim Lambers MAT 772 Fall Semester 21-11 Lecture 21 Notes These notes correspond to Sections 12.6, 12.7 and 12.8 in the text. Multistep Methods All of the numerical methods that we have developed for solving

### DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science

DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI097 TU) Thursday January 014, 18:0-1:0

### The Initial Value Problem for Ordinary Differential Equations

Chapter 5 The Initial Value Problem for Ordinary Differential Equations In this chapter we begin a study of time-dependent differential equations, beginning with the initial value problem (IVP) for a time-dependent

### 2.29 Numerical Fluid Mechanics Fall 2011 Lecture 20

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 20 REVIEW Lecture 19: Finite Volume Methods Review: Basic elements of a FV scheme and steps to step-up a FV scheme One Dimensional examples d x j x j 1/2

### CHAPTER 1. First-Order Differential Equations and Their Applications. 1.1 Introduction to Ordinary Differential Equations

CHAPTER 1 First-Order Differential Equations and Their Applications 1.1 Introduction to Ordinary Differential Equations Differential equations are found in many areas of mathematics, science, and engineering.

### 1. Method 1: bisection. The bisection methods starts from two points a 0 and b 0 such that

Chapter 4 Nonlinear equations 4.1 Root finding Consider the problem of solving any nonlinear relation g(x) = h(x) in the real variable x. We rephrase this problem as one of finding the zero (root) of a

### Differential 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

### 1 Error Analysis for Solving IVP

cs412: introduction to numerical analysis 12/9/10 Lecture 25: Numerical Solution of Differential Equations Error Analysis Instructor: Professor Amos Ron Scribes: Yunpeng Li, Mark Cowlishaw, Nathanael Fillmore

### z x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.

Chapter 2 First order PDE 2.1 How and Why First order PDE appear? 2.1.1 Physical origins Conservation laws form one of the two fundamental parts of any mathematical model of Continuum Mechanics. These

### Chapter 3 Numerical Methods

Chapter 3 Numerical Methods Part 3 3.4 Differential Algebraic Systems 3.5 Integration of Differential Equations 1 Outline 3.4 Differential Algebraic Systems 3.4.1 Constrained Dynamics 3.4.2 First and Second

### Chapter 9 Implicit integration, incompressible flows

Chapter 9 Implicit integration, incompressible flows The methods we discussed so far work well for problems of hydrodynamics in which the flow speeds of interest are not orders of magnitude smaller than

### Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme

### Finite Difference and Finite Element Methods

Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element

### Quasi-Newton Methods

Newton s Method Pros and Cons Quasi-Newton Methods MA 348 Kurt Bryan Newton s method has some very nice properties: It s extremely fast, at least once it gets near the minimum, and with the simple modifications

### Ordinary Differential equations

Chapter 1 Ordinary Differential equations 1.1 Introduction to Ordinary Differential equation In mathematics, an ordinary differential equation (ODE) is an equation which contains functions of only one

### u n 2 4 u n 36 u n 1, n 1.

Exercise 1 Let (u n ) be the sequence defined by Set v n = u n 1 x+ u n and f (x) = 4 x. 1. Solve the equations f (x) = 1 and f (x) =. u 0 = 0, n Z +, u n+1 = u n + 4 u n.. Prove that if u n < 1, then

### 1 Ordinary Differential Equations

Ordinary Differential Equations.0 Mathematical Background.0. Smoothness Definition. A function f defined on [a, b] is continuous at ξ [a, b] if lim x ξ f(x) = f(ξ). Remark Note that this implies existence

### Ordinary Differential Equations: Initial Value problems (IVP)

Chapter Ordinary Differential Equations: Initial Value problems (IVP) Many engineering applications can be modeled as differential equations (DE) In this book, our emphasis is about how to use computer

### Euler s Method, Taylor Series Method, Runge Kutta Methods, Multi-Step Methods and Stability.

Euler s Method, Taylor Series Method, Runge Kutta Methods, Multi-Step Methods and Stability. REVIEW: We start with the differential equation dy(t) dt = f (t, y(t)) (1.1) y(0) = y 0 This equation can be

### Opleiding Wiskunde & Informatica

Opleiding Wiskunde & Informatica Numerics and continuation for Reaction-Diffusion equations Renzo Baasdam s1524054 Supervisors: Martina Chirilus-Bruckner (MI) & Michael Emmerich (LIACS) BACHELOR THESIS

### 8.1 Bifurcations of Equilibria

1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations

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

### Ordinary 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

### Lecture 4: Training a Classifier

Lecture 4: Training a Classifier Roger Grosse 1 Introduction Now that we ve defined what binary classification is, let s actually train a classifier. We ll approach this problem in much the same way as

### DIFFERENTIAL EQUATIONS COURSE NOTES, LECTURE 2: TYPES OF DIFFERENTIAL EQUATIONS, SOLVING SEPARABLE ODES.

DIFFERENTIAL EQUATIONS COURSE NOTES, LECTURE 2: TYPES OF DIFFERENTIAL EQUATIONS, SOLVING SEPARABLE ODES. ANDREW SALCH. PDEs and ODEs, order, and linearity. Differential equations come in so many different

### CSC321 Lecture 15: Exploding and Vanishing Gradients

CSC321 Lecture 15: Exploding and Vanishing Gradients Roger Grosse Roger Grosse CSC321 Lecture 15: Exploding and Vanishing Gradients 1 / 23 Overview Yesterday, we saw how to compute the gradient descent

### An Analysis of Five Numerical Methods for Approximating Certain Hypergeometric Functions in Domains within their Radii of Convergence

An Analysis of Five Numerical Methods for Approximating Certain Hypergeometric Functions in Domains within their Radii of Convergence John Pearson MSc Special Topic Abstract Numerical approximations of

### Optimization Tutorial 1. Basic Gradient Descent

E0 270 Machine Learning Jan 16, 2015 Optimization Tutorial 1 Basic Gradient Descent Lecture by Harikrishna Narasimhan Note: This tutorial shall assume background in elementary calculus and linear algebra.

### CPSC 540: Machine Learning

CPSC 540: Machine Learning Proximal-Gradient Mark Schmidt University of British Columbia Winter 2018 Admin Auditting/registration forms: Pick up after class today. Assignment 1: 2 late days to hand in

### MATHEMATICAL METHODS INTERPOLATION

MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU

### Solutions 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

### 12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis

ATH 337, by T. Lakoba, University of Vermont 113 12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis 12.1 Formulation of the IBVP and the minimax property of its

### Modeling and Solving Constraints. Erin Catto Blizzard Entertainment

Modeling and Solving Constraints Erin Catto Blizzard Entertainment Basic Idea Constraints are used to simulate joints, contact, and collision. We need to solve the constraints to stack boxes and to keep

### 7.2 Steepest Descent and Preconditioning

7.2 Steepest Descent and Preconditioning Descent methods are a broad class of iterative methods for finding solutions of the linear system Ax = b for symmetric positive definite matrix A R n n. Consider

### Runge-Kutta and Collocation Methods Florian Landis

Runge-Kutta and Collocation Methods Florian Landis Geometrical Numetric Integration p.1 Overview Define Runge-Kutta methods. Introduce collocation methods. Identify collocation methods as Runge-Kutta methods.

### Self-Concordant Barrier Functions for Convex Optimization

Appendix F Self-Concordant Barrier Functions for Convex Optimization F.1 Introduction In this Appendix we present a framework for developing polynomial-time algorithms for the solution of convex optimization

### Implicit Functions, Curves and Surfaces

Chapter 11 Implicit Functions, Curves and Surfaces 11.1 Implicit Function Theorem Motivation. In many problems, objects or quantities of interest can only be described indirectly or implicitly. It is then

### w T 1 w T 2. w T n 0 if i j 1 if i = j

Lyapunov Operator Let A F n n be given, and define a linear operator L A : C n n C n n as L A (X) := A X + XA Suppose A is diagonalizable (what follows can be generalized even if this is not possible -

### Numerical Algorithms for ODEs/DAEs (Transient Analysis)

Numerical Algorithms for ODEs/DAEs (Transient Analysis) Slide 1 Solving Differential Equation Systems d q ( x(t)) + f (x(t)) + b(t) = 0 dt DAEs: many types of solutions useful DC steady state: state no

### An Overview of Fluid Animation. Christopher Batty March 11, 2014

An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.

### Lie Groups for 2D and 3D Transformations

Lie Groups for 2D and 3D Transformations Ethan Eade Updated May 20, 2017 * 1 Introduction This document derives useful formulae for working with the Lie groups that represent transformations in 2D and

### Analog Signals and Systems and their properties

Analog Signals and Systems and their properties Main Course Objective: Recall course objectives Understand the fundamentals of systems/signals interaction (know how systems can transform or filter signals)

### Answer Key 1973 BC 1969 BC 24. A 14. A 24. C 25. A 26. C 27. C 28. D 29. C 30. D 31. C 13. C 12. D 12. E 3. A 32. B 27. E 34. C 14. D 25. B 26.

Answer Key 969 BC 97 BC. C. E. B. D 5. E 6. B 7. D 8. C 9. D. A. B. E. C. D 5. B 6. B 7. B 8. E 9. C. A. B. E. D. C 5. A 6. C 7. C 8. D 9. C. D. C. B. A. D 5. A 6. B 7. D 8. A 9. D. E. D. B. E. E 5. E.

### EAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science

EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Multidimensional Unconstrained Optimization Suppose we have a function f() of more than one

### 13 Numerical Solution of ODE s

13 NUMERICAL SOLUTION OF ODE S 28 13 Numerical Solution of ODE s In simulating dynamical systems, we frequently solve ordinary differential equations. These are of the form dx = f(t, x), dt where the function

### PARTIAL DIFFERENTIAL EQUATIONS

MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL

### Numerical Integration (Quadrature) Another application for our interpolation tools!

Numerical Integration (Quadrature) Another application for our interpolation tools! Integration: Area under a curve Curve = data or function Integrating data Finite number of data points spacing specified

### Strong Stability-Preserving (SSP) High-Order Time Discretization Methods

Strong Stability-Preserving (SSP) High-Order Time Discretization Methods Xinghui Zhong 12/09/ 2009 Outline 1 Introduction Why SSP methods Idea History/main reference 2 Explicit SSP Runge-Kutta Methods

### Computational project: Modelling a simple quadrupole mass spectrometer

Computational project: Modelling a simple quadrupole mass spectrometer Martin Duy Tat a, Anders Hagen Jarmund a a Norges Teknisk-Naturvitenskapelige Universitet, Trondheim, Norway. Abstract In this project

### 1 Lyapunov theory of stability

M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability

### ODE Runge-Kutta methods

ODE Runge-Kutta methods The theory (very short excerpts from lectures) First-order initial value problem We want to approximate the solution Y(x) of a system of first-order ordinary differential equations

### Lecture 6: Backpropagation

Lecture 6: Backpropagation Roger Grosse 1 Introduction So far, we ve seen how to train shallow models, where the predictions are computed as a linear function of the inputs. We ve also observed that deeper

### Solutions to Math 53 First Exam April 20, 2010

Solutions to Math 53 First Exam April 0, 00. (5 points) Match the direction fields below with their differential equations. Also indicate which two equations do not have matches. No justification is necessary.

### Robotics. Dynamics. Marc Toussaint U Stuttgart

Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, Euler-Lagrange equation, Newton-Euler recursion, general robot dynamics, joint space control, reference trajectory

### 5.2 - Euler s Method

5. - Euler s Method Consider solving the initial-value problem for ordinary differential equation: (*) y t f t, y, a t b, y a. Let y t be the unique solution of the initial-value problem. In the previous

### FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland 4 May 2012 Because the presentation of this material

### Nonlinear equations and optimization

Notes for 2017-03-29 Nonlinear equations and optimization For the next month or so, we will be discussing methods for solving nonlinear systems of equations and multivariate optimization problems. We will

### Infinite series, improper integrals, and Taylor series

Chapter 2 Infinite series, improper integrals, and Taylor series 2. Introduction to series In studying calculus, we have explored a variety of functions. Among the most basic are polynomials, i.e. functions

### Engineering Mathematics

Thoroughly Revised and Updated Engineering Mathematics For GATE 2017 and ESE 2017 Prelims Note: ESE Mains Electrical Engineering also covered Publications Publications MADE EASY Publications Corporate

### ACCELERATION OF RUNGE-KUTTA INTEGRATION SCHEMES

ACCELERATION OF RUNGE-KUTTA INTEGRATION SCHEMES PHAILAUNG PHOHOMSIRI AND FIRDAUS E. UDWADIA Received 24 November 2003 A simple accelerated third-order Runge-Kutta-type, fixed time step, integration scheme

### COSC 3361 Numerical Analysis I Ordinary Differential Equations (II) - Multistep methods

COSC 336 Numerical Analysis I Ordinary Differential Equations (II) - Multistep methods Fall 2005 Repetition from the last lecture (I) Initial value problems: dy = f ( t, y) dt y ( a) = y 0 a t b Goal:

### Graded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo

2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical

### THEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)

4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M

### Theory, Solution Techniques and Applications of Singular Boundary Value Problems

Theory, Solution Techniques and Applications of Singular Boundary Value Problems W. Auzinger O. Koch E. Weinmüller Vienna University of Technology, Austria Problem Class z (t) = M(t) z(t) + f(t, z(t)),

### Nonlinear Optimization

Nonlinear Optimization (Com S 477/577 Notes) Yan-Bin Jia Nov 7, 2017 1 Introduction Given a single function f that depends on one or more independent variable, we want to find the values of those variables

### Chapter 10. Initial value Ordinary Differential Equations

Chapter 10 Initial value Ordinary Differential Equations Consider the problem of finding a function y(t) that satisfies the following ordinary differential equation (ODE): dy dt = f(t, y), a t b. The function

### AMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.

Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error

### Numerical Methods in Physics and Astrophysics

Kostas Kokkotas 2 October 17, 2017 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation

### 7 Planar systems of linear ODE

7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution

### Chapter 3 Numerical Methods

Chapter 3 Numerical Methods Part 2 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization 1 Outline 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization Summary 2 Outline 3.2

### Damped harmonic motion

Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,

### Variable Step Runge-Kutta-Nyström Methods for the Numerical Solution of Reversible Systems

Variable Step Runge-Kutta-Nyström Methods for the Numerical Solution of Reversible Systems J. R. Cash and S. Girdlestone Department of Mathematics, Imperial College London South Kensington London SW7 2AZ,

### Math 128A Spring 2003 Week 12 Solutions

Math 128A Spring 2003 Week 12 Solutions Burden & Faires 5.9: 1b, 2b, 3, 5, 6, 7 Burden & Faires 5.10: 4, 5, 8 Burden & Faires 5.11: 1c, 2, 5, 6, 8 Burden & Faires 5.9. Higher-Order Equations and Systems