Numerical Derivatives by Symbolic Tools in MATLAB
|
|
- Theresa Anderson
- 5 years ago
- Views:
Transcription
1 Numerical Derivatives by Symbolic ools i MALAB Mig-Gog Lee * ad Rei-Wei Sog ad Hsuag-Ci Cag mglee@cu.edu.tw * Departmet o Applied Matematics Cug Hua Uiversity Hsicu, aiwa Abstract: Numerical approaces o ordiary dieretial equatios (ODEs usually require Jacobia evaluatios ad sometimes iger order derivatives o system o equatios. e evaluatio o tese umerical derivatives ca eiter use algebraic derivatio or by te computer algebra system (CAS, suc as MAPLE V, but tese derivatios may eiter cause debuggig diiculty i programmig or memory swell i computatio. Fiite dierece ormulas icludig orward dierece, cetral dierece, ad backward dierece are used to derive te Jacobia wit dieret degree o accuracy. Symbolic evaluatio capability o MALAB is used to eace ast derivatio o umerical derivatives give i tis paper. e accuracy o tese umerical derivatives sows tat te evaluatio by symbolic computatio rom MALAB is possible to be equipped ito a umerical sceme o ordiary equatios.. Itroductio e calculatio o derivatives is ecessary i te umerical solutio o dieretial equatios, suc as te Jacobia matrix o Newto iteratios i most o te ordiary equatios solvers; ad iger order coeiciets i aylor series metods, etc. A variety o sotware icludig symbolic programs [5] ad automatic dieretiatio [,,4] ave bee itroduced ad used successully i te applicatio problems. I tis paper, we will use te ideas give i [7] ad to implemet i te MALAB eviromet [6]. We use te symbolic tool give i MALAB ad evaluate umerical derivatives o give uctios wit lexible order o accuracy. ese umerical derivatives ca urter be used i te solutios o ODEs by aylor series metods or oter umerical scemes. oug, Miletics ad Molárka [7] ave sow ow to use umerical derivatives i a sti ordiary equatio by aylor series metods, but we demostrated ow to implemet wit dieret combiatio to give accurate iger order derivatives. I orgaizatio o tis paper, Sectio gives some otatio o te dieretial equatios, Sectio gives some o te iger order umerical derivatives i matrix orms, Sectio 4 gives umerical results by implemet aylor series metod to a sti equatio. Coclusio is give at Sectio 5.. Matematical Formulatios Give a ordiary dieretial equatio y ( x = ( x, x, y ( x = y, x [ x, x + ], (. were x = [ x, x,..., y ( x] : R R. Give p suc tat uctios p + yi ( x C ([ x, x + ], or i =,...,, ad ( x, x = [ ( x, x, ( x, x,..., ( x, x]. e variable x may or may ot be explicitly expressed i te rigt-ad side uctio. I te case tat variable x is expressed explicitly i te Equatio., we ca deie a ew variable Y as ollows: LetY( x = [ x, y ( x, y ( x,..., y + ( x]: R R, te
2 dy ( x Y ( x = = [, ( Y( x,..., ( Y( x] = F( Y( x, (. dx werey ( x = x, Y ( x = y,..., Y( x = y. Assume i ( x, x, i =,...,, is a (p+ times dieretiable uctio, te te classical aylor series expressio [] at x + s about x ca be writte as, x + s = x + y ( x s + y ( x s + y ( x s !! p ( p p+ ( p+ s y ( x + s + s y ( x s, p! ( p +! ( p+ uctio y ( x s is evaluated at some ukow poit betwee x ad x s. (.. Numerical Derivative by Fiite Dierece We implemetig aylor series metod, some derivatives o rigt ad side uctio, i.e., Eq.. must be evaluated. I tis sectio, we will describe ow to use iite dierece to express iger order derivatives o Eq.... Fiite Dierece Expressios o Secod Order Derivatives From Eq.. or Eq.., a secod order derivative o y ca be expressed as: y ( x = ( x = [ ( Y( x,..., ( Y( x] = ( Y( x dy( x ( Y( x dy( x = [(,...,( ] = F. (.. Y dx Y dx Y Sice te irst compoet o uctio F is, te iger order derivatives o te irst compoet o tis uctio are all zero. As a result, we will oly use te otatio o Eq.. or simplicity. By implemetig orward dierecig, backward dierecig, ad cetral dierecig to approximate k ( x te partial derivatives,, or k =,,...,, various ormulas ca be obtaied. Matrix A is Y obtaied by orward dierece, matrix C is obtaied by cetral dierece, matrix D is obtaied by a extrapolatio ormula, ad tese partial derivatives are evaluated at x = x. A {[ ( x, y ( x,... y ( x,..., ( x, y ( x,... y ( x [ ( x, y ( x,... y ( x,..., ( x, y ( x,... y ( x ]],..., = [ + + ] [ [ ( x, y ( x,... y ( x,..., ( x, y ( x,... y ( x + + ] ( x, y ( x,... y ( x,..., ( x, y ( x,... y ( x ]} ] {[ [ (, (,... (,..., (, (,... ( C = x + y x y x x + y x y x ] [ ( x, y ( x,... y ( x,..., ( x, y ( x,... y ( x ],..., [ ( x, y ( x,... y ( x,..., ( x, y ( x,... y ( x [ + + [ ( x, y ( x,... y ( x,..., ( x, y ( x,... y ( x ] ]} ] (.. (..
3 D {[ ( x, y ( x,... y ( x,..., ( x, y ( x,... y ( x 6 6 ( x, y ( x,... y ( x,..., ( x, y ( x,... y ( x = [ ] [ ] + ( x, y ( x,... y ( x,..., ( x, y ( x,... y ( x + [ ] + [ ( x, y ( x,... y ( x,..., ( x, y ( x,... y ( x ]],..., [[ ( x, y ( x,... y ( x,..., ( x, y ( x,... y ( x ] 6 (, (,... (,..., (, (,... ( [ xy x y x x y x y x ] + + [ ( x, y ( x,... y ( x,..., ( x, y ( x,... y ( x ] + (..4 + [ ( x, x,... x +,..., ( x, x,... x + ]]} Because o te complex ormulatio o te expressios give at Equatios..-..4, we use otatios give at [7] to simply, let ki, F ( x, x = [ ( x, y ( x,..., y ( x + i,...,..., ( x, y ( x,..., y ( x + i,...], (..5 k k k =,...,, i =, ±, ±,..., ±, we k=, y ( x = x. By tis otatio, we ca rewrite Eqs...~..4 as ollows:,,,,,, A = [( F F, ( F F,..., ( F F ], (..6 [(,,, (,,,..., (,, C = F F F F F F ], (..7 [(, 6,,, D = F F + F + F,..., 6 (..8,,, ( F 6F + F, + F ]. e ollowig matrix E is obtaied by a special type o iite dierece, ad is give as, [(, 8,,, E = F F + 8 F F,..., (..9,, ( F 8F,, + 8 F F ].. ird Order Derivatives As i Sectio., we ca derive tird order derivative o variable y or te secod derivative o Eq.. as: Y = F = ( F = ( FF + ( F. (.. Y Y Y Y Similarly we also deie a ew uctio so as to simpliy otatios o te derived derivative uctios: kl, i F = (..., yk + i,..., yl +,..., kl, =,,...,, i, =, ±, ±..., (.. kk, ii F = (..., yk + i,..., kl, =,,...,,, i =, ±, ±.... Higer order derivatives were obtaied similarly as te secod order derivatives. e ollowig matrix G is ormulated by irst applyig te cetral dierece sceme to ave a secod order m
4 derivative uctio, ad te applyig a orward dierece sceme to ave a tird order derivative uctio. e iduced ormula is give as:,,,, [( m m m m Gm = F F F + F,..., mm, mm, mm, ( F F + F,..., (.. m, m, m, m, ( F F F + F ], m =,,...,. Aoter tird order derivative ormula H m ca be obtaied by implemetig cetral dierece approximatio cosecutively as:,,,, [( m m m m H m = F F F + F,..., 4 mm, mm, mm, 4( F F + F,..., (..4 m, ( F m, F m, F m, + F ], m =,,..., By usig ormulas G m or H m, te matrix o Y ca be approximated as: [ G, G,..., G ] [ H, H,..., H ] Y. (..5. e ourt Order Derivative By Eq..., te ourt derivative o uctio y or te tird order derivative matrix o Eq.. ca be writte as: IV Y = F = ( FF + ( ( F = FFF + ( F F + Y Y Y Y Y Y (.. +( F( F + ( FF + ( F Y Y Y Y Y I Eq..., is derived similarly as i sectio. by lower order derivative uctios, ad its Y ormula ca be writte as, i,, k, l =,,...,. YiY Yk Similarly, we also deie ew uctios to simpliy te otatio [7], ( Y = (..., y + m,... y + e,..., y + z,..., il, mez k k i l iil, mmz k ( Y = k(..., yi + m,..., yl + z,..., (.. iii, mmm ( Y = (..., y + m,..., k, l, i, =,,...,, me,, z=, ±, ±... k k i We te partial derivatives were implemeted wit tree dieret directios, a ormula is give as:
5 ,,,, [( i i i i i = F F F + F i, i, i, F + F F ],..., (..,, [( i i F F i, i, i, i, i, F + F F + F F ], i, =,,... I te partial derivatives were implemeted wit respect to two same directios, a ormula is give as: i, i, i, i, i, i, ( i = [F F F F + F + F ], i, =,,...,. (..4 I te partial derivative were implemeted wit respect to oly oe directio, a ormula is give as:,,,, ( = [ F + F + F F ], =,,...,. (..5 As or tese derivatives matrices or, it ca be obtaied by combig dieret ormulas, suc Y as orward ad cetral dierece ad te our-poit ormula E. We omit te procedures ere. 4. Numerical Derivatives ad Its Accuracy I tis sectio we will implemet umerical derivatives by iite dierece give i Sec. i a aylor series metod to a sti ordiary equatio. I Sectio 4., we will sow ow to acquire tese umerical derivatives i a lexible way, i.e., we ca calculate wit dieret combiatio rom te cetral dierecig, orward dierecig, ad oter ormulas to obtai more accurate derivatives. Some errors betwee exact derivatives ad umerical derivatives are give to sow teir accuracy. I Sectio 4., we will compare te umerical solutios obtaied by implemetig umerical derivative ad by exact derivatives to a sti ordiary dieretial equatio. 4. Accuracy o te Numerical Derivatives e secod order derivatives was give at Equatio.., te umerical derivative ca be obtaied by applyig Eqs...6~..9. For te tird order umerical derivative o Equatio., it will be give as ollows. Durig te implemetatio, a cetral dierece is irst used to get te secod order derivative, ad a orward dierece sceme is applied urter to get te tird order derivatives. e order o applyig iite dierece scemes is ot te critical issue, sice we ca still cage reely wit tese orderig to get accurate umerical derivatives. e above procedure was give oly to demostrate oe possibility o our procedure. As a result, or te above procedure, te tird order umerical derivative ca be writte as: m, m, m, m, G m = [( F F F + F,..., mm, mm, mm, mm, F F F + F,..., (4.. m, m, m, m, ( F F F F ] +, m =,,,.
6 ere is a mior dierece betwee G m ad G m wic is give at Eq.... Matrix G m is used i te implemetatio o tis paper. Similarly, i te ourt derivative o Equatio., we use cetral dierecig sceme ad backward dierecig sceme or easy implemetatio i programmig, a ormula is give we tree directios o partial dieretiatio are all dieret ad is give as, i, i, i, i, i = [( F F F + F i, i, i, i, i, i, F + F + F F ],...,( F F (4.. i, i, i, i, i, i, F + F F + F + F F ], i, =,,.... We tere are two same directios o partial dieretiatio ad a ormula is give as, i, i, i, i, ( i = [ F F F + F (4.. i, i, i, i, F + F + F F ]. We te directios o partial dieretiatio are all te same, a ormula is give as,,,,, ( = [ F F + F F ]. (4..4 Equatios 4..~4..4 are te ormulas or sow i te ourt derivatives. ese Y ormulas ca be implemeted wit te secod order ormulas A, B, C, D ad E, i tis paper we oly adapt matrices A ad C or coveiece. A sti dieretial equatio [6] is give to test te accuracy o te umerical derivatives o dieret orders. We will compare te accuracy o umerical derivatives o rigt ad side uctio wit te exact derivatives. d x = ( y ( x + y ( x y ( x y ( x qy ( x dt α dy ( x = m x x x x (4..5 dt d x = ( y ( x y ( x, = a, = b, = d, dt γ were α, m, q, ad γ are parameters, a, b ad d are iitial values. Letα =., m=.5, γ =, q=., a=5, b=7, ad d=8. Evaluate Eqs...6~..9, te ollowig error estimates o te irst derivative are estimated ad are give at able. Matrices A, C, D, ad E are all ave good accuracy. i able. We =.5, error o matrices A, C, D, E max( A Y max( C Y max( D Y max( E Y
7 Evaluate matrices o Eqs.4.., te ollowig error estimates o te secod order derivative are estimated ad are give at able. It sows tat matrix G as good accuracy. max ( ( G k k, YYk able. We =.5, error i calculatig G max ( ( G k k, YYk max ( ( G k k, YYk Evaluate matrices o Eqs. 4..~4..4, te ollowig error estimates o te tird order derivatives are estimated ad are give at able.. It sows tat matrix as good accuracy. max ( ( i k i,, k YiYkY able. We =.5, error i calculatig max ( ( i k i,, k YiYkY max ( ( i k i,, k YiYkY ese results sow tat te accuracy o te umerical derivatives all satisies te teoretical results give i [7]. 4. e Numerical Solutio by aylor Series Metod Figure. Numerical Solutios o Equatio 4..5 give by aylor Series Metod wit Numerical Derivatives
8 Followig te same algoritm give at [7], let te error cotrolε =. ad te iterval o itegratio is [, ]. e umerical solutios o aylor series metod by implemetig umerical derivatives give i Sec. 4 are sow i Figure above. e traces are solutios o tree compoets o te system o Equatio able 4. Numerical Values by aylor Series Metod wit ree ways to Obtai Derivatives Subiterval aylor Series Metod umerical values MALAB MAPLE Exact Derivative [.5767, ] [ ,.59885] [.59885, ] ŷ ŷ ŷ ŷ ŷ ŷ ŷ ŷ ŷ Durig te itegratio, we also examie te umerical solutios obtaied by applyig umerical derivatives ad exact derivatives i several subitervals to sow tat tese umerical derivatives do beave similarly wit te solutios obtaied by exact derivatives. ese umerical values eve exact at least up to te ourt digits ad are give at able 4. I te ollowig, we examie te umerical solutios by applyig aylor series metod ad by subroutie give i MALAB, te error betwee tese solutios matc well ad are give at able 5. able 5. Numerical Values by aylor Series Metod wit ree ways to Obtai Derivatives (Compared wit solutios give by implemetig Matlab Solutios Error betwee Subiterval MALAB ad aylor Series Metods [.5767, ] [ ,.59885] [.59885, ] x + s yˆ ( x + s k y 8.744e-7 y.54568e-7 y e-9 y e-7 y e-7 y e-8 y e-6 y e-6 y e-8 k
9 From ables 4 ad 5, at several subitervals o itegratio, umerical solutios obtaied by umerical derivative give i tis paper ad te exact derivative bot matc well. It meas tat te umerical derivatives derived rom applyig iite dierece scemes ca be applied durig umerical solutios o dieretial equatio eve it is sti. 5. Coclusio I tis paper, we ave demostrated ow to calculate umerical derivatives wit iite dierece sceme ad use te umeric to a aylor series metod to solve a sti ordiary equatio. ese umerical derivatives also ave good accuracy compared wit exact derivatives. oug alterative way o idig derivatives are also possible, by results rom tis paper, umerical derivatives by iite dierecig is still capable o rederig derivatives wit good accuracy. aylor series metod may ot be suitable or ast computatioal aspect, it requires may iger orders derivatives, but a possibility o parallel computatio may be a resolutio or te speed cocer. I additio, some results sow tat te aylor series metod is ot suitable or very sti problems, sice it is basically a explicit metod. I order to solve tis stability requiremet, a implicit type o aylor series metod may be approaced i te uture. Reereces [] Biso Cristia, Carle Ala, Havlad Paul, Kademi Peyvad, ad Mauer Adrew, ADIFOR. User s Guide (Revisio D, ecical Report ANL/MCS-M-9, Matematics ad Computer Sciece Divisio, Argoe Natioal Laboratory, USA, (994. [] Griewak Adreas, Evaluatig Derivatives-Priciples ad eciques o Algoritmic Dieretiatio, SIAM, Piladelpia, (. [] Kedall E. Atkiso, Itroductio to Numerical Aalysis, Jo Wiley & Sos Ltd., New York (978. [4] Lee, M.G., Applicatio o Automatic Dieretiatio i Numerical Solutio o a Flexible Mecaism, ICCMSE Iteratioal Coerece o Computatioal Metods i Scieces ad Egieerig, Kastoria, Greece, pp5-59. [5] MAPLE V, User Maual, e Waterloo Maple Ic., Caada (997 [6] MALAB, User Maual, e Matworks Ic., USA (6. [7] Miletics, E. ad Molárka G., aylor Series Metod wit Numerical Derivatives or Iitial Value Problems, Jouaal o Computatioal Metods i Scieces ad Egieerig,vol.,o. (, pp. 9-9.
Solving third order boundary value problem with fifth order block method
Matematical Metods i Egieerig ad Ecoomics Solvig tird order boudary value problem wit it order bloc metod A. S. Abdulla, Z. A. Majid, ad N. Seu Abstract We develop a it order two poit bloc metod or te
More informationAn Improved Self-Starting Implicit Hybrid Method
IOSR Joural o Matematics (IOSR-JM e-issn: 78-78, p-issn:9-76x. Volume 0, Issue Ver. II (Mar-Apr. 04, PP 8-6 www.iosrourals.org A Improved Sel-Startig Implicit Hbrid Metod E. O. Adeea Departmet o Matematics/Statistics,
More informationd y f f dy Numerical Solution of Ordinary Differential Equations Consider the 1 st order ordinary differential equation (ODE) . dx
umerical Solutio o Ordiar Dieretial Equatios Cosider te st order ordiar dieretial equatio ODE d. d Te iitial coditio ca be tae as. Te we could use a Talor series about ad obtai te complete solutio or......!!!
More information1. Introduction. 2. Numerical Methods
America Joural o Computatioal ad Applied Matematics, (5: 9- DOI:.59/j.ajcam.5. A Stud o Numerical Solutios o Secod Order Iitial Value Problems (IVP or Ordiar Dieretial Equatios wit Fourt Order ad Butcer
More informationME 501A Seminar in Engineering Analysis Page 1
Accurac, Stabilit ad Sstems of Equatios November 0, 07 Numerical Solutios of Ordiar Differetial Equatios Accurac, Stabilit ad Sstems of Equatios Larr Caretto Mecaical Egieerig 0AB Semiar i Egieerig Aalsis
More informationNUMERICAL DIFFERENTIAL 1
NUMERICAL DIFFERENTIAL Ruge-Kutta Metods Ruge-Kutta metods are ver popular ecause o teir good eiciec; ad are used i most computer programs or dieretial equatios. Te are sigle-step metods as te Euler metods.
More informationSolving Third Order Boundary Value Problem Using. Fourth Order Block Method
Applied Matematical Scieces, Vol. 7,, o. 5, 69-65 HIKARI Ltd, www.m-ikari.com Solvig Tird Order Boudar Value Problem Usig Fourt Order Block Metod Amad Sa Abdulla, *Zaaria Abdul Maid, ad Norazak Seu, Istitute
More informationThe Advection-Diffusion equation!
ttp://www.d.edu/~gtryggva/cf-course/! Te Advectio-iffusio equatio! Grétar Tryggvaso! Sprig 3! Navier-Stokes equatios! Summary! u t + u u x + v u y = P ρ x + µ u + u ρ y Hyperbolic part! u x + v y = Elliptic
More informationMore Elementary Aspects of Numerical Solutions of PDEs!
ttp://www.d.edu/~gtryggva/cfd-course/ Outlie More Elemetary Aspects o Numerical Solutios o PDEs I tis lecture we cotiue to examie te elemetary aspects o umerical solutios o partial dieretial equatios.
More informationA Pseudo Spline Methods for Solving an Initial Value Problem of Ordinary Differential Equation
Joural of Matematics ad Statistics 4 (: 7-, 008 ISSN 549-3644 008 Sciece Publicatios A Pseudo Splie Metods for Solvig a Iitial Value Problem of Ordiary Differetial Equatio B.S. Ogudare ad G.E. Okeca Departmet
More informationA NUMERICAL METHOD OF SOLVING CAUCHY PROBLEM FOR DIFFERENTIAL EQUATIONS BASED ON A LINEAR APPROXIMATION
U.P.B. Sci. Bull., Series A, Vol. 79, Iss. 4, 7 ISSN -77 A NUMERICAL METHOD OF SOLVING CAUCHY PROBLEM FOR DIFFERENTIAL EQUATIONS BASED ON A LINEAR APPROXIMATION Cristia ŞERBĂNESCU, Marius BREBENEL A alterate
More informationx x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula
NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS. If g( is cotiuous i [a,b], te uder wat coditio te iterative (or iteratio metod = g( as a uique solutio i [a,b]? '( i [a,b].. Wat
More informationALLOCATING SAMPLE TO STRATA PROPORTIONAL TO AGGREGATE MEASURE OF SIZE WITH BOTH UPPER AND LOWER BOUNDS ON THE NUMBER OF UNITS IN EACH STRATUM
ALLOCATING SAPLE TO STRATA PROPORTIONAL TO AGGREGATE EASURE OF SIZE WIT BOT UPPER AND LOWER BOUNDS ON TE NUBER OF UNITS IN EAC STRATU Lawrece R. Erst ad Cristoper J. Guciardo Erst_L@bls.gov, Guciardo_C@bls.gov
More informationStability analysis of numerical methods for stochastic systems with additive noise
Stability aalysis of umerical metods for stoctic systems wit additive oise Yosiiro SAITO Abstract Stoctic differetial equatios (SDEs) represet pysical peomea domiated by stoctic processes As for determiistic
More informationSome Variants of Newton's Method with Fifth-Order and Fourth-Order Convergence for Solving Nonlinear Equations
Copyright, Darbose Iteratioal Joural o Applied Mathematics ad Computatio Volume (), pp -6, 9 http//: ijamc.darbose.com Some Variats o Newto's Method with Fith-Order ad Fourth-Order Covergece or Solvig
More informationNumerical Integration Formulas
Numerical Itegratio Formulas Berli Che Departmet o Computer Sciece & Iormatio Egieerig Natioal Taiwa Normal Uiversity Reerece: 1. Applied Numerical Methods with MATLAB or Egieers, Chapter 19 & Teachig
More informationApplication of a Two-Step Third-Derivative Block Method for Starting Numerov Method
Iteratioal Joural of eoretical ad Applied Matematics 7; (: -5 ttp://wwwsciecepublisiroupcom//itam doi: 648/itam75 Applicatio of a wo-step ird-derivative Block Metod for Starti Numerov Metod Oluwaseu Adeyeye
More informationCHAPTER 6c. NUMERICAL INTERPOLATION
CHAPTER 6c. NUMERICAL INTERPOLATION A. J. Clark School o Egieerig Departmet o Civil ad Evirometal Egieerig y Dr. Irahim A. Assakka Sprig ENCE - Computatio Methods i Civil Egieerig II Departmet o Civil
More informationOn Exact Finite-Difference Scheme for Numerical Solution of Initial Value Problems in Ordinary Differential Equations.
O Exact Fiite-Differece Sceme for Numerical Solutio of Iitial Value Problems i Ordiar Differetial Equatios. Josua Suda, M.Sc. Departmet of Matematical Scieces, Adamawa State Uiversit, Mubi, Nigeria. E-mail:
More informationAccuracy. Computational Fluid Dynamics. Computational Fluid Dynamics. Computational Fluid Dynamics
http://www.d.edu/~gtryggva/cfd-course/ Computatioal Fluid Dyamics Lecture Jauary 3, 7 Grétar Tryggvaso It is clear that although the umerical solutio is qualitatively similar to the aalytical solutio,
More informationCS537. Numerical Analysis and Computing
CS57 Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 Jauary 9 9 What is the Root May physical system ca be
More informationPhys. 201 Mathematical Physics 1 Dr. Nidal M. Ershaidat Doc. 12
Physics Departmet, Yarmouk Uiversity, Irbid Jorda Phys. Mathematical Physics Dr. Nidal M. Ershaidat Doc. Fourier Series Deiitio A Fourier series is a expasio o a periodic uctio (x) i terms o a iiite sum
More informationNumerical Methods for Ordinary Differential Equations
Numerical Methods for Ordiary Differetial Equatios Braislav K. Nikolić Departmet of Physics ad Astroomy, Uiversity of Delaware, U.S.A. PHYS 460/660: Computatioal Methods of Physics http://www.physics.udel.edu/~bikolic/teachig/phys660/phys660.html
More information5. Fast NLMS-OCF Algorithm
5. Fast LMS-OCF Algorithm The LMS-OCF algorithm preseted i Chapter, which relies o Gram-Schmidt orthogoalizatio, has a compleity O ( M ). The square-law depedece o computatioal requiremets o the umber
More informationCS321. Numerical Analysis and Computing
CS Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 September 8 5 What is the Root May physical system ca
More informationLinear Complementarity Problem Simulation of The Meter Stick Trick
Liear Complemetarity Problem Simulatio o Te Meter Stick Trick Yoke Peg Leog, Ze Lv Abstract I tis paper, we ocus o aturally icorporatig cotact modes estimatio i a dyamical simulatio o a object wit multiple
More informationMathematical Expressions for Estimation of Errors in the Formulas which are used to obtaining intermediate values of Biological Activity in QSAR
Iteratioal Joural o Iovatio ad Applied Studies ISSN 08-934 Vol. No. 3 Mar. 03, pp. 7-79 03 Iovative Space o Scietiic Researc Jourals ttp//www.issr-jourals.org/ijias/ Matematical Expressios or Estimatio
More informationTHE ENERGY BALANCE ERROR FOR CIRCUIT TRANSIENT ANALYSIS
THE ENERGY BALANCE ERROR FOR CIRCUIT TRANSIENT ANALYSIS FLORIN CONSTANTINESCU, ALEXANDRU GABRIEL GHEORGHE, MIRUNA NIŢESCU Key words: Trasiet aalysis, Eergy balace error, Time step coice. Two algoritms
More informationComputation Sessional. Numerical Differentiation and Integration
CE 6: Egieerig Computatio Sessioal Numerical Dieretiatio ad Itegratio ti di ad gradiet commad di() Returs te dierece betwee adjacet elemets i. Typically used or uequally spaced itervals = gradiet(, ) Determies
More informationWhere do eigenvalues/eigenvectors/eigenfunctions come from, and why are they important anyway?
Where do eigevalues/eigevectors/eigeuctios come rom, ad why are they importat ayway? I. Bacgroud (rom Ordiary Dieretial Equatios} Cosider the simplest example o a harmoic oscillator (thi o a vibratig strig)
More informationComputation of Hahn Moments for Large Size Images
Joural of Computer Sciece 6 (9): 37-4, ISSN 549-3636 Sciece Publicatios Computatio of Ha Momets for Large Size Images A. Vekataramaa ad P. Aat Raj Departmet of Electroics ad Commuicatio Egieerig, Quli
More informationModification of Weerakoon-Fernando s Method with Fourth-Order of Convergence for Solving Nonlinear Equation
ISSN: 50-08 Iteratioal Joural o AdvacedResearch i Sciece, Egieerig ad Techology Vol. 5, Issue 8, August 018 Modiicatio o Weerakoo-Ferado s Method with Fourth-Order o Covergece or Solvig Noliear Equatio
More informationDetermination of Energy Involved In a Stepwise Size Reduction of Maize, Using a Numerical Approach.
IOSR Joural of Matematics (IOSR-JM) e-issn: 78-578, p-issn: 9-765X. Volume, Issue Ver. V (Ja - Feb. 5), PP 6-65 www.iosrourals.org Determiatio of ergy Ivolved I a Stepwise Size Reductio of Maize, Usig
More informationCHAPTER 6d. NUMERICAL INTERPOLATION
CHAPER 6d. NUMERICAL INERPOLAION A. J. Clark School o Egieerig Departmet o Civil ad Evirometal Egieerig by Dr. Ibrahim A. Assakka Sprig ENCE - Computatio Methods i Civil Egieerig II Departmet o Civil ad
More informationTaylor Polynomials and Approximations - Classwork
Taylor Polyomials ad Approimatios - Classwork Suppose you were asked to id si 37 o. You have o calculator other tha oe that ca do simple additio, subtractio, multiplicatio, or divisio. Fareched\ Not really.
More informationThis chapter describes different methods to discretize the diffusion equation. f z 2 = 0. y ) x f
Chapter 8 Diusio Equatio This chapter describes dieret methods to discretize the diusio equatio 2 t α x 2 + 2 y 2 + 2 z 2 = 0 which represets a combied boudary ad iitial value problem, i.e., requires to
More informationAn Extension of the Szász-Mirakjan Operators
A. Şt. Uiv. Ovidius Costaţa Vol. 7(), 009, 37 44 A Extesio o the Szász-Mirakja Operators C. MORTICI Abstract The paper is devoted to deiig a ew class o liear ad positive operators depedig o a certai uctio
More informationENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 4 Solutions [Numerical Methods]
ENGI 3 Advaced Calculus or Egieerig Facult o Egieerig ad Applied Sciece Problem Set Solutios [Numerical Methods]. Use Simpso s rule with our itervals to estimate I si d a, b, h a si si.889 si 3 si.889
More informationg () n = g () n () f, f n = f () n () x ( n =1,2,3, ) j 1 + j 2 + +nj n = n +2j j n = r & j 1 j 1, j 2, j 3, j 4 = ( 4, 0, 0, 0) f 4 f 3 3!
Higher Derivative o Compositio. Formulas o Higher Derivative o Compositio.. Faà di Bruo's Formula About the ormula o the higher derivative o compositio, the oe by a mathematicia Faà di Bruo i Italy o about
More informationA Self-Starting Hybrid Linear Multistep Method for a Direct Solution of the General Second-Order Initial Value Problem
IOS Joural o Matematics (IOS-JM) ISSN: 78-578. Volume 4 Issue 6 (Ja. - eb. ) PP 7- www.iosrourals.org Sel-Startig Hbrid Liear Multistep Metod or a Direct Solutio o te Geeral Secod-Order Iitial Value Problem
More informationA NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS
Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet
More informationNumerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION
alaysia Joural of athematical Scieces 3(1): 83-93 (9) umerical Coformal appig via a Fredholm Itegral Equatio usig Fourier ethod 1 Ali Hassa ohamed urid ad Teh Yua Yig 1, Departmet of athematics, Faculty
More informationCelestin Chameni Nembua University of YaoundéII, Cameroon. Abstract
A ote o te decompositio of te coefficiet of variatio squared: comparig etropy ad Dagum's metods Celesti Camei Nembua Uiversity of YaoudéII, Cameroo Abstract Te aim of tis paper is to propose a ew decompositio
More informationNumerical Solutions of Second Order Boundary Value Problems by Galerkin Residual Method on Using Legendre Polynomials
IOSR Joural of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 11, Issue 6 Ver. IV (Nov. - Dec. 15), PP 1-11 www.iosrjourals.org Numerical Solutios of Secod Order Boudary Value Problems
More informationOn the convergence, consistence and stability of a standard finite difference scheme
AMERICAN JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH 2, Sciece Huβ, ttp://www.sciub.org/ajsir ISSN: 253-649X, doi:.525/ajsir.2.2.2.74.78 O te covergece, cosistece ad stabilit of a stadard fiite differece
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationComputational Methods CMSC/AMSC/MAPL 460. Quadrature: Integration
Computatioal Metods CMSC/AMSC/MAPL 6 Quadrature: Itegratio Ramai Duraiswami, Dept. o Computer Siee Some material adapted rom te olie slides o Eri Sadt ad Diae O Leary Numerial Itegratio Idea is to do itegral
More informationOn One Justification on the Use of Hybrids for the Solution of First Order Initial Value Problems of Ordinary Differential Equations
Pure and Applied Matematics Journal 7; 6(5: 74 ttp://wwwsciencepublisinggroupcom/j/pamj doi: 648/jpamj765 ISSN: 6979 (Print; ISSN: 698 (Online On One Justiication on te Use o Hybrids or te Solution o First
More informationCaputo s Implicit Solution of Time-Fractional Diffusion Equation Using Half-Sweep AOR Iteration
Global Joural o Pure ad Applied Mathematics. ISSN 0973-768 Volume Number 4 (06) pp. 3469-3479 Research Idia Publicatios http://www.ripublicatio.com/gpam.htm Caputo s Implicit Solutio o Time-Fractioal Diusio
More informationMath 257: Finite difference methods
Math 257: Fiite differece methods 1 Fiite Differeces Remember the defiitio of a derivative f f(x + ) f(x) (x) = lim 0 Also recall Taylor s formula: (1) f(x + ) = f(x) + f (x) + 2 f (x) + 3 f (3) (x) +...
More informationA NOTE ON THE TOTAL LEAST SQUARES FIT TO COPLANAR POINTS
A NOTE ON THE TOTAL LEAST SQUARES FIT TO COPLANAR POINTS STEVEN L. LEE Abstract. The Total Least Squares (TLS) fit to the poits (x,y ), =1,,, miimizes the sum of the squares of the perpedicular distaces
More informationFIR Filter Design: Part I
EEL3: Discrete-Time Sigals ad Systems FIR Filter Desig: Part I. Itroductio FIR Filter Desig: Part I I this set o otes, we cotiue our exploratio o the requecy respose o FIR ilters. First, we cosider some
More informationSolution of Differential Equation from the Transform Technique
Iteratioal Joural of Computatioal Sciece ad Mathematics ISSN 0974-3189 Volume 3, Number 1 (2011), pp 121-125 Iteratioal Research Publicatio House http://wwwirphousecom Solutio of Differetial Equatio from
More informationANSWER KEY WITH SOLUTION PAPER - 2 MATHEMATICS SECTION A 1. B 2. B 3. D 4. C 5. B 6. C 7. C 8. B 9. B 10. D 11. C 12. C 13. A 14. B 15.
TARGET IIT-JEE t [ACCELERATION] V0 to V BATCH ADVANCED TEST DATE : - 09-06 ANSWER KEY WITH SOLUTION PAPER - MATHEMATICS SECTION A. B. B. D. C 5. B 6. C 7. C 8. B 9. B 0. D. C. C. A. B 5. C 6. D 7. A 8.
More informationHyperbolic Heat Equation in Bar and Finite Difference Schemes of Exact Spectrum
atest reds o heoretical ad Applied Mechaics Fluid Mechaics ad Heat & Mass raser Hyperbolic Heat Equatio i Bar ad Fiite Dierece Schemes o Exact Spectrum A BUIKIS H KAIS Abstract he solutios o correspodig
More informationPartial Differential Equations
EE 84 Matematical Metods i Egieerig Partial Differetial Eqatios Followig are some classical partial differetial eqatios were is assmed to be a fctio of two or more variables t (time) ad y (spatial coordiates).
More informationTopic 9 - Taylor and MacLaurin Series
Topic 9 - Taylor ad MacLauri Series A. Taylors Theorem. The use o power series is very commo i uctioal aalysis i act may useul ad commoly used uctios ca be writte as a power series ad this remarkable result
More informationUsing An Accelerating Method With The Trapezoidal And Mid-Point Rules To Evaluate The Double Integrals With Continuous Integrands Numerically
ISSN -50 (Paper) ISSN 5-05 (Olie) Vol.7, No., 017 Usig A Acceleratig Method With The Trapezoidal Ad Mid-Poit Rules To Evaluate The Double Itegrals With Cotiuous Itegrads Numerically Azal Taha Abdul Wahab
More informationA New Hybrid in the Nonlinear Part of Adomian Decomposition Method for Initial Value Problem of Ordinary Differential Equation
Joural of Matematics Researc; Vol No ; ISSN - E-ISSN - Publised b Caadia Ceter of Sciece ad Educatio A New Hbrid i te Noliear Part of Adomia Decompositio Metod for Iitial Value Problem of Ordiar Differetial
More informationFinite Difference Method for the Estimation of a Heat Source Dependent on Time Variable ABSTRACT
Malaysia Joural of Matematical Scieces 6(S): 39-5 () Special Editio of Iteratioal Worsop o Matematical Aalysis (IWOMA) Fiite Differece Metod for te Estimatio of a Heat Source Depedet o Time Variable, Allabere
More informationNUMERICAL SOLUTIONS OF THE FRACTIONAL KdV-BURGERS-KURAMOTO EQUATION
S5 NUMERICAL SOLUTIONS OF THE FRACTIONAL KdV-BURGERS-KURAMOTO EQUATION by Doga KAYA a*, Sema GULBAHAR a, ad Asif YOKUS b a Departmet of Matematics, Istabul Commerce Uiversity, Uskudar, Istabul, Turkey
More informationFurther Methods for Advanced Mathematics (FP2) WEDNESDAY 9 JANUARY 2008
ADVANCED GCE 7/ MATHEMATICS (MEI) Furter Metods for Advaced Matematics (F) WEDNESDAY 9 JANUARY 8 Additioal materials: Aswer Booklet (8 pages) Grap paper MEI Eamiatio Formulae ad Tables (MF) Afteroo Time:
More informationPrinciple Of Superposition
ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give
More informationFuzzy n-normed Space and Fuzzy n-inner Product Space
Global Joural o Pure ad Applied Matheatics. ISSN 0973-768 Volue 3, Nuber 9 (07), pp. 4795-48 Research Idia Publicatios http://www.ripublicatio.co Fuzzy -Nored Space ad Fuzzy -Ier Product Space Mashadi
More informationLecture 21. Numerical differentiation. f ( x+h) f ( x) h h
Lecture Numerical differentiation Introduction We can analytically calculate te derivative of any elementary function, so tere migt seem to be no motivation for calculating derivatives numerically. However
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More informationTMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods
TMA4205 Numerical Liear Algebra The Poisso problem i R 2 : diagoalizatio methods September 3, 2007 c Eiar M Røquist Departmet of Mathematical Scieces NTNU, N-749 Trodheim, Norway All rights reserved A
More informationComputer Derivations of Numerical Differentiation Formulae. Int. J. of Math. Education in Sci. and Tech., V 34, No 2 (March-April 2003), pp
Computer Derivations o Numerical Dierentiation Formulae By Jon H. Matews Department o Matematics Caliornia State University Fullerton USA Int. J. o Mat. Education in Sci. and Tec. V No (Marc-April ) pp.8-87.
More informationThe Discrete-Time Fourier Transform (DTFT)
EEL: Discrete-Time Sigals ad Systems The Discrete-Time Fourier Trasorm (DTFT) The Discrete-Time Fourier Trasorm (DTFT). Itroductio I these otes, we itroduce the discrete-time Fourier trasorm (DTFT) ad
More informationGeneralized H-codes and type II codes over GF(4)
THE JOURNA OF CHINA UNIVERSITIES OF POSTS AND TEECOMMUNICATIONS Volume 15, Issue 2, Jue 28 IN i-qi, WEN ia-mi, ZHENG Wei Geeralized H-codes ad type II codes over GF(4) CC umber O1574 Documet A Article
More informationAN EXTENSION OF A RESULT ABOUT THE ORDER OF CONVERGENCE
Bulleti o Mathematical Aalysis ad Applicatios ISSN: 8-9, URL: http://www.bmathaa.or Volume 3 Issue 3), Paes 5-34. AN EXTENSION OF A RESULT ABOUT THE ORDER OF CONVERGENCE COMMUNICATED BY HAJRUDIN FEJZIC)
More informationNumerical Solution of the Two Point Boundary Value Problems By Using Wavelet Bases of Hermite Cubic Spline Wavelets
Australia Joural of Basic ad Applied Scieces, 5(): 98-5, ISSN 99-878 Numerical Solutio of the Two Poit Boudary Value Problems By Usig Wavelet Bases of Hermite Cubic Splie Wavelets Mehdi Yousefi, Hesam-Aldie
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationAtomic Physics 4. Name: Date: 1. The de Broglie wavelength associated with a car moving with a speed of 20 m s 1 is of the order of. A m.
Name: Date: Atomic Pysics 4 1. Te de Broglie wavelegt associated wit a car movig wit a speed of 0 m s 1 is of te order of A. 10 38 m. B. 10 4 m. C. 10 4 m. D. 10 38 m.. Te diagram below sows tree eergy
More informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
More information1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations
. Liearizatio of a oliear system give i the form of a system of ordiary differetial equatios We ow show how to determie a liear model which approximates the behavior of a time-ivariat oliear system i a
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationPARALEL PREDICTOR-CORRECTOR SCHEMES FOR PARABOLIC PROBLEMS ON GRAPHS
COMPUTATIONAL METHODS IN APPLIED MATHEMATICS, Vol. 10(2010, No. 3, pp. 275 282 c 2010 Istitute of Matematics of te Natioal Academy of Scieces of Belarus PARALEL PREDICTOR-CORRECTOR SCHEMES FOR PARABOLIC
More informationSolutions Numerical Simulation - Homework, 4/15/ The DuFort-Frankel method is given by. f n. x 2. (a) Truncation error: Taylor expansion
Solutios Numerical Simulatio - Homework, 4/15/211 36. The DuFort-Frakel method is give by +1 1 = αl = α ( 2 t 2 1 1 +1 + +1 (a Trucatio error: Taylor epasio i +1 t = m [ m ] m= m! t m. (1 i Sice we cosider
More informationLecture 11. Solution of Nonlinear Equations - III
Eiciecy o a ethod Lecture Solutio o Noliear Equatios - III The eiciecy ide o a iterative ethod is deied by / E r r: rate o covergece o the ethod : total uber o uctios ad derivative evaluatios at each step
More informationApproximate solutions for an acoustic plane wave propagation in a layer with high sound speed gradient
Proceedigs o Acoustics Victor Harbor 7- ovember, Victor Harbor, Australia Approximate solutios or a acoustic plae wave propagatio i a layer with high soud speed gradiet Alex Zioviev ad Adria D. Joes Maritime
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationPAijpam.eu ON TENSOR PRODUCT DECOMPOSITION
Iteratioal Joural of Pure ad Applied Mathematics Volume 103 No 3 2015, 537-545 ISSN: 1311-8080 (prited versio); ISSN: 1314-3395 (o-lie versio) url: http://wwwijpameu doi: http://dxdoiorg/1012732/ijpamv103i314
More informationLECTURE 2 LEAST SQUARES CROSS-VALIDATION FOR KERNEL DENSITY ESTIMATION
Jauary 3 07 LECTURE LEAST SQUARES CROSS-VALIDATION FOR ERNEL DENSITY ESTIMATION Noparametric kerel estimatio is extremely sesitive to te coice of badwidt as larger values of result i averagig over more
More informationTIME-PERIODIC SOLUTIONS OF A PROBLEM OF PHASE TRANSITIONS
Far East Joural o Mathematical Scieces (FJMS) 6 Pushpa Publishig House, Allahabad, Idia Published Olie: Jue 6 http://dx.doi.org/.7654/ms99947 Volume 99, umber, 6, Pages 947-953 ISS: 97-87 Proceedigs o
More informationPolynomial Generalizations and Combinatorial Interpretations for Sequences Including the Fibonacci and Pell Numbers
Ope Joural o Discrete Mathematics,,, - http://dxdoiorg/46/odm6 Published Olie Jauary (http://wwwscirporg/oural/odm) Polyomial Geeralizatios ad Combiatorial Iterpretatios or Seueces Icludig the Fiboacci
More informationThe Numerical Solution of Singular Fredholm Integral Equations of the Second Kind
WDS' Proceedigs of Cotributed Papers, Part I, 57 64, 2. ISBN 978-8-7378-39-2 MATFYZPRESS The Numerical Solutio of Sigular Fredholm Itegral Equatios of the Secod Kid J. Rak Charles Uiversity, Faculty of
More informationA New Method to Order Functions by Asymptotic Growth Rates Charlie Obimbo Dept. of Computing and Information Science University of Guelph
A New Method to Order Fuctios by Asymptotic Growth Rates Charlie Obimbo Dept. of Computig ad Iformatio Sciece Uiversity of Guelph ABSTRACT A ew method is described to determie the complexity classes of
More informationThe picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled
1 Lecture : Area Area ad distace traveled Approximatig area by rectagles Summatio The area uder a parabola 1.1 Area ad distace Suppose we have the followig iformatio about the velocity of a particle, how
More informationFourier Techniques lecture by Håkon Hoel
CSC Nada DN55 Sprig - Differetial Equatios II JOp p (7) Fourier Teciques lecture by Håo Hoel Fourier series... Applicatios to liear PDE... 3 Numerical metods...3 3. Vo Neuma aalysis...4 3. Spectral metods...6
More informationCS475 Parallel Programming
CS475 Parallel Programmig Dieretiatio ad Itegratio Wim Bohm Colorado State Uiversity Ecept as otherwise oted, the cotet o this presetatio is licesed uder the Creative Commos Attributio.5 licese. Pheomea
More informationSTAT 516 Answers Homework 6 April 2, 2008 Solutions by Mark Daniel Ward PROBLEMS
STAT 56 Aswers Homework 6 April 2, 28 Solutios by Mark Daiel Ward PROBLEMS Chapter 6 Problems 2a. The mass p(, correspods to either o the irst two balls beig white, so p(, 8 7 4/39. The mass p(, correspods
More informationAn Insight into Differentiation and Integration
Differetiatio A Isigt ito Differetiatio a Itegratio Differetiatio is basically a task to fi out ow oe variable is cagig i relatio to aoter variable, te latter is usually take as a cause of te cage. For
More informationME NUMERICAL METHODS Fall 2007
ME - 310 NUMERICAL METHODS Fall 2007 Group 02 Istructor: Prof. Dr. Eres Söylemez (Rm C205, email:eres@metu.edu.tr ) Class Hours ad Room: Moday 13:40-15:30 Rm: B101 Wedesday 12:40-13:30 Rm: B103 Course
More informationNumerical Methods in Geophysics: Implicit Methods
Numerical Methods i Geophysics: What is a implicit scheme? Explicit vs. implicit scheme for Newtoia oolig rak-nicholso Scheme (mixed explicit-implicit Explicit vs. implicit for the diffusio equatio Relaxatio
More informationADVANCED SOFTWARE ENGINEERING
ADVANCED SOFTWARE ENGINEERING COMP 3705 Exercise Usage-based Testig ad Reliability Versio 1.0-040406 Departmet of Computer Ssciece Sada Narayaappa, Aeliese Adrews Versio 1.1-050405 Departmet of Commuicatio
More informationLecture 2: Monte Carlo Simulation
STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?
More informationThe total error in numerical differentiation
AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and
More informationNumerical Methods! for! Elliptic Equations-II! Multigrid! Methods! f Analytic Solution of! g = bk. da k dt. = α. ( t)
http://www.d.edu/~gtryggva/cfd-course/! umerical Methods! or! Elliptic Equatios-II! Grétar Tryggvaso! Sprig 13! Examples o elliptic equatios! Direct Methods or 1D problems! Elemetary Iterative Methods!
More informationLIMITS AND DERIVATIVES
Capter LIMITS AND DERIVATIVES. Overview.. Limits of a fuctio Let f be a fuctio defied i a domai wic we take to be a iterval, say, I. We sall study te cocept of it of f at a poit a i I. We say f ( ) is
More information