Numerical Analysis Formulae Booklet

Similar documents
Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations

Lecture 12 APPROXIMATION OF FIRST ORDER DERIVATIVES

Chapter 9 Jordan Block Matrices

MOLECULAR VIBRATIONS

Beam Warming Second-Order Upwind Method

BERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS. Aysegul Akyuz Dascioglu and Nese Isler

Assignment 5/MATH 247/Winter Due: Friday, February 19 in class (!) (answers will be posted right after class)

A conic cutting surface method for linear-quadraticsemidefinite

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

UNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS

Lecture 07: Poles and Zeros

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II

Basic Concepts in Numerical Analysis November 6, 2017

Assignment 7/MATH 247/Winter, 2010 Due: Friday, March 19. Powers of a square matrix

1 Lyapunov Stability Theory

A Conventional Approach for the Solution of the Fifth Order Boundary Value Problems Using Sixth Degree Spline Functions

General Method for Calculating Chemical Equilibrium Composition

4 Inner Product Spaces

AN UPPER BOUND FOR THE PERMANENT VERSUS DETERMINANT PROBLEM BRUNO GRENET

Fourth Order Four-Stage Diagonally Implicit Runge-Kutta Method for Linear Ordinary Differential Equations ABSTRACT INTRODUCTION

COV. Violation of constant variance of ε i s but they are still independent. The error term (ε) is said to be heteroscedastic.

u(x, t) = u 0 (x ct). This Riemann invariant u is constant along characteristics λ with x = x 0 +ct (u(x, t) = u 0 (x 0 )):

means the first term, a2 means the term, etc. Infinite Sequences: follow the same pattern forever.

TESTS BASED ON MAXIMUM LIKELIHOOD

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory

EVALUATION OF FUNCTIONAL INTEGRALS BY MEANS OF A SERIES AND THE METHOD OF BOREL TRANSFORM

Initial-Value Problems for ODEs. numerical errors (round-off and truncation errors) Consider a perturbed system: dz dt

F. Inequalities. HKAL Pure Mathematics. 進佳數學團隊 Dr. Herbert Lam 林康榮博士. [Solution] Example Basic properties

Analysis of Lagrange Interpolation Formula

Numerical Solution of Linear Second Order Ordinary Differential Equations with Mixed Boundary Conditions by Galerkin Method

G S Power Flow Solution

Dr. Shalabh. Indian Institute of Technology Kanpur

Bivariate Vieta-Fibonacci and Bivariate Vieta-Lucas Polynomials

18.413: Error Correcting Codes Lab March 2, Lecture 8

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

h-analogue of Fibonacci Numbers

( q Modal Analysis. Eigenvectors = Mode Shapes? Eigenproblem (cont) = x x 2 u 2. u 1. x 1 (4.55) vector and M and K are matrices.

MATH 247/Winter Notes on the adjoint and on normal operators.

CHAPTER 4 RADICAL EXPRESSIONS

QR Factorization and Singular Value Decomposition COS 323

Mu Sequences/Series Solutions National Convention 2014

d dt d d dt dt Also recall that by Taylor series, / 2 (enables use of sin instead of cos-see p.27 of A&F) dsin

( ) ( ) A number of the form x+iy, where x & y are integers and i = 1 is called a complex number.

Stability For a stable numerical scheme, the errors in the initial condition will not grow unboundedly with time.

α1 α2 Simplex and Rectangle Elements Multi-index Notation of polynomials of degree Definition: The set P k will be the set of all functions:

Supplementary Material for Limits on Sparse Support Recovery via Linear Sketching with Random Expander Matrices

1 0, x? x x. 1 Root finding. 1.1 Introduction. Solve[x^2-1 0,x] {{x -1},{x 1}} Plot[x^2-1,{x,-2,2}] 3

Solution of General Dual Fuzzy Linear Systems. Using ABS Algorithm

2006 Jamie Trahan, Autar Kaw, Kevin Martin University of South Florida United States of America

16 Homework lecture 16

Integral Equation Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, Xin Wang and Karen Veroy

Fast Iterative Method-FIM: Application to the Convection- Diffusion Equations

Research Article Gauss-Lobatto Formulae and Extremal Problems

Singular Value Decomposition. Linear Algebra (3) Singular Value Decomposition. SVD and Eigenvectors. Solving LEs with SVD

LINEAR RECURRENT SEQUENCES AND POWERS OF A SQUARE MATRIX

13. Parametric and Non-Parametric Uncertainties, Radial Basis Functions and Neural Network Approximations

( ) 2 2. Multi-Layer Refraction Problem Rafael Espericueta, Bakersfield College, November, 2006

1 Convergence of the Arnoldi method for eigenvalue problems

Introduction to Numerical Differentiation and Interpolation March 10, !=1 1!=1 2!=2 3!=6 4!=24 5!= 120

X ε ) = 0, or equivalently, lim

The number of observed cases The number of parameters. ith case of the dichotomous dependent variable. the ith case of the jth parameter

