Simple shooting-projection method for numerical solution of two-point Boundary Value Problems
|
|
- Norman Chambers
- 6 years ago
- Views:
Transcription
1 Simple shooting-projection method for numerical solution of two-point Boundary Value Problems Stefan M. Filipo Ivan D. Gospodino a Department of Programming and Computer System Application, University of Chemical Technology and Metallurgy, Kl. Ohridski 8 Blvd., Sofia 756, Bulgaria Corresponding author idg@cornell.edu
2 Keywords Ordinary differential equations; two-point boundary value problem; second-derivative preservation; simple shooting-projection method Abstract This paper presents a novel shooting algorithm for solving two-point Boundary Value Problems (BVPs) for ordinary differential equations. The algorithm includes the following steps: First, a value for the initial condition at the first boundary is guessed and a forward numerical integration of the differential equation is performed so that an Initial Value Problem (IVP) solution, called a shooting trajectory, is obtained. The shooting trajectory starts from the first boundary constraint but typically does not end at the second boundary constraint. Next, the shooting trajectory is transformed into a projection trajectory that is an approximate BVP solution. The projection trajectory satisfies both boundary constraints and has the same second derivative as the shooting trajectory. Finally, the projection trajectory is used to correct the value of the initial condition and the procedure is repeated until convergence.. Introduction This paper studies numerical solutions of two-point BVPs of second-order ordinary differential equations. Probably the most intuitive numerical method for solving two-point BVPs is the simple (single) shooting method [], []. In the simple shooting method one makes a guess for the value of the initial condition at the first boundary and performs numerical integration of the differential equation to obtain a shooting trajectory. The shooting trajectory is an IVP solution but is not a BVP solution. The end of the shooting trajectory differs from the second boundary constraint by some distance, called error. Next, one corrects the initial condition at the first boundary in an effort to decrease the absolute value of the error in the next shooting, and repeats this procedure until an initial condition that produces a zero error is reached. Thus, any simple shooting method for solving two-point BVPs is, in fact, an iterative root finding method for finding the zero of the error as a function of the initial condition. One of the most prominent root finding methods is the Newton s method [3], [4], which at each iteration evaluates the derivative of the error with respect to the initial condition and uses this derivative to produce an estimate for the zero of the error. A variant of this method is the constant-slope Newton s method [5], which evaluates the derivative only at the first iteration and then uses this value of the derivative throughout. There are also methods that do not use derivatives. Instead, they utilize the value of the error at two or more initial conditions in order to estimate the zero of the error. Prominent examples of such methods are the bisection method [6], the secant method [6], and the Muller s method [7]. The proposed simple shooting-projection method for solving two-point BVPs differs from the above-mentioned methods in the way the zero error estimate is sought. At each iteration the method starts from some initial condition, finds the corresponding IVP solution (the shooting trajectory), transforms this IVP solution into an approximate BVP solution, called a projection trajectory, and uses this projection trajectory to obtain a zero error estimate used as an initial condition in the next iteration. The projection trajectory is an approximate BVP solution because by design it has the same second derivative as the shooting trajectory and meets both boundary constraints. As it will become evident in the analysis below, this method is akin to the constant-slope Newton s method but it is derived from different premises and has important differences. The paper is organized as follows: first, the transformation of the IVP solution into an approximate BVP solution is described for functions of one independent variable; then the simple shootingprojection algorithm is formulated; then, the convergence analysis of the algorithm is performed; and finally the algorithm is tested by solving several types of two-point BVPs.
3 . Finding an approximate BVP solution using an IVP solution Let u be a real valued function of a real independent variable t [a, b]. If u satisfies the second-order ordinary differential equation u (t) = f (t, u(t), u (t)), t (a, b), () together with two initial conditions u(a) = u a, u (a) =, () then u is an IVP solution. If u satisfies () together with two boundary constraints u(a) = u a, u(b) = u b, (3) then u is a two-point BVP solution. Since a BVP solution is much harder to find than an IVP solution, we propose a simple way of finding an approximate BVP solution u using a given IVP solution u. This u satisfies the differential equation () approximately and satisfies the boundary constraints (3) exactly. As it will be shown below, this approximate BVP solution could be easily incorporated into an iterative shooting procedure for finding the exact BVP solution. A natural choice of u, which leads to simple mathematical derivations, is one that satisfies the following differential equation: u (t) = u (t), t (a, b), (4) together with the two boundary constraints u(a) = u a, u(b) = u b. (5) In (4) u is the sought unknown function, while u is the given IVP solution. Since u satisfies () we can write: u = f (t, u, u ) = f (t, u, u ) + R, t (a, b), (6) where f (t, u, u ) is expanded around u, and R stands for higher order terms. quation (6) tells us that if R is small, then u is an approximate solution to the differential eqn. (). Since, on the other hand, u satisfies the two boundary constraints (5), it follows that u is an approximate BVP solution. If f is constant or depends only on t, then u is an exact BVP solution. Since the approximate BVP solution (4), (5) is used in an iterative procedure it does not necessarily have to be a close approximation. The iterative procedure will, in general, work even if R in eqn. (6) is not small. quation (4), complemented by eqn. (5), defines a second-derivative preserving transformation of the IVP solution u into an approximate BVP solution u. It is easy to see that this transformation is a projection of u into u. Indeed, if we denote this transformation by P, then u=pu. On the other hand, since u satisfies the boundary conditions, u=pu. Thus, Pu= u=pu=p u, or, since this is true for any u, P=P, which is a definition of a projection operator. We call the function u (or the graph {(t, u), t [a, b]}) a shooting trajectory and the function u (or the graph {(t, u), t [a, b]}) a projection trajectory. An interesting property of the second-derivative preserving transformation is that u minimizes the square of the H seminorm of the difference u-u: 3
4 b S = ( u u ) dt. (7) a In (7), u is not a variable but the given IVP solution. If the function u, with two fixed boundary values (5), is to minimize (7) then the uler-lagrange equation must hold: d dt L L =, u u (8) where L=(u -u ). From where: d ( u u ) =, dt (9) which leads to eqn. (4). quation (7) tells us that the first derivative of u is minimally deviated from the first derivative of u, i.e. the two functions have similar shapes in this sense. 3. Algorithm for solving two-point BVPs using second-derivative preserving transformation of trajectories The following simple shooting-projection algorithm for solving two-point BVPs, summarized graphically in fig. below, is proposed: Step : Make a guess for the initial condition Step : Obtain a shooting trajectory We introduce v = u and numerically integrate (), using the uler s method, as a system of two first-order differential equations. We start from the first boundary constraint u a and the initial condition to obtain the following shooting trajectory (see fig. ): u i = u i- + h v i-, v i = v i- + h f(t i, u i, v i- ), i =,, N, () with initial values u = u a, and v =, where u i and v i are the values of u and t t i =a+h(i-); h =(b-a)/(n-) is the size of the integration step; and (N-) is the number of integration points. The shooting trajectory u is a numerical IVP solution that starts from the first boundary constraint (i.e. u =u a ) but does not meet the second boundary constrain (i.e. u N u b ). The distance u N -u b is denoted by and is called error. Step 3: Obtain a projection trajectory In order to obtain the projection trajectory we discretize eqn. (4) using the central difference scheme u i+ - u i + u i- = u i+ - u i + u i-, i =,, N-, () and impose the two boundary constraints (5): u = u a, u N = u b. () The solution of () together with () yields the projection trajectory (see fig. ): u i = u i - (u N - u b )(i-)/(n-), i =,, N. (3) 4
5 The projection trajectory u is an approximate BVP solution that is obtained from the shooting trajectory u, provided by Step. quation (3) is the sought second-derivative preserving projection transformation u=pu. Note, that the derivations in Step 3 are independent of the integration scheme used to obtain the IVP solution u in Step. Step 4: Obtain a new, corrected initial condition v a According to eqn. (3), the first two points of the projection trajectory are: u = u, u = u - (u N - u b )/(N-). (4) Subtracting the first equation in (4) from the second, dividing by h, and taking into account that h(n-)=b-a, a simple final expression for the new initial condition v a is obtained: v a = - ( )/k, (5) where v a =(u - u )/h, =(u - u )/h, ( )=u N ( ) - u b, and k=b-a (see fig.). Step 5: Replace in Step with v a The proposed algorithm can simply be described as: guess an initial condition; shoot to obtain an IVP solution; obtain an approximate BVP solution using this IVP solution; shoot again in the direction suggested by the approximate BVP solution, etc. Note, that eqn. (3) is used only to derive eqn. (5), which, in fact, does not contain any projection trajectory points. It contains only the endpoint u N of the shooting trajectory. Therefore, while carrying out the numerical procedure the actual construction of the projection trajectory is not necessary, i.e. Step 3 is virtual. This step was presented only to derive the mathematical equations and for conceptual clarity. Thus, the presented simple shooting-projection algorithm repeats steps, 4, and 5 until the shooting trajectory satisfies the second boundary constraint within a prescribed tolerance, i.e. u N u b. 3 5 shooting trajectory projection trajectory u N u 5 u b 5 slope u u 3 u N u a -5 u =u u u a slope t t 3 t Fig.. Second-derivative preserving transformation of the shooting trajectory u into a projection trajectory u. In the iterative algorithm is replaced by v a in the next iteration. Most shooting methods use at least two shooting trajectories per iteration in order to obtain an estimate for the zero of the error. quation (5) shows that the proposed simple shooting-projection method uses only one shooting trajectory to achieve the same. Instead of a second shooting trajectory the method makes use of a (virtual) projection trajectory that is an approximate BVP solution and therefore it is suitable for obtaining an estimate for the root of ( )=. b 5
6 Step 3 of the algorithm is equivalent to a Picard s iteration step [8], but only if u is an IVP solution. The Picard s iteration method is a relaxation method [6] that, if started from the IVP solution u, would produce the same u in the first iteration as the simple shooting-projection method, but will proceed to relax u without resorting to shooting. Thus, the Picard s method will not preserve the second derivative in the following iterations. Like other relaxation methods, and unlike the shooting methods, the Picard s method has convergence issues when the two-point BVP solution is oscillatory. In the Results section an example will be shown where the Picard s and other relaxation methods fail to converge whereas the proposed simple shooting-projection method converges. The form of eqn. (5) is similar in structure to the form of the corresponding equations for the secant, the Newton s shooting, and especially so to the constant-slope Newton s shooting methods. The properties of eqn. (5) are discussed in the next section. 4. Convergence analysis of the simple shooting-projection algorithm A plot of the error as a function of the initial condition is shown on fig. a below with a thick line. The point A in the figure corresponds to the starting initial condition. As can be seen from eqn. (5), the new initial condition v a is obtained by drawing a line with a slope k=b-a through point A. The next initial condition is obtained by drawing a line with the same slope k through point B. ach next initial condition is obtained analogously. Thus, the root of ( )= is reached after successive construction of parallel lines of slope k, shown in fig. a with thin lines. A main feature of the simple shooting-projection method is that the slope k is constant. In this sense the method is akin to the constant-slope Newton s method, where the slope is also constant, but has a different value the value of d/d evaluated at the starting initial condition. Since the starting is arbitrary, this slope will also be arbitrary and may not lead to convergence unless the starting happens to be close to the root. In contrast, the slope k does not depend on the starting. Its value b-a comes from the boundaries of the BVP. It is the most meaningful value to use under the assumptions, used to derive eqn. (5). Other shooting methods, such as the secant and the Newton s shooting methods, differ from the proposed simple shooting-projection method in that they readjust k at each iteration. Readjusting k could be an advantage, but sometimes it could be a drawback, as will be demonstrated in the Results section. Let m be the slope of ( ) at the root. Analogously to the constant-slope Newton s method, if close to the root the simple shooting-projection method converges [7] when - m/k <. (6) Otherwise the method diverges. The method converges linearly except for the trivial case m=k where the convergence is quadratic. Since k is always positive the convergence criterion above leads to the four cases, graphically shown in fig.. A A B k B k m m (a) (b) 6
7 k A m m k B A B (c) (d) Fig.. Convergence/divergence sequences of the simple shooting-projection method for the following four possible cases: (a) convergence when k > m >, (b) convergence when m/ < k < m, (c) divergence when k < m/, and (d) divergence when m <. 5. Results The simple shooting-projection algorithm for solving two-point BVPs was tested on several types of second-order ordinary differential equations. xample. Consider the following two-point BVP u = exp(u), t (, ), (7) u() =ln[(tan +)/], u() =ln[(tan (5/)+)/], with an analytical solution: u(t) =ln[(tan (t / +3/)+)/]. (8) Figure 3 traces the evolution of the shooting and the corresponding projection trajectories obtained after iterations,, 3, and 4 using the simple shooting-projection method. As discussed previously, the projection trajectories are not needed for the iteration process but they are shown in the figure to illustrate the algorithm and to demonstrate the goodness of the approximation for the BVP solution. The first shooting trajectory, shown in fig. 3a, corresponds to a starting initial condition =. The first projection trajectory predicts a negative new initial condition v a. The shooting trajectory in fig. b is obtained using this initial condition. The second projection trajectory is already very close to the analytical solution, shown by the solid curve. Convergence within a prescribed tolerance <. is reached at iteration 6. 7
8 3.5.5 shooting trajectory projection trajectory analytical solution u.5 u t t (a) (b) u.5 u t t (c) (d) Fig. 3. volution of the shooting and projection trajectories when solving eqn. (7) using the simple shooting-projection method for N=3 and starting = at: (a) iteration ; (b) iteration ; (c) iteration 3; (d) iteration 4. Figure 4 plots the graph of ( ) (which is found numerically) and the sequence of parallel lines with slope k=b-a that lead to the root of ( )=. The figure shows that the value of k is close to the value of the slope of at the root Fig. 4. Convergence sequence for the iterative solution of eqn. (7) using the simple shootingprojection method for a starting initial condition = and N=3. The root of ( )= is =
9 xample. Consider the following two-point BVP u = -3u u /t, t (, ), (9) u() = /, u() = / 5, with an analytical solution: u ( t) = t / + t. () The shape of ( ) for eqn. (9) is shown in fig. 5. This is a typical case for which the Newton s shooting method and the constant-slope Newton s shooting method will converge only if the starting initial condition is close to the root. For starting =5 the simple shooting-projection method converges within a prescribed tolerance <. in 7 iterations in the way shown in fig. 5a. For the same starting =5 the Newton s shooting method diverges in the way shown in fig. 5b. Starting from the same =5 the constant-slope Newton s shooting method oscillates around the root perpetually without being able to converge. The first few iterations are shown in fig. 5c. simple shooting-projection method Newton shooting method modified Newton shooting method (a) (b) (c) Fig. 5. Convergence/divergence sequences for the iterative solution of eqn. (9) for a starting initial condition = and N= using: (a) the simple shooting-projection method; (b) the Newton s shooting method; and (c) the constant-slope Newton s shooting method. The root of ( )= is = xample 3. Consider the following two-point BVP u = -u cosh(tu)/, t (, ), () u() =, u() =. This equation has no explicit analytical solution. The function ( ) for eqn. () has multiple roots, four of which are shown in fig. 7 below. We set out to find the solution that corresponds to root 3. Starting from = the simple shooting-projection method finds this solution in iterations within a prescribed tolerance <.. The evolution of the shooting trajectory at each iteration (numbered through ) is shown in fig. 6. 9
10 5 4 3 u t Fig. 6. volution of the shooting trajectory when solving eqn. () using the simple shootingprojection method for N= and starting = (corresponding to shooting trajectory ). Convergence is reached after iterations (corresponding to the shooting trajectory ). Figure 7 shows how root 3 of ( )= is reached. The consecutive values of at each iteration are indicated with,, 3, in the figure. They correspond to the shooting trajectories, shown in fig. 6. The function ( ) has three local extrema between root and root 3. Unlike other methods, the simple shooting-projection method does not have convergence issues around such extrema. It simply marches through them and finds root 3. In fact, the convergence of the method is not influenced by the shape of around the root. It is easy to see that the simple shooting-projection method will find root 3 for any starting that lies between root and root 4. simple shooting-projection method 3 9 root root root 3 root Fig. 7. Convergence sequence for the iterative solution of eqn. () using the simple shootingprojection method for a starting initial condition = (point ) and N=. The root 3 of ( )= is = Figure 8 demonstrates the divergence of the secant shooting method when solving eqn. () starting with the two initial conditions = (point ) and =- (point ). The figure shows that at iteration 6 the secant line becomes nearly horizontal and therefore the secant method practically fails. Although some other pairs of initial conditions may happen to find root 3 this would require a considerable trial-anderror effort.
11 secant shooting method root root root 3 root Fig. 8. Divergence sequence for the iterative solution of eqn. () using the secant shooting method for a starting initial conditions = (point ) and =- (point ), and N=. When Newton s shooting method is used to solve eqn. () for starting = it fails to find root 3. Instead, it finds root 4. The constant-slope Newton s method, started from =, finds root 3 but, as can be seen in fig. 8, if the method is started from any between root and root 3, corresponding to a negative slope of or to a positive slope of that is close (or equal) to zero, the method will fail to find root 3. Figure 6 shows that the solution of the two-point BVP () is oscillatory. As mentioned previously, the relaxation methods, like the Jacobi [9], the Picard s, or the Newton s relaxation [6] methods, have convergence issues with oscillatory solutions, while the shooting methods, in general, do not. All three relaxation methods were tried for solving eqn. (). The relaxation methods require a starting trajectory guess that does not necessarily have to be an IVP solution. In order to meaningfully compare these methods to the simple shooting-projection method, as a starting trajectory guess we used the shooting trajectory, corresponding to = and N= (corresponding to trajectory in fig. 6 and to point in fig. 7). Then all three relaxation methods failed to converge. Figures 4, 5a, and 7 show that the value of k is close to the value of the slope of at the root for all examples, considered in this work. From the presented results one can see that the majority of the well-known existing shooting and relaxation methods for solving two-point BVPs fail to find the solution of eqn. () corresponding to root 3. As discussed, the proposed simple shooting-projection method will find this solution, regardless of the starting, as long as it lies between root and root Conclusion This paper described a novel simple shooting method for numerical solution of two-point BVPs. The method starts with a guess for the initial condition, obtains a shooting trajectory, and then transforms this shooting trajectory into a projection trajectory that is an approximate two-point BVP solution. From the projection trajectory a new initial condition is obtained, which is used for the next iteration. This idea leads to a simple, easy to analyze iteration formula for correcting the initial condition. It was demonstrated that at certain cases the method outperforms prominent shooting and relaxation methods. The proposed simple shooting-projection method has three important features: (i) it uses only one shooting trajectory per iteration in order to correct the initial condition; (ii) the slope k in the iterative correction formula (5) is a constant with a value that comes from the boundaries and thus it does not depend on the value of the initial condition ; and (iii) the projection step is a relaxation step.
12 A recommendation for future work is to perform an analysis that shows how the value of the slope k=b-a relates to the value of the slope of at the root of ( )=. This would shed additional light on the convergence properties of the method. Acknowledgements We would like to thank Dr. Stephen Vavasis, Dr. Grigor Georgiev, Dr. Stefka Dimova, Dr. Geno Nikolov, Dr. Jordanka Angelova, and Ana Avdzieva for the valuable discussions. References [] Herbert B. Keller. (968) Numerical Methods For Two-Point Boundary-Value Problems. Blaisdell. [] Herbert B. Keller. (976) Numerical Solution of Two Point Boundary Value Problems. SIAM. [3] Newton, I. (664-67) Methodus fluxionum et serierum infinitarum. [4] Raphson, J. (69) Analysis aequationum universalis. London. [5] Steward. G. W. (996) Afternotes on Numerical Anlysis. SIAM. [6] Press, W. H., Flannery, B. P., Teukolsky, S. A., Vetterling, W. T. (99) Numerical Recipes in FORTRAN: The Art of Scientific Computing, nd ed. Cambridge, ngland: Cambridge University Press. [7] Muller, David. (956) A Method for Solving Algebraic quations Using an Automatic Computer, Mathematical Tables and Other Aids to Computation,, 8-5. [8] Picard. (893) Sur lapplication des metodes dapproximations succesives a letude de certains equations differentielles ordinaires, Journal de Mathmatiques Pures et Appliques, 9, 7 7. [9] Bronshtein, I. N. and Semendyayev, K. A. (997) Handbook of Mathematics, 3rd ed. New York: Springer-Verlag.
Parameter Estimation in Differential Equations: A Numerical Study of Shooting Methods
Parameter Estimation in Differential Equations: A Numerical Study of Shooting Methods Franz Hamilton Faculty Advisor: Dr Timothy Sauer January 5, 2011 Abstract Differential equation modeling is central
More informationToday s class. Roots of equation Finish up incremental search Open methods. Numerical Methods, Fall 2011 Lecture 5. Prof. Jinbo Bi CSE, UConn
Today s class Roots of equation Finish up incremental search Open methods 1 False Position Method Although the interval [a,b] where the root becomes iteratively closer with the false position method, unlike
More informationCOURSE: NUMERICAL ANALYSIS. LESSON: Methods for Solving Non-Linear Equations
COURSE: NUMERICAL ANALYSIS LESSON: Methods for Solving Non-Linear Equations Lesson Developer: RAJNI ARORA COLLEGE/DEPARTMENT: Department of Mathematics, University of Delhi Page No. 1 Contents 1. LEARNING
More informationModule 1 Lecture Notes 2. Optimization Problem and Model Formulation
Optimization Methods: Introduction and Basic concepts 1 Module 1 Lecture Notes 2 Optimization Problem and Model Formulation Introduction In the previous lecture we studied the evolution of optimization
More information5.5 Newton s Approximation Method
498CHAPTER 5. USING DERIVATIVES TO ANALYZE FUNCTIONS; FURTHER APPLICATIONS 4 3 y = x 4 3 f(x) = x cosx y = cosx 3 3 x = cosx x cosx = 0 Figure 5.: Figure showing the existence of a solution of x = cos
More information99 International Journal of Engineering, Science and Mathematics
Journal Homepage: Applications of cubic splines in the numerical solution of polynomials Najmuddin Ahmad 1 and Khan Farah Deeba 2 Department of Mathematics Integral University Lucknow Abstract: In this
More informationInternational Association of Scientific Innovation and Research (IASIR) (An Association Unifying the Sciences, Engineering, and Applied Research)
International Association of Scientific Innovation and Research (IASIR) (An Association Unifying the Sciences, Engineering, and Applied Research) International Journal of Emerging Technologies in Computational
More informationReals 1. Floating-point numbers and their properties. Pitfalls of numeric computation. Horner's method. Bisection. Newton's method.
Reals 1 13 Reals Floating-point numbers and their properties. Pitfalls of numeric computation. Horner's method. Bisection. Newton's method. 13.1 Floating-point numbers Real numbers, those declared to be
More information(Part - 1) P. Sam Johnson. April 14, Numerical Solution of. Ordinary Differential Equations. (Part - 1) Sam Johnson. NIT Karnataka.
P. April 14, 2015 1/51 Overview We discuss the following important methods of solving ordinary differential equations of first / second order. Picard s method of successive approximations (Method of successive
More informationThe Bisection Method versus Newton s Method in Maple (Classic Version for Windows)
The Bisection Method versus (Classic Version for Windows) Author: Barbara Forrest Contact: baforres@uwaterloo.ca Copyrighted/NOT FOR RESALE version 1.1 Contents 1 Objectives for this Lab i 2 Approximate
More informationCS 450 Numerical Analysis. Chapter 7: Interpolation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationVARIANCE REDUCTION TECHNIQUES IN MONTE CARLO SIMULATIONS K. Ming Leung
POLYTECHNIC UNIVERSITY Department of Computer and Information Science VARIANCE REDUCTION TECHNIQUES IN MONTE CARLO SIMULATIONS K. Ming Leung Abstract: Techniques for reducing the variance in Monte Carlo
More informationResearch Report AI A Numerical Equation Solver in Prolog Michael A. Covington Artificial Intelligence Programs The University of Georgia
Research Report AI 1989 02 A Numerical Equation Solver in Prolog Michael A. Covington Artificial Intelligence Programs The University of Georgia Athens, Georgia 30602 U.S.A. A Numerical Equation Solver
More informationP063 ADAPTIVE TRAVEL-TIME AND RAY- PARAMETER INVERSION OF DENSELY SAMPLED 2-D SEISMIC DATA
1 P063 ADAPTIVE TRAVEL-TIME AND RAY- PARAMETER INVERSION OF DENSELY SAMPLED 2-D SEISMIC DATA I. TRINKS 1, S. SINGH 2, C. CHAPMAN 3, P. BARTON 1, M. BOSCH 4, A. CHERRETT 5 Addresses 1 Bullard Laboratories,
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 36
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 36 In last class, we have derived element equations for two d elasticity problems
More informationIn this chapter, we will discuss several general algorithms that find a solution of the following canonical equation. f(x) = 0 (1.
Chapter 1 Root finding algorithms 1.1 Root finding problem In this chapter, we will discuss several general algorithms that find a solution of the following canonical equation f(x) = 0 (1.1) Here f(x)
More informationModeling with Uncertainty Interval Computations Using Fuzzy Sets
Modeling with Uncertainty Interval Computations Using Fuzzy Sets J. Honda, R. Tankelevich Department of Mathematical and Computer Sciences, Colorado School of Mines, Golden, CO, U.S.A. Abstract A new method
More informationHandout 2 - Root Finding using MATLAB
Handout 2 - Root Finding using MATLAB Middle East Technical University MATLAB has couple of built-in root finding functions. In this handout we ll have a look at fzero, roots and solve functions. It is
More informationf( x ), or a solution to the equation f( x) 0. You are already familiar with ways of solving
The Bisection Method and Newton s Method. If f( x ) a function, then a number r for which f( r) 0 is called a zero or a root of the function f( x ), or a solution to the equation f( x) 0. You are already
More informationMATHEMATICAL ANALYSIS, MODELING AND OPTIMIZATION OF COMPLEX HEAT TRANSFER PROCESSES
MATHEMATICAL ANALYSIS, MODELING AND OPTIMIZATION OF COMPLEX HEAT TRANSFER PROCESSES Goals of research Dr. Uldis Raitums, Dr. Kārlis Birģelis To develop and investigate mathematical properties of algorithms
More informationIntroduction to Exploratory Data Analysis
Introduction to Exploratory Data Analysis Ref: NIST/SEMATECH e-handbook of Statistical Methods http://www.itl.nist.gov/div898/handbook/index.htm The original work in Exploratory Data Analysis (EDA) was
More informationAn Efficient Method for Solving the Direct Kinematics of Parallel Manipulators Following a Trajectory
An Efficient Method for Solving the Direct Kinematics of Parallel Manipulators Following a Trajectory Roshdy Foaad Abo-Shanab Kafr Elsheikh University/Department of Mechanical Engineering, Kafr Elsheikh,
More informationToday. Golden section, discussion of error Newton s method. Newton s method, steepest descent, conjugate gradient
Optimization Last time Root finding: definition, motivation Algorithms: Bisection, false position, secant, Newton-Raphson Convergence & tradeoffs Example applications of Newton s method Root finding in
More informationNUMERICAL METHODS, NM (4776) AS
NUMERICAL METHODS, NM (4776) AS Objectives To provide students with an understanding that many mathematical problems cannot be solved analytically but require numerical methods. To develop a repertoire
More informationAn explicit feature control approach in structural topology optimization
th World Congress on Structural and Multidisciplinary Optimisation 07 th -2 th, June 205, Sydney Australia An explicit feature control approach in structural topology optimization Weisheng Zhang, Xu Guo
More informationREVIEW, pages
REVIEW, pages 69 697 8.. Sketch a graph of each absolute function. Identif the intercepts, domain, and range. a) = ƒ - + ƒ b) = ƒ ( + )( - ) ƒ 8 ( )( ) Draw the graph of. It has -intercept.. Reflect, in
More informationAn interesting related problem is Buffon s Needle which was first proposed in the mid-1700 s.
Using Monte Carlo to Estimate π using Buffon s Needle Problem An interesting related problem is Buffon s Needle which was first proposed in the mid-1700 s. Here s the problem (in a simplified form). Suppose
More information1 2 (3 + x 3) x 2 = 1 3 (3 + x 1 2x 3 ) 1. 3 ( 1 x 2) (3 + x(0) 3 ) = 1 2 (3 + 0) = 3. 2 (3 + x(0) 1 2x (0) ( ) = 1 ( 1 x(0) 2 ) = 1 3 ) = 1 3
6 Iterative Solvers Lab Objective: Many real-world problems of the form Ax = b have tens of thousands of parameters Solving such systems with Gaussian elimination or matrix factorizations could require
More informationVisualization and Analysis of Inverse Kinematics Algorithms Using Performance Metric Maps
Visualization and Analysis of Inverse Kinematics Algorithms Using Performance Metric Maps Oliver Cardwell, Ramakrishnan Mukundan Department of Computer Science and Software Engineering University of Canterbury
More informationThe Level Set Method applied to Structural Topology Optimization
The Level Set Method applied to Structural Topology Optimization Dr Peter Dunning 22-Jan-2013 Structural Optimization Sizing Optimization Shape Optimization Increasing: No. design variables Opportunity
More informationBeams. Lesson Objectives:
Beams Lesson Objectives: 1) Derive the member local stiffness values for two-dimensional beam members. 2) Assemble the local stiffness matrix into global coordinates. 3) Assemble the structural stiffness
More informationThe New Bisection Plus Algorithm by Namir Shammas
The New Bisection Plus Algorithm 1 The New Bisection Plus Algorithm by Namir Shammas Introduction This article presents a new variant for the root-bracketing Bisection algorithm. This new version injects
More informationError in Numerical Methods
Error in Numerical Methods By Brian D. Storey This section will describe two types of error that are common in numerical calculations: roundoff and truncation error. Roundoff error is due to the fact that
More informationLacunary Interpolation Using Quartic B-Spline
General Letters in Mathematic, Vol. 2, No. 3, June 2017, pp. 129-137 e-issn 2519-9277, p-issn 2519-9269 Available online at http:\\ www.refaad.com Lacunary Interpolation Using Quartic B-Spline 1 Karwan
More information1
Zeros&asymptotes Example 1 In an early version of this activity I began with a sequence of simple examples (parabolas and cubics) working gradually up to the main idea. But now I think the best strategy
More informationDescription Syntax Remarks and examples Conformability Diagnostics References Also see
Title stata.com solvenl( ) Solve systems of nonlinear equations Description Syntax Remarks and examples Conformability Diagnostics References Also see Description The solvenl() suite of functions finds
More information4 Visualization and. Approximation
4 Visualization and Approximation b A slope field for the differential equation y tan(x + y) tan(x) tan(y). It is not always possible to write down an explicit formula for the solution to a differential
More informationMulti Layer Perceptron trained by Quasi Newton learning rule
Multi Layer Perceptron trained by Quasi Newton learning rule Feed-forward neural networks provide a general framework for representing nonlinear functional mappings between a set of input variables and
More informationHandout 4 - Interpolation Examples
Handout 4 - Interpolation Examples Middle East Technical University Example 1: Obtaining the n th Degree Newton s Interpolating Polynomial Passing through (n+1) Data Points Obtain the 4 th degree Newton
More informationFour equations are necessary to evaluate these coefficients. Eqn
1.2 Splines 11 A spline function is a piecewise defined function with certain smoothness conditions [Cheney]. A wide variety of functions is potentially possible; polynomial functions are almost exclusively
More informationApproximating Square Roots
Math 560 Fall 04 Approximating Square Roots Dallas Foster University of Utah u0809485 December 04 The history of approximating square roots is one of the oldest examples of numerical approximations found.
More informationPRIMAL-DUAL INTERIOR POINT METHOD FOR LINEAR PROGRAMMING. 1. Introduction
PRIMAL-DUAL INTERIOR POINT METHOD FOR LINEAR PROGRAMMING KELLER VANDEBOGERT AND CHARLES LANNING 1. Introduction Interior point methods are, put simply, a technique of optimization where, given a problem
More informationMETHOD OF LEAST SQUARES FOR COMPUTING NORMALS AND CURVATURES FROM 2D VOLUME FRACTIONS
XVIII International Conference on Water Resources CMWR 2010 J. Carrera (Ed) CIMNE, Barcelona 2010 METHOD OF LEAST SQUARES FOR COMPUTING NORMALS AND CURVATURES FROM 2D VOLUME FRACTIONS Poorya A. Ferdowsi
More informationCourse of study- Algebra Introduction: Algebra 1-2 is a course offered in the Mathematics Department. The course will be primarily taken by
Course of study- Algebra 1-2 1. Introduction: Algebra 1-2 is a course offered in the Mathematics Department. The course will be primarily taken by students in Grades 9 and 10, but since all students must
More informationFunctions 3.6. Fall Math (Math 1010) M / 13
Functions 3.6 Fall 2013 - Math 1010 (Math 1010) M 1010 3.6 1 / 13 Roadmap 3.6 - Functions: Relations, Functions 3.6 - Evaluating Functions, Finding Domains and Ranges (Math 1010) M 1010 3.6 2 / 13 3.6
More informationUsing Modified Euler Method (MEM) for the Solution of some First Order Differential Equations with Initial Value Problems (IVPs).
Using Modified Euler Method (MEM) for the Solution of some First Order Differential Equations with Initial Value Problems (IVPs). D.I. Lanlege, Ph.D. * ; U.M. Garba, B.Sc.; and A. Aluebho, B.Sc. Department
More informationA Modified Spline Interpolation Method for Function Reconstruction from Its Zero-Crossings
Scientific Papers, University of Latvia, 2010. Vol. 756 Computer Science and Information Technologies 207 220 P. A Modified Spline Interpolation Method for Function Reconstruction from Its Zero-Crossings
More informationMath 225 Scientific Computing II Outline of Lectures
Math 225 Scientific Computing II Outline of Lectures Spring Semester 2003 I. Interpolating polynomials Lagrange formulation of interpolating polynomial Uniqueness of interpolating polynomial of degree
More informationNumerical methods in physics
Numerical methods in physics Marc Wagner Goethe-Universität Frankfurt am Main summer semester 2018 Version: April 9, 2018 1 Contents 1 Introduction 3 2 Representation of numbers in computers, roundoff
More informationSYNTHESIS OF PLANAR MECHANISMS FOR PICK AND PLACE TASKS WITH GUIDING LOCATIONS
Proceedings of the ASME 2013 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference IDETC/CIE 2013 August 4-7, 2013, Portland, Oregon, USA DETC2013-12021
More informationThemes in the Texas CCRS - Mathematics
1. Compare real numbers. a. Classify numbers as natural, whole, integers, rational, irrational, real, imaginary, &/or complex. b. Use and apply the relative magnitude of real numbers by using inequality
More informationJournal of Engineering Research and Studies E-ISSN
Journal of Engineering Research and Studies E-ISS 0976-79 Research Article SPECTRAL SOLUTIO OF STEADY STATE CODUCTIO I ARBITRARY QUADRILATERAL DOMAIS Alavani Chitra R 1*, Joshi Pallavi A 1, S Pavitran
More informationSolve Non-Linear Parabolic Partial Differential Equation by Spline Collocation Method
Solve Non-Linear Parabolic Partial Differential Equation by Spline Collocation Method P.B. Choksi 1 and A.K. Pathak 2 1 Research Scholar, Rai University,Ahemdabad. Email:pinalchoksey@gmail.com 2 H.O.D.
More informationNonparametric Regression
Nonparametric Regression John Fox Department of Sociology McMaster University 1280 Main Street West Hamilton, Ontario Canada L8S 4M4 jfox@mcmaster.ca February 2004 Abstract Nonparametric regression analysis
More information2.4. A LIBRARY OF PARENT FUNCTIONS
2.4. A LIBRARY OF PARENT FUNCTIONS 1 What You Should Learn Identify and graph linear and squaring functions. Identify and graph cubic, square root, and reciprocal function. Identify and graph step and
More informationComputers in Engineering Root Finding Michael A. Hawker
Computers in Engineering COMP 208 Root Finding Michael A. Hawker Root Finding Many applications involve finding the roots of a function f(x). That is, we want to find a value or values for x such that
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationAn Efficient, Geometric Multigrid Solver for the Anisotropic Diffusion Equation in Two and Three Dimensions
1 n Efficient, Geometric Multigrid Solver for the nisotropic Diffusion Equation in Two and Three Dimensions Tolga Tasdizen, Ross Whitaker UUSCI-2004-002 Scientific Computing and Imaging Institute University
More informationEuler s Methods (a family of Runge- Ku9a methods)
Euler s Methods (a family of Runge- Ku9a methods) ODE IVP An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable: The equation is coupled with an
More informationImportant Things to Remember on the SOL
Notes Important Things to Remember on the SOL Evaluating Expressions *To evaluate an expression, replace all of the variables in the given problem with the replacement values and use (order of operations)
More informationAn introduction to interpolation and splines
An introduction to interpolation and splines Kenneth H. Carpenter, EECE KSU November 22, 1999 revised November 20, 2001, April 24, 2002, April 14, 2004 1 Introduction Suppose one wishes to draw a curve
More information1. Assumptions. 1. Introduction. 2. Terminology
4. Process Modeling 4. Process Modeling The goal for this chapter is to present the background and specific analysis techniques needed to construct a statistical model that describes a particular scientific
More informationC3 Numerical methods
Verulam School C3 Numerical methods 138 min 108 marks 1. (a) The diagram shows the curve y =. The region R, shaded in the diagram, is bounded by the curve and by the lines x = 1, x = 5 and y = 0. The region
More informationMiddle School Math Course 3
Middle School Math Course 3 Correlation of the ALEKS course Middle School Math Course 3 to the Texas Essential Knowledge and Skills (TEKS) for Mathematics Grade 8 (2012) (1) Mathematical process standards.
More informationIntegrated Math I. IM1.1.3 Understand and use the distributive, associative, and commutative properties.
Standard 1: Number Sense and Computation Students simplify and compare expressions. They use rational exponents and simplify square roots. IM1.1.1 Compare real number expressions. IM1.1.2 Simplify square
More informationOhio Tutorials are designed specifically for the Ohio Learning Standards to prepare students for the Ohio State Tests and end-ofcourse
Tutorial Outline Ohio Tutorials are designed specifically for the Ohio Learning Standards to prepare students for the Ohio State Tests and end-ofcourse exams. Math Tutorials offer targeted instruction,
More informationTruncation Errors. Applied Numerical Methods with MATLAB for Engineers and Scientists, 2nd ed., Steven C. Chapra, McGraw Hill, 2008, Ch. 4.
Chapter 4: Roundoff and Truncation Errors Applied Numerical Methods with MATLAB for Engineers and Scientists, 2nd ed., Steven C. Chapra, McGraw Hill, 2008, Ch. 4. 1 Outline Errors Accuracy and Precision
More informationSecond Triangular Hermite Spline Curves and Its Application
Progress in Applied Mathematics Vol. 4, No. 1, 1, pp. [3 36] DOI: 1.3968/j.pam.19558141.1533 ISSN 195-51X [Print] ISSN 195-58 [Online] www.cscanada.net www.cscanada.org Second Triangular Hermite Spline
More informationExperimental Design and Graphical Analysis of Data
Experimental Design and Graphical Analysis of Data A. Designing a controlled experiment When scientists set up experiments they often attempt to determine how a given variable affects another variable.
More information0607 CAMBRIDGE INTERNATIONAL MATHEMATICS
UNIVERSITY OF CAMBRIDGE INTERNATIONAL EXAMINATIONS International General Certificate of Secondary Education MARK SCHEME for the May/June question paper for the guidance of teachers 7 CAMBRIDGE INTERNATIONAL
More information9.1 Solutions, Slope Fields, and Euler s Method
Math 181 www.timetodare.com 9.1 Solutions, Slope Fields, and Euler s Method In-class work: Model of Population Growth One model for the growth of a population is based on the assumption that the population
More informationStudents are required to have a graphing calculator. Instructors of the course support topics using a TI-84 or TI-89.
AP Calculus AB Course Design and Philosophy Students conceptualize calculus when it is approached in a variety of methods. The course is taught using multiple strategies including algebraic, numerical,
More informationFinite element solution of multi-scale transport problems using the least squares based bubble function enrichment
Finite element solution of multi-scale transport problems using the least squares based bubble function enrichment A. Yazdani a, V. Nassehi b1 a Cranfield University, School of Applied Sciences, Cranfield,
More informationGraphing Functions. 0, < x < 0 1, 0 x < is defined everywhere on R but has a jump discontinuity at x = 0. h(x) =
Graphing Functions Section. of your tetbook is devoted to reviewing a series of steps that you can use to develop a reasonable graph of a function. Here is my version of a list of things to check. You
More informationDepartment of Computer Science, University of Otago
Department of Computer Science, University of Otago Technical Report OUCS-2002-02 Integer root finding, a functional perspective Author: Michael H. Albert Status: submitted to the Journal of Functional
More information3.1 Maxima/Minima Values
3.1 Maxima/Minima Values Ex 1: Find all critical points for the curve given by f (x)=x 5 25 3 x3 +20x 1 on the interval [-3, 2]. Identify the min and max values. We're guaranteed max and min points if
More informationMaximizing an interpolating quadratic
Week 11: Monday, Apr 9 Maximizing an interpolating quadratic Suppose that a function f is evaluated on a reasonably fine, uniform mesh {x i } n i=0 with spacing h = x i+1 x i. How can we find any local
More informationCHAPTER 1. Introduction
ME 475: Computer-Aided Design of Structures 1-1 CHAPTER 1 Introduction 1.1 Analysis versus Design 1.2 Basic Steps in Analysis 1.3 What is the Finite Element Method? 1.4 Geometrical Representation, Discretization
More informationCarnegie Learning Math Series Course 1, A Florida Standards Program. Chapter 1: Factors, Multiples, Primes, and Composites
. Factors and Multiples Carnegie Learning Math Series Course, Chapter : Factors, Multiples, Primes, and Composites This chapter reviews factors, multiples, primes, composites, and divisibility rules. List
More informationLecture VIII. Global Approximation Methods: I
Lecture VIII Global Approximation Methods: I Gianluca Violante New York University Quantitative Macroeconomics G. Violante, Global Methods p. 1 /29 Global function approximation Global methods: function
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationCecil Jones Academy Mathematics Fundamentals
Year 10 Fundamentals Core Knowledge Unit 1 Unit 2 Estimate with powers and roots Calculate with powers and roots Explore the impact of rounding Investigate similar triangles Explore trigonometry in right-angled
More informationHOUGH TRANSFORM CS 6350 C V
HOUGH TRANSFORM CS 6350 C V HOUGH TRANSFORM The problem: Given a set of points in 2-D, find if a sub-set of these points, fall on a LINE. Hough Transform One powerful global method for detecting edges
More informationMetaheuristic Optimization with Evolver, Genocop and OptQuest
Metaheuristic Optimization with Evolver, Genocop and OptQuest MANUEL LAGUNA Graduate School of Business Administration University of Colorado, Boulder, CO 80309-0419 Manuel.Laguna@Colorado.EDU Last revision:
More informationIsn t it Saturday? IMO Trainning Camp NUK, 2004_06_19 2
2004 6 19 Isn t it Saturday? IMO Trainning Camp NUK, 2004_06_19 2 But, I guess you want to be the IMO Trainning Camp NUK, 2004_06_19 3 So, let s go! IMO Trainning Camp NUK, 2004_06_19 4 Outline Brief introduction
More informationMontana City School GRADE 5
Montana City School GRADE 5 Montana Standard 1: Students engage in the mathematical processes of problem solving and reasoning, estimation, communication, connections and applications, and using appropriate
More informationFormula for the asymmetric diffraction peak profiles based on double Soller slit geometry
REVIEW OF SCIENTIFIC INSTRUMENTS VOLUME 69, NUMBER 6 JUNE 1998 Formula for the asymmetric diffraction peak profiles based on double Soller slit geometry Takashi Ida Department of Material Science, Faculty
More informationAttractors and basins of dynamical systems
Electronic Journal of Qualitative Theory of Differential Equations 2011, No. 20, 1-11; http://www.math.u-szeged.hu/ejqtde/ Attractors and basins of dynamical systems Attila Dénes, Géza Makay Abstract There
More informationBezier Curves : Simple Smoothers of Noisy Data
Developments in Statistics and Methodology A. Ferligoj and A. Kramberger (Editors) Metodološki zvezki, 9, Ljubljana : FDV 1993 Bezier Curves : Simple Smoothers of Noisy Data Andrej Blejec1 Abstract In
More informationAn Adaptive Stencil Linear Deviation Method for Wave Equations
211 An Adaptive Stencil Linear Deviation Method for Wave Equations Kelly Hasler Faculty Sponsor: Robert H. Hoar, Department of Mathematics ABSTRACT Wave Equations are partial differential equations (PDEs)
More informationResearch Article Data Visualization Using Rational Trigonometric Spline
Applied Mathematics Volume Article ID 97 pages http://dx.doi.org/.//97 Research Article Data Visualization Using Rational Trigonometric Spline Uzma Bashir and Jamaludin Md. Ali School of Mathematical Sciences
More informationODE IVP. An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable:
Euler s Methods ODE IVP An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable: The equation is coupled with an initial value/condition (i.e., value
More informationALGEBRA II A CURRICULUM OUTLINE
ALGEBRA II A CURRICULUM OUTLINE 2013-2014 OVERVIEW: 1. Linear Equations and Inequalities 2. Polynomial Expressions and Equations 3. Rational Expressions and Equations 4. Radical Expressions and Equations
More informationAIRFOIL SHAPE OPTIMIZATION USING EVOLUTIONARY ALGORITHMS
AIRFOIL SHAPE OPTIMIZATION USING EVOLUTIONARY ALGORITHMS Emre Alpman Graduate Research Assistant Aerospace Engineering Department Pennstate University University Park, PA, 6802 Abstract A new methodology
More informationBig Mathematical Ideas and Understandings
Big Mathematical Ideas and Understandings A Big Idea is a statement of an idea that is central to the learning of mathematics, one that links numerous mathematical understandings into a coherent whole.
More informationFortran 90 Two Commonly Used Statements
Fortran 90 Two Commonly Used Statements 1. DO Loops (Compiled primarily from Hahn [1994]) Lab 6B BSYSE 512 Research and Teaching Methods The DO loop (or its equivalent) is one of the most powerful statements
More informationSecond-order shape optimization of a steel bridge
Computer Aided Optimum Design of Structures 67 Second-order shape optimization of a steel bridge A.F.M. Azevedo, A. Adao da Fonseca Faculty of Engineering, University of Porto, Porto, Portugal Email: alvaro@fe.up.pt,
More informationChapter 23. Linear Motion Motion of a Bug
Chapter 23 Linear Motion The simplest example of a parametrized curve arises when studying the motion of an object along a straight line in the plane We will start by studying this kind of motion when
More informationAnalysis of Replication Control Protocols
Analysis of Replication Control Protocols Darrell D. E. Long University of California, Santa Cruz darrell@cis.ucsc.edu June 22, 2003 Abstract In recent years many replication control protocols have been
More informationIntroduction to the Trigonometric Shooting Reconstruction Method
Technical Note Introduction to the Trigonometric Shooting Reconstruction Method David Rose 1 Todd Ekleberry 2 Gary Wilgus 3 Abstract: This article is an introduction to the use of trigonometry for the
More information