Stepsize Control for Cubic Spline Interpolation
|
|
- Kristian Blake
- 5 years ago
- Views:
Transcription
1 Nonme mnuscript No. (will be inserted by te editor) Stepsize Control for Cubic Spline Interpoltion Dvod Kojste Slkuye Received: dte / Accepted: dte Abstrct Numericl lgoritms wic include stepsize re ffected by globl error, wic consists of bot trunction error nd round-off error. In tese lgoritms s te stepsize decreses, te trunction error lso decreses, but te round-off error my increse. Te problem is now to find te stepsize wic minimizes te globl error. In tis pper, by using te CESTAC metod nd te CADNA librry procedure is proposed to control te stepsize for determining interpolting cubic spline functions. Some numericl experiments re given to sow te effectiveness of te proposed procedure. Keywords Cubic spline stepsize CESTAC metod DSA CADNA librry. Mtemtics Subject Clssifiction (2000) 65D07 65G50. 1 Introduction Finite precision computtions my ffect te stbility of lgoritms nd te ccurcy of computed solutions. In fct, te numericl result provided by n lgoritm is ffected by globl error, wic consists of bot trunction error nd round-off error. Computtion of te solution of mny problems in scientific computing involves stepsize. As decreses, te trunction error lso decreses, but te round-off error in te metod my increse. Here, te problem is to find te optiml stepsize opt wic minimizes te globl error. In generl, it is difficult to estimte te optiml stepsize in n lgoritm. In [4], Cesneux nd Jézéquel sowed tt by using te CESTAC (Contrôle et Estimtion Stocstique des Arrondis de Clcul) metod [20, 19] nd te CADNA (Control of Accurcy nd Debugging for Numericl Appliction) librry [2, 3, 17], one cn estimte te optiml stepsize for te numericl computtion of integrls using te trpezoidl nd Simpson s rules. Ten, Abbsbndy nd Argi [1] developed tis metod to generl closed Newton-Cotes integrtion rules. Te development of te metod to multiple integrls cn be found in [12]. In [14], Slkuye et l. proposed procedure wit stepsize control for solving n one-dimensionl initil vlue problems. In tis pper, procedure is presented to control te stepsize for te cubic spline interpoltions wit eqully spced knots. Tis pper is orgnized s follows. In section 2, some new results on te cubic splines re given. Section 3 is devoted to min results for definition of te procedure. A brief review of te CESTAC D. K. Slkuye Fculty of Mtemticl Sciences, University of Guiln, P.O. Box 1914, Rst, Irn Tel.: Fx: E-mil: kojste@guiln.c.ir, slkuye@gmil.com
2 2 Dvod Kojste Slkuye metod, te DSA nd te CADNA librry is given in Section 4. In section 5, some numericl experiments re given. Section 6 is devoted to some concluding remrks. 2 Some results on te cubic splines Let f be rel-vlued function defined on [, b] nd T n = { = x 0 < x 1 < x 2 <... < x n 1 < x n = b} be prtition of [, b] wit eqully spced knots nd stepsize. Let S n be te clmped interpolting cubic spline function on te prtition T n stisfying S n (x i ) = f(x i ) = f i, i = 0, 1,..., n, nd S n(x 0 ) = f (x 0 ), S n(x n ) = f (x n ). (1) It cn be seen tt (see [15,9]) were S n (x) = α j + β j (x x j ) + γ j (x x j ) 2 + δ j (x x j ) 3, x [x j, x j+1 ], (2) α j = f j, in wic M j s stisfy suc tt β j = f j+1 f j + 2M j + M j+1, γ j = M j 6 2, δ j = M j+1 M j M 0 d 0 M 1 d 1 M 2 = d 2, (3).. M n d n, d j = 6 2 (f j+1 2f j + f j 1 ), j = 1, 2,..., n 1, d 0 = 6 (f 1 f 0 f (x 0 )), d n = 6 (f (x n ) f n f n 1 ). It is well-known tt, if f C 2 [, b], ten S n(x) tends to f (x) s tends to zero for ll x [, b], nd E n := (f (x) n(x)) 2 dx = f (x) 2 dx n(x) 2 dx 0, (4) were equlity olds if nd only if f = S n (see [8,13,15]). Terefore, te evlution of E n would be very useful, becuse it my elp us to stte stopping criterion for coosing suitble n. Hence, we re going to investigte ow tis difference depends on te stepsize. We begin wit te following teorem. Teorem 1 Let S n be te clmped cubic spline interpolting f. Ten Proof From (2), one cn see tt n(x) 2 dx = n 1 3 (M n(x) = M j x j+1 x j=1 n 1 Mj 2 + Mn 2 + j=0 M j M j+1 ). x x j + M j+1, x [x j, x j+1 ]. By tis reltion it is esy to see tt xj+1 x j n(x) 2 dx = 3 (M 2 j + M j M j+1 + M 2 j+1).
3 Stepsize Control for Cubic Spline Interpoltion 3 Terefore n(x) 2 dx = n 1 xj+1 i=0 x j n(x) 2 dx = n 1 3 (M n Mj 2 + Mn 2 + M j M j+1 ). j=1 j=0 Tis completes te proof. Teorem 2 Let f be sufficiently smoot function nd S n be te clmped cubic spline interpolting f. Ten f (x) 2 dx = S n(x) 2 dx + b 720 (f (4) (c)) O( 6 ), (5) for some c in [, b]. Proof It is esy to see tt (see [15, pge 140]) (f (x) n(x)) 2 dx = (f(x) S n (x))f (4) (x)dx. (6) For te ske of simplicity, let g(x) = (f(x) S(x))f (4) (x). By te Euler-Mclurin formul [15, pge 156] we ve xi+1 x i g(x)dx = 2 (g(x i) + g(x i+1 )) 2 12 (g (x i+1 ) g (x i )) (g(3) (x i+1 ) g(3) (x + i )) (g(5) (x i+1 ) g(5) (x + i )) + O(8 ). (7) Obviously, g(x i ) = g(x i+1 ) = 0. We now recll tt (see [16, pge 121] nd [9, pge 487]) f (x j ) = S n(x j ) f (5) (x j ) + O( 6 ), j = i, i + 1, (8) f (x j ) = n(x j ) f (4) (x j ) + O( 4 ), j = i, i + 1, (9) f (3) (x i ) = S n (3) (x + i ) 2 f (4) (x i ) 2 12 f (5) (x i ) + O( 4 ), (10) f (3) (x i+1 ) = S n (3) (x i+1 ) + 2 f (4) (x i+1 ) 2 12 f (5) (x i+1 ) + O( 4 ). (11) Differentiting g nd using Eq. (8) we get g (x i+1 ) g (x i ) = (f (5) (x i+1 )f (4) (x i+1 ) f (5) (x i )f (4) (x i )) 4 + O( 6 ). (12) On te oter nd, by tree times differentiting of g nd using Eqs. (8), (9) nd (10), we obtin g (3) (x + i ) = (f (3) S (3) n )(x + i )f (4) (x i ) + 3(f n)(x + i )f (5) (x i ) +3(f S n)(x + i )f (6) (x i ) + (f S n )(x + i )f (7) (x i ) = 2 (f (4) (x i )) f (5) (x i )f (4) (x i ) + O( 4 ). (13) In te sme mnner, from Eqs. (8), (9) nd (11), we deduce g (3) (x i+1 ) = 2 (f (4) (x i+1 )) f (5) (x i+1 )f (4) (x i+1 ) + O( 4 ). (14) Terefore, from Eqs. (13) nd (14), we see tt
4 4 Dvod Kojste Slkuye g (3) (x i+1 ) g(3) (x + i ) = 2 ((f (4) (x i )) 2 + (f (4) (x i+1 )) 2 ) (f (5) (x i+1 )f (4) (x i+1 ) f (5) (x i )f (4) (x i )) + O( 4 ) = (f (4) (ξ i )) (f (5) (x i+1 )f (4) (x i+1 ) f (5) (x i )f (4) (x i )) + O( 4 ), for some ξ i [x i, x i+1 ]. It is mentioned tt te ltter equlity s been obtined by te intermedite vlue teorem. Substituting (12) nd (15) in Eq.(11), yields xi+1 x i g(x)dx = (f (4) (ξ i )) (f (5) (x i+1 )f (4) (x i+1 ) f (5) (x i )f (4) (x i )) (g(5) (x i+1 ) g(5) (x + i )) + O(8 ). = (f (4) (ξ i )) (f (5) (x i+1 )f (4) (x i+1 ) f (5) (x i )f (4) (x i )) g(6) (η i ) + O( 8 ), for some η i (x i, x i+1 ). Here, it is noted tt te ltter equlity is obtined by te men vlue teorem. Tis togeter wit (6) results in (f (x) n(x)) 2 dx = n 1 xi+1 i=0 x i (f(x) S n (x))f (4) (x)dx = 1 n (f (4) (ξ i )) (f (5) (b)f (4) (b) f (5) ()f (4) ()) 6 + O( 6 ) i=0 = b 720 (f (4) (c)) O( 6 ), for some c in [, b] nd te proof is completed. Tis teorem sows tt S n(x) 2 dx pproximtes f (x) 2 dx nd te pproximtion error is of order 4. In te next section, we present teorem wic elps us to propose procedure to control te stepsize for determining interpolting cubic spline functions. 3 Definition of te procedure We first recll te following definition [4, 10]. Definition 1 Let r nd s be two rel numbers. Te number of exct significnt digits tt re common to r nd s cn be defined in (, + ) by 1. for r s, 2. C r,r = +. C r,s = log 10 r + s 2(r s).
5 Stepsize Control for Cubic Spline Interpoltion 5 Teorem 3 (Jézéquel [10]). Let I() be n pproximtion of order p of n exct vlue, i.e., I I() = C p + O( q ) wit 1 p < q nd C R. If I m is n pproximtion computed wit te stepsize 0 /2 m, ten C Im,I m+1 = C Im,I + log 10 ( 2 p 1 ) + O(2m(p q) ). Teorem 4 Let f be sufficiently smoot function nd S m be te clmped cubic spline interpolting f wit stepsize (b )/2 m. Moreover, let Ten Proof Let I = In tis cse, by Teorem 2 we ve 2 p f (x) 2 dx, nd I m = m(x) 2 dx. C Im,I m+1 = C Im,I + log 10 ( ) + O( 1 ). (15) 4m I() = n(x) 2 dx, = b n. I I() = C 4 + O( 6 ), were C = b 720 (f (4) (c)) 2. Now, letting p = 4, q = 6 nd 0 = b, from Teorem 3 te desired reltion is obtined. Tis teorem sows tt, if te convergence zone is reced, i.e., if te term O( 1 4 ) becomes m negligible, ten te significnt digits common to two successive pproximtions I m nd I m+1 re lso in common wit I, up to one. Note tt te term log 10 ( ) corresponds t most to one bit. Tese teoreticl results ve been obtined by tking into ccount only te trunction error on two successive pproximtions I m nd I m+1. However computed results re lso ffected by roundoff error propgtion. Next, we describe ow round-off errors cn be estimted wit probbilistic pproc in order to determine te exct significnt digits of ny computed result. 4 Te CESTAC metod, te DSA nd te CADNA librry Te Discrete Stocstic Aritmetic (DSA) [18] is probbilistic pproc for nlyzing te round-off error propgtion wic is bsed on te syncronous implementtion of te CESTAC metod nd te stocstic order reltions. In te CESTAC metod te floting-point ritmetic is replced by rndom ritmetic [5,7,18,20]. In fct, computer progrm is run N times in prllel nd in tis wy new ritmetic clled stocstic ritmetic is defined (for more detils see [6]). Te min ide is to use rndom rounding to produce severl smples of ec result of ny ritmetic opertion nd te number of common digits in tese smples estimtes te number of exctly significnt digits in te floting-point result. From tis point of view, Vignes [21] introduced te concept of computtionl zero wic is denoted A computtionl zero is computed result wic s no significnt digit or wic is te mtemticl zero. Te CADNA softwre [2,3,17] is librry wic implements utomticlly te DSA in ny code written in Fortrn. Using te CADNA librry, ec stndrd FP types ve teir corresponding stocstic types. Every intrinsic function nd opertor re overloded for tose types. Wen stocstic vrible is printed, only its significnt digits re displyed to point out its ccurcy. If number s no significnt digit (i.e., computtionl zero), te is displyed. In te next section, we describe ow te use of te CADNA librry llows us to control te stepsize for te clmped interpolting cubic spline functions.
6 6 Dvod Kojste Slkuye Tble 1 Results for Exmple 1. m I m I m I m E E E E E E E E E E E E E E E E E E E E E E E E E E E Numericl exmples According to te previous section nd Teorem 4, by using te CESTAC metod nd te CADNA librry, we propose to use I m I m 1 s stopping criterion. Te use of tis stopping criterion llows us to estimte te optiml stepsize s soon s te difference between I m nd I m 1 is equl to te computtionl zero. In tis cse, by te teoreticl results presented in Section 3, te common significnt digits between I m nd I m 1 re lso te common significnt digits between I m nd I, up to one digit. By te proposed metod, te best pproximtion of I provided by te computer is cosen. In tis cse, S m is te best pproximtion of f provided by te computer. It is necessry to mention tt in ll te presented exmples I m s ve been computed by Teorem 1. Let us now present te exmples nd te results wic we obtined by te Fortrn code of te clmped interpolting cubic spline functions combined wit te CADNA librry, version BETA [17]. All te numericl experiments were computed in double precision. Te dt st function of te CADNA librry s been used to tke into ccount te ssignment error of some dt of te exmples wic re not integer numbers nd cnnot be exctly coded on computer. Exmple 1 Consider te Runge s function f(x) = x 2 on [ 3, 3]. Numericl results re given in Tble 1. As we observe I 15 I 14 Terefore te optiml stepsize is estimted s = 6/2 15 = E 004. To sow te effectiveness of te proposed metod te vlues of S 15 (x) nd f(x) for midpoint of ten subintervls (denoted by x j ) of [ 3, 3], for te obtined optiml stepsize, re given in Tble 2. In column 1 of tis tble te index of te subintervl is given. Here we mention tt te vlues of f re computed by MATLAB code wit 20 digits of ccurcy nd re presented by 15 digits wit rounding. As we see, ll of te digits of S 15 ( x j ), up to one, coincide wit tt of te f( x j ) for ll j. Exmple 2 In tis exmple, we consider f(x) = 3 + x + 5 sin x on [0, 2]. Numericl results re given in Tble 3. Here, we ve I 10 I 9 Terefore te optiml stepsize is estimted s = 2/2 10 = E 003. Similr to Exmple 1, te vlues of S 10 (x) nd f(x) for midpoint of ten subintervls of [0, 2] for te obtined optiml stepsize re given in Tble 4. Oter ssumptions re s before. As we see, ll of te digits of S 10 ( x j ), up to one, coincide wit tt of te f( x j ). Exmple 3 Tis exmple is devoted to te function f(x) = tn x on te intervl [0, 1]. Numericl results re given in Tble 5. In tis exmple I 11 I 10 Terefore te optiml stepsize is
7 Stepsize Control for Cubic Spline Interpoltion 7 Tble 2 Comprison between exct nd numericl results for Exmple 1 t some points. j x j = (x j + x j+1 )/2 f( x j ) S 15 ( x j ) E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E-002 Tble 3 Results for Exmple 2. m I m I m I m E E E E E E E E E E E E E E E E E Tble 4 Comprison between exct nd numericl results for Exmple 2 t some points. j x j = (x j + x j+1 )/2 f( x j ) S 15 ( x j ) E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E+001 estimted s = 1/2 11 = E 004. Te vlues of S 11 (x) nd f(x) for some points re given in Tble 6. Oter ssumptions re similr to previous exmples. As we see, ll of te digits of S 11 ( x j ), up to one, coincide wit tt of f( x j ). 6 Conclusion In tis pper, teorem s been stted to provide stopping criterion to control te stepsize in te clmped interpolting cubic spline functions. We observed tt te use of te CESTAC metod nd te CADNA librry llows us to estimte te optiml stepsize. Numericl exmples sow tt te proposed metod is effective. It is esy to sow tt te results of tis pper old for oter type of splines. Acknowledgements Te utor would like to tnk to te nonymous referee for is/er comments. Tis work is prtilly supported by University of Guiln.
8 8 Dvod Kojste Slkuye Tble 5 Results for Exmple 3. m I m I m I m E E E E E E E E E E E E E E E E E E E Tble 6 Comprison between exct nd numericl results for Exmple 3 t some points. j x j = (x j + x j+1 )/2 f( x j ) S 15 ( x j ) E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E+001 References 1. S. Abbsbndy nd M. A. Friborzi Argi, A stocstic sceme for solving definite integrls, Appl. Numer. Mt., 55 (2005) J. M. Cesneux, L ritmétique stocstique et le logiciel CADNA, Hbilittion à diriger des recerces, Université Pierre et Mrie Curie, Pris, J. M. Cesneux, CADNA: An ADA tool for round-off errors nlysis nd for numericl debugging, In ADA in Aerospce, Brcelone, Spin, December J. M. Cesneux nd F. Jézéquel, Dynmicl Control of Computtions Using te Trpezodil nd Simpson s rules, J. Universl Comput. Sci., 4 (1998) J. M. Cesneux, Study of te computing ccurcy by using probbilistic pproc, in: C. Ullric (Ed.), Contribution to Computer Aritmetic nd Self-Vlidting Numericl Metods, IMACS, New Brunswick, NJ, J. M. Cesneux, Stocstic ritmetic properties, IMACS Comput. Appl. Mt. (1992) J. M. Cesneux nd J. Vignes, Les fondements de l ritmétique stocstique, C. R. Acd. Sci. Pris, Sér. I Mt., 315 (1992) C. de Boor, A prcticl guide to splines, Revised edition, Springer-Verlg F. B. Hildebrnd, Introduction to numericl nlysis, second edition, Dover Publictions, Inc., New York, F. Jézéquel, A dynmicl strtegy for pproximtion metods, C. R. Mécnique, 334 (2006) F. Jézéquel nd J. M. Cesneux, CADNA: librry for estimting round-off error propgtion, Computer Pysics Communictions, 178 (2008) F. Jézéquel, Contrôle dynmique de métodes d pproximtion, Hbilittion á diriger des recerces, Université Pierre et Mrie Curie, Pris, Februry 2005, Avilble t: ttp://www-np.lip6.fr/ jezequel/hdr.tml. 13. A. Qurteroni, R. Scco nd F. Sleri, Numericl mtemtics, Springer-Verlg, New York, D. K. Slkuye, F. Toutounin nd H. S. Yzdi, A procedure wit stepsize control for solving n one-dimensionl IVPs, Mtemtics nd Computers in Simultion, 79 (2008) J. Stoer nd R. Bulirsc, Introduction to Numericl Anlysis, Springer-Verlg, Tird Edition, S. S. Sstry, Introductory metods of numericl nlysis, Fourt edition, Prentice-Hll of Indi, Te CADNA librry, URL ddress: ttp:// 18. J. Vignes, Discrete stocstic ritmetic for vlidting results of numericl softwre, Numericl Algoritms, 37 (2004) J. Vignes, New metods for evluting te vlidity of te results of mtemticl computtions, Mt. Comp. Simul., 20 (1978) J. Vignes, A stocstic ritmetic for relible scientific computtion, Mt. Comp. Simul., 35 (1993) J. Vignes, Zéro mtémtique et zéro informtique, C. R. Acd. Sci. Pris Sér. I Mt., 303 (1986)
Lecture Note 4: Numerical differentiation and integration. Xiaoqun Zhang Shanghai Jiao Tong University
Lecture Note 4: Numericl differentition nd integrtion Xioqun Zng Sngi Jio Tong University Lst updted: November, 0 Numericl Anlysis. Numericl differentition.. Introduction Find n pproximtion of f (x 0 ),
More informationKey words. Numerical quadrature, piecewise polynomial, convergence rate, trapezoidal rule, midpoint rule, Simpson s rule, spectral accuracy.
O SPECTRA ACCURACY OF QUADRATURE FORMUAE BASED O PIECEWISE POYOMIA ITERPOATIO A KURGAOV AD S TSYKOV Abstrct It is well-known tt te trpezoidl rule, wile being only second-order ccurte in generl, improves
More informationMath Week 5 concepts and homework, due Friday February 10
Mt 2280-00 Week 5 concepts nd omework, due Fridy Februry 0 Recll tt ll problems re good for seeing if you cn work wit te underlying concepts; tt te underlined problems re to be nded in; nd tt te Fridy
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationRGMIA Research Report Collection, Vol. 1, No. 1, SOME OSTROWSKI TYPE INEQUALITIES FOR N-TIME DIFFERENTIA
ttp//sci.vut.edu.u/rgmi/reports.tml SOME OSTROWSKI TYPE INEQUALITIES FOR N-TIME DIFFERENTIABLE MAPPINGS AND APPLICATIONS P. CERONE, S.S. DRAGOMIR AND J. ROUMELIOTIS Astrct. Some generliztions of te Ostrowski
More informationUniversity of Houston, Department of Mathematics Numerical Analysis II
University of Houston, Deprtment of Mtemtics Numericl Anlysis II 6 Glerkin metod, finite differences nd colloction 6.1 Glerkin metod Consider sclr 2nd order ordinry differentil eqution in selfdjoint form
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More information1 The Riemann Integral
The Riemnn Integrl. An exmple leding to the notion of integrl (res) We know how to find (i.e. define) the re of rectngle (bse height), tringle ( (sum of res of tringles). But how do we find/define n re
More informationNumerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
More information63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1
3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =
More informationTopic 6b Finite Difference Approximations
/8/8 Course Instructor Dr. Rymond C. Rump Oice: A 7 Pone: (95) 747 6958 E Mil: rcrump@utep.edu Topic 6b Finite Dierence Approximtions EE 486/5 Computtionl Metods in EE Outline Wt re inite dierence pproximtions?
More informationCMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature
CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy
More informationDOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES
DOIG PHYSICS WITH MATLAB MATHEMATICAL ROUTIES COMPUTATIO OF OE-DIMESIOAL ITEGRALS In Cooper School of Physics, University of Sydney in.cooper@sydney.edu.u DOWLOAD DIRECTORY FOR MATLAB SCRIPTS mth_integrtion_1d.m
More informationExam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1
Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution
More informationChapter 2. Numerical Integration also called quadrature. 2.2 Trapezoidal Rule. 2.1 A basic principle Extending the Trapezoidal Rule DRAWINGS
S Cpter Numericl Integrtion lso clled qudrture Te gol of numericl integrtion is to pproximte numericlly. f(x)dx Tis is useful for difficult integrls like sin(x) ; sin(x ); x + x 4 Or worse still for multiple-dimensionl
More informationNUMERICAL INTEGRATION
NUMERICAL INTEGRATION How do we evlute I = f (x) dx By the fundmentl theorem of clculus, if F (x) is n ntiderivtive of f (x), then I = f (x) dx = F (x) b = F (b) F () However, in prctice most integrls
More informationIII. Lecture on Numerical Integration. File faclib/dattab/lecture-notes/numerical-inter03.tex /by EC, 3/14/2008 at 15:11, version 9
III Lecture on Numericl Integrtion File fclib/dttb/lecture-notes/numerical-inter03.tex /by EC, 3/14/008 t 15:11, version 9 1 Sttement of the Numericl Integrtion Problem In this lecture we consider the
More informationMath& 152 Section Integration by Parts
Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible
More information3.4 Numerical integration
3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,
More informationDerivatives, Finding Roots, Interpolating and Integration
Derivtives, Finding Roots, Interpolting nd Integrtion Doing Numericl Derivtives Altoug tking derivtives is very esy opertion for ll of you, tere re situtions were it must e done numericlly. For exmple,
More informationThe Riemann Integral
Deprtment of Mthemtics King Sud University 2017-2018 Tble of contents 1 Anti-derivtive Function nd Indefinite Integrls 2 3 4 5 Indefinite Integrls & Anti-derivtive Function Definition Let f : I R be function
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More information1 The Lagrange interpolation formula
Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x
More informationGENERALIZATIONS OF WEIGHTED TRAPEZOIDAL INEQUALITY FOR MONOTONIC MAPPINGS AND ITS APPLICATIONS. (b a)3 [f(a) + f(b)] f x (a,b)
GENERALIZATIONS OF WEIGHTED TRAPEZOIDAL INEQUALITY FOR MONOTONIC MAPPINGS AND ITS APPLICATIONS KUEI-LIN TSENG, GOU-SHENG YANG, AND SEVER S. DRAGOMIR Abstrct. In this pper, we estblish some generliztions
More informationAN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION
Applied Mthemtics E-Notes, 5(005), 53-60 c ISSN 1607-510 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ AN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION
More informationNumerical Integration
Chpter 1 Numericl Integrtion Numericl differentition methods compute pproximtions to the derivtive of function from known vlues of the function. Numericl integrtion uses the sme informtion to compute numericl
More informationLecture 12: Numerical Quadrature
Lecture 12: Numericl Qudrture J.K. Ryn@tudelft.nl WI3097TU Delft Institute of Applied Mthemtics Delft University of Technology 5 December 2012 () Numericl Qudrture 5 December 2012 1 / 46 Outline 1 Review
More informationQUADRATURE is an old-fashioned word that refers to
World Acdemy of Science Engineering nd Technology Interntionl Journl of Mthemticl nd Computtionl Sciences Vol:5 No:7 011 A New Qudrture Rule Derived from Spline Interpoltion with Error Anlysis Hdi Tghvfrd
More informationFuzzy transform to approximate solution of boundary value problems via optimal coefficients
217 Interntionl Conference on Hig Performnce Computing & Simultion Fuzzy trnsform to pproximte solution of boundry vlue problems vi optiml coefficients Zr Alijni Mtemtics nd Sttistics Dept University of
More informationCzechoslovak Mathematical Journal, 55 (130) (2005), , Abbotsford. 1. Introduction
Czechoslovk Mthemticl Journl, 55 (130) (2005), 933 940 ESTIMATES OF THE REMAINDER IN TAYLOR S THEOREM USING THE HENSTOCK-KURZWEIL INTEGRAL, Abbotsford (Received Jnury 22, 2003) Abstrct. When rel-vlued
More informationLecture 14 Numerical integration: advanced topics
Lecture 14 Numericl integrtion: dvnced topics Weinn E 1,2 nd Tiejun Li 2 1 Deprtment of Mthemtics, Princeton University, weinn@princeton.edu 2 School of Mthemticl Sciences, Peking University, tieli@pku.edu.cn
More informationS. S. Dragomir. 2, we have the inequality. b a
Bull Koren Mth Soc 005 No pp 3 30 SOME COMPANIONS OF OSTROWSKI S INEQUALITY FOR ABSOLUTELY CONTINUOUS FUNCTIONS AND APPLICATIONS S S Drgomir Abstrct Compnions of Ostrowski s integrl ineulity for bsolutely
More informationNumerical Analysis. 10th ed. R L Burden, J D Faires, and A M Burden
Numericl Anlysis 10th ed R L Burden, J D Fires, nd A M Burden Bemer Presenttion Slides Prepred by Dr. Annette M. Burden Youngstown Stte University July 9, 2015 Chpter 4.1: Numericl Differentition 1 Three-Point
More informationArithmetic Mean Derivative Based Midpoint Rule
Applied Mthemticl Sciences, Vol. 1, 018, no. 13, 65-633 HIKARI Ltd www.m-hikri.com https://doi.org/10.1988/ms.018.858 Arithmetic Men Derivtive Bsed Midpoint Rule Rike Mrjulis 1, M. Imrn, Symsudhuh Numericl
More informationKeywords : Generalized Ostrowski s inequality, generalized midpoint inequality, Taylor s formula.
Generliztions of the Ostrowski s inequlity K. S. Anstsiou Aristides I. Kechriniotis B. A. Kotsos Technologicl Eductionl Institute T.E.I.) of Lmi 3rd Km. O.N.R. Lmi-Athens Lmi 3500 Greece Abstrct Using
More informationAdvanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004
Advnced Clculus: MATH 410 Notes on Integrls nd Integrbility Professor Dvid Levermore 17 October 2004 1. Definite Integrls In this section we revisit the definite integrl tht you were introduced to when
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationNumerical Integration - (4.3)
Numericl Itegrtio - (.). Te Degree of Accurcy of Qudrture Formul: Te degree of ccurcy of qudrture formul Qf is te lrgest positive iteger suc tt x k dx Qx k, k,,,...,. Exmple fxdx 9 f f,,. Fid te degree
More informationLecture 17. Integration: Gauss Quadrature. David Semeraro. University of Illinois at Urbana-Champaign. March 20, 2014
Lecture 17 Integrtion: Guss Qudrture Dvid Semerro University of Illinois t Urbn-Chmpign Mrch 0, 014 Dvid Semerro (NCSA) CS 57 Mrch 0, 014 1 / 9 Tody: Objectives identify the most widely used qudrture method
More informationSome estimates on the Hermite-Hadamard inequality through quasi-convex functions
Annls of University of Criov, Mth. Comp. Sci. Ser. Volume 3, 7, Pges 8 87 ISSN: 13-693 Some estimtes on the Hermite-Hdmrd inequlity through qusi-convex functions Dniel Alexndru Ion Abstrct. In this pper
More informationNumerical quadrature based on interpolating functions: A MATLAB implementation
SEMINAR REPORT Numericl qudrture bsed on interpolting functions: A MATLAB implementtion by Venkt Ayylsomyjul A seminr report submitted in prtil fulfillment for the degree of Mster of Science (M.Sc) in
More informationEuler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )
Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s
More informationNumerical integration
2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter
More informationS. S. Dragomir. 1. Introduction. In [1], Guessab and Schmeisser have proved among others, the following companion of Ostrowski s inequality:
FACTA UNIVERSITATIS NIŠ) Ser Mth Inform 9 00) 6 SOME COMPANIONS OF OSTROWSKI S INEQUALITY FOR ABSOLUTELY CONTINUOUS FUNCTIONS AND APPLICATIONS S S Drgomir Dedicted to Prof G Mstroinni for his 65th birthdy
More informationThe Riemann-Lebesgue Lemma
Physics 215 Winter 218 The Riemnn-Lebesgue Lemm The Riemnn Lebesgue Lemm is one of the most importnt results of Fourier nlysis nd symptotic nlysis. It hs mny physics pplictions, especilly in studies of
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationSections 5.2: The Definite Integral
Sections 5.2: The Definite Integrl In this section we shll formlize the ides from the lst section to functions in generl. We strt with forml definition.. The Definite Integrl Definition.. Suppose f(x)
More informationLECTURE 19. Numerical Integration. Z b. is generally thought of as representing the area under the graph of fèxè between the points x = a and
LECTURE 9 Numericl Integrtion Recll from Clculus I tht denite integrl is generlly thought of s representing the re under the grph of fèxè between the points x = nd x = b, even though this is ctully only
More information情報処理学会研究報告 IPSJ SIG Technical Report Vol.2012-ARC-202 No.20 Vol.2012-HPC-137 No /12/14 Euler-Maclaurin 1 Euler-Maclaurin 1) Taylor Taylor Euler
Euler-Mclurin Euler-Mclurin Tylor Tylor Euler-Mclurin Tylor Performnce of te Numericl Integrtion by Euler-Mclurin summtion formul Hirosi Hirym Euler-Mclurin summtion formul cn be regrded s error evlution
More informationDefinite integral. Mathematics FRDIS MENDELU
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the
More informationMidpoint Approximation
Midpoint Approximtion Sometimes, we need to pproximte n integrl of the form R b f (x)dx nd we cnnot find n ntiderivtive in order to evlute the integrl. Also we my need to evlute R b f (x)dx where we do
More information6.5 Numerical Approximations of Definite Integrals
Arknss Tech University MATH 94: Clculus II Dr. Mrcel B. Finn 6.5 Numericl Approximtions of Definite Integrls Sometimes the integrl of function cnnot be expressed with elementry functions, i.e., polynomil,
More informationUndergraduate Research
Undergrdute Reserch A Trigonometric Simpson s Rule By Ctherine Cusimno Kirby nd Sony Stnley Biogrphicl Sketch Ctherine Cusimno Kirby is the dughter of Donn nd Sm Cusimno. Originlly from Vestvi Hills, Albm,
More informationOrthogonal Polynomials
Mth 4401 Gussin Qudrture Pge 1 Orthogonl Polynomils Orthogonl polynomils rise from series solutions to differentil equtions, lthough they cn be rrived t in vriety of different mnners. Orthogonl polynomils
More informationFUNCTIONS OF α-slow INCREASE
Bulletin of Mthemticl Anlysis nd Applictions ISSN: 1821-1291, URL: http://www.bmth.org Volume 4 Issue 1 (2012), Pges 226-230. FUNCTIONS OF α-slow INCREASE (COMMUNICATED BY HÜSEYIN BOR) YILUN SHANG Abstrct.
More informationDefinite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová (Mendel University) Definite integrl MENDELU / Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function
More informationUnit #9 : Definite Integral Properties; Fundamental Theorem of Calculus
Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl
More information12 Basic Integration in R
14.102, Mt for Economists Fll 2004 Lecture Notes, 10/14/2004 Tese notes re primrily bsed on tose written by Andrei Bremzen for 14.102 in 2002/3, nd by Mrek Pyci for te MIT Mt Cmp in 2003/4. I ve mde only
More informationZ b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but...
Chpter 7 Numericl Methods 7. Introduction In mny cses the integrl f(x)dx cn be found by finding function F (x) such tht F 0 (x) =f(x), nd using f(x)dx = F (b) F () which is known s the nlyticl (exct) solution.
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationProperties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives
Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn
More informationf(x)dx . Show that there 1, 0 < x 1 does not exist a differentiable function g : [ 1, 1] R such that g (x) = f(x) for all
3 Definite Integrl 3.1 Introduction In school one comes cross the definition of the integrl of rel vlued function defined on closed nd bounded intervl [, b] between the limits nd b, i.e., f(x)dx s the
More informationAN INEQUALITY OF OSTROWSKI TYPE AND ITS APPLICATIONS FOR SIMPSON S RULE AND SPECIAL MEANS. I. Fedotov and S. S. Dragomir
RGMIA Reserch Report Collection, Vol., No., 999 http://sci.vu.edu.u/ rgmi AN INEQUALITY OF OSTROWSKI TYPE AND ITS APPLICATIONS FOR SIMPSON S RULE AND SPECIAL MEANS I. Fedotov nd S. S. Drgomir Astrct. An
More informationNumerical Integration. Newton Cotes Formulas. Quadrature. Newton Cotes Formulas. To approximate the integral b
Numericl Integrtion Newton Cotes Formuls Given function f : R R nd two rel numbers, b R, < b, we clculte (pproximtely) the integrl I(f,, b) = f (x) dx K. Frischmuth (IfM UR) Numerics for CSE 08/09 8 /
More informationLecture 1. Functional series. Pointwise and uniform convergence.
1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is
More informationLecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature
Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 6: Singulr Integrls, Open Qudrture rules, nd Guss Qudrture (Compiled 6 August 7) In this lecture we discuss the
More information2. General formula for Runge-Kutta methods
NOLTA, IEICE Pper Equivlent circuits for implicit unge-kutt metods in circuit simultors for nonliner circuits Ysuiko Toym 1), Ttusy Kuwzki 1, nd Jun Sirtki 1 Mkiko Okumur 1 1 Kngw Institute of Tecnology,
More informationa < a+ x < a+2 x < < a+n x = b, n A i n f(x i ) x. i=1 i=1
Mth 33 Volume Stewrt 5.2 Geometry of integrls. In this section, we will lern how to compute volumes using integrls defined by slice nlysis. First, we recll from Clculus I how to compute res. Given the
More informationAPPROXIMATE INTEGRATION
APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be
More informationThe Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve.
Clculus Li Vs The Fundmentl Theorem of Clculus. The Totl Chnge Theorem nd the Are Under Curve. Recll the following fct from Clculus course. If continuous function f(x) represents the rte of chnge of F
More informationIntegration Techniques
Integrtion Techniques. Integrtion of Trigonometric Functions Exmple. Evlute cos x. Recll tht cos x = cos x. Hence, cos x Exmple. Evlute = ( + cos x) = (x + sin x) + C = x + 4 sin x + C. cos 3 x. Let u
More informationMath 426: Probability Final Exam Practice
Mth 46: Probbility Finl Exm Prctice. Computtionl problems 4. Let T k (n) denote the number of prtitions of the set {,..., n} into k nonempty subsets, where k n. Argue tht T k (n) kt k (n ) + T k (n ) by
More informationLecture 23: Interpolatory Quadrature
Lecture 3: Interpoltory Qudrture. Qudrture. The computtion of continuous lest squres pproximtions to f C[, b] required evlutions of the inner product f, φ j = fxφ jx dx, where φ j is polynomil bsis function
More informationLecture 20: Numerical Integration III
cs4: introduction to numericl nlysis /8/0 Lecture 0: Numericl Integrtion III Instructor: Professor Amos Ron Scribes: Mrk Cowlishw, Yunpeng Li, Nthnel Fillmore For the lst few lectures we hve discussed
More informationAn iterative method for solving nonlinear functional equations
J. Mth. Anl. Appl. 316 (26) 753 763 www.elsevier.com/locte/jm An itertive method for solving nonliner functionl equtions Vrsh Dftrdr-Gejji, Hossein Jfri Deprtment of Mthemtics, University of Pune, Gneshkhind,
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More informationTheoretical foundations of Gaussian quadrature
Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of
More informationNumerical Methods I Orthogonal Polynomials
Numericl Methods I Orthogonl Polynomils Aleksndr Donev Cournt Institute, NYU 1 donev@cournt.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fll 2014 Nov 6th, 2014 A. Donev (Cournt Institute) Lecture IX
More informationNEW INEQUALITIES OF SIMPSON S TYPE FOR s CONVEX FUNCTIONS WITH APPLICATIONS. := f (4) (x) <. The following inequality. 2 b a
NEW INEQUALITIES OF SIMPSON S TYPE FOR s CONVEX FUNCTIONS WITH APPLICATIONS MOHAMMAD ALOMARI A MASLINA DARUS A AND SEVER S DRAGOMIR B Abstrct In terms of the first derivtive some ineulities of Simpson
More informationB.Sc. in Mathematics (Ordinary)
R48/0 DUBLIN INSTITUTE OF TECHNOLOGY KEVIN STREET, DUBLIN 8 B.Sc. in Mthemtics (Ordinry) SUPPLEMENTAL EXAMINATIONS 01 Numericl Methods Dr. D. Mckey Dr. C. Hills Dr. E.A. Cox Full mrks for complete nswers
More informationThe final exam will take place on Friday May 11th from 8am 11am in Evans room 60.
Mth 104: finl informtion The finl exm will tke plce on Fridy My 11th from 8m 11m in Evns room 60. The exm will cover ll prts of the course with equl weighting. It will cover Chpters 1 5, 7 15, 17 21, 23
More informationDefinite Integrals. The area under a curve can be approximated by adding up the areas of rectangles = 1 1 +
Definite Integrls --5 The re under curve cn e pproximted y dding up the res of rectngles. Exmple. Approximte the re under y = from x = to x = using equl suintervls nd + x evluting the function t the left-hnd
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More information4.4 Areas, Integrals and Antiderivatives
. res, integrls nd ntiderivtives 333. Ares, Integrls nd Antiderivtives This section explores properties of functions defined s res nd exmines some connections mong res, integrls nd ntiderivtives. In order
More informationHarmonic Mean Derivative - Based Closed Newton Cotes Quadrature
IOSR Journl of Mthemtics (IOSR-JM) e-issn: - p-issn: 9-X. Volume Issue Ver. IV (My. - Jun. 0) PP - www.iosrjournls.org Hrmonic Men Derivtive - Bsed Closed Newton Cotes Qudrture T. Rmchndrn D.Udykumr nd
More informationTHE HANKEL MATRIX METHOD FOR GAUSSIAN QUADRATURE IN 1 AND 2 DIMENSIONS
THE HANKEL MATRIX METHOD FOR GAUSSIAN QUADRATURE IN 1 AND 2 DIMENSIONS CARLOS SUERO, MAURICIO ALMANZAR CONTENTS 1 Introduction 1 2 Proof of Gussin Qudrture 6 3 Iterted 2-Dimensionl Gussin Qudrture 20 4
More informationFundamental Theorem of Calculus
Funmentl Teorem of Clculus Liming Png 1 Sttement of te Teorem Te funmentl Teorem of Clculus is one of te most importnt teorems in te istory of mtemtics, wic ws first iscovere by Newton n Leibniz inepenently.
More informationMath 20C Multivariable Calculus Lecture 5 1. Lines and planes. Equations of lines (Vector, parametric, and symmetric eqs.). Equations of lines
Mt 2C Multivrible Clculus Lecture 5 1 Lines nd plnes Slide 1 Equtions of lines (Vector, prmetric, nd symmetric eqs.). Equtions of plnes. Distnce from point to plne. Equtions of lines Slide 2 Definition
More informationOptimal Iterate of the Power and Inverse Iteration Methods
Optimal Iterate of the Power and Inverse Iteration Methods Davod Khojasteh Salkuyeh and Faezeh Toutounian Department of Mathematics, University of Mohaghegh Ardabili, P.O. Box. 79, Ardabil, Iran E-mail:
More informationf (a) + f (b) f (λx + (1 λ)y) max {f (x),f (y)}, x, y [a, b]. (1.1)
TAMKANG JOURNAL OF MATHEMATICS Volume 41, Number 4, 353-359, Winter 1 NEW INEQUALITIES OF HERMITE-HADAMARD TYPE FOR FUNCTIONS WHOSE SECOND DERIVATIVES ABSOLUTE VALUES ARE QUASI-CONVEX M. ALOMARI, M. DARUS
More informationCOT4501 Spring Homework VII
COT451 Spring 1 Homework VII The ssignment is due in clss on Thursdy, April 19, 1. There re five regulr problems nd one computer problem (using MATLAB). For written problems, you need to show your work
More informationSolutions to Assignment #8
Mth 1 Numericl Anlysis (Bueler) December 9, 29 Solutions to Assignment #8 Problems 64, exercise 14: The nswer turns out to be yes, which mens tht I hve to be orgnized in writing it up There re lot of fcts
More information38 Riemann sums and existence of the definite integral.
38 Riemnn sums nd existence of the definite integrl. In the clcultion of the re of the region X bounded by the grph of g(x) = x 2, the x-xis nd 0 x b, two sums ppered: ( n (k 1) 2) b 3 n 3 re(x) ( n These
More informationINTRODUCTION TO INTEGRATION
INTRODUCTION TO INTEGRATION 5.1 Ares nd Distnces Assume f(x) 0 on the intervl [, b]. Let A be the re under the grph of f(x). b We will obtin n pproximtion of A in the following three steps. STEP 1: Divide
More informationDiscrete Least-squares Approximations
Discrete Lest-squres Approximtions Given set of dt points (x, y ), (x, y ),, (x m, y m ), norml nd useful prctice in mny pplictions in sttistics, engineering nd other pplied sciences is to construct curve
More informationMath 554 Integration
Mth 554 Integrtion Hndout #9 4/12/96 Defn. A collection of n + 1 distinct points of the intervl [, b] P := {x 0 = < x 1 < < x i 1 < x i < < b =: x n } is clled prtition of the intervl. In this cse, we
More informationTaylor Polynomial Inequalities
Tylor Polynomil Inequlities Ben Glin September 17, 24 Abstrct There re instnces where we my wish to pproximte the vlue of complicted function round given point by constructing simpler function such s polynomil
More informationSolution to Fredholm Fuzzy Integral Equations with Degenerate Kernel
Int. J. Contemp. Mth. Sciences, Vol. 6, 2011, no. 11, 535-543 Solution to Fredholm Fuzzy Integrl Equtions with Degenerte Kernel M. M. Shmivnd, A. Shhsvrn nd S. M. Tri Fculty of Science, Islmic Azd University
More informationAbstract inner product spaces
WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the
More information