Numerical integration methods
|
|
- MargaretMargaret Beasley
- 6 years ago
- Views:
Transcription
1 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 crgo, nd rigged. This is becuse it ws impossible to clculte the wter displced volume, i.e. the integrl of the hull-describing function, with the rudimentry mth known t tht time. Consequently, there ws no wy to properly estimte the buoynt force. Thus, the loction of the wterline ws unknown until the ship ws fully flot. Additionlly, not every integrl cn be computed nlyticlly, even for reltively simple functions. Exmple The Guss error function defined s erf(y) = 2 π y 0 e x2 dx cnnot be clculted using only elementry functions. 1.1 Integrtion problem sttement At first, we consider the evlution of one dimensionl integrl, lso clled qudrture, since this opertion is equivlent to finding of the re under the curve of given function. b f(x)dx Surprisingly, one does not need to know ny high level mth to do so. All you need is precision scle nd scissors. The recipe goes like this: plot the function on some preferbly 1
2 dense nd hevy mteril, cut the re of interest with scissors, mesure the mss of the cutout, divide the obtined mss by the density nd thickness of the mteril nd you re done. Well, this is, of course, not very precise nd sounds so low tech. We will employ more modern numericl methods. Words of wisdom Once you become proficient with computers, there is tendency to do every problem numericlly. Resist it! If you cn find n nlyticl solution for problem, do so. It is usully fster to clculte, nd more importntly, will provide you some physicl intuition for the problem overll. Numericl solutions usully do not possess predictive power. 1.2 The Rectngle method The definition of the integrl vi Riemnn s sum: b N 1 f(x)dx = lim f(x i )(x i+1 x i ) (1.1) N i=1 where x i b nd N is the number of points, f(x) f 1 f 2 f 3 f 4 f 5 f 6 f i f N-1 f N X 1 X 2 X 3 X 4 X 5 X 6 X i X N-1 X N h b x Figure 1.1: The rectngle method illustrtion. Shded boxes show the pproximtions of the re under the curve. Here, f 1 = f(x 1 = ), f 2 = f(x 2 ), f 3 = f(x 3 ),..., f N = f(x N = b). Riemnn s definition gives us directions for the rectngle method: pproximte the re under the curve by set of boxes or rectngles (see fig. 1.1). For simplicity, we evlute our function t N points eqully seprted by the distnce h on the intervl (, b). 2
3 Rectngle method integrl estimte where b f(x)dx N 1 i=1 f(x i )h, (1.2) h = b N 1, x i = + (i 1)h, x 1 = nd x N = b (1.3) The MATLAB implementtion of this method is quite simple: Listing 1.1: integrte_in_1d.m (vilble t with_matlab_book/ch_integrtion/code/integrte_in_1d.m) function integrl1d = integrte_in_1d(f,, b) % integrtion with simple rectngle/box method % int_^b f(x) dx N=100; % number of points in the sum x=linspce(,b,n); s=0; for xi=x(1:end 1) % we need to exclude x(end)=b s = s + f(xi); end %% now we clculte the integrl integrl1d = s*(b )/(N 1); To demonstrte how to use it, let s check 1 0 x2 dx = 1/3 f x.^2; integrte_in_1d(f,0,1) ns = Well, is quite close to the expected vlue of 1/3. The smll devition from the exct result is due to the reltively smll number of points; we used N = 100. If you hve previous experience in low level lnguges (from the rry functions implementtion point of view) like C or Jv, the bove implementtion is the first thing tht comes to your mind. While it is correct, it does not use MATLAB s bility to use mtrices s rguments of function. A better wy, which voids using loop, is shown below. Listing 1.2: integrte_in_1d_mtlb_wy.m (vilble t progrmming_with_matlab_book/ch_integrtion/code/integrte_in_1d_mtlb_wy. 3
4 m) function integrl1d = integrte_in_1d_mtlb_wy(f,, b) % integrtion with simple rectngle/box method % int_^b f(x) dx N=100; % number of points in the sum x=linspce(,b,n); % if function f cn work with vector rgument then we cn do integrl1d = (b )/(N 1)*sum( f( x(1:end 1) ) ); % we exclude x(end)=b Words of wisdom In MATLAB, loops re generlly slower compred to the equivlent evlution of function with vector or mtrix rgument. Try to void loops which iterte evlution over mtrix elements Rectngle method lgorithmic error While the rectngle method is very simple to understnd nd implement, it hs wful lgorithmic error nd slow convergence. Let s see why this is so. A closer look t fig. 1.1 revels tht box often underestimtes (see, for exmple, the box between x 1 nd x 2 ) or overestimtes (for exmple, the box between x 6 nd x 7 ) the re. This is the lgorithmic error which is proportionl to f h 2 /2 to the first order, i.e. the re of tringle between box nd the curve. Since we hve N 1 intervls, this error ccumultes in the worst cse scenrio. So we cn sy tht the lgorithmic error (E) is Rectngle method lgorithmic error estimte ) ( ) E = O ((N 1) h2 (b ) 2 f 2 = O 2N f (1.4) In the lst term, we replced N 1 with N under the ssumption tht N is lrge. The O symbol represents the big O nottion, i.e. there is n unknown proportionlity coefficient. 1.3 Trpezoidl method The pproximtion of ech intervl with trpezoid (see fig. 1.2) is n ttempt to circumvent the wekness of the rectngle method. In other words, we do liner pproximtion of the 4
5 f(x) f 1 f 2 f 3 f 4 f 5 f 6 f i f N-1 f N X 1 X 2 X 3 X 4 X 5 X 6 X i X N-1 X N h b x Figure 1.2: The illustrtion of the trpezoidl method. Shded trpezoids show the pproximtions of the re under the curve. integrted function. Reclling the formul for the re of trpezoid, we pproximte the re of ech intervl s h(f i+1 + f i )/2, nd then we sum them ll. Trpezoidl method integrl estimte b f(x)dx h ( 1 2 f 1 + f 2 + f f N 2 + f N f N) (1.5) The 1/2 coefficient disppers for inner points, since they re counted twice: once for the left nd once for the right trpezoid re. Trpezoidl method lgorithmic error To evlute the lgorithmic error, one should note tht we re using liner pproximtion for the function nd ignoring the second order term. Reclling the Tylor expnsion, this mens tht to the first order, we re ignoring contribution of the second derivtive (f ) terms. With bit of ptience, it cn be shown tht Trpezoidl method lgorithmic error estimte ( ) (b ) 3 E = O 12N f 2 (1.6) Let s compre the integrl estimte with rectngle (eq. (1.2)) nd trpezoidl (eq. (1.5)) 5
6 methods. b b f(x)dx h (f 2 + f f N 2 + f N 1 ) + h (f 1 ), (1.7) f(x)dx h (f 2 + f f N 2 + f N 1 ) + h 1 2 (f 1 + f N ) (1.8) It is esy to see tht they re lmost identicl, with the only difference in the second term of either h (f 1 ) for the rectngle method or h 1 2 (f 1 + f N ) for the trpezoidl one. While this might seem like minor chnge in the underlying formul, it results in the lgorithmic error decresing s 1/N 2 for the trpezoidl method, which is wy better thn the 1/N dependence for the rectngle method. 1.4 Simpson s method f(x) f 1 f 2 f 3 f 4 f 5 f 6 f N-2 f i f N-1 f N X 1 X 2 X 3 X 4 X 5 X 6 X i X N-2 X N-1 h b X N x Figure 1.3: Simpson s method illustrtion The next logicl step is to pproximte the function with second order curves, i.e. prbols (see fig. 1.3). This leds us to the Simpson method. Let s consider triplet of consequent points (x i 1, f i 1 ), (x i, f i ), nd (x i+1, f i+1 ). The re under the prbol pssing through ll of these points is h/3 (f i 1 + 4f i + f i+1 ). Then we sum the res of ll triplets nd obtin Simpson s method integrl estimte b f(x)dx h 1 3 (f 1 + 4f 2 + 2f 3 + 4f f N 2 + 4f N 1 + f N ) (1.9) N must be in specil form N = 2k + 1, i.e. odd, nd 3 Yet gin, the first (f 1 ) nd lst (f N ) points re specil, since they re counted only once, while the edge points f 3, f 5,... re counted twice s members of the left nd right triplets. 6
7 Simpson s method lgorithmic error Since we re using more terms of the function s Tylor expnsion, the convergence of the method is improving nd the error of the method decreses with N even fster. Simpson method lgorithmic error estimte ( ) (b ) 5 E = O 180N f (4) 4 (1.10) Perhps one would be surprised to see f (4) nd N 4. This is becuse when we integrte triplet re, we go by the h to the left nd to the right of the centrl point, so the terms proportionl to the x 3 f (3) re cnceled out. 1.5 Generlized formul for integrtion A creful reder my hve lredy noticed tht the integrtion formuls for ll bove methods cn be written in the sme generl form. Generlized formul for numericl integrtion b where w i is the weight coefficient. f(x)dx h i=1 N f(x i )w i, (1.11) So, one hs no excuse to use the rectngle nd trpezoid method over Simpson s, since the clcultionl burden is exctly the sme, but the clcultion error for Simpson s method drops drsticlly s the number of points increses. Strictly speking, even Simpson s method is not so superior in comprison to others which use even higher order polynomil pproximtion for the function. These higher order methods would hve exctly the sme form s eq. (1.11). The difference will be in the weight coefficients. One cn see more detiled discussion of this issue in [1]. 7
8 y b y y x x b x Figure 1.4: The pond re estimte vi the Monte Crlo method. 1.6 Monte Crlo Integrtion Toy exmple: finding the re of pond Suppose we wnt to estimte the re of pond. Nively, we might rush to grb mesuring tpe. However, we could just wlk round nd throw stones in every possible direction, nd drw n imginry box round the pond to see how mny stones lnded inside (N inside ) the pond out of their totl number (N totl ). As illustrted in fig. 1.4, the rtio of bove numbers should be proportionl to the rtio of the pond re to the box re. So, we conclude tht the estimted re of the pond (A pond ) A pond = N inside N totl A box (1.12) where A box = (b x x )(b y y ) is the box re. The bove estimte requires tht the thrown stones re distributed rndomly nd uniformly 1. It is lso good ide to keep the re of the surrounding box s tight s possible to increse the probbility of hitting the pond. Imgine wht would hppen if the box ws huge compred to the pond size: we would need lot of stones (while we hve only finite mount of them) to hit the pond even once, nd until then the pond re estimte would be zero in ccordnce with eq. (1.12). 1 This is not trivil tsk. We will tlk bout this in chpter 11. For now, we will use rnd function provided by MATLAB. 8
9 1.6.1 Nive Monte Crlo integrtion y H f(x) b x Figure 1.5: Estimte of the integrl by counting points under the curve nd in the surrounding box. Well, the clcultion of the re under the curve is not much different from the mesurement of the pond re s shown in fig So bx x f(x)dx = N inside N totl A box Monte Crlo integrtion derived The bove method is not optiml. Let s focus our ttention on nrrow strip round the point x b (see fig. 1.6). For this strip N inside N totl H f(x b ) (1.13) So, there is no need to wste ll of these resources if we cn get the f(x b ) right wy. Thus, the improved estimte of the integrl by the Monte Crlo method looks like 9
10 H y f(x ) b x b b x Figure 1.6: Explntion for the improved Monte Crlo method The Monte Crlo method integrl estimte Choose N uniformly nd rndomly distributed points x i inside [, b] b f(x)dx b N N f(x i ) (1.14) i=1 The Monte Crlo method lgorithmic error A creful look t eq. (1.14) shows tht 1/N N i=1 f(x i) =< f > is ctully sttisticl estimte of the function men. Sttisticlly, we re only sure in this estimte within the stndrd devition of the men estimte. This brings us to the following expression. 10
11 Monte Crlo method lgorithmic error estimte ( b ) E = O < f 2 > < f > 2 N (1.15) where < f > = 1 N < f 2 > = 1 N N f(x i ) i=1 N f 2 (x i ) i=1 The right most term under the squre root is the estimte of the function stndrd devition σ = < f 2 > < f > 2. Looking t the bove expression, you might conclude tht you hve wsted precious minutes of your life reding bout very mediocre method, which reduces its error proportionlly to 1/ N. This is worse even compred to the otherwise wful rectngle method. Hold your horses. In the following section, we will see how the Monte Crlo method outshines ll others for the cse of multidimensionl integrtion. 1.7 Multidimensionl integrtion If someone sks us to clculte multidimensionl integrl, we just need to pply our knowledge of how to del with single dimension integrls. For exmple, in the two dimensionl cse, we simply rerrnge terms: bx by x y f(x, y) dx dy = bx The lst single dimension integrl is the function of only x: x by dx dy f(x, y) (1.16) y by y dy f(x, y) = F (x) (1.17) So, the two-dimensionl integrl boils down to the chin of two one-dimension ones, which we re lredy fit to process: bx by x y f(x, y) dx dy = bx x dx F (x) (1.18) 11
12 Miniml exmple for integrtion in 2D The listing below shows how to do 2D integrls by chining the one-dimension ones; note tht it piggybcks on the single dimensionl integrl in listing 1.1 (but feel free to use ny other method). Listing 1.3: integrte_in_2d.m (vilble t with_matlab_book/ch_integrtion/code/integrte_in_2d.m) function integrl2d=integrte_in_2d(f, xrnge, yrnge) % Integrtes function f in 2D spce % f is hndle to function of x, y i.e. f(x,y) should be vlid % xrnge is vector contining lower nd upper limits of integrtion % long the first dimension. % xrnge = [x_lower x_upper] % yrnge is similr but for the second dimension % We will define (Integrl f(x,y) dy) s Fx(x) Fx f(x,y), yrnge(1), yrnge(2) ); % ^^^^^ we fix 'x', ^^^^^^^^^^^here we reuse this lredy fixed x % so it reds s Fy(y) % This is quite cumbersome. % It is probbly impossible to do generl D dimensionl cse. % Notice tht mtlb folks implemented integrl, integrl2, integrl3 % but they did not do ny for higher thn 3 dimensions. integrl2d = integrte_in_1d(fx, xrnge(1), xrnge(2) ); end Let s clculte 2 dx (2x 2 + y 2 ) dy (1.19) f 2*x.^2 + y.^2; integrte_in_2d( f, [0,2], [0,1] ) ns = It is esy to see tht the exct nswer is 6. The observed devition from the nlyticl result is due to the smll number of points used in the clcultion. 12
13 1.8 Multidimensionl integrtion with Monte Crlo The bove "chin" method cn be expnded to ny number of dimensions. Cn we rest now? Not so fst. Note tht if we would like to split the integrtion region by N points in every of the D dimensions, then the number of evlutions, nd thus the clcultion time, grows N D. This renders the rectngle, trpezoidl, Simpson s, nd like methods useless for high-dimension integrls. The Monte Crlo method is notble exception; it looks very simple even for multidimensionl cse, mintins the sme N evlution time, nd its error is still 1/ N. A 3D cse, for exmple, would look like this bx x by bz dx dy dz f(x, y, z) (b x x )(b y y )(b z z ) y z N nd the generl form is shown below Monte Crlo method in D-spce N f(x i, y i, z i ) (1.20) i=1 dv D f( x) = dx 1 dx 2 dx 3...dx D f( x) V D V D V D N N f( x i ) (1.21) i=1 where V D is the D-dimensionl volume, x i rndomly nd uniformly distributed points in the volume V D Monte Crlo method demonstrtion To see how elegnt nd simple the implementtion of the Monte Crlo method cn be, we will evlute the integrl in eq. (1.19). f 2*x.^2 + y.^2; bx=2; x=0; by=1; y=0; % first we prepre x nd y components of rndom points % rnd provides uniformly distributed points in the (0,1) intervl N=1000; % Number of rndom points x=x+(bx x)*rnd(1,n); % 1 row, N columns y=y+(by y)*rnd(1,n); % 1 row, N columns % finlly integrl evlution integrl2d = (bx x)*(by y)/n * sum( f(x,y) ) 13
14 integrl2d = We used only 1000 points nd the result is quite close to the nlyticl vlue of Numericl integrtion gotchs Using very lrge number of points Since the lgorithmic error of numericl integrtion methods drops with n increse of N, it is very tempting to increse it. But we must remember bout round-off errors nd resist this tempttion. So, h should not be too smll or equivlently N should not be too big. N definitely should not be infinite s the Riemnn s sum in eq. (1.1) prescribes, since our lifetime is finite. Using too few points sin(x) π x, rdins Figure 1.7: The exmple of bdly chosen points ( ) for integrtion of the sin(x) function ( ). The smples give the impression tht the function is horizontl line. 14
15 There is dnger of under smpling if the integrted function chnges very quickly nd we put points very sprsely (see for exmple fig. 1.7). We should first plot the function (if it is not too computtionlly txing), then choose pproprite mount of points: there should be t lest 2 nd idelly more points for every vlley nd hill of the function 2. To see if you hit sweet spot, try to double the mount of points nd see if the estimte of the integrl stops chnging drsticlly. This is the foundtion for the so clled dptive integrtion method, where n lgorithm utomticlly decides on the required number of points MATLAB functions for integrtion Below is short summry for MATLAB s built-in functions for integrtion: 1D integrtion integrl trpz employs modified trpezoidl method qud employs modified Simpson s method Here is how to use qud to clculte 2 0 3x3 dx >> f 3*x.^2; >> qud(f, 0, 2) ns = As expected, the nswer is 8. Note tht qud expects your function to be vector friendly. Multidimensionl integrtion integrl2 2D cse integrl3 3D cse Let s clculte the integrl 1 0 >> f x.^2 + y.^2; >> integrl2( f, 0, 1, 0, 1) ns = (x2 + y 2 )dxdy which equls to 2/3. There re mny others built-ins. See MATLAB s numericl integrtion documenttion to lern more. 2 We will discuss it in detil in the section devoted to the discrete Fourier trnsform (chpter 15). 15
16 MATLAB s implementtions re more powerful thn those which we discussed, but deep inside, they use similr methods Self-Study Generl comments: Do not forget to run some test cses. MATLAB hs built-in numericl integrtion methods, such s qud. You might check the vlidity of your implementtions with nswers produced by this MATLAB built-in function. qud requires your function to be ble to work with n rry of x points, otherwise it will fil. Of course, it is lwys better to compre to the exct nlyticlly clculted vlue. Problem 1: Implement the trpezoidl numericl integrtion method. Cll your function trpezint(f,,b,n), where nd b re left nd right limits of integrtion, N is the number of points, nd f is hndle to the function. Problem 2: Implement the Simpson numericl integrtion method. simpsonint(f,,b,n). Remember bout specil form of N=2k+1. Cll your function Problem 3: Implement the Monte-Crlo numericl integrtion method. Cll your function montecrloint(f,,b,n). Problem 4: For your tests clculte 10 0 [exp( x) + (x/1000) 3 ]dx Plot the integrl bsolute error from its true vlue for the bove methods (include rectngulr method s well) vs. different number of points N. Try to do it from smll N=3 to N=10 6. Use loglog plotting function for better representtion (mke sure tht you hve enough points in ll res of the plot). Cn you relte the trends to eqs. (1.4), (1.6), (1.9) nd (1.15). Why does the error strt to grow with lrger N? Does it grow for ll methods? Why? Problem 5: Clculte π/2 0 sin(401x)dx Compre your result with the exct nswer 1/401. number of points to clculte this integrl. Provide discussion bout required 16
17 Problem 6: Clculte 1 1 dx 1 0 dy(x 4 + y 2 + x y) Compre results of the Monte-Crlo numericl integrtion method nd MATLAB s built-in integrl2. Problem 7: Implement method to find the volume of the N-dimensionl sphere for the rbitrry dimension N nd the sphere rdius R V (N, R) = dx 1 dx 2 dx 3 dx N x 2 1 +x2 2 +x x2 N R2 Clculte the volume of the sphere with R = 1 nd N =
18 Bibliogrphy [1] W. H. Press, S. A. Teukolsky, W. T. Vetterling, nd B. P. Flnnery. Numericl Recipes 3rd Edition: The Art of Scientific Computing. Cmbridge University Press, New York, NY, USA, 3 edition,
P(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 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 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 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 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 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 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 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 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 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 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 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 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 informationFig.1. Let a source of monochromatic light be incident on a slit of finite width a, as shown in Fig. 1.
Answer on Question #5692, Physics, Optics Stte slient fetures of single slit Frunhofer diffrction pttern. The slit is verticl nd illuminted by point source. Also, obtin n expression for intensity distribution
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 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 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 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 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 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 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 informationApproximate computations
Living with floting-point numers Stndrd normlized representtion (sign + frction + exponent): Approximte computtions Rnges of vlues: Representtions for:, +, +0, 0, NN (not numer) Jordi Cortdell Deprtment
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 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 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 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 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 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 information1 Quad-Edge Construction Operators
CS48: Computer Grphics Hndout # Geometric Modeling Originl Hndout #5 Stnford University Tuesdy, 8 December 99 Originl Lecture #5: 9 November 99 Topics: Mnipultions with Qud-Edge Dt Structures Scribe: Mike
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 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 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 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 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 informationINTRODUCTION TO SIMPLICIAL COMPLEXES
INTRODUCTION TO SIMPLICIAL COMPLEXES CASEY KELLEHER AND ALESSANDRA PANTANO 0.1. Introduction. In this ctivity set we re going to introduce notion from Algebric Topology clled simplicil homology. The min
More informationApproximation by NURBS with free knots
pproximtion by NURBS with free knots M Rndrinrivony G Brunnett echnicl University of Chemnitz Fculty of Computer Science Computer Grphics nd Visuliztion Strße der Ntionen 6 97 Chemnitz Germny Emil: mhrvo@informtiktu-chemnitzde
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 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 informationArea & Volume. Chapter 6.1 & 6.2 September 25, y = 1! x 2. Back to Area:
Bck to Are: Are & Volume Chpter 6. & 6. Septemer 5, 6 We cn clculte the re etween the x-xis nd continuous function f on the intervl [,] using the definite integrl:! f x = lim$ f x * i )%x n i= Where fx
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 informationSIMPLIFYING ALGEBRA PASSPORT.
SIMPLIFYING ALGEBRA PASSPORT www.mthletics.com.u This booklet is ll bout turning complex problems into something simple. You will be ble to do something like this! ( 9- # + 4 ' ) ' ( 9- + 7-) ' ' Give
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 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 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 informationToday. Search Problems. Uninformed Search Methods. Depth-First Search Breadth-First Search Uniform-Cost Search
Uninformed Serch [These slides were creted by Dn Klein nd Pieter Abbeel for CS188 Intro to AI t UC Berkeley. All CS188 mterils re vilble t http://i.berkeley.edu.] Tody Serch Problems Uninformed Serch Methods
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 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 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 informationCHAPTER III IMAGE DEWARPING (CALIBRATION) PROCEDURE
CHAPTER III IMAGE DEWARPING (CALIBRATION) PROCEDURE 3.1 Scheimpflug Configurtion nd Perspective Distortion Scheimpflug criterion were found out to be the best lyout configurtion for Stereoscopic PIV, becuse
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 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 informationIterated Integrals. f (x; y) dy dx. p(x) To evaluate a type I integral, we rst evaluate the inner integral Z q(x) f (x; y) dy.
Iterted Integrls Type I Integrls In this section, we begin the study of integrls over regions in the plne. To do so, however, requires tht we exmine the importnt ide of iterted integrls, in which inde
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 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 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 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 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 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 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 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 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 informationMid-term exam. Scores. Fall term 2012 KAIST EE209 Programming Structures for EE. Thursday Oct 25, Student's name: Student ID:
Fll term 2012 KAIST EE209 Progrmming Structures for EE Mid-term exm Thursdy Oct 25, 2012 Student's nme: Student ID: The exm is closed book nd notes. Red the questions crefully nd focus your nswers on wht
More informationAgenda & Reading. Class Exercise. COMPSCI 105 SS 2012 Principles of Computer Science. Arrays
COMPSCI 5 SS Principles of Computer Science Arrys & Multidimensionl Arrys Agend & Reding Agend Arrys Creting & Using Primitive & Reference Types Assignments & Equlity Pss y Vlue & Pss y Reference Copying
More informationA New Learning Algorithm for the MAXQ Hierarchical Reinforcement Learning Method
A New Lerning Algorithm for the MAXQ Hierrchicl Reinforcement Lerning Method Frzneh Mirzzdeh 1, Bbk Behsz 2, nd Hmid Beigy 1 1 Deprtment of Computer Engineering, Shrif University of Technology, Tehrn,
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 information50 AMC LECTURES Lecture 2 Analytic Geometry Distance and Lines. can be calculated by the following formula:
5 AMC LECTURES Lecture Anlytic Geometry Distnce nd Lines BASIC KNOWLEDGE. Distnce formul The distnce (d) between two points P ( x, y) nd P ( x, y) cn be clculted by the following formul: d ( x y () x )
More informationRay surface intersections
Ry surfce intersections Some primitives Finite primitives: polygons spheres, cylinders, cones prts of generl qudrics Infinite primitives: plnes infinite cylinders nd cones generl qudrics A finite primitive
More informationInternational Journal of Scientific & Engineering Research, Volume 4, Issue 12, December ISSN
Interntionl Journl of Scientific & Engineering Reserch, Volume 4, Issue 1, December-1 ISSN 9-18 Generlised Gussin Qudrture over Sphere K. T. Shivrm Abstrct This pper presents Generlised Gussin qudrture
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 informationCSCI 446: Artificial Intelligence
CSCI 446: Artificil Intelligence Serch Instructor: Michele Vn Dyne [These slides were creted by Dn Klein nd Pieter Abbeel for CS188 Intro to AI t UC Berkeley. All CS188 mterils re vilble t http://i.berkeley.edu.]
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 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 informationSimplifying Algebra. Simplifying Algebra. Curriculum Ready.
Simplifying Alger Curriculum Redy www.mthletics.com This ooklet is ll out turning complex prolems into something simple. You will e le to do something like this! ( 9- # + 4 ' ) ' ( 9- + 7-) ' ' Give this
More informationIntroduction Transformation formulae Polar graphs Standard curves Polar equations Test GRAPHS INU0114/514 (MATHS 1)
POLAR EQUATIONS AND GRAPHS GEOMETRY INU4/54 (MATHS ) Dr Adrin Jnnett MIMA CMth FRAS Polr equtions nd grphs / 6 Adrin Jnnett Objectives The purpose of this presenttion is to cover the following topics:
More informationCS 221: Artificial Intelligence Fall 2011
CS 221: Artificil Intelligence Fll 2011 Lecture 2: Serch (Slides from Dn Klein, with help from Sturt Russell, Andrew Moore, Teg Grenger, Peter Norvig) Problem types! Fully observble, deterministic! single-belief-stte
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 informationStained Glass Design. Teaching Goals:
Stined Glss Design Time required 45-90 minutes Teching Gols: 1. Students pply grphic methods to design vrious shpes on the plne.. Students pply geometric trnsformtions of grphs of functions in order to
More informationAgilent Mass Hunter Software
Agilent Mss Hunter Softwre Quick Strt Guide Use this guide to get strted with the Mss Hunter softwre. Wht is Mss Hunter Softwre? Mss Hunter is n integrl prt of Agilent TOF softwre (version A.02.00). Mss
More informationAnswer Key Lesson 6: Workshop: Angles and Lines
nswer Key esson 6: tudent Guide ngles nd ines Questions 1 3 (G p. 406) 1. 120 ; 360 2. hey re the sme. 3. 360 Here re four different ptterns tht re used to mke quilts. Work with your group. se your Power
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 information3.5.1 Single slit diffraction
3..1 Single slit diffrction ves pssing through single slit will lso diffrct nd produce n interference pttern. The reson for this is to do with the finite width of the slit. e will consider this lter. Tke
More informationAngle Properties in Polygons. Part 1 Interior Angles
2.4 Angle Properties in Polygons YOU WILL NEED dynmic geometry softwre OR protrctor nd ruler EXPLORE A pentgon hs three right ngles nd four sides of equl length, s shown. Wht is the sum of the mesures
More informationAI Adjacent Fields. This slide deck courtesy of Dan Klein at UC Berkeley
AI Adjcent Fields Philosophy: Logic, methods of resoning Mind s physicl system Foundtions of lerning, lnguge, rtionlity Mthemtics Forml representtion nd proof Algorithms, computtion, (un)decidility, (in)trctility
More informationLecture 10 Evolutionary Computation: Evolution strategies and genetic programming
Lecture 10 Evolutionry Computtion: Evolution strtegies nd genetic progrmming Evolution strtegies Genetic progrmming Summry Negnevitsky, Person Eduction, 2011 1 Evolution Strtegies Another pproch to simulting
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 informationClass-XI Mathematics Conic Sections Chapter-11 Chapter Notes Key Concepts
Clss-XI Mthemtics Conic Sections Chpter-11 Chpter Notes Key Concepts 1. Let be fixed verticl line nd m be nother line intersecting it t fixed point V nd inclined to it t nd ngle On rotting the line m round
More information4.4 Improper Integrals
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
More informationCS143 Handout 07 Summer 2011 June 24 th, 2011 Written Set 1: Lexical Analysis
CS143 Hndout 07 Summer 2011 June 24 th, 2011 Written Set 1: Lexicl Anlysis In this first written ssignment, you'll get the chnce to ply round with the vrious constructions tht come up when doing lexicl
More informationEssential Question What are some of the characteristics of the graph of a rational function?
8. TEXAS ESSENTIAL KNOWLEDGE AND SKILLS A..A A..G A..H A..K Grphing Rtionl Functions Essentil Question Wht re some of the chrcteristics of the grph of rtionl function? The prent function for rtionl functions
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 informationLists in Lisp and Scheme
Lists in Lisp nd Scheme Lists in Lisp nd Scheme Lists re Lisp s fundmentl dt structures, ut there re others Arrys, chrcters, strings, etc. Common Lisp hs moved on from eing merely LISt Processor However,
More informationAlgorithm Design (5) Text Search
Algorithm Design (5) Text Serch Tkshi Chikym School of Engineering The University of Tokyo Text Serch Find sustring tht mtches the given key string in text dt of lrge mount Key string: chr x[m] Text Dt:
More informationMisrepresentation of Preferences
Misrepresenttion of Preferences Gicomo Bonnno Deprtment of Economics, University of Cliforni, Dvis, USA gfbonnno@ucdvis.edu Socil choice functions Arrow s theorem sys tht it is not possible to extrct from
More informationCSEP 573 Artificial Intelligence Winter 2016
CSEP 573 Artificil Intelligence Winter 2016 Luke Zettlemoyer Problem Spces nd Serch slides from Dn Klein, Sturt Russell, Andrew Moore, Dn Weld, Pieter Abbeel, Ali Frhdi Outline Agents tht Pln Ahed Serch
More information1 Drawing 3D Objects in Adobe Illustrator
Drwing 3D Objects in Adobe Illustrtor 1 1 Drwing 3D Objects in Adobe Illustrtor This Tutoril will show you how to drw simple objects with three-dimensionl ppernce. At first we will drw rrows indicting
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 informationA Heuristic Approach for Discovering Reference Models by Mining Process Model Variants
A Heuristic Approch for Discovering Reference Models by Mining Process Model Vrints Chen Li 1, Mnfred Reichert 2, nd Andres Wombcher 3 1 Informtion System Group, University of Twente, The Netherlnds lic@cs.utwente.nl
More informationGeometric transformations
Geometric trnsformtions Computer Grphics Some slides re bsed on Shy Shlom slides from TAU mn n n m m T A,,,,,, 2 1 2 22 12 1 21 11 Rows become columns nd columns become rows nm n n m m A,,,,,, 1 1 2 22
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 information