Given a table of data poins of an unknown or complicated function f : we want to find a (simpler) function p s.t. px (

Maps on Triangular Matrix Algebras

13. Artificial Neural Networks for Function Approximation

Non-degenerate Perturbation Theory

Asymptotic Formulas Composite Numbers II

Chapter 4 Multiple Random Variables

A scalar t is an eigenvalue of A if and only if t satisfies the characteristic equation of A: det (A ti) =0

Lecture Note to Rice Chapter 8

Point Estimation: definition of estimators

L5 Polynomial / Spline Curves

Application of Legendre Bernstein basis transformations to degree elevation and degree reduction

Generalized One-Step Third Derivative Implicit Hybrid Block Method for the Direct Solution of Second Order Ordinary Differential Equation

u(x, t) = u 0 (x ct). This Riemann invariant u is constant along characteristics λ with x = x 0 +ct (u(x, t) = u 0 (x 0 )):

On the construction of symmetric nonnegative matrix with prescribed Ritz values

Generalized Linear Regression with Regularization

Taylor s Series and Interpolation. Interpolation & Curve-fitting. CIS Interpolation. Basic Scenario. Taylor Series interpolates at a specific

On the convergence of derivatives of Bernstein approximation

ANALYSIS ON THE NATURE OF THE BASIC EQUATIONS IN SYNERGETIC INTER-REPRESENTATION NETWORK

7.0 Equality Contraints: Lagrange Multipliers

ON THE ELEMENTARY SYMMETRIC FUNCTIONS OF A SUM OF MATRICES

( ) ( ) ( ( )) ( ) ( ) ( ) ( ) ( ) = ( ) ( ) + ( ) ( ) = ( ( )) ( ) + ( ( )) ( ) Review. Second Derivatives for f : y R. Let A be an m n matrix.

Exercises for Square-Congruence Modulo n ver 11

End of Finite Volume Methods Cartesian grids. Solution of the Navier-Stokes Equations. REVIEW Lecture 17: Higher order (interpolation) schemes

ECE606: Solid State Devices Lecture 13 Solutions of the Continuity Eqs. Analytical & Numerical

A GMRES Method for Solving Fuzzy Linear Equations

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

Log1 Contest Round 2 Theta Complex Numbers. 4 points each. 5 points each

Algorithms Theory, Solution for Assignment 2

3D Geometry for Computer Graphics. Lesson 2: PCA & SVD

Fractional Order Finite Difference Scheme For Soil Moisture Diffusion Equation And Its Applications

5 Short Proofs of Simplified Stirling s Approximation

A Remark on the Uniform Convergence of Some Sequences of Functions

Discrete Adomian Decomposition Method for. Solving Burger s-huxley Equation

A unified matrix representation for degree reduction of Bézier curves

ρ < 1 be five real numbers. The

Non-linear Matrix Equations: Equilibrium Analysis of Markov Chains

FREQUENCY ANALYSIS OF A DOUBLE-WALLED NANOTUBES SYSTEM

The Mathematical Appendix

Ideal multigrades with trigonometric coefficients

Transcription:

Numercal Aalyss Formulae Booklet. Iteratve Scemes for Systems of Lear Algebrac Equatos:.... Taylor Seres... 3. Fte Dfferece Approxmatos... 3 4. Egevalues ad Egevectors of Matrces.... 3 5. Vector ad Matrx Norms... 5 6. Ital Value Problem solvers... 6 7. Polyomal Iterpolato... 7 8. Secod Order Dfferece Equatos... 7 9. Euler s Formula ad oter Trgoometrc Idettes... 7 0. Tomas Algortm (TDMA)... 8. Newto Rapso Sceme... 9. Lax-Wedroff Sceme... 9 3. Crak-Ncolso sceme.....9 b.03.05 UFMFX9-30-3 Formula Booklet Page of 9

. Iteratve Scemes for Systems of Lear Algebrac Equatos:. Successve Over-Relaxato (SOR) terato sceme for solvg were s a matrx: ( ). Te Jacob terato matrx for solvg s gve by (F.) were, wt beg a dagoal matrx, ad are strctly lower ad strctly upper tragular matrces, respectvely.. Taylor Seres (a) Te Taylor seres expaso of a fucto f ( x ) about te pot 0 as te form f ( x ) f ( x) f ( x) f ( x) f ( x)! 3!! were x p x. 3 f ( ) ( x) f ( )! ( ) ( p) (F.) (b) Te Taylor seres expaso of a fucto f ( x, y k) about te pot defed by 0, k 0 as te form up to ad cludg secod order terms of f ( x, y k) f f x kf y! f xx kf xy k! f yy (F.) UFMFX9-30-3 Formula Booklet Page of 9

3. Fte Dfferece Approxmatos (a) Regular grd: u x u u cetral dfferece (F3.) u x u u forward dfferece (F3.) u u u u x cetral dfferece (F3.3) (b) Irregular grd: 0 x x u x x x u ( ) u u ( ) x (F3.4) u x u ( ) u u ( ) (F3.5) 4. Egevalues ad Egevectors of Matrces. (a) Egevalues of commo trdagoal system Te egevalues of a commo m m trdagoal matrx: a b c a b c a b c a b c a as m dstct egevalues gve by: s s a bc cos s,, m m (F4.) UFMFX9-30-3 Formula Booklet Page 3 of 9

(b) () Commo egevectors If matrx A ad B sare commo egevectors te tese egevectors are also egevectors k A k B for ay scalars k, k. of () If matrx A s vertble, te A ad A sare te same egevectors. (c) Gersgor s Crcle Teorem Gve a matrx A, te all ts egevalues le te rego D D were D z : z a R ad R j j j a. (F4.) (d) Oter propertes regardg egevalues of matrces () All te egevalues of a real symmetrc matrx are real. () If matrx A s a real trdagoal matrx ad all ts o-dagoal etres are of te same sg (eter all egatve or all postve), te all ts egevalues are real. UFMFX9-30-3 Formula Booklet Page 4 of 9

5. Vector ad Matrx Norms A vector orm of a vector deoted by, s a real-valued fucto of tat satsfes te followg propertes: s ay scalar (te dot deotg multplcato) te tragle equalty Some commo vector orms are gve below for a vector : u u u u m (F5.) u u u u m (F5.) max,,, m u u u u (F5.3) Gve a matrx ad gve a vector orm, te te duced orm of te matrx, deoted by, s defed as: Usg te defto of te duced matrx orm (F5.4), te followg matrx orms (of matrces wt real valued etres) result from te vector orms defed (F5.)-(F5.3), represetg te spectral radus (F5.5) UFMFX9-30-3 Formula Booklet Page 5 of 9

6. Ital Value Problem solvers For te frst order tal value problem y f ( t,, y(0) y 0 6. Euler s Metod y y f 6. Te Modfed Euler Metod y y f P y y f f P 6.3 Te Classcal d Order Ruge-Kutta Metod k f t, y k f t, y k y y k 6.4 Te Classcal 4 t Order Ruge-Kutta Metod k f t, y 4 3 k f t, y k k3 f t, y k k f t, y k y y k k k k 6 3 4 I te above expressos f f t, y ). ( UFMFX9-30-3 Formula Booklet Page 6 of 9

7. Polyomal Iterpolato Te polyomal of degree k passg troug k pots ( t, f), ( t, f ), ( t, f ),... s gve by: p ( t) f rf r( r ) f r( r )( r ) f... r( r )...( r k ) k f (F7.) 3 k! 3! ( k)! were t t r, f f f, k k f t, y( t) pk ( t) r( r )( r )...( r k ) f ( c), t c t. (F7.) k! ad 8. Secod Order Dfferece Equatos Te costat coeffcet secod order dfferece equato a y b y c y 0 were a, b ad c are costats, as a geeral soluto, provded p p, gve by were p ad p satsfy te auxlary equato y Ap Bp ap bp c 0 Teorem: Te roots of te quadratc equato are bouded by oe modulus f ad oly f all te followg codtos old: Q() 0, Q( ) 0, c a. 9. Euler s Formula ad oter Trgoometrc Idettes Euler s Formula: e cos s (F9.) cos s cos (F9.) UFMFX9-30-3 Formula Booklet Page 7 of 9

0. Tomas Algortm (TDMA) Cosder te followg trdagoal system of lear equatos, were te v s are te ukows: b c 0 0 0 v d a b c 0 0 v d 0 a3 b3. 0 v 3 d 3 0 0.. cm.. 0 0 0 am b m v m d m (F0.) Te TDMA states tat te soluto to ts trdagoal system s also a system to te followg bdagoal system of lear equatos: ˆ cˆ 0 0 0 v d 0 cˆ 0 0 v ˆ d 0 0. 0 v ˆ 3 d3 (F0.) 0 0. cˆ m.. 0 0 0 0 v m dˆ m were te forward step s gve by te followg: c cˆ, b c cˆ for,.., m b cˆ a ˆ, b d d d dˆ a d m ˆ for,.., b ˆ c a Te backward step, to obta te values of v, s to perform back substtuto o (F0.). UFMFX9-30-3 Formula Booklet Page 8 of 9

. Newto Rapso Sceme Cosder te followg system of m omogeeous olear equatos: f U, U,..., U 0,,..., m. (F.) m Te te Newto Rapso metod for solvg te system (F.) s a teratve sceme as follows: U U ( p) ( p) ( p) were te colum matrx A ( p ) f were A s te mm matrx gve by: A f ( p) s a soluto to te followg lear system of equatos: j, wt te etres evaluated at te p t terato stage U j ad te colum matrx f also as ts etres evaluated at te p t terato stage.. Lax-Wedroff Sceme Te Lax-Wedroff metod for te umercal soluto of equatos of te form: s gve by: were (F.) ( ) (F.) were, for example,. Wt ( ) ( ) (F.3) ( ) ( ) (F.4) 3. Crak Ncolso Sceme Te Crak-Ncolso sceme for te approxmato to te partal dfferetal equato gve by: s ( ) (( ) ( )) UFMFX9-30-3 Formula Booklet Page 9 of 9