Onesided finitedifference approximations suitable for use with Richardson extrapolation


 Lester Logan
 4 years ago
 Views:
Transcription
1 Journal of Computatonal Physcs 219 (2006) Short note Onesded fntedfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng, Indan Insttute of Technology, Kharagpur , West Bengal, Inda Receved 30 March 2006; receved n revsed form 12 May 2006; accepted 29 May 2006 Avalable onlne 10 July Abstract New expressons for onesded fntedfference approxmatons are proposed. In these approxmatons the oddorder error terms are elmnated whle the evenorder terms are left to be taken care of by Rchardson extrapolaton. The effectve local truncaton error s shown to be less than for hgherorder onesded fntedfference approxmatons but the solutons for a test problem are shown to have comparable accuracy for both approxmatons. Ó 2006 Elsever Inc. All rghts reserved. Keywords: Fntedfference method; Onesded dfference approxmaton; Order of accuracy; Rchardson extrapolaton 1. Introducton A fntedfference approxmaton s one of the commonly used methods for numercal soluton of ordnary and partal dfferental equatons. The approxmatons most often used have secondorder accuracy. The order of accuracy can be ncreased ether by usng hgherorder fnte dfference approxmatons or by Rchardson extrapolaton [2,3]. For nternal grd ponts, f central dfference approxmatons are used, the truncaton error has the form A 2 h 2 + A 4 h 4 + A 6 h 6 +, where h s the dstance between adjacent grd ponts and A 2,A 4,A 6,... are constants. Dong the computaton wth two or three dfferent grd szes and carryng out one or two extrapolatons to elmnate the error terms of second order or of second and fourth order we obtan a soluton wth fourth or sxthorder accuracy. Snce each of the computatons uses a secondorderaccurate fntedfference approxmaton there s no dffculty wth the boundary condtons for a Drchlet problem. However, a problem arses when the boundary condton nvolves a dervatve and we use a onesded fntedfference approxmaton for the dervatve. By usng three grd ponts, ncludng one on the boundary, we can obtan a onesded approxmaton whch s secondorderaccurate [1]. The truncaton error has the form B 2 h 2 + B 3 h 3 + B 4 h 4 +. So f we carry out one extrapolaton to elmnate the secondorder error term the result would have thrdorder rather than fourthorder accuracy. One way of obtanng fourthorder accuracy * Correspondng author. Tel.: ; fax: Emal address: (S.N. Bhattacharyya) $  see front matter Ó 2006 Elsever Inc. All rghts reserved. do: j.jcp
2 14 K. Rahul, S.N. Bhattacharyya Journal of Computatonal Physcs 219 (2006) s by usng a onesded dfference approxmaton whch has a hgher order of accuracy. Such hgherorderaccurate onesded fntedfference approxmatons are known [1]. To obtan fourthorder accuracy, one could use a onesded fntedfference approxmaton whch has fourthorder accuracy. Ths would nvolve fve grd ponts, ncludng one on the boundary. The addtonal two ponts are requred n order to elmnate the terms of the order of h 2 and h 3 n the truncaton error. However, f we are gong to use Rchardson extrapolaton t s not necessary to elmnate the order h 2 term n the truncaton error for the onesded fntedfference approxmatons. So we can use just one more pont to obtan a onesded dfference approxmaton whch has a truncaton error of the form B 2 h 2 + B 4 h 4 + B 5 h 5 +. After one extrapolaton, ths would gve a result whch has fourthorder accuracy. Snce ths approxmaton nvolves smaller number of grd ponts compared to the fourthorder approxmaton whch nvolves ponts farther away from the boundary, although both the methods have fourthorder accuracy, we expect the coeffcent of h 4 n the truncaton error to be smaller n the frst case. Contnung n ths manner, f we requre sxthorder accuracy, we can use one more pont to obtan a onesded fntedfference approxmaton whch has a truncaton error of the form of B 2 h 2 + B 4 h 4 + B 6 h 6 + B 7 h 7 +. In ths study, we derve onesded fntedfference approxmatons wth oddorder terms n the truncaton error elmnated and the evenorder terms left to be taken care of by Rchardson extrapolaton. We compare the local truncaton error wth that of hgherorder onesded fntedfference approxmatons. We then apply these approxmatons to a test problem and check the accuracy of the solutons obtaned. 2. Onesded fntedfference approxmatons Let us consder a functon (x) and defne a set of unformly spaced grd ponts x, =0,1,...,N, wth the grd spacng h = x +1 x. We represent (x )by. Smlarly 0 (x ), 00 (x ),..., are represented by 0 ; 00 ;..., where the prmes denote dfferentaton. We can wrte a Taylor seres expanson for around x k ¼ kh 0 þ k2 h k3 h þ k4 h 4 24 ðvþ k5 h ðvþ þ k6 h ðvþ : ð1þ From Eq. (1) we readly obtan frstorderaccurate onesded fntedfference approxmatons [1] 0 ¼ ð 1 Þ þ OðhÞ; ð2þ h where the upper and lower sgns are to be used to the rght and to the left of x. Wrtng Eq. (1) for k = 1 and 2 and elmnatng 00 between them, we obtan secondorderaccurate onesded fntedfference approxmatons [1] 0 ¼ ð3 4 1 þ 2 Þ þ Oðh 2 Þ: ð3þ 2h We now derve an approxmaton for use wth Rchardson extrapolaton whch has truncaton error of the form B 2 h 2 + B 4 h 4 + B 5 h 5 +. Wrtng Eq. (1) for k = 1, 2 and 3 and elmnatng 00 and ðvþ between them we obtan the onesded approxmatons 0 ¼ ð þ Þ þ h 2 1 6h h ðvþ þ: ð4þ Alternatvely, we can use a fourthorderaccurate onesded fntedfference approxmaton obtaned by wrtng Eq. (1) for k = 1, 2, 3, and 4 and elmnatng 00, 000 and ðvþ, to obtan 0 ¼ ð þ þ 3 4 Þ þ h h 120 ðvþ þ: ð5þ When usng the approxmatons gven by Eq. (4), f we use Rchardson extrapolaton to elmnate the secondorder error term then the leadng order error s due to the fourthorder term. Effectvely the local truncaton error for ths approxmaton becomes ð11=120þh 4 ðvþ. If nstead of Eq. (4) we use Eq. (5) the local truncaton error s ð24=120þh 4 ðvþ. Thus the onesded fntedfference approxmatons we have proposed, gven by Eq. (4), when used together wth Rchardson extrapolaton, effectvely gve a local truncaton error whch s
3 K. Rahul, S.N. Bhattacharyya Journal of Computatonal Physcs 219 (2006) φ =0 2 φ =1 2 φ =0 φ =0 φ =0 ( φ x ) x=1 = r( φ x ) x=1+ φ =0 Fg. 1. The test problem. slghtly less than half the local truncaton error when usng hgherorder onesded approxmatons gven by Eq. (5). Contnung the logc used n dervng Eq. (4) we can derve onesded fntedfference approxmatons whch have a truncaton error of the form of B 2 h 2 + B 4 h 4 + B 6 h 6 + B 7 h 7 +. These are 0 ¼ ð þ þ 4 Þ : ð6þ 20h When used together wth central dfference approxmatons for nternal grd ponts, f we carry out two Rchardson extrapolatons to elmnate the error terms of order h 2 and h 4, we can obtan a soluton wth sxthorder accuracy. Table 1 Soluton of the test problem usng a frstorderaccurate onesded fntedfference approxmaton at the nterface N N N N=2 N=4 N2 N N=
4 16 K. Rahul, S.N. Bhattacharyya Journal of Computatonal Physcs 219 (2006) Applcaton to a test problem In ths secton, we use the fntedfference approxmatons derved n the prevous secton to compute the numercal soluton for a test problem and check the order of accuracy obtaned. The test problem s defned as r 2 ¼ 1000 for 0 < x < 1; 0 < y < 1; r 2 ¼ 0 for 1 < x < 2; 0 < y < 1; o ox x¼1 ¼ r o ox x¼1þ for 0 < y < 1; ¼ 0 for x ¼ 0or2; 0 6 y 6 1; ¼ 0 for 0 6 x 6 2; y ¼ 0or1: Ths s shown n Fg. 1. For most of our numercal calculatons we use r = 2. We defne a set of grd ponts (x j,y k ) where x j = jh for 0 6 x 6 2N, y k = kh for 0 6 x 6 N and h =1N. For the nternal grd ponts, except the grd ponts at the nterface x = 1, we use a secondorder central dfference approxmaton ðr 2 Þ j;k ¼ j;k 1 þ j 1;k 4 j;k þ jþ1;k þ j;kþ1 : ð8þ h 2 The truncaton error for ths approxmaton has the form A 2 h 2 + A 4 h 4 + A 6 h 6 +. The boundary condton = 0 on the doman boundary s easly taken nto account. For the condton at the nterface, x = 1, we study the effect of usng the dfferent onesded fntedfference approxmatons derved n the prevous secton. To ð7þ Table 2 Soluton of the test problem usng a secondorderaccurate onesded fntedfference approxmaton at the nterface N N N N=2 N=4 N2 N N,N2 N,N2 N=2;N=4 N=4;N=8 N2,N4 N=2 N;N=2 N=2;N=
5 begn wth we use the frstorderaccurate onesded approxmatons gven by Eq. (2). The nterface condton n fntedfference approxmaton becomes N;j N 1;j h K. Rahul, S.N. Bhattacharyya Journal of Computatonal Physcs 219 (2006) ¼ r Nþ1;j N;j : ð9þ h The system of equatons gven by Eqs. (8) and (9) are solved numercally. The solutons at x = 0.25, 0.75, 1.25 and 1.75 and y = for dfferent number of grd ponts are shown n Table 1, where N represents the soluton obtaned usng N grd ntervals along each coordnate axs on each unt square. We observe that at x = 0.25 and y = for small number of grd ponts t almost appears that we have secondorder accuracy but as we ncrease the number of grd ponts the order of accuracy s seen to decrease. Near the nterface x = 0.75 and 1.25, y = we observe close to frstorder accuracy. Further at we observe that the accuracy s even worse but, wth ncrease n number of grd ponts, appears to tend to frst order. Ths probably occurs because the soluton n the square on the rght s drven by the source term n the square on the left through the nterface condton. So an naccurate modellng of the nterface condton leads to naccurate soluton n the entre square on the rght, though why the order of accuracy at x = 1.75 seems worse than at x = 1.25 s not clear. Thus we fnd that although the approxmatons for the nternal grd ponts s secondorderaccurate, snce we use a frstorderaccurate approxmaton for the condton at the nterface the soluton s only frstorderaccurate. We repeat the calculaton usng secondorderaccurate onesded approxmatons gven by Eq. (3) for the condton at the nterface. The results are shown n Table 2 and we clearly see secondorder accuracy. We next carry out Rchardson extrapolaton to elmnate the h 2 term n the error. The extrapolated values are shown n Table 3 Soluton of the test problem usng at the nterface a secondorderaccurate onesded fntedfference approxmaton wth the thrdorder error term elmnated N N N,N2 N,N2 N=2;N=4 N=4;N=8 N2,N4 N;N=2 N=2;N=
6 18 K. Rahul, S.N. Bhattacharyya Journal of Computatonal Physcs 219 (2006) Table 4 Soluton of the test problem usng a fourthorderaccurate onesded fntedfference approxmaton at the nterface N N N,N2 N,N2 N=2;N=4 N=4;N=8 N2,N4 N;N=2 N=2;N= the table as N,N2. Here N,N2 denotes the value obtaned by the extrapolaton usng N and N2. The rato of the dfferences appear close to 8 ndcatng thrdorder accuracy. If the nterface condton were not present we should have got fourthorder accuracy after one extrapolaton. In order to obtan fourthorder accuracy after one extrapolaton we use onesded fntedfference approxmatons gven by Eq. (4). The results are shown n Table 3. We observe that after one extrapolaton we clearly obtan fourthorder accuracy. Next we use the approxmatons gven by Eq. (5) and the results are shown n Table 4. Agan after one extrapolaton we obtan fourthorder accuracy. Comparng Tables 3 and 4 we fnd that although the local truncaton error s smaller for the approxmatons n Eq. (4) the accuracy of the soluton for the test problem s smlar to that obtaned usng Eq. (5). However, Eq. (4) stll has the advantage that t s easer to derve. Next we use the approxmaton gven by Eq. (6) and the results are shown n Table 5. Here N,N2,N4 s the value obtaned by usng Rchardson extrapolaton to elmnate the order h 4 error term between N,N2 and N2,N4. After two extrapolatons the dfferences are very small and consequently, t s dffcult to show that we obtan sxthorder accuracy. However, we do observe that agreement to 9 sgnfcant fgures s obtaned wth reasonable number of grd ponts. Alternatvely f we were to use sxthorderaccurate onesded fnte dfference approxmatons the algebra nvolved n dervng the expressons would be extremely tedous. Thus the advantage of Eq. (6) s that, when used together wth Rchardson extrapolaton, t provdes hgh accuracy wthout requrng very tedous algebra. So far the calculatons have been carred out for r = 2. We now check whether the onesded approxmatons work well for larger values of r. The results for the test problem wth r = 10 computed usng Eq. (4) are shown n Table 6. Agan we fnd that ths provdes fourthorder accuracy after one Rchardson extrapolaton, as expected. The magntude of the errors are also comparable wth those for r =2.
7 K. Rahul, S.N. Bhattacharyya Journal of Computatonal Physcs 219 (2006) Table 5 Soluton of the test problem usng at the nterface a secondorderaccurate onesded fntedfference approxmaton wth the thrd and ffthorder error terms elmnated N N N,N2 N,N2,N Table 6 Soluton of the test problem usng at the nterface a secondorderaccurate onesded fntedfference approxmaton wth the thrdorder error term elmnated and wth r = 10 N N N,N2 N,N2 N=2;N=4 N=4;N=8 N2,N4 N;N=2 N=2;N= (contnued on next page)
8 20 K. Rahul, S.N. Bhattacharyya Journal of Computatonal Physcs 219 (2006) Table 6 (contnued) N N N,N2 N,N2 N=2;N=4 N=4;N=8 N2,N4 N;N=2 N=2;N= Concluson In ths study, we have proposed some onesded fntedfference approxmatons for use wth Rchardson extrapolaton. The essental logc s to use extra grd ponts to elmnate the oddorder terms n the truncaton error but leave the evenorder terms to be elmnated by Rchardson extrapolaton. Usng a test problem we have demonstrated that the computed results have the order of accuracy we would expect. These onesded fntedfference approxmatons, when used together wth Rchardson extrapolaton are shown to have smaller local truncaton error than the conventonal hgherorder onesded fntedfference approxmatons. However, for the test problem for both approxmatons the solutons have comparable accuracy. The onesded dfference approxmatons proposed n ths manuscrpt nvolve smaller number of grd ponts and are easer to derve than the conventonal onesded approxmatons whch provde same order of accuracy. References [1] J.D. Anderson Jr., Computatonal Flud Dynamcs, McGrawHll, New York, [2] R.L. Burden, J.D. Fares, Numercal Analyss, BrooksCole, Pacfc Grove, CA, [3] G. Dahlqust, Å. Björck, Numercal Methods (Translated by N. Anderson), PrentceHall, Englewood Clffs, NJ, 1974.
NUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06FnteDfference Method (Onedmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTERVI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3d rotaton around the xaxs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3d rotaton around the yaxs:
More informationIndeterminate pinjointed frames (trusses)
Indetermnate pnjonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationLecture 2: Numerical Methods for Differentiations and Integrations
Numercal Smulaton of Space Plasmas (I [AP4036] Lecture 2 by LngHsao Lyu March, 2018 Lecture 2: Numercal Methods for Dfferentatons and Integratons As we have dscussed n Lecture 1 that numercal smulaton
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More information2.29 Numerical Fluid Mechanics
REVIEW Lecture 10: Sprng 2015 Lecture 11 Classfcaton of Partal Dfferental Equatons PDEs) and eamples wth fnte dfference dscretzatons Parabolc PDEs Ellptc PDEs Hyperbolc PDEs Error Types and Dscretzaton
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12
REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13
CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationNONCENTRAL 7POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NONCENTRAL 7POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationME 501A Seminar in Engineering Analysis Page 1
umercal Solutons of oundaryvalue Problems n Os ovember 7, 7 umercal Solutons of oundary Value Problems n Os Larry aretto Mechancal ngneerng 5 Semnar n ngneerng nalyss ovember 7, 7 Outlne Revew stff equaton
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, MskolcEgyetemváros,
More informationChapter 3 Differentiation and Integration
MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton
More informationApplication of BSpline to Numerical Solution of a System of Singularly Perturbed Problems
Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of BSplne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More informationALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION
ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LIX, 013, f.1 DOI: 10.478/v101570100y ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION BY ION
More informationAGC Introduction
. Introducton AGC 3 The prmary controller response to a load/generaton mbalance results n generaton adjustment so as to mantan load/generaton balance. However, due to droop, t also results n a nonzero
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationarxiv: v1 [math.ho] 18 May 2008
Recurrence Formulas for Fbonacc Sums Adlson J. V. Brandão, João L. Martns 2 arxv:0805.2707v [math.ho] 8 May 2008 Abstract. In ths artcle we present a new recurrence formula for a fnte sum nvolvng the Fbonacc
More informationSpinrotation coupling of the angularly accelerated rigid body
Spnrotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 Emal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More informationC/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1
C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned
More informationPART 8. Partial Differential Equations PDEs
he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef
More informationNumerical Simulation of OneDimensional Wave Equation by NonPolynomial Quintic Spline
IOSR Journal of Matematcs (IOSRJM) eissn: 78578, pissn: 319765X. Volume 14, Issue 6 Ver. I (Nov  Dec 018), PP 630 www.osrournals.org Numercal Smulaton of OneDmensonal Wave Equaton by NonPolynomal
More informationResearch Article Green s Theorem for Sign Data
Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationBernoulli Numbers and Polynomials
Bernoull Numbers and Polynomals T. Muthukumar tmk@tk.ac.n 17 Jun 2014 The sum of frst n natural numbers 1, 2, 3,..., n s n n(n + 1 S 1 (n := m = = n2 2 2 + n 2. Ths formula can be derved by notng that
More informationFormulas for the Determinant
page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use
More informationThe internal structure of natural numbers and one method for the definition of large prime numbers
The nternal structure of natural numbers and one method for the defnton of large prme numbers Emmanul Manousos APM Insttute for the Advancement of Physcs and Mathematcs 3 Poulou str. 53 Athens Greece Abstract
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 008 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More informationFTCS Solution to the Heat Equation
FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use
More informationLeast squares cubic splines without Bsplines S.K. Lucas
Least squares cubc splnes wthout Bsplnes S.K. Lucas School of Mathematcs and Statstcs, Unversty of South Australa, Mawson Lakes SA 595 emal: stephen.lucas@unsa.edu.au Submtted to the Gazette of the Australan
More informationMaejo International Journal of Science and Technology
Maejo Int. J. Sc. Technol. ()  Full Paper Maejo Internatonal Journal of Scence and Technology ISSN  Avalable onlne at www.mjst.mju.ac.th Fourthorder method for sngularly perturbed sngular boundary value
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons Â and ˆb avalable. Then the best thng we can do s to solve Âˆx ˆb exactly whch
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationSingular Value Decomposition: Theory and Applications
Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real
More informationCALCULUS CLASSROOM CAPSULES
CALCULUS CLASSROOM CAPSULES SESSION S86 Dr. Sham Alfred Rartan Valley Communty College salfred@rartanval.edu 38th AMATYC Annual Conference Jacksonvlle, Florda November 8, 202 2 Calculus Classroom Capsules
More informationMAE140  Linear Circuits  Fall 10 Midterm, October 28
M140  Lnear rcuts  Fall 10 Mdterm, October 28 nstructons () Ths exam s open book. You may use whatever wrtten materals you choose, ncludng your class notes and textbook. You may use a hand calculator
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture  30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have
More informationECONOMICS 351*A MidTerm Exam  Fall Term 2000 Page 1 of 13 pages. QUEEN'S UNIVERSITY AT KINGSTON Department of Economics
ECOOMICS 35*A MdTerm Exam  Fall Term 000 Page of 3 pages QUEE'S UIVERSITY AT KIGSTO Department of Economcs ECOOMICS 35*  Secton A Introductory Econometrcs Fall Term 000 MIDTERM EAM ASWERS MG Abbott
More informationAdvanced Circuits Topics  Part 1 by Dr. Colton (Fall 2017)
Advanced rcuts Topcs  Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed
More information= z 20 z n. (k 20) + 4 z k = 4
Problem Set #7 solutons 7.2.. (a Fnd the coeffcent of z k n (z + z 5 + z 6 + z 7 + 5, k 20. We use the known seres expanson ( n+l ( z l l z n below: (z + z 5 + z 6 + z 7 + 5 (z 5 ( + z + z 2 + z + 5 5
More informationSolution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method
Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another wellknown teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222
More informationMA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials
MA 323 Geometrc Modellng Course Notes: Day 13 Bezer Curves & Bernsten Polynomals Davd L. Fnn Over the past few days, we have looked at de Casteljau s algorthm for generatng a polynomal curve, and we have
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationA New Refinement of Jacobi Method for Solution of Linear System Equations AX=b
Int J Contemp Math Scences, Vol 3, 28, no 17, 819827 A New Refnement of Jacob Method for Soluton of Lnear System Equatons AX=b F Naem Dafchah Department of Mathematcs, Faculty of Scences Unversty of Gulan,
More informationThe optimal delay of the second test is therefore approximately 210 hours earlier than =2.
THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 615086 provdes approxmaton formulas for the PF for smple
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 19329466 Vol. 10, Issue 1 (June 2015), pp. 223229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fntedmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationFundamental loopcurrent method using virtual voltage sources technique for special cases
Fundamental loopcurrent method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,
More informationTemperature. Chapter Heat Engine
Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the
More informationNumerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method
Appled Mathematcs, 6, 7, 54 Publshed Onlne Jul 6 n ScRes. http://www.scrp.org/journal/am http://.do.org/.436/am.6.77 umercal Solutons of a Generalzed th Order Boundar Value Problems Usng Power Seres Approxmaton
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00
ONE IMENSIONAL TRIANGULAR FIN EXPERIMENT Techncal Advsor: r..c. Look, Jr. Verson: /3/ 7. GENERAL OJECTIVES a) To understand a onedmensonal epermental appromaton. b) To understand the art of epermental
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationRestricted divisor sums
ACTA ARITHMETICA 02 2002) Restrcted dvsor sums by Kevn A Broughan Hamlton) Introducton There s a body of work n the lterature on varous restrcted sums of the number of dvsors of an nteger functon ncludng
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationSupplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More informationSolution for singularly perturbed problems via cubic spline in tension
ISSN 76769 England UK Journal of Informaton and Computng Scence Vol. No. 06 pp.669 Soluton for sngularly perturbed problems va cubc splne n tenson K. Aruna A. S. V. Rav Kant Flud Dynamcs Dvson Scool
More informationFinite Differences, Interpolation, and Numerical Differentiation
4 Fnte Dfferences, Interpolaton, and Numercal Dfferentaton 4. INTRODUCTION Lnear nterpolaton s dscussed n the precedng chapter as a method for fndng a partcular root of a polynomal, or, transcendental
More informationComplex Numbers Alpha, Round 1 Test #123
Complex Numbers Alpha, Round Test #3. Wrte your 6dgt ID# n the I.D. NUMBER grd, leftjustfed, and bubble. Check that each column has only one number darkened.. In the EXAM NO. grd, wrte the 3dgt Test
More information1. Estimation, Approximation and Errors Percentages Polynomials and Formulas Identities and Factorization 52
ontents ommonly Used Formulas. Estmaton, pproxmaton and Errors. Percentages. Polynomals and Formulas 8. Identtes and Factorzaton. Equatons and Inequaltes 66 6. Rate and Rato 8 7. Laws of Integral Indces
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationAPPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS
Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent
More informationImportant Instructions to the Examiners:
Summer 0 Examnaton Subject & Code: asc Maths (70) Model Answer Page No: / Important Instructons to the Examners: ) The Answers should be examned by key words and not as wordtoword as gven n the model
More informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationIntegrals and Invariants of EulerLagrange Equations
Lecture 16 Integrals and Invarants of EulerLagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationIrregular vibrations in multimass discretecontinuous systems torsionally deformed
(2) 4 48 Irregular vbratons n multmass dscretecontnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscretecontnuous systems consstng of an arbtrary number rgd bodes connected
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More information18.1 Introduction and Recap
CS787: Advanced Algorthms Scrbe: Pryananda Shenoy and Shjn Kong Lecturer: Shuch Chawla Topc: Streamng Algorthmscontnued) Date: 0/26/2007 We contnue talng about streamng algorthms n ths lecture, ncludng
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 63 ISSN: 97537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationNEW GENERAL FORMULAS FOR PYTHAGOREAN TRIPLES
NEW GENERAL FORMULAS FOR PYTHAGOREAN TRIPLES EDUARDO CALOOY ROQUE Abstract. Ths paper shows the General Formulas for Pythagorean trples that were derved from the dfferences of the sdes of a rght trangle.
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 97457X (Onlne), ISSN 9745718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 1314 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationGrid Generation around a Cylinder by Complex Potential Functions
Research Journal of Appled Scences, Engneerng and Technolog 4(): 53535, 0 ISSN: 0407467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around
More informationThe exponential map of GL(N)
The exponental map of GLN arxv:hepth/9604049v 9 Apr 996 Alexander Laufer Department of physcs Unversty of Konstanz P.O. 5560 M 678 78434 KONSTANZ Aprl 9, 996 Abstract A fnte expanson of the exponental
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationMAE140  Linear Circuits  Winter 16 Midterm, February 5
Instructons ME140  Lnear Crcuts  Wnter 16 Mdterm, February 5 () Ths exam s open book. You may use whatever wrtten materals you choose, ncludng your class notes and textbook. You may use a hand calculator
More informationModelli Clamfim Equazione del Calore Lezione ottobre 2014
CLAMFIM Bologna Modell 1 @ Clamfm Equazone del Calore Lezone 17 15 ottobre 2014 professor Danele Rtell danele.rtell@unbo.t 1/24? Convoluton The convoluton of two functons g(t) and f(t) s the functon (g
More informationConsistency & Convergence
/9/007 CHE 374 Computatonal Methods n Engneerng Ordnary Dfferental Equatons Consstency, Convergence, Stablty, Stffness and Adaptve and Implct Methods ODE s n MATLAB, etc Consstency & Convergence Consstency
More informationProcedia Computer Science
Avalable onlne at www.scencedrect.com Proceda Proceda Computer Computer Scence Scence 1 (01) 00 (009) 589 597 000 000 Proceda Computer Scence www.elsever.com/locate/proceda Internatonal Conference on Computatonal
More informationChapter 4 The Wave Equation
Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a secondorder lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as
More informationA MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS
Journal of Mathematcs and Statstcs 9 (1): 48, 1 ISSN 1549644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF
More informationNumerical Solution of OneDimensional Heat and Wave Equation by NonPolynomial Quintic Spline
Internatonal Journal of Mathematcal Modellng & Computatons Vol. 05, No. 04, Fall 2015, 291305 Numercal Soluton of OneDmensonal Heat and Wave Equaton by NonPolynomal Quntc Splne J. Rashdna a, and M.
More informationA NUMERICAL COMPARISON OF LANGRANGE AND KANE S METHODS OF AN ARM SEGMENT
Internatonal Conference Mathematcal and Computatonal ology 0 Internatonal Journal of Modern Physcs: Conference Seres Vol. 9 0 68 75 World Scentfc Publshng Company DOI: 0.4/S009450059 A NUMERICAL COMPARISON
More information