More Elementary Aspects of Numerical Solutions of PDEs!
|
|
- Brian Bridges
- 5 years ago
- Views:
Transcription
1 ttp:// 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. Here we will examie: Stability o te liear advectio-diusio equatio Multidimesioal liear advectio-diusio equatio Steady state (boudary value problems) Grétar Tryggvaso Sprig 0 As log as accuracy is reasoable, itegratio at larger time steps is more eiciet ad desirable. Ca we icrease te time step ideiitely? Let s repeat te -D advectio-diusio equatio wit larges time step. Use Δt = 0., istead o Δt = 0.05 Evolutio or U=; D=0.05; k= N= Δt=0. Exact Numerical Istead o decayig as it sould, te amplitude o te umerical solutio keeps icreasig. Ideed, i we cotiued te calculatios, we would evetually produce umbers larger ta te computer ca adle. Tis results i a overlow or NaN (Not a Number). Ordiary Dieretial Equatio
2 Take: d dt = Te exact solutio is ( t) = e t Forward Euler ODE Example + = "t wit iitial coditio (0) = + oly i t " + = ( "t) + = ( "t) Δt=0.5 Δt=.5 Δt=.5 ODE Example + oly i t " However Backward Euler + = + "t + ODE Example + = ( "t) = ( "t) Obviously, oscillates uless = = ( "t) or all Δt t " + = ( + t) Stability Aalysis o te Advectio-Diusio Equatio: Vo Neuma Metod Geerally, stability aalysis o te ull oliear system o equatios is too ivolved to be practical, ad we study a model problem tat i some way mimics te ull equatios. Te liear advectio-diusio equatio is oe suc model equatio, ad we will apply vo Neuma's metod to ceck te stability o a simple iite dierece approximatio to tat equatio. Cosider te -D advectio-diusio equatio: + U = D t x x I iite-dierece orm: U = D " t + Look at te evolutio o a small perturbatio =
3 Te evolutio o te perturbatio is govered by: U = D " t Write te error as a wave (expad as a Fourier series): " ikx k e k = " $ = $ ( x ) = $ Droppig te subscript ikx = e Recall: e ikx = coskx + isikx Te error at ode is: = e ikx + Te error at + ad - ca be writte as = e ikx + = e ik (x + ) = e ikx e ik " = e ikx " = e ik (x ") = e ikx e "ik Substitutig = e ikx ito + " t yields + + e ikx " e ikx t = e ikx e ik + U + " " + U (eik e ikx " e " ik e ikx ) = D (eik e ikx " e ikx + e " ik e ikx ) " = e ikx e "ik = D + " + " + " t Te equatio or te error is: + + U (eik " e " ik ) = D Solvig or te ratio o te errors: " t ik ik D" t = ( e e ) + ( e (eik " + e "ik ) U ik ik + e ) + Dividig by te error amplitude at : Usig: e ik + e ik = cosk; ampliicatio actor " t ik ik D" t = ( e e ) + ( e = U"t D"t isik + (cosk ) U ik ik = 4 D"t si k i U"t sik e ik e ik = isik; + e ) si =" cos Te ratio o te error amplitude at + ad is: + =" 4 Dt si k " i Ut sik Stability requires tat + " Sice te ampliicatio actor is a complex umber, ad k, te wave umber o te error, ca be aytig, te determiatio o te stability limit is sligtly ivolved. We will look at two special cases: (a) U = 0 ad (b) D = 0
4 (a) Cosider irst te case we U = 0, so te problem reduces to a pure diusio + = " 4 Dt si k Sice si () te ampliicatio actor is always less ta, ad we id tat it is bigger ta - i " 4 Dt " Dt " (b) Cosider ow te oter limit were D = 0 ad we ave a pure advectio problem. + =" i Ut sik Sice te ampliicatio actor as te orm +i() te absolute value o tis complex umber is always larger ta uity ad te metod is ucoditioally ustable or tis case. + i Ut sik For te geeral case we must ivestigate te stability coditio i more detail. We will ot do so ere, but simply quote te results: For a two-dimesioal problem, assume a error o te orm A stability aalysis gives: i, = i(kxi +ly ) e td " ad U t D " D" t 4 ( U + V ) " t ad 4 D Notice tat ig velocity ad low viscosity lead to istability accordig to te secod restrictio. For a tree-dimesioal problem we get: Dt " 6 ad ( U + V + W ) t " 8 D Stability Now you kow Covergece te solutio to te iite-dierece equatio approaces te true solutio to te PDE avig te same iitial ad boudary coditios as te mes is reied. Lax s Equivalece Teorem Give a properly posed iitial value problem ad a iitedierece approximatio to it tat satisies te cosistecy coditio, stability is te ecessary ad suiciet coditio or covergece. Summary Itroduced a ormal metod to examie weter a give iite dierece approximatio is stable or ot Itroduced Lax s equivalet teorem.
5 Multidimesioal Equatios Two-Dimesioal Advectio- Diusio Equatio We will use te model equatio: +U t x + V y = D " $ x + % y & to demostrate ow to solve a partial equatio (iitial value problem) umerically. Te extesio to two-dimesios is relatively straigt orward, oce te oe-dimesioal problem is ully uderstood. Multidimesioal Equatios For a two-dimesioal low discretize te variables o a two-dimesioal grid + - (x, y) i - i i+ i, + = (x, y + ) i, = (x,y) i +, = (x +, y) + i, i, Δt Multidimesioal Equatios t + U x + V y = D x Te discrete equatio is: = U i+, D i, i+, + y V i, + i, + i, i, + i, i, i, + Solve or i, or + i, = i, + i, = i, Multidimesioal Equatios + Δt U i+, D i+, ΔtU i, + i, V i, + + i, + i, + i, ( i+, i, ) ΔtV ( i, + i, ) ( ) Accuracy: O t, + ΔtD A stability aalysis gives: + 4 i, i+, + i, + i, + + ( i, 4 i, ) td " 4 ad ( U + V ) t " 4 D Example
6 Multidimesioal Equatios t + U x + V y = D x + y + i, = i, "tu Multidimesioal Equatios ( i+, i, ) + "td i +, =NY + i, + i, + + ( i, 4 i, ) x = 0 y = 0 =0.0 =.0 Uiorm low troug te domai U=- V=0 i, stored at eac grid poit or give o te boudary =0.0 = = i= i= i=nx Boudary coditios Multidimesioal Equatios Were is give, we simply speciy its value Were te ormal derivative is speciied, we approximate te value at te boudary by oesided diereces At te i= boudary, or example, y = 0 y i, ad by usig i, = 0 we id tat: i, = i, Multidimesioal Equatios % two-dimesioal usteady diusio by te FTCS sceme % =3;m=3;step=0;D=0.05;legt=.0;=legt/(-); dt=.0*0.5**/d;=zeros(,m);o=zeros(,m);time=0.0; u=-0.0; v=-.0; (:,)=.0; or l=:step,l,time old o;mes(); axis([0 0 m 0.5]);pause; o=; or i=:-, or =:m- (i,)=o(i,)-(0.5*dt*u/)*(o(i+,)-o(i-,))-... (0.5*dt*v/)*(o(i,+)-o(i,-))+... (D*dt/^)*(o(i+,)+ o(i,+)+o(i-,)+o(i, -)-4*o(i,)); ed,ed or i=:, (i,)=(i,);ed;or =:m,(,)=(,);(m,)=(m-,);ed; time=time+dt; ed; Te usteady evolutio o te solutio Multidimesioal Equatios Multidimesioal Boudary Value Problems (Steady-State)
7 Boudary Value Problems Cosider te Poisso Equatio: + = S x y Tis equatio as a solutio i or is speciied o te boudary Use stadard iite diereces to discretize: i +, + i, i, + + i, i, + i, = S i, For uiorm grids: i+, + i, i, ca be writte as Solve or i, : Boudary Value Problems + i, + i, + i, i+, + i, + i, + i, + 4 i, = S = S i, i, = ( i +, i, i, i, + S i, ) Boudary Value Problems Boudary Value Problems Solve or i, ad use te rigt ad side to compute a ew value. Deote te old values by α ad te ew oes wit α+ i, ( ) + = 4 i +, + i", + i, " + i, + " S i, Tis iteratio process Jacobi iteratio is very robust but may iteratios are required to reac a accurate solutio. Te iteratio must be carried out util te solutio is suicietly accurate. To measure te error, deie te residual: R i, = i+, + i, + i, + i, + 4 i, S i, At steady-state te residual sould be zero. Te poitwise residual or te average absolute residual ca be used, depedig o te problem. Ote, simpler criteria, suc as te cage rom oe iteratio to te ext is used Boudary Value Problems Altoug te Jacobi iteratio is a very robust iteratio tecique, it coverges VERY slowly. We tereore seek a way to ACCELERATE te covergece to steady-state, makig use o te act tat it is oly te steady-state tat is o iterest. Here we itroduce te Gauss-Seidler metod ad te Successive Over-Relaxatio (SOR) metod. Boudary Value Problems Te Jacobi iteratio ca be improved somewat by usig ew values as soo as tey become available. + - i- i i+ + i, = 4 ( i +, + i", + + i, + or =:m or i=: iterate ed ed + + i, " " S i, ) From a programmig poit o view, Gauss-Seidler iteratio is eve simpler ta Jacobi iteratio sice oly oe vector wit values is eeded.
8 Boudary Value Problems Te Gauss-Seidler iteratio ca be accelerated eve urter by various acceleratio teciques. Te simplest oe is te Successive Over-Relaxatio (SOR) iteratio + i, = " 4 ( i +, + i, + ( ") i, + + i, i, + S i, ) Example Te SOR iteratio is very simple to program, ust as te Gauss-Seidler iteratio. Te user must select te coeiciet. It must be bouded by <β<. β=.5 is usually a good startig value. Boudary Value Problems Boudary Value Problems x + y = 0 = = % two-dimesioal steady-state problem by SOR =40;m=40;iteratios=5000;legt=.0;=legt/(-); T=zeros(,m);bb=.7; T(0:-0,)=.0; or l=:iteratios, or i=:-, or =:m- T(i,)=bb*0.5*(T(i+,)+... T(i,+)+T(i-,)+T(i,-))+(.0-bb)*T(i,); ed,ed % id residual res=0; or i=:-, or =:m- res=res+abs(t(i+,)+... T(i,+)+T(i-,)+T(i,-)-4*T(i,))/^; ed,ed l,res/((m-)*(-)) % Prit iteratio ad residual i (res/((m-)*(-)) < 0.00), break,ed ed; cotour(t); Boudary Value Problems Te program is easily modiied or te Jacobi ad te Gauss-Seidler iteratio: Average absolute error: 0.00 Number o iteratios Jacobi: 989 Gauss-Seidler: 986 SOR (β =.5): 30 SOR (β =.7): 6 SOR (β =.9): 9 SOR (β =.95): Te coverged solutio: Boudary Value Problems x + y =
9 Summary Stability itroduced te vo Neuma metod. Fairly mecaical process, we will provide more isigt by te iite volume poit o view Multidimesioal advectio-diusio equatio. Essetially te same as te oe-dimesioal problem Iterative metods or boudary value problems. Elemetary approaces to steady state problems
Accuracy. 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 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 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 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 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 informationAnalysis of a Numerical Scheme An Example
http://www.d.edu/~gtryggva/cfd-course/ Computatioal Fluid Dyamics Lecture 3 Jauary 5, 7 Aalysis of a Numerical Scheme A Example Grétar Tryggvaso Numerical Aalysis Example Use the leap-frog method (cetered
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 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 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 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/! Numerical Methods! or! Elliptic Equatios-II! Grétar Tryggvaso! Sprig! Examples o elliptic equatios! Direct Methods or D problems! Elemetary Iterative Methods! Note
More information( ) κ ( f j +1. ( ) 1 + κ ( f j +2. ( ) j +1/ 2 ) = Uf j +1/ 2. = r j +1/ 2. Ψ( r) = r + r 1 + r. Ψ( r) = minmod(r,1) j 1/ 2 ( F n +1/ 2 F n
ttp://users.wpi.edu/~gretar/me6.tml! Numerical Metods! or Hyperbolic! Equatios VIII! Grétar Tryggvaso! Sprig! Upwid or oliear systems:!!flux vector splittig!!goduovʼs metod! Higer Order Metod!!Higer order
More informationComputational Fluid Dynamics. Lecture 3
Computatioal Fluid Dyamics Lecture 3 Discretizatio Cotiued. A fourth order approximatio to f x ca be foud usig Taylor Series. ( + ) + ( + ) + + ( ) + ( ) = a f x x b f x x c f x d f x x e f x x f x 0 0
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 informationSolving 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 informationUNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014
UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 6C Problem Set 4 Bejami Stahl November 6, 4 BOAS, P. 63, PROBLEM.-5 The Laguerre differetial equatio, x y + ( xy + py =, will be solved
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 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 informationNumerical Derivatives by Symbolic Tools in MATLAB
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
More informationNumerical Methods for Elliptic Equations-I. Computational Fluid Dynamics. φ f = S. x + 2 f. x a f x + b f x + cf = s. f y b f.
http://www.d.edu/~gtryggva/cfd-course/ http://www.d.edu/~gtryggva/cfd-course/ Computatioal Fluid Dyamics Lecture 13 March 6, 17 umerical Methods or Elliptic Equatios-I Grétar Tryggvaso Examples o elliptic
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
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 information4 Conditional Distribution Estimation
4 Coditioal Distributio Estimatio 4. Estimators Te coditioal distributio (CDF) of y i give X i = x is F (y j x) = P (y i y j X i = x) = E ( (y i y) j X i = x) : Tis is te coditioal mea of te radom variable
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 informationFinite Dierence Schemes
MATH-459 Numerical Methods for Coservatio Laws by Prof. Ja S. Hesthave Solutio set 2: Fiite Dierece Schemes Exercise 2. Cosistecy A method is cosistet if its local trucatio error T k satises T k (x, t)
More informationTaylor expansion: Show that the TE of f(x)= sin(x) around. sin(x) = x - + 3! 5! L 7 & 8: MHD/ZAH
Taylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. A ay poit i the eighbourhood of 0, the fuctio ƒ() ca be represeted by a power series of the followig form: X 0 f(a) f() f() ( ) f( ) ( )
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 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 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 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 informationb i u x i U a i j u x i u x j
M ath 5 2 7 Fall 2 0 0 9 L ecture 1 9 N ov. 1 6, 2 0 0 9 ) S ecod- Order Elliptic Equatios: Weak S olutios 1. Defiitios. I this ad the followig two lectures we will study the boudary value problem Here
More informationCastiel, Supernatural, Season 6, Episode 18
13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio
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 informationOptimization Methods MIT 2.098/6.255/ Final exam
Optimizatio Methods MIT 2.098/6.255/15.093 Fial exam Date Give: December 19th, 2006 P1. [30 pts] Classify the followig statemets as true or false. All aswers must be well-justified, either through a short
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 informationNumerical Astrophysics: hydrodynamics
Numerical Astrophysics: hydrodyamics Part 1: Numerical solutios to the Euler Equatios Outlie Numerical eperimets are a valuable tool to study astrophysical objects (where we ca rarely do direct eperimets).
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 informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More informationLecture 7 Testing Nonlinear Inequality Restrictions 1
Eco 75 Lecture 7 Testig Noliear Iequality Restrictios I Lecture 6, we discussed te testig problems were te ull ypotesis is de ed by oliear equality restrictios: H : ( ) = versus H : ( ) 6= : () We sowed
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 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 information1 The Euler Forward scheme (schéma d Euler explicite)
TP : Fiite differece metod for a Europea optio M2 Modélisatio aléatoire - Uiversité Deis-Diderot Cours EDP e Fiace et Métodes Numériques December 8, 2017 1 Te Euler Forward sceme (scéma d Euler explicite)
More informationMost text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t
Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said
More informationBeyond simple iteration of a single function, or even a finite sequence of functions, results
A Primer o the Elemetary Theory of Ifiite Compositios of Complex Fuctios Joh Gill Sprig 07 Abstract: Elemetary meas ot requirig the complex fuctios be holomorphic Theorem proofs are fairly simple ad are
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 information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationMath 312 Lecture Notes One Dimensional Maps
Math 312 Lecture Notes Oe Dimesioal Maps Warre Weckesser Departmet of Mathematics Colgate Uiversity 21-23 February 25 A Example We begi with the simplest model of populatio growth. Suppose, for example,
More informationWave Equation! ( ) with! b = 0; a =1; c = c 2. ( ) = det ( ) = 0. α = ±c. α = 1 2a b ± b2 4ac. c 2. u = f. v = f x ; t c v. t u. x t. t x = 2 f.
Compuaioal Fluid Dyamics p://www.d.edu/~gryggva/cfd-course/ Compuaioal Fluid Dyamics Wave equaio Wave Equaio c Firs wrie e equaio as a sysem o irs order equaios Iroduce u ; v ; Gréar Tryggvaso Sprig yieldig
More informationLøsningsførslag i 4M
Norges tekisk aturviteskapelige uiversitet Istitutt for matematiske fag Side 1 av 6 Løsigsførslag i 4M Oppgave 1 a) A sketch of the graph of the give f o the iterval [ 3, 3) is as follows: The Fourier
More informationPAPER : IIT-JAM 2010
MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure
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 informationNumerical Method for Blasius Equation on an infinite Interval
Numerical Method for Blasius Equatio o a ifiite Iterval Alexader I. Zadori Omsk departmet of Sobolev Mathematics Istitute of Siberia Brach of Russia Academy of Scieces, Russia zadori@iitam.omsk.et.ru 1
More informationDEGENERACY AND ALL THAT
DEGENERACY AND ALL THAT Te Nature of Termodyamics, Statistical Mecaics ad Classical Mecaics Termodyamics Te study of te equilibrium bulk properties of matter witi te cotext of four laws or facts of experiece
More informationTR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT
TR/46 OCTOBER 974 THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION by A. TALBOT .. Itroductio. A problem i approximatio theory o which I have recetly worked [] required for its solutio a proof that the
More informationIterative Techniques for Solving Ax b -(3.8). Assume that the system has a unique solution. Let x be the solution. Then x A 1 b.
Iterative Techiques for Solvig Ax b -(8) Cosider solvig liear systems of them form: Ax b where A a ij, x x i, b b i Assume that the system has a uique solutio Let x be the solutio The x A b Jacobi ad Gauss-Seidel
More informationRecursive Algorithms. Recurrences. Recursive Algorithms Analysis
Recursive Algorithms Recurreces Computer Sciece & Egieerig 35: Discrete Mathematics Christopher M Bourke cbourke@cseuledu A recursive algorithm is oe i which objects are defied i terms of other objects
More informationPartial Differential Equations
Partial Differetial Equatios Part 2 Massimo Ricotti ricotti@astro.umd.edu Uiversity of Marylad Partial Differetial Equatios p.1/15 Upwid differecig I additio to amplitude errors (istability or dampig),
More informationChapter 8. Euler s Gamma function
Chapter 8 Euler s Gamma fuctio The Gamma fuctio plays a importat role i the fuctioal equatio for ζ(s that we will derive i the ext chapter. I the preset chapter we have collected some properties of the
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 informationIntroduction to the Multigrid Method
Semiar: Te Iterplay betwee Matematical Modellig ad Numerical Simulatio Itroductio to te Multigrid Metod Bogojeska Jasmia JASS, 005 Abstract Te metods for solvig liear systems of equatios ca be divided
More information(b) What is the probability that a particle reaches the upper boundary n before the lower boundary m?
MATH 529 The Boudary Problem The drukard s walk (or boudary problem) is oe of the most famous problems i the theory of radom walks. Oe versio of the problem is described as follows: Suppose a particle
More informationIntroduction to Optimization Techniques. How to Solve Equations
Itroductio to Optimizatio Techiques How to Solve Equatios Iterative Methods of Optimizatio Iterative methods of optimizatio Solutio of the oliear equatios resultig form a optimizatio problem is usually
More informationwavelet collocation method for solving integro-differential equation.
IOSR Joural of Egieerig (IOSRJEN) ISSN (e): 5-3, ISSN (p): 78-879 Vol. 5, Issue 3 (arch. 5), V3 PP -7 www.iosrje.org wavelet collocatio method for solvig itegro-differetial equatio. Asmaa Abdalelah Abdalrehma
More informationDefinition 4.2. (a) A sequence {x n } in a Banach space X is a basis for X if. unique scalars a n (x) such that x = n. a n (x) x n. (4.
4. BASES I BAACH SPACES 39 4. BASES I BAACH SPACES Sice a Baach space X is a vector space, it must possess a Hamel, or vector space, basis, i.e., a subset {x γ } γ Γ whose fiite liear spa is all of X ad
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 informationNumerical Methods in Fourier Series Applications
Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets
More informationPROBLEM SET 5 SOLUTIONS 126 = , 37 = , 15 = , 7 = 7 1.
Math 7 Sprig 06 PROBLEM SET 5 SOLUTIONS Notatios. Give a real umber x, we will defie sequeces (a k ), (x k ), (p k ), (q k ) as i lecture.. (a) (5 pts) Fid the simple cotiued fractio represetatios of 6
More informationL 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5!
aylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. At ay poit i the eighbourhood of =0, the fuctio ca be represeted as a power series of the followig form: X 0 f(a) f() ƒ() f()= ( ) f( ) (
More informationChapter 4 : Laplace Transform
4. Itroductio Laplace trasform is a alterative to solve the differetial equatio by the complex frequecy domai ( s = σ + jω), istead of the usual time domai. The DE ca be easily trasformed ito a algebraic
More informationOn Strongly Consistent Finite Dierence Approximations
D.Michels et al. (KAUST,JINR,SSU) Strogly Cosistet Approximatios 17 April 2018 1 / 25 O Strogly Cosistet Fiite Dierece Approximatios Domiik Michels 1, Vladimir Gerdt 2, Dmitry Lyakhov 1, ad Yuri Blikov
More informationTHE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.
THE SOLUTION OF NONLINEAR EQUATIONS f( ) = 0. Noliear Equatio Solvers Bracketig. Graphical. Aalytical Ope Methods Bisectio False Positio (Regula-Falsi) Fied poit iteratio Newto Raphso Secat The root of
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 informationBoundary layer problem on conveyor belt. Gabriella Bognár University of Miskolc 3515 Miskolc-Egyetemváros, Hungary
Boudary layer problem o coveyor belt Gabriella Bogár Uiversity of Miskolc 355 Miskolc-Egyetemváros, Hugary e-mail: matvbg@ui-miskolc.hu Abstract: A techologically importat source of the boudary layer pheomeo
More information1 Approximating Integrals using Taylor Polynomials
Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................
More informationSequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet
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 informationSOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.
SOLUTIONS TO EXAM 3 Problem Fid the sum of the followig series 2 + ( ) 5 5 2 5 3 25 2 2 This series diverges Solutio: Note that this defies two coverget geometric series with respective radii r 2/5 < ad
More informationA sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as
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 informationPolynomials with Rational Roots that Differ by a Non-zero Constant. Generalities
Polyomials with Ratioal Roots that Differ by a No-zero Costat Philip Gibbs The problem of fidig two polyomials P(x) ad Q(x) of a give degree i a sigle variable x that have all ratioal roots ad differ by
More informationChapter 8. Euler s Gamma function
Chapter 8 Euler s Gamma fuctio The Gamma fuctio plays a importat role i the fuctioal equatio for ζ(s) that we will derive i the ext chapter. I the preset chapter we have collected some properties of the
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 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 informationThe Phi Power Series
The Phi Power Series I did this work i about 0 years while poderig the relatioship betwee the golde mea ad the Madelbrot set. I have fially decided to make it available from my blog at http://semresearch.wordpress.com/.
More informationThe Growth of Functions. Theoretical Supplement
The Growth of Fuctios Theoretical Supplemet The Triagle Iequality The triagle iequality is a algebraic tool that is ofte useful i maipulatig absolute values of fuctios. The triagle iequality says that
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 informationLecture 2: Finite Difference Methods in Heat Transfer
Lecture 2: Fiite Differece Methods i Heat Trasfer V.Vuorie Aalto Uiversity School of Egieerig Heat ad Mass Trasfer Course, Autum 2016 November 1 st 2017, Otaiemi ville.vuorie@aalto.fi Overview Part 1 (
More informationSimilarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall
Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig
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 informationENGI 9420 Lecture Notes 3 - Numerical Methods Page 3.01
ENGI 940 Lecture Notes 3 - Numerical Methods Page 3.0 3. Numerical Methods The majority of equatios of iterest i actual practice do ot admit ay aalytic solutio. Eve equatios as simple as = e ad I = e d
More informationSeries III. Chapter Alternating Series
Chapter 9 Series III With the exceptio of the Null Sequece Test, all the tests for series covergece ad divergece that we have cosidered so far have dealt oly with series of oegative terms. Series with
More informationNotes on iteration and Newton s method. Iteration
Notes o iteratio ad Newto s method Iteratio Iteratio meas doig somethig over ad over. I our cotet, a iteratio is a sequece of umbers, vectors, fuctios, etc. geerated by a iteratio rule of the type 1 f
More informationn 3 ln n n ln n is convergent by p-series for p = 2 > 1. n2 Therefore we can apply Limit Comparison Test to determine lutely convergent.
06 微甲 0-04 06-0 班期中考解答和評分標準. ( poits) Determie whether the series is absolutely coverget, coditioally coverget, or diverget. Please state the tests which you use. (a) ( poits) (b) ( poits) (c) ( poits)
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 informationCHAPTER 11 Limits and an Introduction to Calculus
CHAPTER Limits ad a Itroductio to Calculus Sectio. Itroductio to Limits................... 50 Sectio. Teciques for Evaluatig Limits............. 5 Sectio. Te Taget Lie Problem................. 50 Sectio.
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 informationf(x) dx as we do. 2x dx x also diverges. Solution: We compute 2x dx lim
Math 3, Sectio 2. (25 poits) Why we defie f(x) dx as we do. (a) Show that the improper itegral diverges. Hece the improper itegral x 2 + x 2 + b also diverges. Solutio: We compute x 2 + = lim b x 2 + =
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 information-ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION
NEW NEWTON-TYPE METHOD WITH k -ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION R. Thukral Padé Research Cetre, 39 Deaswood Hill, Leeds West Yorkshire, LS7 JS, ENGLAND ABSTRACT The objective
More informationNotes for Lecture 11
U.C. Berkeley CS78: Computatioal Complexity Hadout N Professor Luca Trevisa 3/4/008 Notes for Lecture Eigevalues, Expasio, ad Radom Walks As usual by ow, let G = (V, E) be a udirected d-regular graph with
More information