A Class of Linearly Implicit Numerical Methods for Solving Stiff Ordinary Differential Equations
|
|
- Agnes Lynch
- 6 years ago
- Views:
Transcription
1 The Open Numerical Method Journal, 2010, 2, Open Acce A Cla o Linearl Implicit Numerical Method or Solving Sti Ordinar Dierential Equation S.S. Filippov * and A.V. Tglian Keldh Intitute o Applied Mathematic, Ruian Academ o Science, Mocow , Ruia Abtract: We introduce ABC-cheme, a new cla o linearl implicit one-tep method or numerical integration o ti ordinar dierential equation tem. Formula o ABC-cheme invoke the Jacobian o dierential tem imilar to the method o Roenbrock tpe, but unlike the latter the include alo the quare o the Jacobian matrix. Keword: Solving ti ordinar dierential equation, Linearl implicit method, ABC-cheme, A-tabilit, L-tabilit. 1. INTRODUCTION We propoe a new cla o one-tep numerical method or olving ti ordinar dierential equation. Thee method emplo the Jacobian o a dierential tem and, in ditinction rom Roenbrock method [1], the quare o Jacobian i alo involved in their ormula. The irt two onetage method o thi kind were reported b S.S. Filippov and M.V. Bulatov (Conerence on Scientiic Computation, Geneva, Switzerland, June 26-29, 2002, p. 26); ee alo [2]. The term ABC-cheme or uch method wa uggeted later in [3]. In Section 2 one-tage ABC-cheme are deined and ome reult obtained or them are preented. Section 3 contain everal example o one-tage ABC-cheme. Multitage ABC-cheme are introduced in Section 4. Two example o two-tage ABC-cheme are given in Section 5. Some reult o a numerical experiment with ABC-cheme compared with thoe obtained b the ue o implicit Runge- Kutta method are preented in Section ONE-STAGE ABC-SCHEMES Deinition 1 A one-tage ABC-cheme or numerical integration o a Cauch problem or an autonomou tem o n ordinar dierential equation (x) = ((x)), (x 0 ) = 0 (1) i deined a ollow: (I + Ah + Bh 2 2 )[ 1 (h) 0 ] = (I + Ch )h (2) Here, A, B, and C are the coeicient that determine a particular method, 1 (h) i the deired numerical olution ater one tep o integration with the tep ize h, (x) and *Addre correpondence to thi author at the Keldh Intitute o Applied Mathematic o the Ruian Academ o Science, Miukaa q. 4, Mocow , Ruia; Tel: +7(495) ; Fax: +7(499) ; ilippov@ keldh.ru () are n-dimenional vector unction, i the Jacobian matrix, and I i the identit matrix. We conider the irt tep o integration a a repreentative one or the ubequent tep and write,,... without argument or ( 0 ), ( 0 ),.... The ollowing tatement or one-tage ABC-cheme can be eail proved in tandard wa (ee e.g. [4] and [5]). Theorem 1 The convergence order o method (2) i not le then one at an choice o real coeicient A, B, and C. Theorem 2 The order o method (2) equal two, i C = A In thi cae, the principal error term i equal to (x 0 + h) 1 (h) = h3 3! ( + ), where = 1+ 3A + 6B (3) and n 2 = i j k, j, k=1 j k n 2 = i j j, k=1 j k or i = 1, K,n. k Theorem 3 The tabilit unction o ABC-cheme (2) i given b 1+ (1 + A)z + (B + C)z2 1+ Az+ Bz / Bentham Open
2 2 The Open Numerical Method Journal, 2010, Volume 2 Filippov and Tglian Theorem 4 The ABC-cheme (2) o order two are A-table, i A 1 2, B A Theorem 5 The ABC-cheme (2) o order two are L-table, i B = A 1 2. Furthermore, ome important reult or linear autonomou tem ollow immediatel rom the above theorem. Corollar 1 ABC-cheme (2) approximate olution to linear tem (1) having contant coeicient with order three, i B = A 2 1 6, A 1 2. In thi cae, Eq. (3) ield = 0, and we have a amil o method depending on the ingle parameter A: I + Ah A h 2 2 ( (h) ) = h + A h 2 (4) with the principal error term (x 0 + h) 1 (h) = h4 4! (1 + 2A) 3 and tabilit unction 1+ (1 + A)z + A z 2 1+ Az A z 2 Thee method are A-table or A 12, and at A = 23 the method i alo L-table. Corollar 2 With A = 1 2 method (4) take the orm I 1 2 h h2 2 ( 1 (h) 0 ) = h. It give ourth order approximation or the olution o linear autonomou tem o dierential equation. It principal term o local error i then equal to (x 0 + h) 1 (h) = h5 1 5! 6 4. Thi method i A-table with the tabilit unction z z z z2 Remark 1 Solving linear tem o algebraic equation (2) in it original orm eem to be rather expenive. Indeed, in addition to ~ n 3 /3 multiplication and diviion that are needed or the LU- decompoition o the matrix in the let-hand ide 3 o Eq. (2), extra n multiplication are required or quaring a ull matrix o nth order, i.e. totall ~4n 3 /3 multiplicative operation. However, it i poible to avoid matrix multiplication b decompoing the matrix in the let-hand ide o Eq. (2) a ollow: I + Ah + Bh 2 2 = (I + Fh )(I + Gh ) where the new coeicient F and G are real or complex number depending on the value o A and B. In thi cae, onl two LU- decompoition o the matrice (I + Fh ) and (I + Gh ) are needed, i.e. totall ~2n 3 /3 multiplicative operation. Moreover, the amount o arithmetical operation can be once more halved, i we conine ourelve with the choice B = A 2 /4, becaue in thi cae F = G = A 2, and onl one LU- decompoition i needed ( cheap ABCcheme). Fig. (1) illutrate ome eential reult obtained or the one-tage ABC-cheme o order 2. Region A contain all pair o coeicient (A, B) correponding to A-table method (Theorem 4). The thick line L indicate L-table method (Theorem 5). The dahed line i the locu o all the method (4) with = 0 (Corollar 1). The dotted parabola B = A 2 /4 repreent the cheap ABC-cheme (Remark 1). The example rom the next ection are indicated b mall circle with correponding number. 3 EXAMPLES OF ONE-STAGE ABC-SCHEMES Each o the example given below i indicated in Fig. (1) b a mall circle with the number o the correponding example. Example 1 The choice A = 1 2 and B = C = 0 give an A-table method o the orm I h 2 [ 1 (h) 0 ] = h with = 1/ 2 in the principal term o the local error and the tabilit unction z z Actuall, thi i a Roenbrock tpe method, though it i not mentioned in [1]. Example 2 Now let u et A = 1, B = C = 12. In thi cae, we get an L-table method
3 A Cla o Linearl Implicit Numerical Method or Solving The Open Numerical Method Journal, 2010, Volume 2 3 Fig. (1). A-table and L-table one-tage econd order ABC-cheme. I h + h2 2 2 [ 1(h) 0 ] = h h2 2 Eq. (3) give = 1 or thi method, and the tabilit unction o it i given b 1 1 z z2 The method wa derived in other wa and dicued in [2]. Example 3 The choice A = 23, B = C = 16 give the L-table method I 2h 3 + h2 6 2 [ (h) h2 1 0 ] = h 6 with = 0 in the principal term o local error. Thereore, thi method i a member o the amil deined b Eq. (4). It tabilit unction ha the orm z z z2 It wa alo mentioned in [3]. Example 4 The choice A = 1 2, B = 112, and C = 0 give an A- table method decribed in Corollar 2 (ee above). The correponding value o rom Eq. (3) i equal to zero. Thi method i alo a member o the amil decribed b Eq. (4). It give ourth order approximation or the olution o linear autonomou tem o dierential equation. Example 5 With the choice B = A 2 /4, we get cheap ABCcheme that minimize the cot o olving the tem o linear algebraic equation (2) (ee Remark 1 above). Then Eq. (2) take the ollowing orm or the econd order method: 2 I Ah [ (h) ] = h + A h 2 Setting A = give an L-table method. The tabilit unction o thi method i 1+ ( 2 1)z z 2 Eq. (3) ield the value o equal to Example 6 Another example o a cheap ABC-cheme give the choice A = 1 3 1/ Thi time the value o i equal to zero. Thi mean that the method i alo a member o the amil decribed b Eq. (4), and it give third order approximation or the olution o linear autonomou tem o dierential equation. 4. MULTISTAGE ABC-SCHEMES Deinition 2 A multitage ABC-cheme or numerical integration o a Cauch problem or an autonomou tem o n ordinar dierential equation (1) i deined a ollow:
4 4 The Open Numerical Method Journal, 2010, Volume 2 Filippov and Tglian (I + A i h + B i h 2 2 )[u i (h) 0 ] = ( i I + C i h )h (u i1 (h))(i = 1,K.) 1 (h) = i = 1 i u i (h) (5) Here, A i, B i, C i, i, and i are the coeicient that determine a particular method, i the number o tage ( 1), 1 (h) i the deired numerical olution ater one tep o integration with the tep ize h ( 1 (h) i the weighted um o partial olution u i (h) obtained on the ith tage, u 0 (h) 0 ); (x) and () are n-dimenional vector unction, i the Jacobian matrix, and I i the identit matrix. We conider the irt tep o integration a a repreentative one or the ubequent tep and write,,... without argument or ( 0 ), ( 0 ),.... Note that the number o coeicient that deine a particular multitage ABC-cheme in the cae 2 i ubtantiall more then or one-tage ABC-cheme. Thi act enable one to contruct method o order higher then 2, but it i alo the caue o diicultie that encounter in the anali o order condition and tabilit unction. Theorem 6 Stabilit unction o multitage ABC-cheme can be written in the ollowing orm: i R i (z) The tabilit unction R i (z) o equential tage are evaluated recurivel: R i (z) = 1+ P i (z) Q i (z) R i1 (z) (i = 1, K, ) where R 0 (z) = 1, P i (z) = i z + C i z 2, Q i (z) = 1+ A i z + B i z 2 The proo o thi theorem i traightorward. One ha to appl the ormula rom Deinition 2 to Dahlquit tet equation ' =, where i a complex number (ee e.g. [5]), and then put h = z. 5. EXAMPLES OF TWO-STAGE ABC-SCHEMES The gain o uing cheap ABC-cheme (ee Remark 1 above) become till more in the cae o multitage ABCcheme. I we put A i = A, B i = B or all tage and aume B = A 2 4, then onl a ingle LU-decompoition will be needed on each tep o integration, ince the matrix in the let-hand ide o Eq. (5) i the ame or all tage. For thi reaon, we conine ourelve with two example o cheap two-tage ABC-cheme. Example 1 The choice 1 = 2 = 1, 1 = 2/3, 2 = 1/3, A A = A B 1 = B 2 = A 2 /4 give a amil o two-tage third order method depending on a ingle parameter A. In thi cae, C 1 = 3 4 A A C 2 = 3 2 A2 + 2A The tabilit unction at z take the orm 1 = 2, R() = A A 3 Thee method are A-table at the value o A between approximatel 0.75 and 0.4. The value A 0.59 correpond to an L-table method. Example 2 The choice 1 = 1/ 3, 2 = 1, 1 = 0, 2 = 1, A 1 = A 2 = A, B 1 = B 2 = A 2 /4 give again a amil o twotage third order method depending on a ingle parameter A, but in thi cae C 1 = 1 4 A A and the tabilit unction at C 2 = A z now take the orm R() = A A A A 4 Further reult or thee method will be preented elewhere. 6. NUMERICAL EXPERIMENT For our numerical experiment, we have choen a particular cae o the ingularl perturbed tet problem uggeted b Kap [6], namel, the initial value problem 1 (x) = (2 + 1 ) 1 (x) (x), 2 (x) = 1 (x) 2 (x) 2 2 (x), 1 (0) = 2 (0) = 1, 0 x 1. The exact olution o thi problem 1 (x) = e 2 x, 2 (x) = e x doe not depend on. However, the problem become ver ti, a 0. We compare the reult o numerical integration perormed with the ue o our method: - method 1 i the one-tage ABC-cheme rom Example 3 o Section 3; - method 2 i the implicit midpoint rule (one-tage Gau method [4, 5]); - method 3 i the two-tage cheap ABC-cheme rom Example 1 o Section 5 with A = 0.59 ; - method 4 i the two-tage Gau method [5].
5 A Cla o Linearl Implicit Numerical Method or Solving The Open Numerical Method Journal, 2010, Volume 2 5 Table 1. The Comparion o Reult Obtained with the Ue o ABC-Scheme and Gau Method Method 1 Method 2 Method 3 Method Method 1, 2, 3, and 4 have claical order 2, 2, 3, and 4, repectivel. Thee method were ued or the numerical integration o the above problem with everal diminihing value o and with two contant value o tep ize, h = 1/40 and h = 1/80. Table 1 contain the ollowing value: i the Euclidean norm o the abolute value o error or h = 1/80at the endpoint o the integration interval; = log 2 ( e 40 / e 80 ) i the actual order o accurac etimated uing the reult o integration with h = 1/ 40 and h = 1/80. In the cae o method 1 and 3 the computation wa perormed uing double preciion. In the cae o method 2 and 4 the data (evaluated with comparable preciion) are taken rom Table in [7]. Oberve that, or mall value o, ABC-cheme give better reult than implicit Runge-Kutta method. Note that the implicit Runge-Kutta method emplo Newton iteration, i.e. the are more expenive then the ABCcheme. One can clearl ee the phenomenon o lowering o the actual order o accurac at mall value o or the method 3 and 4, which i in accordance with the theor o B-convergence [5, 7]. ACKNOWLEGEMENTS The author are grateul to Proeor J.G. Verwer or hi encouraging criticim and to Proeor E. Hairer or drawing our attention to a poible wa to eliminate matrix multiplication. We are alo thankul to Dr M.V. Bulatov or a ueul dicuion o one-tage ABC-cheme and to Dr G.Yu. Kulikov or hi valuable obervation. REFERENCES [1] Roenbrock HH. Some general implicit procee or the numerical olution o dierential equation. Comput J 1962/63; 5: [2] Bulatov MV. Contruction o a one-tage L-table econd-order method. Dierential Equation 2003; 39: [3] Filippov SS. ABC-cheme or ti tem o ordinar dierential equation. Doklad Mathematic 2004; 70: [4] Hairer E, Nørett SP, Wanner G. Solving ordinar dierential equation, I: Nonti problem. Springer-Verlag: Berlin [5] Hairer E, Wanner G. Solving ordinar dierential equation, II: Sti and dierential-algebraic problem. Springer-Verlag: Berlin [6] Kap P. Roenbrock-tpe method. In: Dahlquit G, Jeltch R. Ed., Numerical method or olving ti initial value problem. Int. ür Geometrie und praktiche Math. der RWTH Aachen 1981; Bericht No. 9. [7] Dekker K, Verwer JG. Stabilit o Runge-Kutta method or ti nonlinear dierential equation. North-Holland: Amterdam Received: October 30, 2009 Revied: Januar 28, 2010 Accepted: Februar 02, 2010 Filippov and Tglian; Licenee Bentham Open. Thi i an open acce article licened under the term o the Creative Common Attribution Non-Commercial Licene ( which permit unretricted, non-commercial ue, ditribution and reproduction in an medium, provided the work i properl cited.
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS
CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3
More information7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281
72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition
More informationPosition. If the particle is at point (x, y, z) on the curved path s shown in Fig a,then its location is defined by the position vector
34 C HAPTER 1 KINEMATICS OF A PARTICLE 1 1.5 Curvilinear Motion: Rectangular Component Occaionall the motion of a particle can bet be decribed along a path that can be epreed in term of it,, coordinate.
More informationSocial Studies 201 Notes for November 14, 2003
1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationSocial Studies 201 Notes for March 18, 2005
1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the
More informationTarzan s Dilemma for Elliptic and Cycloidal Motion
Tarzan Dilemma or Elliptic and Cycloidal Motion Yuji Kajiyama National Intitute o Technology, Yuge College, Shimo-Yuge 000, Yuge, Kamijima, Ehime, 794-593, Japan kajiyama@gen.yuge.ac.jp btract-in thi paper,
More informationOBSERVER-BASED REDUCED ORDER CONTROLLER DESIGN FOR THE STABILIZATION OF LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS
International Journal o Computer Science, Engineering and Inormation Technology (IJCSEIT, Vol.1, No.5, December 2011 OBSERVER-BASED REDUCED ORDER CONTROLLER DESIGN FOR THE STABILIZATION OF LARGE SCALE
More informationKey Mathematical Backgrounds
Ke Mathematical Background Dierential Equation Ordinar Linear Partial Nonlinear: mooth, nonmooth,, piecewie linear Map Linear Nonlinear Equilibrium/Stead-State State Solution Linearization Traner Function
More informationMATEMATIK Datum: Tid: eftermiddag. A.Heintz Telefonvakt: Anders Martinsson Tel.:
MATEMATIK Datum: 20-08-25 Tid: eftermiddag GU, Chalmer Hjälpmedel: inga A.Heintz Telefonvakt: Ander Martinon Tel.: 073-07926. Löningar till tenta i ODE och matematik modellering, MMG5, MVE6. Define what
More informationOne Class of Splitting Iterative Schemes
One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi
More informationOne Dimensional Modeling of the Shape Memory Effect
odeling and Numerical Simulation o aterial Science, 3, 3, 4-8 http://dx.doi.org/.436/mnm.3.347 Publihed Online October 3 (http://www.cirp.org/journal/mnm) One Dimenional odeling o the Shape emory Eect
More informationChapter 4. The Laplace Transform Method
Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination
More informationA Simplified Methodology for the Synthesis of Adaptive Flight Control Systems
A Simplified Methodology for the Synthei of Adaptive Flight Control Sytem J.ROUSHANIAN, F.NADJAFI Department of Mechanical Engineering KNT Univerity of Technology 3Mirdamad St. Tehran IRAN Abtract- A implified
More informationAlternate Dispersion Measures in Replicated Factorial Experiments
Alternate Diperion Meaure in Replicated Factorial Experiment Neal A. Mackertich The Raytheon Company, Sudbury MA 02421 Jame C. Benneyan Northeatern Univerity, Boton MA 02115 Peter D. Krau The Raytheon
More informationApproximate Analytical Solution for Quadratic Riccati Differential Equation
Iranian J. of Numerical Analyi and Optimization Vol 3, No. 2, 2013), pp 21-31 Approximate Analytical Solution for Quadratic Riccati Differential Equation H. Aminikhah Abtract In thi paper, we introduce
More informationImplicit Second Derivative Hybrid Linear Multistep Method with Nested Predictors for Ordinary Differential Equations
American Scientiic Research Journal or Engineering, Technolog, and Sciences (ASRJETS) ISSN (Print) -44, ISSN (Online) -44 Global Societ o Scientiic Research and Researchers http://asretsournal.org/ Implicit
More informationTransfer Function Approach to the Model Matching Problem of Nonlinear Systems
Tranfer Function Approach to the Model Matching Problem of Nonlinear Stem Mirolav Halá Ülle Kotta Claude H. Moog Intitute of Control and Indutrial Informatic, Fac. of Electrical Engineering and IT, Slovak
More informationNumerical algorithm for the analysis of linear and nonlinear microstructure fibres
Numerical algorithm for the anali of linear and nonlinear microtructure fibre Mariuz Zdanowicz *, Marian Marciniak, Marek Jaworki, Igor A. Goncharenko National Intitute of Telecommunication, Department
More informationOptimal Coordination of Samples in Business Surveys
Paper preented at the ICES-III, June 8-, 007, Montreal, Quebec, Canada Optimal Coordination of Sample in Buine Survey enka Mach, Ioana Şchiopu-Kratina, Philip T Rei, Jean-Marc Fillion Statitic Canada New
More informationClustering Methods without Given Number of Clusters
Clutering Method without Given Number of Cluter Peng Xu, Fei Liu Introduction A we now, mean method i a very effective algorithm of clutering. It mot powerful feature i the calability and implicity. However,
More informationCHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL
98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i
More information1. The F-test for Equality of Two Variances
. The F-tet for Equality of Two Variance Previouly we've learned how to tet whether two population mean are equal, uing data from two independent ample. We can alo tet whether two population variance are
More informationAn Approach for Solving Multi-Objective Linear Fractional Programming Problem and It s Comparison with Other Techniques
International Journal of Scientific and Innovative Mathematical Reearch (IJSIMR) Volume 5, Iue 11, 2017, PP 1-5 ISSN 2347-307X (Print) & ISSN 2347-3142 (Online) DOI: http://dxdoiorg/1020431/2347-31420511001
More informationMoment of Inertia of an Equilateral Triangle with Pivot at one Vertex
oment of nertia of an Equilateral Triangle with Pivot at one Vertex There are two wa (at leat) to derive the expreion f an equilateral triangle that i rotated about one vertex, and ll how ou both here.
More informationFinite Difference Formulae for Unequal Sub- Intervals Using Lagrange s Interpolation Formula
Int. Journal o Mat. Analyi, Vol., 9, no. 7, 85-87 Finite Dierence Formulae or Unequal Sub- Interval Uing Lagrange Interpolation Formula Aok K. Sing a and B. S. Badauria b Department o Matematic, Faculty
More informationSuggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall
Beyond Significance Teting ( nd Edition), Rex B. Kline Suggeted Anwer To Exercie Chapter. The tatitic meaure variability among core at the cae level. In a normal ditribution, about 68% of the core fall
More informationIntroduction to Laplace Transform Techniques in Circuit Analysis
Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found
More informationDynamical Behavior Analysis and Control of a Fractional-order Discretized Tumor Model
6 International Conference on Information Engineering and Communication Technolog (IECT 6) ISBN: 978--6595-375-5 Dnamical Behavior Anali and Control of a Fractional-order Dicretied Tumor Model Yaling Zhang,a,
More informationSource slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis
Source lideplayer.com/fundamental of Analytical Chemitry, F.J. Holler, S.R.Crouch Chapter 6: Random Error in Chemical Analyi Random error are preent in every meaurement no matter how careful the experimenter.
More informationA Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems
A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement
More informationSTRAIN LIMITS FOR PLASTIC HINGE REGIONS OF CONCRETE REINFORCED COLUMNS
13 th World Conerence on Earthquake Engineering Vancouver, B.C., Canada Augut 1-6, 004 Paper No. 589 STRAIN LIMITS FOR PLASTIC HINGE REGIONS OF CONCRETE REINFORCED COLUMNS Rebeccah RUSSELL 1, Adolo MATAMOROS,
More informationApplication of Extended Scaling Law to the Surface Tension of Fluids of Wide Range of Molecular Shapes
Application o Extended caling Law to the urace enion o Fluid o Wide Range o Molecular hape Mohammad Hadi Ghatee, Ali oorghali (Department o Chemitry, College o cience, hiraz Univerity, hiraz 71454, Iran)
More informationEvolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis
Proceeding of 01 4th International Conference on Machine Learning and Computing IPCSIT vol. 5 (01) (01) IACSIT Pre, Singapore Evolutionary Algorithm Baed Fixed Order Robut Controller Deign and Robutne
More informationUnified Design Method for Flexure and Debonding in FRP Retrofitted RC Beams
Unified Deign Method for Flexure and Debonding in FRP Retrofitted RC Beam G.X. Guan, Ph.D. 1 ; and C.J. Burgoyne 2 Abtract Flexural retrofitting of reinforced concrete (RC) beam uing fibre reinforced polymer
More informationLecture 7: Testing Distributions
CSE 5: Sublinear (and Streaming) Algorithm Spring 014 Lecture 7: Teting Ditribution April 1, 014 Lecturer: Paul Beame Scribe: Paul Beame 1 Teting Uniformity of Ditribution We return today to property teting
More informationON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS. Volker Ziegler Technische Universität Graz, Austria
GLASNIK MATEMATIČKI Vol. 1(61)(006), 9 30 ON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS Volker Ziegler Techniche Univerität Graz, Autria Abtract. We conider the parameterized Thue
More informationLaplace Adomian Decomposition Method for Solving the Nonlinear Volterra Integral Equation with Weakly Kernels
Studie in Nonlinear Science (4): 9-4, ISSN -9 IDOSI Publication, Laplace Adomian Decompoition Method for Solving the Nonlinear Volterra Integral Equation with Weakly Kernel F.A. Hendi Department of Mathematic
More informationON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS
ON A CERTAIN FAMILY OF QUARTIC THUE EQUATIONS WITH THREE PARAMETERS VOLKER ZIEGLER Abtract We conider the parameterized Thue equation X X 3 Y (ab + (a + bx Y abxy 3 + a b Y = ±1, where a, b 1 Z uch that
More informationReformulation of Block Implicit Linear Multistep Method into Runge Kutta Type Method for Initial Value Problem
International Journal of Science and Technology Volume 4 No. 4, April, 05 Reformulation of Block Implicit Linear Multitep Method into Runge Kutta Type Method for Initial Value Problem Muhammad R., Y. A
More informationAn Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem
An Inequality for Nonnegative Matrice and the Invere Eigenvalue Problem Robert Ream Program in Mathematical Science The Univerity of Texa at Dalla Box 83688, Richardon, Texa 7583-688 Abtract We preent
More informationIEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation
IEOR 316: Fall 213, Profeor Whitt Topic for Dicuion: Tueday, November 19 Alternating Renewal Procee and The Renewal Equation 1 Alternating Renewal Procee An alternating renewal proce alternate between
More informationCHAPTER 6. Estimation
CHAPTER 6 Etimation Definition. Statitical inference i the procedure by which we reach a concluion about a population on the bai of information contained in a ample drawn from that population. Definition.
More informationJump condition at the boundary between a porous catalyst and a homogeneous fluid
From the SelectedWork of Francico J. Valde-Parada 2005 Jump condition at the boundary between a porou catalyt and a homogeneou fluid Francico J. Valde-Parada J. Alberto Ochoa-Tapia Available at: http://work.bepre.com/francico_j_valde_parada/12/
More informationSTOCHASTIC GENERALIZED TRANSPORTATION PROBLEM WITH DISCRETE DISTRIBUTION OF DEMAND
OPERATIONS RESEARCH AND DECISIONS No. 4 203 DOI: 0.5277/ord30402 Marcin ANHOLCER STOCHASTIC GENERALIZED TRANSPORTATION PROBLEM WITH DISCRETE DISTRIBUTION OF DEMAND The generalized tranportation problem
More informationAn Efficient Class of Estimators for the Finite Population Mean in Ranked Set Sampling
Open Journal of Statitic 06 6 46-435 Publihed Online June 06 in SciRe http://cirporg/journal/oj http://ddoiorg/0436/oj0663038 An Efficient Cla of Etimator for the Finite Population Mean in Ranked Set Sampling
More informationComparing Means: t-tests for Two Independent Samples
Comparing ean: t-tet for Two Independent Sample Independent-eaure Deign t-tet for Two Independent Sample Allow reearcher to evaluate the mean difference between two population uing data from two eparate
More informationTheoretical study of the dual harmonic system and its application on the CSNS/RCS
Theoretical tudy of the dual harmonic ytem and it application on the CSNS/RCS Yao-Shuo Yuan, Na Wang, Shou-Yan Xu, Yue Yuan, and Sheng Wang Dongguan branch, Intitute of High Energy Phyic, CAS, Guangdong
More information1. Preliminaries. In [8] the following odd looking integral evaluation is obtained.
June, 5. Revied Augut 8th, 5. VA DER POL EXPASIOS OF L-SERIES David Borwein* and Jonathan Borwein Abtract. We provide concie erie repreentation for variou L-erie integral. Different technique are needed
More informationON ASYMPTOTIC FORMULA OF THE PARTITION FUNCTION p A (n)
#A2 INTEGERS 15 (2015) ON ASYMPTOTIC FORMULA OF THE PARTITION FUNCTION p A (n) A David Chritopher Department of Mathematic, The American College, Tamilnadu, India davchrame@yahoocoin M Davamani Chritober
More informationμ + = σ = D 4 σ = D 3 σ = σ = All units in parts (a) and (b) are in V. (1) x chart: Center = μ = 0.75 UCL =
Our online Tutor are available 4*7 to provide Help with Proce control ytem Homework/Aignment or a long term Graduate/Undergraduate Proce control ytem Project. Our Tutor being experienced and proficient
More informationYoram Gat. Technical report No. 548, March Abstract. A classier is said to have good generalization ability if it performs on
A bound concerning the generalization ability of a certain cla of learning algorithm Yoram Gat Univerity of California, Berkeley Technical report No. 548, March 999 Abtract A claier i aid to have good
More informationProblem Set 8 Solutions
Deign and Analyi of Algorithm April 29, 2015 Maachuett Intitute of Technology 6.046J/18.410J Prof. Erik Demaine, Srini Devada, and Nancy Lynch Problem Set 8 Solution Problem Set 8 Solution Thi problem
More information4-4 E-field Calculations using Coulomb s Law
1/21/24 ection 4_4 -field calculation uing Coulomb Law blank.doc 1/1 4-4 -field Calculation uing Coulomb Law Reading Aignment: pp. 9-98 1. xample: The Uniform, Infinite Line Charge 2. xample: The Uniform
More informationSolving Differential Equations by the Laplace Transform and by Numerical Methods
36CH_PHCalter_TechMath_95099 3//007 :8 PM Page Solving Differential Equation by the Laplace Tranform and by Numerical Method OBJECTIVES When you have completed thi chapter, you hould be able to: Find the
More informationOnline supplementary information
Electronic Supplementary Material (ESI) for Soft Matter. Thi journal i The Royal Society of Chemitry 15 Online upplementary information Governing Equation For the vicou flow, we aume that the liquid thickne
More informationSimple Observer Based Synchronization of Lorenz System with Parametric Uncertainty
IOSR Journal of Electrical and Electronic Engineering (IOSR-JEEE) ISSN: 78-676Volume, Iue 6 (Nov. - Dec. 0), PP 4-0 Simple Oberver Baed Synchronization of Lorenz Sytem with Parametric Uncertainty Manih
More informationImproving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization
1976 MONTHLY WEATHER REVIEW VOLUME 15 Improving the Efficiency of a Digital Filtering Scheme for Diabatic Initialization PETER LYNCH Met Éireann, Dublin, Ireland DOMINIQUE GIARD CNRM/GMAP, Météo-France,
More informationCake ltration analysis the eect of the relationship between the pore liquid pressure and the cake compressive stress
Chemical Engineering Science 56 (21) 5361 5369 www.elevier.com/locate/ce Cake ltration analyi the eect of the relationhip between the pore liquid preure and the cake compreive tre C. Tien, S. K. Teoh,
More informationSimulation Study on the Shock Properties of the Double-Degree-of-Freedom Cushioning Packaging System
Proceeding of the 7th IAPRI World Conference on Packaging Simulation Study on the Shock Propertie of the Double-Degree-of-Freedom Cuhioning Packaging Sytem Xia Zhu, Qiaoqiao Yan, Xiaoling Yao, Junbin Chen,
More informationSingular Value Inequalities for Compact Normal Operators
dvance in Linear lgebra & Matrix Theory, 3, 3, 34-38 Publihed Online December 3 (http://www.cirp.org/ournal/alamt) http://dx.doi.org/.436/alamt.3.347 Singular Value Inequalitie for Compact Normal Operator
More informationWhite Rose Research Online URL for this paper: Version: Accepted Version
Thi i a repoitory copy of Identification of nonlinear ytem with non-peritent excitation uing an iterative forward orthogonal leat quare regreion algorithm. White Roe Reearch Online URL for thi paper: http://eprint.whiteroe.ac.uk/107314/
More informationON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang
Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang
More informationPractice Problems - Week #7 Laplace - Step Functions, DE Solutions Solutions
For Quetion -6, rewrite the piecewie function uing tep function, ketch their graph, and find F () = Lf(t). 0 0 < t < 2. f(t) = (t 2 4) 2 < t In tep-function form, f(t) = u 2 (t 2 4) The graph i the olid
More informationA SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES. Sanghyun Cho
A SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES Sanghyun Cho Abtract. We prove a implified verion of the Nah-Moer implicit function theorem in weighted Banach pace. We relax the
More informationIf Y is normally Distributed, then and 2 Y Y 10. σ σ
ull Hypothei Significance Teting V. APS 50 Lecture ote. B. Dudek. ot for General Ditribution. Cla Member Uage Only. Chi-Square and F-Ditribution, and Diperion Tet Recall from Chapter 4 material on: ( )
More informationLecture 4 Topic 3: General linear models (GLMs), the fundamentals of the analysis of variance (ANOVA), and completely randomized designs (CRDs)
Lecture 4 Topic 3: General linear model (GLM), the fundamental of the analyi of variance (ANOVA), and completely randomized deign (CRD) The general linear model One population: An obervation i explained
More informationin a circular cylindrical cavity K. Kakazu Department of Physics, University of the Ryukyus, Okinawa , Japan Y. S. Kim
Quantization of electromagnetic eld in a circular cylindrical cavity K. Kakazu Department of Phyic, Univerity of the Ryukyu, Okinawa 903-0, Japan Y. S. Kim Department of Phyic, Univerity of Maryland, College
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world leading publiher of Open Acce book Built by cientit, for cientit 3,5 8,.7 M Open acce book available International author and editor Download Our author are among the 5 Countrie
More informationTRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL
GLASNIK MATEMATIČKI Vol. 38583, 73 84 TRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL p-laplacian Haihen Lü, Donal O Regan and Ravi P. Agarwal Academy of Mathematic and Sytem Science, Beijing, China, National
More informationDimensional Analysis A Tool for Guiding Mathematical Calculations
Dimenional Analyi A Tool for Guiding Mathematical Calculation Dougla A. Kerr Iue 1 February 6, 2010 ABSTRACT AND INTRODUCTION In converting quantitie from one unit to another, we may know the applicable
More informationColorado School of Mines. Computer Vision. Professor William Hoff Dept of Electrical Engineering &Computer Science.
Proeor William Ho Dept o Electrical Engineering &Computer Science http://inide.mine.edu/~who/ Uncertaint Uncertaint Let a that we have computed a reult (uch a poe o an object), rom image data How do we
More informationCS 170: Midterm Exam II University of California at Berkeley Department of Electrical Engineering and Computer Sciences Computer Science Division
1 1 April 000 Demmel / Shewchuk CS 170: Midterm Exam II Univerity of California at Berkeley Department of Electrical Engineering and Computer Science Computer Science Diviion hi i a cloed book, cloed calculator,
More informationConvex Hulls of Curves Sam Burton
Convex Hull of Curve Sam Burton 1 Introduction Thi paper will primarily be concerned with determining the face of convex hull of curve of the form C = {(t, t a, t b ) t [ 1, 1]}, a < b N in R 3. We hall
More informationCodes Correcting Two Deletions
1 Code Correcting Two Deletion Ryan Gabry and Frederic Sala Spawar Sytem Center Univerity of California, Lo Angele ryan.gabry@navy.mil fredala@ucla.edu Abtract In thi work, we invetigate the problem of
More informationEfficiency Analysis of a Multisectoral Economic System
Efficienc Anali of a Multiectoral Economic Stem Mikuláš uptáčik Vienna Univerit of Economic and Buine Vienna, Autria Bernhard Böhm Vienna Univerit of Technolog Vienna,Autria Paper prepared for the 7 th
More informationTHE DIVERGENCE-FREE JACOBIAN CONJECTURE IN DIMENSION TWO
THE DIVERGENCE-FREE JACOBIAN CONJECTURE IN DIMENSION TWO J. W. NEUBERGER Abtract. A pecial cae, called the divergence-free cae, of the Jacobian Conjecture in dimenion two i proved. Thi note outline an
More informationMath 273 Solutions to Review Problems for Exam 1
Math 7 Solution to Review Problem for Exam True or Fale? Circle ONE anwer for each Hint: For effective tudy, explain why if true and give a counterexample if fale (a) T or F : If a b and b c, then a c
More informationA CATEGORICAL CONSTRUCTION OF MINIMAL MODEL
A ATEGORIAL ONSTRUTION OF MINIMAL MODEL A. Behera, S. B. houdhury M. Routaray Department of Mathematic National Intitute of Technology ROURKELA - 769008 (India) abehera@nitrkl.ac.in 512ma6009@nitrkl.ac.in
More information(b) Is the game below solvable by iterated strict dominance? Does it have a unique Nash equilibrium?
14.1 Final Exam Anwer all quetion. You have 3 hour in which to complete the exam. 1. (60 Minute 40 Point) Anwer each of the following ubquetion briefly. Pleae how your calculation and provide rough explanation
More informationJan Purczyński, Kamila Bednarz-Okrzyńska Estimation of the shape parameter of GED distribution for a small sample size
Jan Purczyńki, Kamila Bednarz-Okrzyńka Etimation of the hape parameter of GED ditribution for a mall ample ize Folia Oeconomica Stetinenia 4()/, 35-46 04 Folia Oeconomica Stetinenia DOI: 0.478/foli-04-003
More informationCalculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient
Ŕ periodica polytechnica Mechanical Engineering 54/1 21 15 2 doi: 1.3311/pp.me.21-1.3 web: http:// www.pp.bme.hu/ me c Periodica Polytechnica 21 RESERCH RTICLE Calculation of the temperature of boundary
More informationManagement, Nakhon Pathom Rajabhat University, 85 Malaiman Road, Muang, Nakhon Pathom 73000, Kingdom of Thailand
The Importance o Denity Dependent Flow and Solute Tranport Modeling to imulate Seawater Intruion into a Coatal Aquier Sytem Phatcharaak Arlai 1 and Manred Koch 2 1) Dr.-Ing., Head, Reearch Unit or Sutainable
More informationPreemptive scheduling on a small number of hierarchical machines
Available online at www.ciencedirect.com Information and Computation 06 (008) 60 619 www.elevier.com/locate/ic Preemptive cheduling on a mall number of hierarchical machine György Dóa a, Leah Eptein b,
More informationBUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL
J. Autral. Math. Soc. Ser. B 4(999), 332 349 BUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL J. LEE and J.-M. VANDEN-BROECK 2 (Received 22 April 995; revied 23 April 996)
More informationAvoiding Forbidden Submatrices by Row Deletions
Avoiding Forbidden Submatrice by Row Deletion Sebatian Wernicke, Jochen Alber, Jen Gramm, Jiong Guo, and Rolf Niedermeier Wilhelm-Schickard-Intitut für Informatik, niverität Tübingen, Sand 13, D-72076
More informationBogoliubov Transformation in Classical Mechanics
Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How
More informationinto a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get
Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Laplace Tranform Paul Dawkin Table of Content Preface... Laplace Tranform... Introduction... The Definition... 5 Laplace Tranform... 9 Invere Laplace Tranform... Step Function...4
More informationStandard Guide for Conducting Ruggedness Tests 1
Deignation: E 69 89 (Reapproved 996) Standard Guide for Conducting Ruggedne Tet AMERICA SOCIETY FOR TESTIG AD MATERIALS 00 Barr Harbor Dr., Wet Conhohocken, PA 948 Reprinted from the Annual Book of ASTM
More informationFinite-Dimensional Control of Parabolic PDE Systems Using Approximate Inertial Manifolds
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 16, 39840 1997 ARTICLE NO AY975649 Finite-Dimenional Control o Parabolic PDE Sytem Uing Approximate Inertial Maniold Panagioti D Chritoide and Prodromo
More informationStochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions
Stochatic Optimization with Inequality Contraint Uing Simultaneou Perturbation and Penalty Function I-Jeng Wang* and Jame C. Spall** The John Hopkin Univerity Applied Phyic Laboratory 11100 John Hopkin
More informationLecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004
18.997 Topic in Combinatorial Optimization April 29th, 2004 Lecture 21 Lecturer: Michel X. Goeman Scribe: Mohammad Mahdian 1 The Lovaz plitting-off lemma Lovaz plitting-off lemma tate the following. Theorem
More informationInference for Two Stage Cluster Sampling: Equal SSU per PSU. Projections of SSU Random Variables on Each SSU selection.
Inference for Two Stage Cluter Sampling: Equal SSU per PSU Projection of SSU andom Variable on Eac SSU election By Ed Stanek Introduction We review etimating equation for PSU mean in a two tage cluter
More informationEfficient Methods of Doppler Processing for Coexisting Land and Weather Clutter
Efficient Method of Doppler Proceing for Coexiting Land and Weather Clutter Ça gatay Candan and A Özgür Yılmaz Middle Eat Technical Univerity METU) Ankara, Turkey ccandan@metuedutr, aoyilmaz@metuedutr
More informationMolecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions
Original Paper orma, 5, 9 7, Molecular Dynamic Simulation of Nonequilibrium Effect ociated with Thermally ctivated Exothermic Reaction Jerzy GORECKI and Joanna Natalia GORECK Intitute of Phyical Chemitry,
More informationNumerical Approach for Solving Stiff Differential Equations Through the Extended Trapezoidal Rule Formulae
Aerican Journal o Matheatical and Coputer Modelling 7; (4): -6 http://wwwciencepublihinggroupco//ac doi: 648/ac74 Nuerical Approach or Solving Sti Dierential Equation Through the Extended Trapezoidal Rule
More informationList coloring hypergraphs
Lit coloring hypergraph Penny Haxell Jacque Vertraete Department of Combinatoric and Optimization Univerity of Waterloo Waterloo, Ontario, Canada pehaxell@uwaterloo.ca Department of Mathematic Univerity
More informationAdvanced methods for ODEs and DAEs
Lecture : Implicit Runge Kutta method Bojana Roić, 9. April 7 What you need to know before thi lecture numerical integration: Lecture from ODE iterative olver: Lecture 5-8 from ODE 9. April 7 Bojana Roić
More informationCOHOMOLOGY AS A LOCAL-TO-GLOBAL BRIDGE
COHOMOLOGY AS A LOCAL-TO-GLOBAL BRIDGE LIVIU I. NICOLAESCU ABSTRACT. I dicu low dimenional incarnation of cohomology and illutrate how baic cohomological principle lead to a proof of Sperner lemma. CONTENTS.
More informationOn The Approximate Solution of Linear Fuzzy Volterra-Integro Differential Equations of the Second Kind
AUSTRALIAN JOURNAL OF BASIC AND APPLIED SCIENCES ISSN:99-878 EISSN: 39-8 Journal ome page: www.ajbaweb.com On Te Approimate Solution of Linear Fuzzy Volterra-Integro Differential Equation of te Second
More information