Spectral Methods for Time-dependent Variable-coefficient PDE Based on Block Gaussian Quadrature
|
|
- Erick McCoy
- 5 years ago
- Views:
Transcription
1 Spectrl Methods for Time-dependent Vrible-coefficient PDE Bsed on Block Gussin Qudrture Jmes V. Lmbers Abstrct Block Krylov subspce spectrl (KSS) methods re best-of-both-worlds compromise between explicit nd implicit time-stepping methods for vrible-coefficient PDE, in tht they combine the efficiency of explicit methods nd the stbility of implicit methods, while lso chieving spectrl ccurcy in spce nd high-order ccurcy in time. Block KSS methods compute ech Fourier coefficient of the solution using techniques developed by Gene Golub nd Gérrd Meurnt for pproximting elements of functions of mtrices by block Gussin qudrture in the spectrl, rther thn physicl, domin. This pper describes how block KSS methods cn be pplied to vriety of equtions, nd lso demonstrtes their superiority, in terms of ccurcy nd efficiency, to other Krylov subspce methods in the literture. Keywords: spectrl methods, Gussin qudrture, block Lnczos method, Mxwell s eqution, het eqution 1 Introduction In [11 clss of methods, clled block Krylov subspce spectrl (KSS) methods, ws introduced for the purpose of solving prbolic vrible-coefficient PDE. These methods re bsed on techniques developed by Golub nd Meurnt in [4 for pproximting elements of function of mtrix by Gussin qudrture in the spectrl domin. In [12, these methods were generlized to the secondorder wve eqution, for which these methods hve exhibited even higher-order ccurcy. It hs been shown in these references tht KSS methods, by employing different pproximtions of the solution opertor for ech Fourier coefficient of the solution, chieve higher-order ccurcy in time thn other Krylov subspce methods (see, for exmple, [9) for stiff systems of ODE, nd they re lso quite stble, considering tht they re Submitted Februry 19, The University of Southern Mississippi, Deprtment of Mthemtics, Httiesburg, MS USA Tel/Fx: /5818 Emil: Jmes.Lmbers@usm.edu explicit methods. They re lso effective for solving systems of coupled equtions, such s Mxwell s equtions [16. In this pper, we review block KSS methods, s pplied to vrious types of PDE, nd compre their performnce to other Krylov subspce methods from the literture. Section 2 reviews the min properties of block KSS methods, s pplied to the prbolic problems for which they were designed. Section 3 discusses implementtion detils, nd demonstrtes why KSS methods need to explicitly generte only one Krylov subspce, lthough informtion from severl is used. In Section 4, we discuss modifictions tht must be mde to block KSS methods in order to pply them to systems of coupled wve equtions, such s Mxwell s equtions. Numericl results re presented in Section 5, nd conclusions re stted in Section 6. 2 Krylov Subspce Spectrl Methods We first review block KSS methods, which re esier to describe for prbolic problems. Let S(t) = exp[ Lt represent the exct solution opertor of the problem u t + Lu = 0, t > 0, (1) with pproprite initil conditions nd periodic boundry conditions. The opertor L is second-order, self-djoint, positive definite differentil opertor. Let, denote the stndrd inner product of functions defined on [0, 2π. Block Krylov subspce spectrl methods, introduced in [11, use Gussin qudrture on the spectrl domin to compute the Fourier coefficients of the solution. These methods re time-stepping lgorithms tht compute the solution t time t 1, t 2,..., where t n = nδt for some choice of Δt. Given the computed solution u(x, t n ) t time t n, the solution t time t n+1 is computed by pproximting the Fourier coefficients tht would be obtined by pplying the exct solution opertor to u(x, t n ), 1 ˆu(ω, t n+1 ) = e iωx, S(Δt) u(x, t n ). (2) 2π
2 In [4 Golub nd Meurnt describe method for computing quntities of the form u T f(a)v, (3) where u nd v re N-vectors, A is n N N symmetric positive definite mtrix, nd f is smooth function. Our gol is to pply this method with A = L N where L N is spectrl discretiztion of L, f(λ) = exp( λt) for some t, nd the vectors u nd v re obtined from ê ω nd u n 1, where ê ω is discretiztion of 2π e iωx nd u n is the pproximte solution t time t n, evluted on n N-point uniform grid. The bsic ide is s follows: since the mtrix A is symmetric positive definite, it hs rel eigenvlues b = λ 1 λ 2 λ N = > 0, (4) nd corresponding orthogonl eigenvectors q j, j = 1,..., N. Therefore, the quntity (3) cn be rewritten s N u T f(a)v = f(λ j )u T q j q T j v. (5) which cn lso be viewed s Riemnn-Stieltjes integrl u T f(a)v = I[f = f(λ) dα(λ). (6) As discussed in [4, the integrl I[f cn be pproximted using Gussin qudrture rules, which yields n pproximtion of the form I[f = K w j f(λ j ) + R[f, (7) where the nodes λ j, j = 1,..., K, s well s the weights w j, j = 1,..., K, cn be obtined using the symmetric Lnczos lgorithm if u = v, nd the unsymmetric Lnczos lgorithm if u = v (see [7). In the cse u = v, there is possibility tht the weights my not be positive, which destbilizes the qudrture rule (see [1 for detils). Insted, we consider [ u v T f(a) [ u v, (8) which results in the 2 2 mtrix [ u f(λ) dμ(λ) = T f(a)u u T f(a)v v T f(a)u v T f(a)v, (9) where μ(λ) is 2 2 mtrix function of λ, ech entry of which is mesure of the form α(λ) from (6). In [4 Golub nd Meurnt showed how block method cn be used to generte qudrture formuls. We will describe this process here in more detil. The integrl f(λ) dμ(λ) is now 2 2 symmetric mtrix nd the most generl K-node qudrture formul is of the form f(λ) dμ(λ) = K W j f(t j )W j + error, (10) with T j nd W j being symmetric 2 2 mtrices. By digonlizing ech T j, we obtin the simpler formul f(λ) dμ(λ) = 2K f(λ j )v j v T j + error, (11) where, for ech j, λ j is sclr nd v j is 2-vector. Ech node λ j is n eigenvlue of the mtrix M 1 B1 T B 1 M 2 B T 2 T K = , (12) B K 2 M K 1 BK 1 T B K 1 M K which is block-tringulr mtrix of order 2K. The vector v j consists of the first two elements of the corresponding normlized eigenvector. To compute the mtrices M j nd B j, we use the block Lnczos lgorithm, which ws proposed by Golub nd Underwood in [6. We re now redy to describe block KSS methods. For ech wve number ω = N/2 + 1,..., N/2, we define R 0 (ω) = [ ê ω u n nd compute the QR fctoriztion R 0 (ω) = X 1 (ω)b 0 (ω). We then crry out block Lnczos itertion, pplied to the discretized opertor L N, to obtin block tridigonl mtrix T K (ω) of the form (12), where ech entry is function of ω. Then, we cn express ech Fourier coefficient of the pproximte solution t the next time step s [û n+1 ω = [ B0 H E12 H exp[ T K (ω)δte 12 B 0 (13) where E 12 = [ e 1 e 2. The computtion of (13) consists of computing the eigenvlues nd eigenvectors of T K (ω) in order to obtin the nodes nd weights for Gussin qudrture, s described erlier. This lgorithm hs locl temporl ccurcy O(Δt 2K 1 ) [11. Furthermore, block KSS methods re more ccurte thn the originl KSS methods described in [14, even though they hve the sme order of ccurcy, becuse the solution u n plys greter role in the determintion of the qudrture nodes. They re lso more effective 12
3 for problems with oscilltory or discontinuous coefficients [11. Block KSS methods re even more ccurte for the second-order wve eqution, for which block Lnczos itertion is used to compute both the solution nd its time derivtive. In [12, Theorem 6, it is shown tht when the leding coefficient is constnt nd the coefficient q(x) is bndlimited, the 1-node KSS method, which hs secondorder ccurcy in time, is lso unconditionlly stble. In generl, s shown in [12, the locl temporl error is O(Δt 4K 2 ) when K block Gussin nodes re used. 3 Implementtion KSS methods compute Jcobi mtrix corresponding to ech Fourier coefficient, in contrst to trditionl Krylov subspce methods tht normlly use only single Krylov subspce generted by the initil dt or the solution from the previous time step. While it would pper tht KSS methods incur substntil mount of dditionl computtionl expense, tht is not ctully the cse, becuse nerly ll of the Krylov subspces tht they compute re closely relted by the wve number ω, in the 1-D cse, or ω = (ω 1, ω 2,..., ω n ) in the n-d cse. In fct, the only Krylov subspce tht is explicitly computed is the one generted by the solution from the previous time step, of dimension (K + 1), where K is the number of block Gussin qudrture nodes. In ddition, the verges of the coefficients of L j, for j = 0, 1, 2,..., 2K 1, re required, where L is the sptil differentil opertor. When the coefficients of L re independent of time, these cn be computed once, during preprocessing step. This computtion cn be crried out in O(N log N) opertions using symbolic clculus [13, 15. With these considertions, the lgorithm for single time step of 1-node block KSS method for solving (1), where Lu = pu xx + q(x)u, with pproprite initil conditions nd periodic boundry conditions, is s follows. We denote the verge of function f(x) on [0, 2π by f, nd the computed solution t time t n by u n. ˆu n = fft(u n ), v = Lu n, ˆv = fft(v) for ech ω do α 1 = pω 2 + q (in preprocessing step) β 1 = ˆv(ω) α 1ˆu n (ω) α 2 = u n, v 2 Re [ˆu n (ω)v(ω) + α 1 u n (ω) 2 e ω = [ u [ n, u n ˆu n (ω) 2 1/2 α 1 β 1 /e ω T ω = β 1 /e ω α 2 /e 2 ω ˆu n+1 (ω) = [e TωΔt 11ˆu n (ω) + [e TωΔt 12 e ω end u n+1 = ifft(ˆu n+1 ) It should be noted tht for prbolic problem such s (1), the loop over ω only needs to ccount for nonnegligible Fourier coefficients of the solution, which re reltively few due to the smoothness of solutions to such problems. 4 Appliction to Mxwell s Equtions We consider Mxwell s eqution on the cube [0, 2π 3, with periodic boundry conditions. Assuming nonconductive mteril with no losses, we hve div Ê = 0, div Ĥ = 0, (14) curl Ê = μ Ĥ t, curl Ĥ = ε Ê t, (15) where Ê, Ĥ re the vectors of the electric nd mgnetic fields, nd ε, μ re the electric permittivity nd mgnetic permebility, respectively. Tking the curl of both sides of (15) yields με 2 Ê t 2 = ΔÊ + μ 1 curl Ê μ, (16) με 2 Ĥ t 2 = ΔĤ + ε 1 curl Ĥ ε. (17) In this section, we discuss generliztions tht must be mde to block KSS methods in order to pply them to non-self-djoint system of coupled equtions such s (16). Additionl detils re given in [16. First, we consider the following 1-D problem, 2 u + Lu = 0, t > 0, (18) t2 with pproprite initil conditions, nd periodic boundry conditions, where u : [0, 2π [0, ) R n for n > 1, nd L(x, D) is n n n mtrix where the (i, j) entry is n differentil opertor L ij (x, D) of the form m ij L ij (x, D)u(x) = ij μ (x)d μ u, μ=0 D = d dx, (19) with sptilly vrying coefficients ij μ, μ = 0, 1,..., m ij. Generliztion of KSS methods to system of the form (18) cn proceed s follows. For i, j = 1,..., n, let L ij (D) be the constnt-coefficient opertor obtined by verging the coefficients of L ij (x, D) over [0, 2π. Then, for ech wve number ω, we define L(ω) be the mtrix with entries L ij (ω), i.e., the symbols of L ij (D) evluted t ω.
4 Next, we compute the spectrl decomposition of L(ω) for ech ω. For j = 1,..., n, let q j (ω) be the Schur vectors of L(ω). Then, we define our test nd tril functions by φ j,ω (x) = q j (ω) e iωx. For Mxwell s equtions, the mtrix A N tht discretizes the opertor AÊ = 1 ) (ΔÊ με + μ 1 curl Ê μ is not symmetric, nd for ech coefficient of the solution, the resulting qudrture nodes λ j, j = 1,..., 2K, from (11) re now complex nd must be obtined by strightforwrd modifiction of block Lnczos itertion for unsymmetric mtrices. 5 Numericl Results In this section, we compre the performnce of block KSS methods with vrious methods bsed on exponentil integrtors [8, 10, Prbolic Problems We first consider 1-D prbolic problem of the form (1), where the differentil opertor L is defined by Lu(x) = pu (x) + q(x)u(x), where p 0.4 nd q(x) cos x 0.02 sin x cos 2x sin 2x cos 3x is constructed so s to hve the smoothness of function with three continuous derivtives, s is the initil dt u(x, 0). Periodic boundry conditions re imposed. Since the exct solution is not vilble, the error is estimted by tking the l 2 -norm of the reltive difference between ech solution, nd tht of solution computed using smller time step Δt = 1/64 nd the mximum number of grid points. The results re shown in Figures 1 nd 2. As the number of grid points is doubled, only the block KSS method shows n improvement in ccurcy; the preconditioned Lnczos method exhibits slight degrdtion in performnce, while the explicit fourth-order exponentil integrtor-bsed method requires tht the time step be reduced by fctor of 4 before it cn deliver the expected order of convergence; similr behvior ws demonstrted for n explicit 3rd-order method from [9 in [14. The preconditioned Lnczos method requires 8 Lnczos itertions to mtch the ccurcy of block KSS method tht uses only 2. On the other hnd, the block KSS method incurs dditionl expense due to (1) the computtion of the moments of L, for ech Fourier coefficient, nd (2) the exponentition of seprte Jcobi mtrices for ech Fourier coefficient. These expenses re mitigted by the fct tht the first tkes plce once, during preprocessing stge, nd both tsks require n mount of work tht is proportionl not to the number of grid points, but to the number of non-negligible Fourier coefficients of the solution. We solve this problem using the following methods: A 2-node block KSS method. Ech time step requires construction of Krylov subspce of dimension 3 generted by the solution, nd the coefficients of L 2 nd L 3 re computed during preprocessing step. A preconditioned Lnczos itertion for pproximting e τa v, introduced in [17 for pproximting the mtrix exponentil of sectoril opertors, nd dpted in [19 for efficient ppliction to the solution of prbolic PDE. In this pproch, Lnczos itertion is pplied to (I+hA) 1, where h is prmeter, in order to obtin restricted rtionl pproximtion of the mtrix exponentil. We use m = 4 nd m = 8 Lnczos itertions, nd choose h = Δt/10, s in [19. A method bsed on exponentil integrtors, from [8, tht is of order 3 when the Jcobin is pproximted to within O(Δt). We use m = 8 Lnczos itertions. Figure 1: Estimtes of reltive error t t = 0.1 in solutions of (1) computed using preconditioned exponentil integrtor [19 with 4 nd 8 Lnczos itertions, nd 2- node block KSS method. All methods compute solutions on n N-point grid, with time step Δt, for vrious vlues of N nd Δt.
5 in these experiments, we use m = 2 Lnczos itertions to pproximte the Jcobin. It should be noted tht when m is incresed, even to substntil degree, the results re negligibly ffected. Figure 3 demonstrtes the convergence behvior for both methods. At both sptil resolutions, the block KSS method exhibits pproximtely 6th-order ccurcy in time s Δt decreses, except tht for N = 16, the sptil error rising from trunction of Fourier series is significnt enough tht the overll error fils to decrese below the level chieved t Δt = 1/8. For N = 32, the solution is sufficiently resolved in spce, nd the order of overgence s Δt 0 is pproximtely 6.1. Figure 2: Estimtes of reltive error t t = 0.1 in solutions of (1) computed using 4th-order method bsed on n exponentil integrtor [10 with 8 Lnczos itertions, nd 2-node block KSS method. Both methods compute solutions on n N-point grid, with time step Δt, for vrious vlues of N nd Δt. 5.2 Mxwell s Equtions We now pply 2-node block KSS method to (16), with initil conditions Ê Ê(x, y, z, 0) = F(x, y, z), (x, y, z, 0) = G(x, y, z), t (20) with periodic boundry conditions. The coefficients μ nd ε re given by μ(x, y, z) = cos z cos y sin y cos(y + z) cos(y z) cos x cos(x z) sin(x z) cos(x + y) cos(x y), ε(x, y, z) = cos z cos y cos(y + z) cos x cos(x z) cos(x + y) cos(x y). The components of F nd G re generted in similr fshion, except tht the x- nd z-components re zero. We use block KSS method tht uses K = 2 block qudrture nodes per coefficient in the bsis described in Section 4, tht is 6th-order ccurte in time, nd cosine method bsed on Gutschi-type exponentil integrtor [8, 10. This method is second-order in time, nd We lso note tht incresing the resolution does not pose ny difficulty from stbility point of view. Unlike explicit finite-difference schemes tht re constrined by CFL condition, KSS methods do not require reduction in the time step to offset reduction in the sptil step in order to mintin boundedness of the solution, becuse their domin of dependence includes the entire sptil domin for ny Δt. The Gutschi-type exponentil integrtor method is second-order ccurte, s expected, nd delivers nerly identicl results for both sptil resolutions, but even with Krylov subspce of much higher dimension thn tht used in the block KSS method, it is only ble to chieve t most second-order ccurcy, wheres block KSS method, using Krylov subspce of dimension 3, chieves sixth-order ccurcy. This is due to the incorportion of the moments of the sptil differentil opertor into the computtion, nd the use of Gussin qudrture rules specificlly tilored to ech Fourier coefficient. 6 Summry nd Future Work We hve demonstrted tht block KSS methods cn be pplied to Mxwell s equtions with smoothly vrying coefficients, by pproprite generliztion of their ppliction to the sclr second-order wve eqution, in wy tht preserves the order of ccurcy chieved for the wve eqution. Furthermore, it hs been demonstrted tht while trditionl Krylov subspce methods bsed on exponentil integrtors re most effective for prbolic problems, especilly when ided by preconditioning s in [19, KSS methods perform best when pplied to hyperbolic problems, in view of their much higher order of ccurcy. Future work will extend the pproch described in this pper to more relistic pplictions involving Mxwell s equtions by using symbol modifiction to efficiently implement perfectly mtched lyers (see [2), nd vrious techniques (see [3, 18) to effectively hndle discontinuous coefficients.
6 [9 Hochbruck, M., Lubich, C.: On Krylov Subspce Approximtions to the Mtrix Exponentil Opertor. SIAM Journl of Numericl Anlysis 34 (1997) [10 Hochbruck, M., Lubich, C., Selhofer, H.: Exponentil Integrtors for Lrge Systems of Differentil Equtions. SIAM Journl of Scientific Computing 19 (1998) [11 Lmbers, J. V.: Enhncement of Krylov Subspce Spectrl Methods by Block Lnczos Itertion. Electronic Trnsctions on Numericl Anlysis 31 (2008) Figure 3: Estimtes of reltive error t t = 1 in solutions of (16), (20) computed using cosine method bsed on Gutschi-type exponentil integrtor [8, 10 with 2 Lnczos itertions, nd 2-node block KSS method. Both methods compute solutions on n N 3 -point grid, with time step Δt, for vrious vlues of N nd Δt. References [1 Atkinson, K.: An Introduction to Numericl Anlysis, 2nd Ed. Wiley (1989) [2 Berenger, J,: A perfectly mtched lyer for the bsorption of electromgnetic wves. J. Comp. Phys. 114 (1994) [3 Gelb, A., Tnner, J.: Robust Reprojection Methods for the Resolution of the Gibbs Phenomenon. Appl. Comput. Hrmon. Anl. 20 (2006) [4 Golub, G. H., Meurnt, G.: Mtrices, Moments nd Qudrture. Proceedings of the 15th Dundee Conference, June-July 1993, Griffiths, D. F., Wtson, G. A. (eds.), Longmn Scientific & Technicl (1994) [5 Golub, G. H., Gutknecht, M. H.: Modified Moments for Indefinite Weight Functions. Numerische Mthemtik 57 (1989) [6 Golub, G. H., Underwood, R.: The block Lnczos method for computing eigenvlues. Mthemticl Softwre III, J. Rice Ed., (1977) [7 Golub, G. H, Welsch, J.: Clcultion of Guss Qudrture Rules. Mth. Comp. 23 (1969) [8 Hochbruck, M., Lubich, C.: A Gutschi-type method for oscilltory second-order differentil equtions, Numerische Mthemtik 83 (1999) [12 Lmbers, J. V.: An Explicit, Stble, High-Order Spectrl Method for the Wve Eqution Bsed on Block Gussin Qudrture. IAENG Journl of Applied Mthemtics 38 (2008) [13 Lmbers, J. V.: Krylov Subspce Spectrl Methods for the Time-Dependent Schrödinger Eqution with Non-Smooth Potentils. Numericl Algorithms 51 (2009) [14 Lmbers, J. V.: Krylov Subspce Spectrl Methods for Vrible-Coefficient Initil-Boundry Vlue Problems. Electronic Trnsctions on Numericl Anlysis 20 (2005) [15 Lmbers, J. V.: Prcticl Implementtion of Krylov Subspce Spectrl Methods. Journl of Scientific Computing 32 (2007) [16 Lmbers, J. V.: A Spectrl Time-Domin Method for Computtionl Electrodynmics. Advnces in Applied Mthemtics nd Mechnics 1 (2009) [17 Moret, I., Novti, P.: RD-rtionl pproximtion of the mtrix exponentil opertor. BIT 44 (2004) [18 Vllius, T., Honknen, M.: Reformultion of the Fourier nodl method with dptive sptil resolution: ppliction to multilevel profiles. Opt. Expr. 10(1) (2002) [19 vn den Eshof, J., Hochbruck, M.: Preconditioning Lnczos pproximtions to the mtrix exponentil. SIAM Journl of Scientific Computing 27 (2006)
A Spectral Time-Domain Method for Computational Electrodynamics
A Spectrl Time-Domin Method for Computtionl Electrodynmics Jmes V. Lmbers Abstrct Ever since its introduction by Kne Yee over forty yers go, the finite-difference time-domin (FDTD) method hs been widely-used
More informationMatrices, Moments and Quadrature, cont d
Jim Lmbers MAT 285 Summer Session 2015-16 Lecture 2 Notes Mtrices, Moments nd Qudrture, cont d We hve described how Jcobi mtrices cn be used to compute nodes nd weights for Gussin qudrture rules for generl
More informationA Spectral Time-Domain Method for Computational Electrodynamics
A Spectral Time-Domain Method for Computational Electrodynamics James V. Lambers Abstract Block Krylov subspace spectral (KSS) methods have previously been applied to the variable-coefficient heat equation
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 informationPhysics 116C Solution of inhomogeneous ordinary differential equations using Green s functions
Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner
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 informationFredholm Integral Equations of the First Kind Solved by Using the Homotopy Perturbation Method
Int. Journl of Mth. Anlysis, Vol. 5, 211, no. 19, 935-94 Fredholm Integrl Equtions of the First Kind Solved by Using the Homotopy Perturbtion Method Seyyed Mhmood Mirzei Deprtment of Mthemtics, Fculty
More informationConstruction of Gauss Quadrature Rules
Jim Lmbers MAT 772 Fll Semester 2010-11 Lecture 15 Notes These notes correspond to Sections 10.2 nd 10.3 in the text. Construction of Guss Qudrture Rules Previously, we lerned tht Newton-Cotes qudrture
More informationAn approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95
An pproximtion to the rithmetic-geometric men G.J.O. Jmeson, Mth. Gzette 98 (4), 85 95 Given positive numbers > b, consider the itertion given by =, b = b nd n+ = ( n + b n ), b n+ = ( n b n ) /. At ech
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 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 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 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 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 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 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 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 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 informationMath 270A: Numerical Linear Algebra
Mth 70A: Numericl Liner Algebr Instructor: Michel Holst Fll Qurter 014 Homework Assignment #3 Due Give to TA t lest few dys before finl if you wnt feedbck. Exercise 3.1. (The Bsic Liner Method for Liner
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 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 information1 Linear Least Squares
Lest Squres Pge 1 1 Liner Lest Squres I will try to be consistent in nottion, with n being the number of dt points, nd m < n being the number of prmeters in model function. We re interested in solving
More informationTravelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing
Applied Mthemtics E-Notes 8(8) - c IN 67-5 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ Trvelling Profile olutions For Nonliner Degenerte Prbolic Eqution And Contour Enhncement In Imge
More informationA Nonclassical Collocation Method For Solving Two-Point Boundary Value Problems Over Infinite Intervals
Austrlin Journl of Bsic nd Applied Sciences 59: 45-5 ISS 99-878 A onclssicl Colloction Method For Solving o-oint Boundry Vlue roblems Over Infinite Intervls M Mlei nd M vssoli Kni Deprtment of Mthemtics
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 informationMATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1
MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further
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 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 informationImproper Integrals, and Differential Equations
Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted
More informationModule 11.4: nag quad util Numerical Integration Utilities. Contents
Qudrture Module Contents Module 11.4: ng qud util Numericl Integrtion Utilities ng qud util provides utility procedures for computtion involving integrtion of functions, e.g., the computtion of the weights
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 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 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 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 informationMath 6395: Hyperbolic Conservation Laws and Numerical Methods. Hyperbolic Conservation Laws are a type of PDEs arise from physics: fluid/gas dynamics.
Mth 6395: Hyperbolic Conservtion Lws nd Numericl Methods Course structure Hyperbolic Conservtion Lws re type of PDEs rise from physics: fluid/gs dynmics. Distintive solution structures Other types of PDEs
More informationMath 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of
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 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 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 informationCLOSED EXPRESSIONS FOR COEFFICIENTS IN WEIGHTED NEWTON-COTES QUADRATURES
Filomt 27:4 (2013) 649 658 DOI 10.2298/FIL1304649M Published by Fculty of Sciences nd Mthemtics University of Niš Serbi Avilble t: http://www.pmf.ni.c.rs/filomt CLOSED EXPRESSIONS FOR COEFFICIENTS IN WEIGHTED
More informationPart IB Numerical Analysis
Prt IB Numericl Anlysis Theorems with proof Bsed on lectures by G. Moore Notes tken by Dexter Chu Lent 2016 These notes re not endorsed by the lecturers, nd I hve modified them (often significntly) fter
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 informationBest Approximation. Chapter The General Case
Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given
More informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More informationThe Velocity Factor of an Insulated Two-Wire Transmission Line
The Velocity Fctor of n Insulted Two-Wire Trnsmission Line Problem Kirk T. McDonld Joseph Henry Lbortories, Princeton University, Princeton, NJ 08544 Mrch 7, 008 Estimte the velocity fctor F = v/c nd the
More informationTangent Line and Tangent Plane Approximations of Definite Integral
Rose-Hulmn Undergrdute Mthemtics Journl Volume 16 Issue 2 Article 8 Tngent Line nd Tngent Plne Approximtions of Definite Integrl Meghn Peer Sginw Vlley Stte University Follow this nd dditionl works t:
More informationNumerical Linear Algebra Assignment 008
Numericl Liner Algebr Assignment 008 Nguyen Qun B Hong Students t Fculty of Mth nd Computer Science, Ho Chi Minh University of Science, Vietnm emil. nguyenqunbhong@gmil.com blog. http://hongnguyenqunb.wordpress.com
More informationLecture Note 9: Orthogonal Reduction
MATH : Computtionl Methods of Liner Algebr 1 The Row Echelon Form Lecture Note 9: Orthogonl Reduction Our trget is to solve the norml eution: Xinyi Zeng Deprtment of Mthemticl Sciences, UTEP A t Ax = A
More informationChapter 0. What is the Lebesgue integral about?
Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous
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 informationMonte Carlo method in solving numerical integration and differential equation
Monte Crlo method in solving numericl integrtion nd differentil eqution Ye Jin Chemistry Deprtment Duke University yj66@duke.edu Abstrct: Monte Crlo method is commonly used in rel physics problem. The
More informationSolutions of Klein - Gordan equations, using Finite Fourier Sine Transform
IOSR Journl of Mthemtics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 13, Issue 6 Ver. IV (Nov. - Dec. 2017), PP 19-24 www.iosrjournls.org Solutions of Klein - Gordn equtions, using Finite Fourier
More informationACM 105: Applied Real and Functional Analysis. Solutions to Homework # 2.
ACM 05: Applied Rel nd Functionl Anlysis. Solutions to Homework # 2. Andy Greenberg, Alexei Novikov Problem. Riemnn-Lebesgue Theorem. Theorem (G.F.B. Riemnn, H.L. Lebesgue). If f is n integrble function
More informationDefinition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim
Mth 9 Course Summry/Study Guide Fll, 2005 [1] Limits Definition of Limit: We sy tht L is the limit of f(x) s x pproches if f(x) gets closer nd closer to L s x gets closer nd closer to. We write lim f(x)
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 4 Contravariance, Covariance, and Spacetime Diagrams
Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz
More informationChapter 7 Notes, Stewart 8e. 7.1 Integration by Parts Trigonometric Integrals Evaluating sin m x cos n (x) dx...
Contents 7.1 Integrtion by Prts................................... 2 7.2 Trigonometric Integrls.................................. 8 7.2.1 Evluting sin m x cos n (x)......................... 8 7.2.2 Evluting
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties
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 informationGoals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite
Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct
More informationSTURM-LIOUVILLE BOUNDARY VALUE PROBLEMS
STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS Throughout, we let [, b] be bounded intervl in R. C 2 ([, b]) denotes the spce of functions with derivtives of second order continuous up to the endpoints. Cc 2
More informationOrdinary differential equations
Ordinry differentil equtions Introduction to Synthetic Biology E Nvrro A Montgud P Fernndez de Cordob JF Urchueguí Overview Introduction-Modelling Bsic concepts to understnd n ODE. Description nd properties
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 informationConservation Law. Chapter Goal. 5.2 Theory
Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very
More informationCAAM 453 NUMERICAL ANALYSIS I Examination There are four questions, plus a bonus. Do not look at them until you begin the exam.
Exmintion 1 Posted 23 October 2002. Due no lter thn 5pm on Mondy, 28 October 2002. Instructions: 1. Time limit: 3 uninterrupted hours. 2. There re four questions, plus bonus. Do not look t them until you
More informationRecitation 3: More Applications of the Derivative
Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech
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 informationPartial Derivatives. Limits. For a single variable function f (x), the limit lim
Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the right-hnd side limit equls to the left-hnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles
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 informationMAC-solutions of the nonexistent solutions of mathematical physics
Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE
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 informationSection 4.8. D v(t j 1 ) t. (4.8.1) j=1
Difference Equtions to Differentil Equtions Section.8 Distnce, Position, nd the Length of Curves Although we motivted the definition of the definite integrl with the notion of re, there re mny pplictions
More informationDensity of Energy Stored in the Electric Field
Density of Energy Stored in the Electric Field Deprtment of Physics, Cornell University c Tomás A. Aris October 14, 01 Figure 1: Digrm of Crtesin vortices from René Descrtes Principi philosophie, published
More informationContents. Outline. Structured Rank Matrices Lecture 2: The theorem Proofs Examples related to structured ranks References. Structure Transport
Contents Structured Rnk Mtrices Lecture 2: Mrc Vn Brel nd Rf Vndebril Dept. of Computer Science, K.U.Leuven, Belgium Chemnitz, Germny, 26-30 September 2011 1 Exmples relted to structured rnks 2 2 / 26
More information1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ),
1. Guss-Jcobi qudrture nd Legendre polynomils Simpson s rule for evluting n integrl f(t)dt gives the correct nswer with error of bout O(n 4 ) (with constnt tht depends on f, in prticulr, it depends on
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 informationChapter 3. Vector Spaces
3.4 Liner Trnsformtions 1 Chpter 3. Vector Spces 3.4 Liner Trnsformtions Note. We hve lredy studied liner trnsformtions from R n into R m. Now we look t liner trnsformtions from one generl vector spce
More informationOverview of Calculus I
Overview of Clculus I Prof. Jim Swift Northern Arizon University There re three key concepts in clculus: The limit, the derivtive, nd the integrl. You need to understnd the definitions of these three things,
More informationA Convergence Theorem for the Improper Riemann Integral of Banach Space-valued Functions
Interntionl Journl of Mthemticl Anlysis Vol. 8, 2014, no. 50, 2451-2460 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/10.12988/ijm.2014.49294 A Convergence Theorem for the Improper Riemnn Integrl of Bnch
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 informationEfficient Computation of a Class of Singular Oscillatory Integrals by Steepest Descent Method
Applied Mthemticl Sciences, Vol. 8, 214, no. 31, 1535-1542 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/1.12988/ms.214.43166 Efficient Computtion of Clss of Singulr Oscilltory Integrls by Steepest Descent
More information1 2-D Second Order Equations: Separation of Variables
Chpter 12 PDEs in Rectngles 1 2-D Second Order Equtions: Seprtion of Vribles 1. A second order liner prtil differentil eqution in two vribles x nd y is A 2 u x + B 2 u 2 x y + C 2 u y + D u 2 x + E u +
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 informationQuadrature Rules for Evaluation of Hyper Singular Integrals
Applied Mthemticl Sciences, Vol., 01, no. 117, 539-55 HIKARI Ltd, www.m-hikri.com http://d.doi.org/10.19/ms.01.75 Qudrture Rules or Evlution o Hyper Singulr Integrls Prsnt Kumr Mohnty Deprtment o Mthemtics
More informationChapters 4 & 5 Integrals & Applications
Contents Chpters 4 & 5 Integrls & Applictions Motivtion to Chpters 4 & 5 2 Chpter 4 3 Ares nd Distnces 3. VIDEO - Ares Under Functions............................................ 3.2 VIDEO - Applictions
More informationMath 31S. Rumbos Fall Solutions to Assignment #16
Mth 31S. Rumbos Fll 2016 1 Solutions to Assignment #16 1. Logistic Growth 1. Suppose tht the growth of certin niml popultion is governed by the differentil eqution 1000 dn N dt = 100 N, (1) where N(t)
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =
More informationCBE 291b - Computation And Optimization For Engineers
The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn
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 informationMAT 168: Calculus II with Analytic Geometry. James V. Lambers
MAT 68: Clculus II with Anlytic Geometry Jmes V. Lmbers Februry 7, Contents Integrls 5. Introduction............................ 5.. Differentil Clculus nd Quotient Formuls...... 5.. Integrl Clculus nd
More informationc n φ n (x), 0 < x < L, (1) n=1
SECTION : Fourier Series. MATH4. In section 4, we will study method clled Seprtion of Vribles for finding exct solutions to certin clss of prtil differentil equtions (PDEs. To do this, it will be necessry
More informationJim Lambers MAT 169 Fall Semester Lecture 4 Notes
Jim Lmbers MAT 169 Fll Semester 2009-10 Lecture 4 Notes These notes correspond to Section 8.2 in the text. Series Wht is Series? An infinte series, usully referred to simply s series, is n sum of ll of
More informationMATH 144: Business Calculus Final Review
MATH 144: Business Clculus Finl Review 1 Skills 1. Clculte severl limits. 2. Find verticl nd horizontl symptotes for given rtionl function. 3. Clculte derivtive by definition. 4. Clculte severl derivtives
More informationNumerical Integration. 1 Introduction. 2 Midpoint Rule, Trapezoid Rule, Simpson Rule. AMSC/CMSC 460/466 T. von Petersdorff 1
AMSC/CMSC 46/466 T. von Petersdorff 1 umericl Integrtion 1 Introduction We wnt to pproximte the integrl I := f xdx where we re given, b nd the function f s subroutine. We evlute f t points x 1,...,x n
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More informationThe Algebra (al-jabr) of Matrices
Section : Mtri lgebr nd Clculus Wshkewicz College of Engineering he lgebr (l-jbr) of Mtrices lgebr s brnch of mthemtics is much broder thn elementry lgebr ll of us studied in our high school dys. In sense
More information221B Lecture Notes WKB Method
Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using
More informationImplicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments
Implicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments James V. Lambers September 24, 2008 Abstract This paper presents a reformulation
More informationBest Approximation in the 2-norm
Jim Lmbers MAT 77 Fll Semester 1-11 Lecture 1 Notes These notes correspond to Sections 9. nd 9.3 in the text. Best Approximtion in the -norm Suppose tht we wish to obtin function f n (x) tht is liner combintion
More information