4.4 Improper Integrals
|
|
- Kelley Powers
- 6 years ago
- Views:
Transcription
1 4.4 Improper Integrls 4 which contin no singulrities, nd where the endpoints re lso nonsingulr. qromb, in such circumstnces, tkes mny, mny fewer function evlutions thn either of the routines in 4.2. For exmple, the integrl 2 0 x 4 log(x + x 2 +)dx converges (with prmeters s shown bove) on the very first extrpoltion, fter just 5 clls to trpzd, while qsimp requires 8 clls (8 times s mny evlutions of the integrnd) nd qtrp requires 3 clls (mking 256 times s mny evlutions of the integrnd). CITED REFERENCES AND FURTHER READING: Stoer, J., nd Bulirsch, R. 980, Introduction to Numericl Anlysis (New York: Springer-Verlg), Dhlquist, G., nd Bjorck, A. 974, Numericl Methods (Englewood Cliffs, NJ: Prentice-Hll), Rlston, A., nd Rbinowitz, P. 978, A First Course in Numericl Anlysis, 2nd ed. (New York: McGrw-Hill), Improper Integrls For our present purposes, n integrl will be improper if it hs ny of the following problems: its integrnd goes to finite limiting vlue t finite upper nd lower limits, but cnnot be evluted right on one of those limits (e.g., sin x/x t x =0) its upper limit is, or its lower limit is it hs n integrble singulrity t either limit (e.g., x /2 t x =0) it hs n integrble singulrity t known plce between its upper nd lower limits it hs n integrble singulrity t n unknown plce between its upper nd lower limits If n integrl is infinite (e.g., x dx), or does not exist in limiting sense (e.g., cos xdx), we do not cll it improper; we cll it impossible. No mount of clever lgorithmics will return meningful nswer to n ill-posed problem. In this section we will generlize the techniques of the preceding two sections to cover the first four problems on the bove list. A more dvnced discussion of qudrture with integrble singulrities occurs in Chpter 8, notbly 8.3. The fifth problem, singulrity t unknown loction, cn relly only be hndled by the use of vrible stepsize differentil eqution integrtion routine, s will be given in Chpter 6. We need workhorse like the extended trpezoidl rule (eqution 4..), but one which is n open formul in the sense of 4., i.e., does not require the integrnd to be evluted t the endpoints. Eqution (4..9), the extended midpoint rule, is the best choice. The reson is tht (4..9) shres with (4..) the deep property
2 42 Chpter 4. Integrtion of Functions of hving n error series tht is entirely even in h. Indeed there is formul, not s well known s it ought to be, clled the Second Euler-Mclurin summtion formul, xn h[f 3/2 + f 5/2 + f 7/2 + +f N 3/2 +f N /2 ] x + B 2h 2 4 (f N f )+ + B 2kh 2k ( 2 2k+ )(f (2k ) N f (2k ) )+ (2k)! (4.4.) This eqution cn be derived by writing out (4.2.) with stepsize h, then writing it out gin with stepsize h/2, then subtrcting the first from twice the second. It is not possible to double the number of steps in the extended midpoint rule nd still hve the benefit of previous function evlutions (try it!). However, it is possible to triple the number of steps nd do so. Shll we do this, or double nd ccept the loss? On the verge, tripling does fctor 3 of unnecessry work, since the right number of steps for desired ccurcy criterion my in fct fll nywhere in the logrithmic intervl implied by tripling. For doubling, the fctor is only 2, but we lose n extr fctor of 2 in being unble to use ll the previous evlutions. Since.732 < 2.44, it is better to triple. Here is the resulting routine, which is directly comprble to trpzd. #define FUNC(x) ((*func)(x)) flot midpnt(flot (*func)(flot), flot, flot b, int n) This routine computes the nth stge of refinement of n extended midpoint rule. func is input s pointer to the function to be integrted between limits nd b, lso input. When clled with n=, the routine returns the crudest estimte of f(x)dx. Subsequent clls with n=2,3,... (in tht sequentil order) will improve the ccurcy of s by dding (2/3) 3 n- dditionl interior points. s should not be modified between sequentil clls. flot x,tnm,sum,del,ddel; sttic flot s; if (n == ) return (s=(b-)*func(0.5*(+b))); else for(it=,j=;j<n-;j++) it *= 3; tnm=it; del=(b-)/(3.0*tnm); ddel=del+del; x=+0.5*del; sum=0.0; for (j=;j<=it;j++) x += ddel; x += del; s=(s+(b-)*sum/tnm)/3.0; return s; The dded points lternte in spcing between del nd ddel. The new sum is combined with the old integrl to give refined integrl.
3 4.4 Improper Integrls 43 The routine midpnt cn exctly replce trpzd in driver routine like qtrp ( 4.2); one simply chnges trpzd(func,,b,j) to midpnt(func,,b, j), nd perhps lso decreses the prmeter JMAX since 3 JMAX (from step tripling) is much lrger number thn 2 JMAX (step doubling). The open formul implementtion nlogous to Simpson s rule (qsimp in 4.2) substitutes midpnt for trpzd nd decreses JMAX s bove, but now lso chnges the extrpoltion step to be s=(9.0*st-ost)/8.0; since, when the number of steps is tripled, the error decreses to /9th its size, not /4th s with step doubling. Either the modified qtrp or the modified qsimp will fix the first problem on the list t the beginning of this section. Yet more sophisticted is to generlize Romberg integrtion in like mnner: #include <mth.h> #define EPS.0e-6 #define JMAX 4 #define JMAXP (JMAX+) #define K 5 flot qromo(flot (*func)(flot), flot, flot b, flot (*choose)(flot(*)(flot), flot, flot, int)) Romberg integrtion on n open intervl. Returns the integrl of the function func from to b, using ny specified integrting function choose nd Romberg s method. Normlly choose will be n open formul, not evluting the function t the endpoints. It is ssumed tht choose triples the number of steps on ech cll, nd tht its error series contins only even powers of the number of steps. The routines midpnt, midinf, midsql, midsqu, midexp, re possible choices for choose. The prmeters hve the sme mening s in qromb. void polint(flot x[], flot y[], int n, flot x, flot *y, flot *dy); void nrerror(chr error_text[]); int j; flot ss,dss,h[jmaxp+],s[jmaxp]; h[]=.0; for (j=;j<=jmax;j++) s[j]=(*choose)(func,,b,j); if (j >= K) polint(&h[j-k],&s[j-k],k,0.0,&ss,&dss); if (fbs(dss) <= EPS*fbs(ss)) return ss; h[j+]=h[j]/9.0; This is where the ssumption of step tripling nd n even error series is used. nrerror("too mny steps in routing qromo"); return 0.0; Never get here. Don t be put off by qromo s complicted ANSI declrtion. A typicl invoction (integrting the Bessel function Y 0 (x) from 0 to 2) is simply #include "nr.h" flot nswer;... nswer=qromo(bessy0,0.0,2.0,midpnt);
4 44 Chpter 4. Integrtion of Functions The differences between qromo nd qromb ( 4.3) re so slighttht it is perhps grtuitous to list qromo in full. It, however, is n excellent driver routine for solving ll the other problems of improper integrls in our first list (except the intrctble fifth), s we shll now see. The bsic trick for improper integrls is to mke chnge of vribles to eliminte the singulrity, or to mp n infinite rnge of integrtion to finite one. For exmple, the identity / /b ( ) t 2 f dt b > 0 (4.4.2) t cn be used with either b nd positive, or with nd b negtive, nd works for ny function which decreses towrds infinity fster thn /x 2. You cn mke the chnge of vrible implied by (4.4.2) either nlyticlly nd then use (e.g.) qromo nd midpnt to do the numericl evlution, or you cn let the numericl lgorithm mke the chnge of vrible for you. We prefer the ltter method s being more trnsprent to the user. To implement eqution (4.4.2) we simply write modified version of midpnt, clled midinf, which llows b to be infinite (or, more precisely, very lrge number on your prticulr mchine, such s 0 30 ), or to be negtive nd infinite. #define FUNC(x) ((*funk)(.0/(x))/((x)*(x))) Effects the chnge of vrible. flot midinf(flot (*funk)(flot), flot, flot bb, int n) This routine is n exct replcement for midpnt, i.e., returns the nth stge of refinement of the integrl of funk from to bb, except tht the function is evluted t evenly spced points in /x rther thn in x. This llows the upper limit bb to be s lrge nd positive s the computer llows, or the lower limit to be s lrge nd negtive, but not both. nd bb must hve the sme sign. flot x,tnm,sum,del,ddel,b,; sttic flot s; b=.0/; These two sttements chnge the limits of integrtion. =.0/bb; if (n == ) From this point on, the routine is identicl to midpnt. return (s=(b-)*func(0.5*(+b))); else for(it=,j=;j<n-;j++) it *= 3; tnm=it; del=(b-)/(3.0*tnm); ddel=del+del; x=+0.5*del; sum=0.0; for (j=;j<=it;j++) x += ddel; x += del; return (s=(s+(b-)*sum/tnm)/3.0);
5 4.4 Improper Integrls 45 If you need to integrte from negtive lower limit to positive infinity, you do this by breking the integrl into two pieces t some positive vlue, for exmple, nswer=qromo(funk,-5.0,2.0,midpnt)+qromo(funk,2.0,.0e30,midinf); Where should you choose the brekpoint? At sufficiently lrge positive vlue so tht the function funk is t lest beginning to pproch its symptotic decrese to zero vlue t infinity. The polynomil extrpoltion implicit in the second cll to qromo dels with polynomil in /x, not in x. To del with n integrl tht hs n integrble power-lw singulrity t its lower limit, one lso mkes chnge of vrible. If the integrnd diverges s (x ) γ, 0 γ<, ner x =, use the identity (b ) γ t γ γ f(t γ + )dt (b >) (4.4.3) γ 0 If the singulrity is t the upper limit, use the identity (b ) γ t γ γ f(b t γ )dt (b >) (4.4.4) γ 0 If there is singulrity t both limits, divide the integrl t n interior brekpoint s in the exmple bove. Equtions (4.4.3) nd (4.4.4) re prticulrly simple in the cse of inverse squre-root singulrities, cse tht occurs frequently in prctice: for singulrity t, nd b 0 b 0 2tf( + t 2 )dt (b >) (4.4.5) 2tf(b t 2 )dt (b >) (4.4.6) for singulrity t b. Once gin, we cn implement these chnges of vrible trnsprently to the user by defining substitute routines for midpnt which mke the chnge of vrible utomticlly: #include <mth.h> #define FUNC(x) (2.0*(x)*(*funk)(+(x)*(x))) flot midsql(flot (*funk)(flot), flot, flot bb, int n) This routine is n exct replcement formidpnt, except tht it llows for n inverse squre-root singulrity in the integrnd t the lower limit. flot x,tnm,sum,del,ddel,,b; sttic flot s; b=sqrt(bb-); =0.0; if (n == ) The rest of the routine is exctly like midpnt nd is omitted.
6 46 Chpter 4. Integrtion of Functions Similrly, #include <mth.h> #define FUNC(x) (2.0*(x)*(*funk)(bb-(x)*(x))) flot midsqu(flot (*funk)(flot), flot, flot bb, int n) This routine is n exct replcement formidpnt, except tht it llows for n inverse squre-root singulrity in the integrnd t the upper limit bb. flot x,tnm,sum,del,ddel,,b; sttic flot s; b=sqrt(bb-); =0.0; if (n == ) The rest of the routine is exctly like midpnt nd is omitted. One lst exmple should suffice to show how these formuls re derived in generl. Suppose the upper limit of integrtion is infinite, nd the integrnd flls off exponentilly. Then we wnt chnge of vrible tht mps e x dx into (±)dt (with the sign chosen to keep the upper limit of the new vrible lrger thn the lower limit). Doing the integrtion gives by inspection so tht x= x= t = e x or x = log t (4.4.7) t=e t=0 The user-trnsprent implementtion would be #include <mth.h> #define FUNC(x) ((*funk)(-log(x))/(x)) f( log t) dt t (4.4.8) flot midexp(flot (*funk)(flot), flot, flot bb, int n) This routine is n exct replcement for midpnt, except tht bb is ssumed to be infinite (vlue pssed not ctully used). It is ssumed tht the function funk decreses exponentilly rpidly t infinity. flot x,tnm,sum,del,ddel,,b; sttic flot s; b=exp(-); =0.0; if (n == ) The rest of the routine is exctly like midpnt nd is omitted. CITED REFERENCES AND FURTHER READING: Acton, F.S. 970, Numericl Methods Tht Work; 990, corrected edition (Wshington: Mthemticl Assocition of Americ), Chpter 4.
7 4.5 Gussin Qudrtures nd Orthogonl Polynomils 47 Dhlquist, G., nd Bjorck, A. 974, Numericl Methods (Englewood Cliffs, NJ: Prentice-Hll), 7.4.3, p Stoer, J., nd Bulirsch, R. 980, Introduction to Numericl Anlysis (New York: Springer-Verlg), 3.7, p Gussin Qudrtures nd Orthogonl Polynomils In the formuls of 4., the integrl of function ws pproximted by the sum of its functionl vlues t set of eqully spced points, multiplied by certin ptly chosen weighting coefficients. We sw tht s we llowed ourselves more freedom in choosing the coefficients, we could chieve integrtion formuls of higher nd higher order. The ide of Gussin qudrtures is to give ourselves the freedom to choose not only the weighting coefficients, but lso the loction of the bscisss t which the function is to be evluted: They will no longer be eqully spced. Thus, we will hve twice the number of degrees of freedom t our disposl; it will turn out tht we cn chieve Gussin qudrture formuls whose order is, essentilly, twice tht of the Newton-Cotes formul with the sme number of function evlutions. Does this sound too good to be true? Well, in sense it is. The ctch is fmilir one, which cnnot be overemphsized: High order is not the sme s high ccurcy. High order trnsltes to high ccurcy only when the integrnd is very smooth, in the sense of being well-pproximted by polynomil. There is, however, one dditionl feture of Gussin qudrture formuls tht dds to their usefulness: We cn rrnge the choice of weights nd bscisss to mke the integrl exct for clss of integrnds polynomils times some known function W (x) rther thn for the usul clss of integrnds polynomils. The function W (x) cn then be chosen to remove integrble singulrities from the desired integrl. Given W (x), in other words, nd given n integer N, we cn find set of weights w j nd bscisss x j such tht the pproximtion W (x)f(x)dx N w j f(x j ) (4.5.) j= is exct if f(x) is polynomil. For exmple, to do the integrl exp( cos 2 x) dx (4.5.2) x 2 (not very nturl looking integrl, it must be dmitted), we might well be interested in Gussin qudrture formul bsed on the choice W (x) = x 2 (4.5.3) in the intervl (, ). (This prticulr choice is clled Guss-Chebyshev integrtion, for resons tht will become cler shortly.)
There are also formulas of higher order for this situation, but we will refrain from giving them.
136 Chapter 4. Integration of Functions N = 1 2 3 4 (total after N = 4) Figure 4.2.1. Sequential calls to the routine trapzd incorporate the information from previous calls and evaluate the integrand only
More informationImproper Integrals. October 4, 2017
Improper Integrls October 4, 7 Introduction We hve seen how to clculte definite integrl when the it is rel number. However, there re times when we re interested to compute the integrl sy for emple 3. Here
More informationIntegration. September 28, 2017
Integrtion September 8, 7 Introduction We hve lerned in previous chpter on how to do the differentition. It is conventionl in mthemtics tht we re supposed to lern bout the integrtion s well. As you my
More informationIntegration. October 25, 2016
Integrtion October 5, 6 Introduction We hve lerned in previous chpter on how to do the differentition. It is conventionl in mthemtics tht we re supposed to lern bout the integrtion s well. As you my hve
More informationDouble Integrals. MATH 375 Numerical Analysis. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Double Integrals
Double Integrls MATH 375 Numericl Anlysis J. Robert Buchnn Deprtment of Mthemtics Fll 2013 J. Robert Buchnn Double Integrls Objectives Now tht we hve discussed severl methods for pproximting definite integrls
More informationSection 10.4 Hyperbolas
66 Section 10.4 Hyperbols Objective : Definition of hyperbol & hyperbols centered t (0, 0). The third type of conic we will study is the hyperbol. It is defined in the sme mnner tht we defined the prbol
More informationMath 464 Fall 2012 Notes on Marginal and Conditional Densities October 18, 2012
Mth 464 Fll 2012 Notes on Mrginl nd Conditionl Densities klin@mth.rizon.edu October 18, 2012 Mrginl densities. Suppose you hve 3 continuous rndom vribles X, Y, nd Z, with joint density f(x,y,z. The mrginl
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 informationP(r)dr = probability of generating a random number in the interval dr near r. For this probability idea to make sense we must have
Rndom Numers nd Monte Crlo Methods Rndom Numer Methods The integrtion methods discussed so fr ll re sed upon mking polynomil pproximtions to the integrnd. Another clss of numericl methods relies upon using
More informationMath/CS 467/667 Programming Assignment 01. Adaptive Gauss Quadrature. q(x)p 4 (x) = 0
Adptive Guss Qudrture 1. Find n orthogonl polynomil p 4 of degree 4 such tht 1 1 q(x)p 4 (x) = 0 for every polynomil q(x) of degree 3 or less. You my use Mple nd the Grm Schmidt process s done in clss.
More informationSolutions to Math 41 Final Exam December 12, 2011
Solutions to Mth Finl Em December,. ( points) Find ech of the following its, with justifiction. If there is n infinite it, then eplin whether it is or. ( ) / ln() () (5 points) First we compute the it:
More informationIntroduction to Integration
Introduction to Integrtion Definite integrls of piecewise constnt functions A constnt function is function of the form Integrtion is two things t the sme time: A form of summtion. The opposite of differentition.
More information)
Chpter Five /SOLUTIONS Since the speed ws between nd mph during this five minute period, the fuel efficienc during this period is between 5 mpg nd 8 mpg. So the fuel used during this period is between
More informationMath 142, Exam 1 Information.
Mth 14, Exm 1 Informtion. 9/14/10, LC 41, 9:30-10:45. Exm 1 will be bsed on: Sections 7.1-7.5. The corresponding ssigned homework problems (see http://www.mth.sc.edu/ boyln/sccourses/14f10/14.html) At
More informationIf f(x, y) is a surface that lies above r(t), we can think about the area between the surface and the curve.
Line Integrls The ide of line integrl is very similr to tht of single integrls. If the function f(x) is bove the x-xis on the intervl [, b], then the integrl of f(x) over [, b] is the re under f over the
More information12-B FRACTIONS AND DECIMALS
-B Frctions nd Decimls. () If ll four integers were negtive, their product would be positive, nd so could not equl one of them. If ll four integers were positive, their product would be much greter thn
More information9.1 apply the distance and midpoint formulas
9.1 pply the distnce nd midpoint formuls DISTANCE FORMULA MIDPOINT FORMULA To find the midpoint between two points x, y nd x y 1 1,, we Exmple 1: Find the distnce between the two points. Then, find the
More information1. SEQUENCES INVOLVING EXPONENTIAL GROWTH (GEOMETRIC SEQUENCES)
Numbers nd Opertions, Algebr, nd Functions 45. SEQUENCES INVOLVING EXPONENTIAL GROWTH (GEOMETRIC SEQUENCES) In sequence of terms involving eponentil growth, which the testing service lso clls geometric
More information10.5 Graphing Quadratic Functions
0.5 Grphing Qudrtic Functions Now tht we cn solve qudrtic equtions, we wnt to lern how to grph the function ssocited with the qudrtic eqution. We cll this the qudrtic function. Grphs of Qudrtic Functions
More information4452 Mathematical Modeling Lecture 4: Lagrange Multipliers
Mth Modeling Lecture 4: Lgrnge Multipliers Pge 4452 Mthemticl Modeling Lecture 4: Lgrnge Multipliers Lgrnge multipliers re high powered mthemticl technique to find the mximum nd minimum of multidimensionl
More informationsuch that the S i cover S, or equivalently S
MATH 55 Triple Integrls Fll 16 1. Definition Given solid in spce, prtition of consists of finite set of solis = { 1,, n } such tht the i cover, or equivlently n i. Furthermore, for ech i, intersects i
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 informationUnit 5 Vocabulary. A function is a special relationship where each input has a single output.
MODULE 3 Terms Definition Picture/Exmple/Nottion 1 Function Nottion Function nottion is n efficient nd effective wy to write functions of ll types. This nottion llows you to identify the input vlue with
More informationLecture 7: Integration Techniques
Lecture 7: Integrtion Techniques Antiderivtives nd Indefinite Integrls. In differentil clculus, we were interested in the derivtive of given rel-vlued function, whether it ws lgeric, eponentil or logrithmic.
More informationMA1008. Calculus and Linear Algebra for Engineers. Course Notes for Section B. Stephen Wills. Department of Mathematics. University College Cork
MA1008 Clculus nd Liner Algebr for Engineers Course Notes for Section B Stephen Wills Deprtment of Mthemtics University College Cork s.wills@ucc.ie http://euclid.ucc.ie/pges/stff/wills/teching/m1008/ma1008.html
More informationMatlab s Numerical Integration Commands
Mtlb s Numericl Integrtion Commnds The relevnt commnds we consider re qud nd dblqud, triplequd. See the Mtlb help files for other integrtion commnds. By the wy, qud refers to dptive qudrture. To integrte:
More information9 4. CISC - Curriculum & Instruction Steering Committee. California County Superintendents Educational Services Association
9. CISC - Curriculum & Instruction Steering Committee The Winning EQUATION A HIGH QUALITY MATHEMATICS PROFESSIONAL DEVELOPMENT PROGRAM FOR TEACHERS IN GRADES THROUGH ALGEBRA II STRAND: NUMBER SENSE: Rtionl
More informationRational Numbers---Adding Fractions With Like Denominators.
Rtionl Numbers---Adding Frctions With Like Denomintors. A. In Words: To dd frctions with like denomintors, dd the numertors nd write the sum over the sme denomintor. B. In Symbols: For frctions c nd b
More informationMATH 25 CLASS 5 NOTES, SEP
MATH 25 CLASS 5 NOTES, SEP 30 2011 Contents 1. A brief diversion: reltively prime numbers 1 2. Lest common multiples 3 3. Finding ll solutions to x + by = c 4 Quick links to definitions/theorems Euclid
More informationMidterm 2 Sample solution
Nme: Instructions Midterm 2 Smple solution CMSC 430 Introduction to Compilers Fll 2012 November 28, 2012 This exm contins 9 pges, including this one. Mke sure you hve ll the pges. Write your nme on the
More informationSubtracting Fractions
Lerning Enhncement Tem Model Answers: Adding nd Subtrcting Frctions Adding nd Subtrcting Frctions study guide. When the frctions both hve the sme denomintor (bottom) you cn do them using just simple dding
More information6.2 Volumes of Revolution: The Disk Method
mth ppliction: volumes by disks: volume prt ii 6 6 Volumes of Revolution: The Disk Method One of the simplest pplictions of integrtion (Theorem 6) nd the ccumultion process is to determine so-clled volumes
More informationx )Scales are the reciprocal of each other. e
9. Reciprocls A Complete Slide Rule Mnul - eville W Young Chpter 9 Further Applictions of the LL scles The LL (e x ) scles nd the corresponding LL 0 (e -x or Exmple : 0.244 4.. Set the hir line over 4.
More information18.3 Integral Equations with Singular Kernels
18.3 Integrl Equtions with Singulr Kernels 797 This procedure cn be repeted s with Romberg integrtion. The generl consensus is tht the best of the higher order methods is the block-by-block method (see
More informationSection 3.1: Sequences and Series
Section.: Sequences d Series Sequences Let s strt out with the definition of sequence: sequence: ordered list of numbers, often with definite pttern Recll tht in set, order doesn t mtter so this is one
More information5/9/17. Lesson 51 - FTC PART 2. Review FTC, PART 1. statement as the Integral Evaluation Theorem as it tells us HOW to evaluate the definite integral
Lesson - FTC PART 2 Review! We hve seen definition/formul for definite integrl s n b A() = lim f ( i )Δ = f ()d = F() = F(b) F() n i=! where F () = f() (or F() is the ntiderivtive of f() b! And hve seen
More informationMATH 2530: WORKSHEET 7. x 2 y dz dy dx =
MATH 253: WORKSHT 7 () Wrm-up: () Review: polr coordintes, integrls involving polr coordintes, triple Riemnn sums, triple integrls, the pplictions of triple integrls (especilly to volume), nd cylindricl
More informationThis notebook investigates the properties of non-integer differential operators using Fourier analysis.
Frctionl erivtives.nb Frctionl erivtives by Fourier ecomposition by Eric Thrne 4/9/6 This notebook investigtes the properties of non-integer differentil opertors using Fourier nlysis. In[]:=
More information6.3 Definite Integrals and Antiderivatives
Section 6. Definite Integrls nd Antiderivtives 8 6. Definite Integrls nd Antiderivtives Wht ou will lern out... Properties of Definite Integrls Averge Vlue of Function Men Vlue Theorem for Definite Integrls
More informationFig.25: the Role of LEX
The Lnguge for Specifying Lexicl Anlyzer We shll now study how to uild lexicl nlyzer from specifiction of tokens in the form of list of regulr expressions The discussion centers round the design of n existing
More informationCHAPTER 8 Quasi-interpolation methods
CHAPTER 8 Qusi-interpoltion methods In Chpter 5 we considered number of methods for computing spline pproximtions. The strting point for the pproximtion methods is dt set tht is usully discrete nd in the
More informationIn the last lecture, we discussed how valid tokens may be specified by regular expressions.
LECTURE 5 Scnning SYNTAX ANALYSIS We know from our previous lectures tht the process of verifying the syntx of the progrm is performed in two stges: Scnning: Identifying nd verifying tokens in progrm.
More informationBefore We Begin. Introduction to Spatial Domain Filtering. Introduction to Digital Image Processing. Overview (1): Administrative Details (1):
Overview (): Before We Begin Administrtive detils Review some questions to consider Winter 2006 Imge Enhncement in the Sptil Domin: Bsics of Sptil Filtering, Smoothing Sptil Filters, Order Sttistics Filters
More informationCS201 Discussion 10 DRAWTREE + TRIES
CS201 Discussion 10 DRAWTREE + TRIES DrwTree First instinct: recursion As very generic structure, we could tckle this problem s follows: drw(): Find the root drw(root) drw(root): Write the line for the
More information2 Computing all Intersections of a Set of Segments Line Segment Intersection
15-451/651: Design & Anlysis of Algorithms Novemer 14, 2016 Lecture #21 Sweep-Line nd Segment Intersection lst chnged: Novemer 8, 2017 1 Preliminries The sweep-line prdigm is very powerful lgorithmic design
More information1 The Definite Integral
The Definite Integrl Definition. Let f be function defined on the intervl [, b] where
More informationChapter Spline Method of Interpolation More Examples Electrical Engineering
Chpter. Spline Method of Interpoltion More Exmples Electricl Engineering Exmple Thermistors re used to mesure the temperture of bodies. Thermistors re bsed on mterils chnge in resistnce with temperture.
More informationParallel Square and Cube Computations
Prllel Squre nd Cube Computtions Albert A. Liddicot nd Michel J. Flynn Computer Systems Lbortory, Deprtment of Electricl Engineering Stnford University Gtes Building 5 Serr Mll, Stnford, CA 945, USA liddicot@stnford.edu
More informationThe Reciprocal Function Family. Objectives To graph reciprocal functions To graph translations of reciprocal functions
- The Reciprocl Function Fmil Objectives To grph reciprocl functions To grph trnsltions of reciprocl functions Content Stndrds F.BF.3 Identif the effect on the grph of replcing f () b f() k, kf(), f(k),
More informationPointwise convergence need not behave well with respect to standard properties such as continuity.
Chpter 3 Uniform Convergence Lecture 9 Sequences of functions re of gret importnce in mny res of pure nd pplied mthemtics, nd their properties cn often be studied in the context of metric spces, s in Exmples
More informationMTH 146 Conics Supplement
105- Review of Conics MTH 146 Conics Supplement In this section we review conics If ou ne more detils thn re present in the notes, r through section 105 of the ook Definition: A prol is the set of points
More informationDynamic Programming. Andreas Klappenecker. [partially based on slides by Prof. Welch] Monday, September 24, 2012
Dynmic Progrmming Andres Klppenecker [prtilly bsed on slides by Prof. Welch] 1 Dynmic Progrmming Optiml substructure An optiml solution to the problem contins within it optiml solutions to subproblems.
More informationFall 2018 Midterm 1 October 11, ˆ You may not ask questions about the exam except for language clarifications.
15-112 Fll 2018 Midterm 1 October 11, 2018 Nme: Andrew ID: Recittion Section: ˆ You my not use ny books, notes, extr pper, or electronic devices during this exm. There should be nothing on your desk or
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd processes. Introducing technology
More informationPresentation Martin Randers
Presenttion Mrtin Rnders Outline Introduction Algorithms Implementtion nd experiments Memory consumption Summry Introduction Introduction Evolution of species cn e modelled in trees Trees consist of nodes
More informationCHAPTER 5 Spline Approximation of Functions and Data
CHAPTER 5 Spline Approximtion of Functions nd Dt This chpter introduces number of methods for obtining spline pproximtions to given functions, or more precisely, to dt obtined by smpling function. In Section
More informationThe Basic Properties of the Integral
The Bsic Properties of the Integrl When we compute the derivtive of complicted function, like + sin, we usull use differentition rules, like d [f()+g()] d f()+ d g(), to reduce the computtion d d d to
More informationLecture Overview. Knowledge-based systems in Bioinformatics, 1MB602. Procedural abstraction. The sum procedure. Integration as a procedure
Lecture Overview Knowledge-bsed systems in Bioinformtics, MB6 Scheme lecture Procedurl bstrction Higher order procedures Procedures s rguments Procedures s returned vlues Locl vribles Dt bstrction Compound
More informationMatrices and Systems of Equations
Mtrices Mtrices nd Sstems of Equtions A mtri is rectngulr rr of rel numbers. CHAT Pre-Clculus Section 8. m m m............ n n n mn We will use the double subscript nottion for ech element of the mtri.
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd business. Introducing technology
More informationIntroduction. Chapter 4: Complex Integration. Introduction (Cont d)
Introduction Chpter 4: Complex Integrtion Li, Yongzho Stte Key Lbortory of Integrted Services Networks, Xidin University October 10, 2010 The two-dimensionl nture of the complex plne required us to generlize
More informationThe Fundamental Theorem of Calculus
MATH 6 The Fundmentl Theorem of Clculus The Fundmentl Theorem of Clculus (FTC) gives method of finding the signed re etween the grph of f nd the x-xis on the intervl [, ]. The theorem is: FTC: If f is
More information6.3 Volumes. Just as area is always positive, so is volume and our attitudes towards finding it.
6.3 Volumes Just s re is lwys positive, so is volume nd our ttitudes towrds finding it. Let s review how to find the volume of regulr geometric prism, tht is, 3-dimensionl oject with two regulr fces seprted
More informationRepresentation of Numbers. Number Representation. Representation of Numbers. 32-bit Unsigned Integers 3/24/2014. Fixed point Integer Representation
Representtion of Numbers Number Representtion Computer represent ll numbers, other thn integers nd some frctions with imprecision. Numbers re stored in some pproximtion which cn be represented by fixed
More informationHomework. Context Free Languages III. Languages. Plan for today. Context Free Languages. CFLs and Regular Languages. Homework #5 (due 10/22)
Homework Context Free Lnguges III Prse Trees nd Homework #5 (due 10/22) From textbook 6.4,b 6.5b 6.9b,c 6.13 6.22 Pln for tody Context Free Lnguges Next clss of lnguges in our quest! Lnguges Recll. Wht
More informationGENG2140 Modelling and Computer Analysis for Engineers
GENG4 Moelling n Computer Anlysis or Engineers Letures 9 & : Gussin qurture Crete y Grn Romn Joles, PhD Shool o Mehnil Engineering, UWA GENG4 Content Deinition o Gussin qurture Computtion o weights n points
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd processes. Introducing technology
More informationCOMP 423 lecture 11 Jan. 28, 2008
COMP 423 lecture 11 Jn. 28, 2008 Up to now, we hve looked t how some symols in n lphet occur more frequently thn others nd how we cn sve its y using code such tht the codewords for more frequently occuring
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd business. Introducing technology
More informationNumerical integration methods
Chpter 1 Numericl integrtion methods The bility to clculte integrls is quite importnt. The uthor ws told tht, in the old dys, the gun ports were cut into ship only fter it ws flot, loded with equivlent
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd processes. Introducing technology
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd business. Introducing technology
More informationStudy Guide for Exam 3
Mth 05 Elementry Algebr Fll 00 Study Guide for Em Em is scheduled for Thursdy, November 8 th nd ill cover chpters 5 nd. You my use "5" note crd (both sides) nd scientific clcultor. You re epected to no
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd processes. Introducing technology
More informationIf you are at the university, either physically or via the VPN, you can download the chapters of this book as PDFs.
Lecture 5 Wlks, Trils, Pths nd Connectedness Reding: Some of the mteril in this lecture comes from Section 1.2 of Dieter Jungnickel (2008), Grphs, Networks nd Algorithms, 3rd edition, which is ville online
More information9 Graph Cutting Procedures
9 Grph Cutting Procedures Lst clss we begn looking t how to embed rbitrry metrics into distributions of trees, nd proved the following theorem due to Brtl (1996): Theorem 9.1 (Brtl (1996)) Given metric
More informationSmall Business Networking
Why network is n essentil productivity tool for ny smll business Effective technology is essentil for smll businesses looking to increse the productivity of their people nd business. Introducing technology
More informationSummer Review Packet For Algebra 2 CP/Honors
Summer Review Pcket For Alger CP/Honors Nme Current Course Mth Techer Introduction Alger uilds on topics studied from oth Alger nd Geometr. Certin topics re sufficientl involved tht the cll for some review
More informationWebAssign Lesson 1-3a Substitution Part 1 (Homework)
WeAssign Lesson -3 Sustitution Prt (Homework) Current Score : / 3 Due : Fridy, June 7 04 :00 AM MDT Jimos Skriletz Mth 75, section 3, Summer 04 Instructor: Jimos Skriletz. /.5 points Suppose you hve the
More informationOutline. Tiling, formally. Expression tile as rule. Statement tiles as rules. Function calls. CS 412 Introduction to Compilers
CS 412 Introduction to Compilers Andrew Myers Cornell University Lectur8 Finishing genertion 9 Mr 01 Outline Tiling s syntx-directed trnsltion Implementing function clls Implementing functions Optimizing
More informationa(e, x) = x. Diagrammatically, this is encoded as the following commutative diagrams / X
4. Mon, Sept. 30 Lst time, we defined the quotient topology coming from continuous surjection q : X! Y. Recll tht q is quotient mp (nd Y hs the quotient topology) if V Y is open precisely when q (V ) X
More informationCSCI 104. Rafael Ferreira da Silva. Slides adapted from: Mark Redekopp and David Kempe
CSCI 0 fel Ferreir d Silv rfsilv@isi.edu Slides dpted from: Mrk edekopp nd Dvid Kempe LOG STUCTUED MEGE TEES Series Summtion eview Let n = + + + + k $ = #%& #. Wht is n? n = k+ - Wht is log () + log ()
More informationSlides for Data Mining by I. H. Witten and E. Frank
Slides for Dt Mining y I. H. Witten nd E. Frnk Simplicity first Simple lgorithms often work very well! There re mny kinds of simple structure, eg: One ttriute does ll the work All ttriutes contriute eqully
More informationMath 35 Review Sheet, Spring 2014
Mth 35 Review heet, pring 2014 For the finl exm, do ny 12 of the 15 questions in 3 hours. They re worth 8 points ech, mking 96, with 4 more points for netness! Put ll your work nd nswers in the provided
More informationZZ - Advanced Math Review 2017
ZZ - Advnced Mth Review Mtrix Multipliction Given! nd! find the sum of the elements of the product BA First, rewrite the mtrices in the correct order to multiply The product is BA hs order x since B is
More informationEngineer-to-Engineer Note
Engineer-to-Engineer Note EE-295 Technicl notes on using Anlog Devices DSPs, processors nd development tools Visit our Web resources http://www.nlog.com/ee-notes nd http://www.nlog.com/processors or e-mil
More information8.2 Areas in the Plane
39 Chpter 8 Applictions of Definite Integrls 8. Ares in the Plne Wht ou will lern out... Are Between Curves Are Enclosed Intersecting Curves Boundries with Chnging Functions Integrting with Respect to
More informationf[a] x + f[a + x] x + f[a +2 x] x + + f[b x] x
Bsic Integrtion This chpter contins the fundmentl theory of integrtion. We begin with some problems to motivte the min ide: pproximtion by sum of slices. The chpter confronts this squrely, nd Chpter 3
More informationMA 124 (Calculus II) Lecture 2: January 24, 2019 Section A3. Professor Jennifer Balakrishnan,
Wht is on tody Professor Jennifer Blkrishnn, jbl@bu.edu 1 Velocity nd net chnge 1 2 Regions between curves 3 1 Velocity nd net chnge Briggs-Cochrn-Gillett 6.1 pp. 398-46 Suppose you re driving long stright
More informationChapter 2 Sensitivity Analysis: Differential Calculus of Models
Chpter 2 Sensitivity Anlysis: Differentil Clculus of Models Abstrct Models in remote sensing nd in science nd engineering, in generl re, essentilly, functions of discrete model input prmeters, nd/or functionls
More informationCompilers Spring 2013 PRACTICE Midterm Exam
Compilers Spring 2013 PRACTICE Midterm Exm This is full length prctice midterm exm. If you wnt to tke it t exm pce, give yourself 7 minutes to tke the entire test. Just like the rel exm, ech question hs
More informationCSE 401 Midterm Exam 11/5/10 Sample Solution
Question 1. egulr expressions (20 points) In the Ad Progrmming lnguge n integer constnt contins one or more digits, but it my lso contin embedded underscores. Any underscores must be preceded nd followed
More information2014 Haskell January Test Regular Expressions and Finite Automata
0 Hskell Jnury Test Regulr Expressions nd Finite Automt This test comprises four prts nd the mximum mrk is 5. Prts I, II nd III re worth 3 of the 5 mrks vilble. The 0 Hskell Progrmming Prize will be wrded
More informationFunctor (1A) Young Won Lim 10/5/17
Copyright (c) 2016-2017 Young W. Lim. Permission is grnted to copy, distribute nd/or modify this document under the terms of the GNU Free Documenttion License, Version 1.2 or ny lter version published
More informationEngineer To Engineer Note
Engineer To Engineer Note EE-186 Technicl Notes on using Anlog Devices' DSP components nd development tools Contct our technicl support by phone: (800) ANALOG-D or e-mil: dsp.support@nlog.com Or visit
More informationLU Decomposition. Mechanical Engineering Majors. Authors: Autar Kaw
LU Decomposition Mechnicl Engineering Mjors Authors: Autr Kw Trnsforming Numericl Methods Eduction for STEM Undergrdutes // LU Decomposition LU Decomposition LU Decomposition is nother method to solve
More informationEngineer To Engineer Note
Engineer To Engineer Note EE-169 Technicl Notes on using Anlog Devices' DSP components nd development tools Contct our technicl support by phone: (800) ANALOG-D or e-mil: dsp.support@nlog.com Or visit
More informationDefinition of Regular Expression
Definition of Regulr Expression After the definition of the string nd lnguges, we re redy to descrie regulr expressions, the nottion we shll use to define the clss of lnguges known s regulr sets. Recll
More informationProduct of polynomials. Introduction to Programming (in C++) Numerical algorithms. Product of polynomials. Product of polynomials
Product of polynomils Introduction to Progrmming (in C++) Numericl lgorithms Jordi Cortdell, Ricrd Gvldà, Fernndo Orejs Dept. of Computer Science, UPC Given two polynomils on one vrile nd rel coefficients,
More informationScanner Termination. Multi Character Lookahead. to its physical end. Most parsers require an end of file token. Lex and Jlex automatically create an
Scnner Termintion A scnner reds input chrcters nd prtitions them into tokens. Wht hppens when the end of the input file is reched? It my be useful to crete n Eof pseudo-chrcter when this occurs. In Jv,
More information