Transport Calculations. Tseelmaa Byambaakhuu, Dean Wang*, and Sicong Xiao
|
|
- Damon Howard
- 6 years ago
- Views:
Transcription
1 A Locl hp Adptive Diffusion Synthetic Accelertion Method for Neutron Trnsport Clcultions Tseelm Bymbkhuu, Den Wng*, nd Sicong Xio University of Msschusetts Lowell, 1 University Ave, Lowell, MA USA *Emil: Den_Wng@uml.edu Number of pges: 21 Number of figures: 7 Number of tbles: 1 1
2 ABSTRACT We present locl hp dptive diffusion synthetic ccelertion (DSA) method for neutron trnsport clcultions. This new DSA method, clled DG-DSA, solves the diffusion eqution on corse mesh using the interior penlty discontinuous Glerkin methods. We investigte vrious numericl spects of the DG-DSA method such s convergence, locl hp dpttion, nd penlty size. We demonstrte tht our DG-DSA method cn effectively nd efficiently ccelerte trnsport source itertions. Keywords: Accelertion methods, DSA, DG-DSA 2
3 I. INTRODUCTION There hs been considerble reserch on diffusion synthetic ccelertion (DSA) of trnsport source itertions. 1-4 The DSA method is bsed on the use of diffusion clcultion for pproximting the itertive error of the trnsport source itertion. Most work ws focused on the development of so-clled consistent (or prtilly consistent) diffusion discretiztions in order to obtin the optiml convergence performnce. While the use of consistent discretiztions mkes the DSA method effective, it does not necessrily mke it efficient becuse the numericl solution of discretized elliptic diffusion problem itself cn be costly, prticulrly when the problem size becomes lrge. Adms nd Mrtin developed DSA scheme bsed on discontinuous finite element (DFE) diffusion discretiztion to ccelerte DFE trnsport itertions. 5 They noted tht their DFE bsed DSA scheme ws very effective, however, the efficient solution of the DFE diffusion equtions in 2D geometries remined n open question. In this pper, we present new discontinuous Glerkin (DG) discretiztion of the diffusion eqution, clled DG-DSA, which cn effectively nd efficiently ccelerte the S N trnsport itertions. The novelty of our method is tht the diffusion eqution is solved on corse mesh grid using the DG methods insted on the trnsport fine mesh, nd the DG diffusion discretiztion incorportes locl hp dpttion, i.e., locl dpttion of mesh size nd/or polynomil degree, bsed on locl totl cross section (or opticl thickness). Therefore, the resulting number of degrees of freedom of the DG 3
4 discretiztion is much less thn the conventionl consistent DSA discretiztions, nd thus DG-DSA cn chieve significnt improvement in computtionl efficiency. The remining pper is orgnized s follows. In Sec. II, we present in detil the formultion nd lgorithm of the DG-DSA method. A numericl study of DG-DSA is crried out in Sec. III, focusing on vrious numericl spects such s convergence performnce, locl hp dpttion, nd penlty. Sec. IV concludes the pper with brief summry nd discussion. II. DG-DSA FORMULATION AND ALGORITHM In this pper, we introduce the DG-DSA method bsed on monoenergetic S N neutron trnsport fixed source problems on 2D Crtesin geometry. The scttering nd neutron source re ssumed isotropic. The flowchrt of the DG-DSA lgorithm is shown in Fig. 1. 4
5 Fig. 1. Flowchrt of the DG-DSA lgorithm. The l &' itertion cycle begins with the S N trnsport eqution with itertion indices is expressed s μ ) )* ψ,-. / x, y, μ, η + η ) )5 ψ,-. / x, y, μ, η + Σ 7 ψ,-. / x, y, μ, η = : ; < φ, x, y + > < Q x, y, (1) where φ nd ψ re the sclr flux nd ngulr flux, respectively. Σ & nd re the totl cross section nd scttering cross section. μ nd η re the neutron ngulr directions. x 5
6 nd y re the sptil positions. Q is the externl neutron source. l is the source itertion index nd l + 1/2 is the intermedite step. During ech source itertion, the diffusion is utilized to pproximte the trnsport itertive flux error s ) )* > ) E: F )* δφ,->/h x, y ) )5 > ) E: F )5 δφ,->/h x, y + Σ I δφ,->/h x, y = Σ J [φ,->/h x, y φ, x, y ], (2) where Σ I is the bsorption cross section. We solve the bove diffusion eqution with priml discontinuous Glerkin method (DG) 6, which is commonly used method for solving elliptic nd prbolic problems. D δφ x, y + Σ I δφ x, y = f(x, y), in Ω, (3) Reflective BC: D δφ x, y n = 0, on Γ XYZ, (3b) Vcuum BC: D δφ x, y n = > δφ x, y, H on Γ [I\, (3c) where Ω is polygonl domin in R H, Γ XYZ nd Γ [I\ re the disjoint sets tht prtition the domin boundry. n is the unit norml vector to the boundry exterior to Ω. The functions D nd f re defined s D = > E: F, (4) f x, y = Σ J [φ,->/h x, y φ, x, y ]. (4b) Multiplied by test function υ nd integrted over one element E gives υ D δφ + Σ I δφ υ = f υ, E E, (5) 6
7 where E is the discretiztion of Ω, i.e. the mesh. Appling the Green s theorem on the first term υ D δφ = D δφ υ D δφ n E υ ), E E, (6) where n E denotes n outwrd norml vector to, E, the boundry of element E. Summing over ll the elements gives E (D δφ υ + Σ I δφ υ) Y h ijf D δφ n e υ D δφ Y Y h klm h opq Y n e υ = E fυ, (7) where Γ rs7 is set contining only interior edges nd n e is n outwrd unit norml vector of edge e. Two new opertors re introduced: jump nd verge Interior edge: υ = > H υ ḷ + υ / l, (8) υ = υ ḷ υ / l, e = E > Y E H Y, (8b) Boundry edge: υ = υ = υ ḷ, e = E > Y Ω, (8c) where Ω is the boundry of the domin Ω. For continuous nd second order differentible functions, the jump nd verge opertors re simpler δφ = 0, e Γ rs7, (9) D δφ n e = 0, e Γ rs7, (9b) D δφ n e = D δφ n e, e Γ rs7. (9c) 7
8 Using Eqs. (8) nd (9) on Eq. (7), it is strightforwrd to obtin the vritionl formultion (or wek formultion) of Eq. (3). E (D δφ υ + Σ I δφ υ) D δφ n e υ Y h ijf + Y > H Y h opq [δφ][υ] = fυ Y E. (10) Using the smoothness of the solution δφ expressed in Eq. 9, we cn dd the following two terms to the wek formultion ε Y h ijf D υ n e [δφ] + Y x l Y Y Y h ijf δφ [ υ] = 0. (11) where e is the edge length, σ Y is nonnegtive rel penlty number nd ε is nother prmeter tht my tke the vlue of {-1, 0, 1}. Thus, the vritionl formultion cn be rewritten s { δφ, υ = L(υ), (12) where { δφ, υ is biliner form nd L(υ) is liner form which re defined s: { δφ, υ = E (D δφ υ + Σ I δφ υ) D δφ n e υ Y h ijf + Y > H Y h opq δφ υ + ε D υ n e [δφ] Y Y h ijf + Y x l Y Y Y h ijf δφ [ υ], (13) L υ = E fυ. (13b) Depending on the choice of prmeter ε, the methods cn be clled differently. A detiled discussion of these three types of DG methods cn be in Reference [6]. ε = 1, Symmetric interior penlty Glerkin SIPG +1, Nonsymmetric interior penlty Glerkin NIPG 0, Incomplete interior penlty Glerkin IIPG. (14) 8
9 DG-DSA utilizes finite element spce D E, the spce of discontinuous polynomils. The globl bsis functions of D E hve support contined in ech element nd D E = spn{p r : 1 i N, E E}, (15) with P r x, y = p r x, y, (x, y) E 0, (x, y) E, (15b) where {p r } is set of locl bsis functions tht re chosen to be the monomil bsis functions, trnslted from the intervl (-1,1) for qudrilterl mesh: p r x, y = * * š.// H(* š. * š ) œ 5 5.// H(5. 5 ) ž, I + J = i, 0 i k, (15c) nd (x ->/H, y ->/H ) is the midpoint of the element E bounded by (x, x -> ) (y, y -> ). This yields the locl dimension N = ( ->)( -H) H where k is the highest polynomil degree of n element E., (15d) Finlly, the DG-DSA method is to find δφ,->/h (x, y) in D E such tht { δφ,->/h (x, y), υ = L υ, υ D E. (16) The solution δφ,->/h (x, y) cn be n expnsion of locl bsis functions s δφ,->/h x, y = rª> α r E P r (x, y), (17) where α r is the unknown rel numbers to be solved. Substituting Eq. (17) into Eq. (16) gives Aα = b, (18) 9
10 where α is the vector with components α r, b is the vector with components L P, nd A is sprse mtrix expressed s rª> A = { P r, P E, E E, 1 j N ±. (18b) At the end of the l &' source itertion, the sclr flux cn be updted in the next trnsport itertion s φ,-> x, y = φ,->/h x, y + δφ,->/h (x, y) (19) The trnsport source itertion will continue until the convergence criterion is stisfied. The DG-DSA method is locl dptive, which mens tht both the mesh size nd the polynomil degree re loclly djustble bsed on locl totl cross section. The penlty number, σ Y, is chosen to give the optimum performnce bsed on scoping nlysis, which will be discussed lter. III. NUMERICAL RESULTS III.A. Numericl Convergence Study A numericl study of the DG-DSA ccelertion performnce ws crried out bsed on 2D S N fixed source model problem, which is homogeneous 6cm 6cm squre with the reflective boundry on the left nd bottom sides nd the vcuum boundry on the top nd right sides. The domin is discretized into 5 5 uniform corse-mesh cells. The fine-mesh number in ech corse-mesh cell is The numericl solution for 10
11 the S N trnsport ws obtined on the fine mesh grid (60 60) using the Guss-Legendre S 12 qudrture set for ngulr discretiztion nd the dimond difference (DD) method for sptil discretiztion. The DG-DSA results were obtined on vrious corse mesh grids, which ws determined bsed on the totl cross sections. The corse mesh (5 5) ws used for smll cross sections (Σ 7 1 cm > ), nd the fine mesh (60 60) ws used for lrge cross sections (Σ 7 > 6 cm > ). For medium cross sections (1 < Σ 7 < 6 cm > ), the DG-DSA mesh size ws determined by mintining the opticl thickness (i.e., Σ & Δ, where Δ is the corse-mesh size) round 1.2. The symmetric interior penlty Glerkin method (SIPG) ws used with piecewise liner polynomils nd the optimized penlty number. The S 12 with DG-DSA ws implemented in MATLAB. In order to chrcterize the convergence behvior, we estimte the spectrl rdius numericlly s defined by ρ = lim, º»¼.»¼» ¼» ¼½.. (20) Note tht the convergence is rpid when ρ 1, nd it slows down when ρ increses. When ρ 1 the scheme fils to converge. Fig. 2 presents the numericl spectrl rdius of DG-DSA s function of totl cross section for vrious scttering rtios, c. It shows tht the DG-DSA method is very effective nd stble for wide rnge of totl cross sections (or opticl thickness). The generl trend is tht the convergence rte decreses with the increse in totl cross sections up to Σ cm > (i.e., the opticl thickness of 20), therefter the spectrl 11
12 rdius decreses becuse the S N solution tends to the diffusion limit. In ddition, the spectrl rdius increses with the scttering rtio in generl. 1 x = 1.2 cm x = 1.2/Σ 7 x = 0.1 cm Spectrl Rdius c = 0.6 c = 0.8 c = 0.9 c = S t (cm -1 ) Fig. 2. Spectrl rdius vs. Σ 7. The bove convergence nlysis cn be used to develop the locl mesh refinement strtegy for DG-DSA. For exmple, in typicl light wter rector the totl cross section of wter in fst neutron groups is less thn 1 cm >, nd it is lrger thn 1 cm > in therml groups. The totl cross section of the fuel is typiclly less thn 1 cm >. Therefore, for fst group trnsport clcultions the DG-DSA mesh cn be s lrge s fuel pin size (~1.2 cm) or even lrger. Only for therml groups, it requires reltively fine mesh. It is interesting to note tht in most neutronics trnsport codes the conventionl consistent DSA discretiztion uses the sme mesh s the trnsport discretiztion, which typiclly hs more thn 100 cells in fuel pin. However, our DG-DSA discretiztion cn hve less thn 20 cells to chieve the sme convergence performnce. An exmple of the conventionl DSA nd our DG-DSA mesh structures is shown Fig. 12
13 3. It shows tht the number of degrees of freedom of the DG-DSA method is much less thn those of the conventionl DSA methods. This is why our DG-DSA scheme cn be very efficient. A detiled comprison is given in Sec. III.C. Consistent DSA Mesh DG-DSA Mesh Fig. 3. Consistent DSA mesh vs. DG-DSA mesh. III.B. Locl p Adpttion We study the locl p (polynomil degree) dptivity of the DG-DSA method bsed on 2D monoenergetic trnsport fixed source problem with homogeneous cross section. Similr to the problem in Sec. III.A, the model problem is 6cm 6cm squre with the reflective boundry on the left nd bottom sides nd the vcuum boundry on the top nd right sides. The domin is divided into 5 5 uniform corse-mesh cells. The finemesh number in ech corse-mesh is The numericl S N trnsport solutions were obtined on the fine mesh (60 60) using DD for sptil discretiztion nd the Guss- Legendre S >H qudrture set for ngulr discretiztion. The DG-DSA solutions were obtined on the corse mesh (5 5), using piecewise constnt (P0), liner (P1), nd qudrtic (P2) polynomils, respectively. 13
14 The numericl results re presented in Fig. 4. The converged flux is shown in Fig. 4, nd the flux reltive error s function of trnsport itertion is shown in Fig. 4b. It shows tht the DG-DSA ccelertion scheme is more effective with higher polynomils, but the computtionl sving decreses with incresing polynomil degree (e.g., P1 vs. P2). The use of higher polynomils is more expensive becuse of lrger number of degrees of freedom. For nucler neutron trnsport problems, we recommend P1 for the region of lrge totl cross sections (Σ 7 > 1 cm > ), nd P1 or P0 for smll cross sections (Σ 7 < 1 cm > ). () Converged sclr flux 14
15 Flux Reltive Error 1.0E E E E E-08 SI P0 P1 P2 1.0E Trnsport Sweep # (b) Flux reltive error vs. Itertion number Fig. 4. Numericl results of locl p dpttion. III.C. Locl h Adpttion In this section, we solve the sme problem s bove with inhomogeneous cross sections s shown in Fig. 5. The problem hs 5 5 uniform corse-mesh cells, nd ech corse cell consists of fine-mesh cells. This cse is mimic of mini fuel ssembly, but there is dditionlly highly bsorbing region (in blue) t the loction (2, 2). The numericl solutions for the S N trnsport were obtined on the fine mesh (60 60) using the Guss-Legendre S 12 qudrture set for ngulr discretiztion nd the dimond difference (DD) method for sptil discretiztion. The DG-DSA solutions were obtined on both the fine-mesh (FM) nd corse-mesh (CM) grids. Piecewise liner polynomil functions were used for the DG-DSA solutions. It should be noted tht locl mesh refinement, 6 6, ws pplied to the bsorbing region (in blue), where the opticl thickness is lrge. The resulting mesh is nonconforming mesh with some 15
16 hnging nodes. The totl number of cells in the DG-DSA mesh is 252, while the S N mesh hs 3600 cells. The numericl results re shown in Fig. 6. The converged sclr flux is plotted in Fig. 6. The convergence performnce, i.e., the flux reltive error vs. trnsport sweep number, is illustrted in Fig. 6b. The results of the unccelerted S N source itertion re shown for comprison. It shows tht DG-DSA cn effectively converge the S N itertions on the corse mesh s on the fine mesh. Fig. 5. Specifictions of 2D problem. 16
17 () Converged sclr flux Flux Reltive Error 1.0E E E E E-08 SI FM DG-DSA CM DG-DSA 1.0E Trnsport Sweep # (b) Flux Reltive Error vs. Itertion number Fig. 6. Numericl results of locl h dpttion. The comprison of computing time is summrized in TABLE I for the bove problem. It shows tht DG-DSA is very effective when solved on consistent fine mesh, but computtionlly inefficient becuse the discretized liner diffusion system is very lrge with the number of degrees of freedom of (only 756 for the corse- 17
18 mesh solution). This cse demonstrtes tht the DG-DSA method cn effectively nd efficiently ccelerte the trnsport itertion by using corse mesh grid. In ddition, s compred with continuous Glerkin (CG) methods, n dvntge of the DG methods is tht they cn be discretized on nonconforming mesh which enbles flexible implementtion of locl hp dpttion, while the CG implementtion on the noncomforming mesh is much more involving nd complicted. It should be noted tht in this study the liner system of DG-DSA, i.e., Eq. (18), ws solved simply using the MATLAB built-in bckslsh function. TABLE I. Computtionl Performnce Comprison* SI FM DG-DSA CM DG-DSA Number of degrees of freedom Globl mtrix size b Trnsport itertion number Trnsport time (s) Globl mtrix ssembly time (s) c DG-DSA time (s) Totl clcultion time (s) Speedup *Computtionl results were obtined with MATLAB R2017 on McBook Pro with Processor 2.9 GHz Intel Core i7. SI: totl number of ngulr flux unknowns; DG-DSA: totl number of polynomil coefficients. b Mtrix A of Eq. (18b). c Computing time for clculting the entries of the globl mtrix. III.D. Penlty Number Our DG-DSA method is bsed on the primry DG pproch which employs interior penlty to stbilize the numericl solution. It is found in our model problems tht the 18
19 ccelertion performnce of DG-DSA is sensitive to the penlty number, σ Y, lthough it is numericlly stble for wide rnge of penlty. To chieve the optiml convergence performnce, the suggested penlty number for SIPG s function of totl cross section is depicted in Fig. 7. It shows tht for smll cross section (< 0.8 cm > ), the convergence of DG-DSA is not much sensitive to penlty, for which wide rnge of vlues cn be chosen. For 0.8 < Σ 7 < 1 cm >, the penlty cn be determined bsed on liner interpoltion between 1 nd 0.3. For Σ 7 > 1 cm >, constnt penlty of 0.3 cn be used. Nevertheless, for different implementtion it is suggested tht similr sensitive study be conducted to find the optiml rnge of the penlty number se St Fig. 7. Penlty number vs. totl cross section. IV. CONCLUSIONS In this pper, we hve presented our newly developed diffusion synthetic ccelertion method, DG-DSA, for speeding up the convergence of neutron trnsport clcultions. This new DSA method cn gretly improve the computtionl efficiency of the convectionl DSA methods by using the DG methods. The novelty of DG-DSA is 19
20 tht it reduces the number of degrees of freedom by discretizing the diffusion eqution on corse mesh with locl hp dption. Our numericl results hve demonstrted its rpid convergence performnce nd efficiency. In ddition, it is worth mentioning tht for LWR pplictions DG-DSA cn further employ mesh dpttion for different neutron energy group, i.e., the mesh for fst group clcultions cn be corser thn tht for therml groups (e.g., pin or qurter ssembly size) since the totl cross sections of wter nd fuel in fst groups re smller thn therml group cross sections. Finlly, the ppliction of DG-DSA for other unstructured mesh nd trnsport solvers will be studied in the future. ACKNOWLEDGEMENTS This reserch ws supported by the Deprtment of Energy Nucler Energy University Progrm(NEUP). REFERENCES 1. R. E. ALCOUFFE, Diffusion Synthetic Accelertion Methods for the Dimond- Differenced Discrete-Ordintes Equtions, Nucl. Sci. Eng., 64, 344 (1977). 2. E. W. LARSEN, Unconditionlly Stble Diffusion- Synthetic Accelertion Methods for the Slb Geometry Discrete Ordintes Equtions. Prt I: Theory, Nucl. Sci. Eng., 82, 47 (1982). 3. M. L. ADAMS nd E. W. LARSEN., Fst Itertive Methods for Discrete Ordintes Prticle Trnsport Clcultions, Prog. Nucl. Energy., 40, 3, (2002). 4. J. S. WARSA, T. A. WAREING, nd J. E. MOREL, Krylov Itertive Methods nd 20
21 the Degrded Effectiveness of Diffusion Synthetic Accelertion for Multidimensionl S N Clcultions in Problems with Mteril Discontinuities, Nucl. Sci. Eng., 147, 218 (2004). 5. M. L. ADAMS nd W. R. MARTIN, Diffusion Synthetic Accelertion of Discontinuous Finite Element Trnsport Itertions, Nucl. Sci. Eng., 111, 145 (1992). 6. B. Riviere, Discontinuous Glerkin Methods for Solving Elliptic nd Prbolic Equtions: Theory nd Implementtion, Society for Industril nd Applied Mthemtics, (2008). 21
13: Diffusion in 2 Energy Groups
3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups
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 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 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 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 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 informationNUMERICAL INTEGRATION
NUMERICAL INTEGRATION How do we evlute I = f (x) dx By the fundmentl theorem of clculus, if F (x) is n ntiderivtive of f (x), then I = f (x) dx = F (x) b = F (b) F () However, in prctice most integrls
More 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 informationIntroduction to Finite Element Method
Introduction to Finite Element Method Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pn.pl/ tzielins/ Tble of Contents 1 Introduction 1 1.1 Motivtion nd generl concepts.............
More informationReview of Gaussian Quadrature method
Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge
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 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 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 informationNumerical Methods I Orthogonal Polynomials
Numericl Methods I Orthogonl Polynomils Aleksndr Donev Cournt Institute, NYU 1 donev@cournt.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fll 2014 Nov 6th, 2014 A. Donev (Cournt Institute) Lecture IX
More informationEuler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )
Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s
More 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 informationChapter 3 Solving Nonlinear Equations
Chpter 3 Solving Nonliner Equtions 3.1 Introduction The nonliner function of unknown vrible x is in the form of where n could be non-integer. Root is the numericl vlue of x tht stisfies f ( x) 0. Grphiclly,
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 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 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 informationApplicable Analysis and Discrete Mathematics available online at
Applicble Anlysis nd Discrete Mthemtics vilble online t http://pefmth.etf.rs Appl. Anl. Discrete Mth. 4 (2010), 23 31. doi:10.2298/aadm100201012k NUMERICAL ANALYSIS MEETS NUMBER THEORY: USING ROOTFINDING
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 informationNew data structures to reduce data size and search time
New dt structures to reduce dt size nd serch time Tsuneo Kuwbr Deprtment of Informtion Sciences, Fculty of Science, Kngw University, Hirtsuk-shi, Jpn FIT2018 1D-1, No2, pp1-4 Copyright (c)2018 by The Institute
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 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 information1 Error Analysis of Simple Rules for Numerical Integration
cs41: introduction to numericl nlysis 11/16/10 Lecture 19: Numericl Integrtion II Instructor: Professor Amos Ron Scries: Mrk Cowlishw, Nthnel Fillmore 1 Error Anlysis of Simple Rules for Numericl Integrtion
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 informationLecture 20: Numerical Integration III
cs4: introduction to numericl nlysis /8/0 Lecture 0: Numericl Integrtion III Instructor: Professor Amos Ron Scribes: Mrk Cowlishw, Yunpeng Li, Nthnel Fillmore For the lst few lectures we hve discussed
More 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 informationWe partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.
Mth 255 - Vector lculus II Notes 4.2 Pth nd Line Integrls We begin with discussion of pth integrls (the book clls them sclr line integrls). We will do this for function of two vribles, but these ides cn
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 3 Polynomial Interpolation: Numerical Differentiation and Integration.
Advnced Computtionl Fluid Dynmics AA215A Lecture 3 Polynomil Interpoltion: Numericl Differentition nd Integrtion Antony Jmeson Winter Qurter, 2016, Stnford, CA Lst revised on Jnury 7, 2016 Contents 3 Polynomil
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 informationPENALIZED LEAST SQUARES FITTING. Manfred von Golitschek and Larry L. Schumaker
Serdic Mth. J. 28 2002), 329-348 PENALIZED LEAST SQUARES FITTING Mnfred von Golitschek nd Lrry L. Schumker Communicted by P. P. Petrushev Dedicted to the memory of our collegue Vsil Popov Jnury 4, 942
More informationLecture 13 - Linking E, ϕ, and ρ
Lecture 13 - Linking E, ϕ, nd ρ A Puzzle... Inner-Surfce Chrge Density A positive point chrge q is locted off-center inside neutrl conducting sphericl shell. We know from Guss s lw tht the totl chrge on
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 informationSection 6.1 INTRO to LAPLACE TRANSFORMS
Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform
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 informationPredict Global Earth Temperature using Linier Regression
Predict Globl Erth Temperture using Linier Regression Edwin Swndi Sijbt (23516012) Progrm Studi Mgister Informtik Sekolh Teknik Elektro dn Informtik ITB Jl. Gnesh 10 Bndung 40132, Indonesi 23516012@std.stei.itb.c.id
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 informationMAT 772: Numerical Analysis. James V. Lambers
MAT 772: Numericl Anlysis Jmes V. Lmbers August 23, 2016 2 Contents 1 Solution of Equtions by Itertion 7 1.1 Nonliner Equtions....................... 7 1.1.1 Existence nd Uniqueness................ 7 1.1.2
More informationLecture 12: Numerical Quadrature
Lecture 12: Numericl Qudrture J.K. Ryn@tudelft.nl WI3097TU Delft Institute of Applied Mthemtics Delft University of Technology 5 December 2012 () Numericl Qudrture 5 December 2012 1 / 46 Outline 1 Review
More informationCalculus I-II Review Sheet
Clculus I-II Review Sheet 1 Definitions 1.1 Functions A function is f is incresing on n intervl if x y implies f(x) f(y), nd decresing if x y implies f(x) f(y). It is clled monotonic if it is either incresing
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 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 information5.7 Improper Integrals
458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the
More informationIII. Lecture on Numerical Integration. File faclib/dattab/lecture-notes/numerical-inter03.tex /by EC, 3/14/2008 at 15:11, version 9
III Lecture on Numericl Integrtion File fclib/dttb/lecture-notes/numerical-inter03.tex /by EC, 3/14/008 t 15:11, version 9 1 Sttement of the Numericl Integrtion Problem In this lecture we consider the
More informationMeasuring Electron Work Function in Metal
n experiment of the Electron topic Mesuring Electron Work Function in Metl Instructor: 梁生 Office: 7-318 Emil: shling@bjtu.edu.cn Purposes 1. To understnd the concept of electron work function in metl nd
More informationP 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)
1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this
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 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 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 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 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 informationChapter 3 Polynomials
Dr M DRAIEF As described in the introduction of Chpter 1, pplictions of solving liner equtions rise in number of different settings In prticulr, we will in this chpter focus on the problem of modelling
More informationA Modified ADM for Solving Systems of Linear Fredholm Integral Equations of the Second Kind
Applied Mthemticl Sciences, Vol. 6, 2012, no. 26, 1267-1273 A Modified ADM for Solving Systems of Liner Fredholm Integrl Equtions of the Second Kind A. R. Vhidi nd T. Dmercheli Deprtment of Mthemtics,
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 informationLecture 14 Numerical integration: advanced topics
Lecture 14 Numericl integrtion: dvnced topics Weinn E 1,2 nd Tiejun Li 2 1 Deprtment of Mthemtics, Princeton University, weinn@princeton.edu 2 School of Mthemticl Sciences, Peking University, tieli@pku.edu.cn
More 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 informationDOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES
DOIG PHYSICS WITH MATLAB MATHEMATICAL ROUTIES COMPUTATIO OF OE-DIMESIOAL ITEGRALS In Cooper School of Physics, University of Sydney in.cooper@sydney.edu.u DOWLOAD DIRECTORY FOR MATLAB SCRIPTS mth_integrtion_1d.m
More 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 informationA027 Uncertainties in Local Anisotropy Estimation from Multi-offset VSP Data
A07 Uncertinties in Locl Anisotropy Estimtion from Multi-offset VSP Dt M. Asghrzdeh* (Curtin University), A. Bon (Curtin University), R. Pevzner (Curtin University), M. Urosevic (Curtin University) & B.
More informationNumerical Integration
Chpter 1 Numericl Integrtion Numericl differentition methods compute pproximtions to the derivtive of function from known vlues of the function. Numericl integrtion uses the sme informtion to compute numericl
More 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 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 information3 Mathematics of the Poisson Equation
3 Mthemtics of the Poisson Eqution 3. Green functions nd the Poisson eqution () The Dirichlet Green function stisfies the Poisson eqution with delt-function chrge 2 G D (r, r o ) = δ 3 (r r o ) (3.) nd
More informationSolving the Poisson Partial Differential Equation using Spectral Polynomial Methods
Solving the Poisson Prtil Differentil Eqution using Spectrl Polynomil Methods Seungkeol Choe April 8th, 4 Abstrct In this report, we present spectrl polynomil method for solving the Poisson eqution with
More informationChapter 1. Basic Concepts
Socrtes Dilecticl Process: The Þrst step is the seprtion of subject into its elements. After this, by deþning nd discovering more bout its prts, one better comprehends the entire subject Socrtes (469-399)
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 informationWeek 10: Riemann integral and its properties
Clculus nd Liner Algebr for Biomedicl Engineering Week 10: Riemnn integrl nd its properties H. Führ, Lehrstuhl A für Mthemtik, RWTH Achen, WS 07 Motivtion: Computing flow from flow rtes 1 We observe the
More informationTerminal Velocity and Raindrop Growth
Terminl Velocity nd Rindrop Growth Terminl velocity for rindrop represents blnce in which weight mss times grvity is equl to drg force. F 3 π3 ρ L g in which is drop rdius, g is grvittionl ccelertion,
More informationLinear and Non-linear Feedback Control Strategies for a 4D Hyperchaotic System
Pure nd Applied Mthemtics Journl 017; 6(1): 5-13 http://www.sciencepublishinggroup.com/j/pmj doi: 10.11648/j.pmj.0170601.1 ISSN: 36-9790 (Print); ISSN: 36-981 (Online) Liner nd Non-liner Feedbck Control
More informationModelling of the near infra-red radiation pulse propagation in biological tissues for medical imaging application
JOURNAL OF INTENSE PULSED LASERS AND APPLICATIONS IN ADVANCED PHYSICS Vol. 3, No. 4, p. 4-45 Modelling of the ner infr-red rdition pulse propgtion in biologicl tissues for medicl imging ppliction A. SAOULI
More informationMain topics for the First Midterm
Min topics for the First Midterm The Midterm will cover Section 1.8, Chpters 2-3, Sections 4.1-4.8, nd Sections 5.1-5.3 (essentilly ll of the mteril covered in clss). Be sure to know the results of the
More informationSummary of equations chapters 7. To make current flow you have to push on the charges. For most materials:
Summry of equtions chpters 7. To mke current flow you hve to push on the chrges. For most mterils: J E E [] The resistivity is prmeter tht vries more thn 4 orders of mgnitude between silver (.6E-8 Ohm.m)
More informationArithmetic Mean Derivative Based Midpoint Rule
Applied Mthemticl Sciences, Vol. 1, 018, no. 13, 65-633 HIKARI Ltd www.m-hikri.com https://doi.org/10.1988/ms.018.858 Arithmetic Men Derivtive Bsed Midpoint Rule Rike Mrjulis 1, M. Imrn, Symsudhuh Numericl
More informationNUMERICAL ANALYSIS MEETS NUMBER THEORY: USINGROOTFINDINGMETHODSTOCALCULATE INVERSES MOD p n
Applicble Anlysis nd Discrete Mthemtics vilble online t http://pefmth.etf.bg.c.yu Appl. Anl. Discrete Mth. x (xxxx), xxx xxx. doi:10.2298/aadmxxxxxxxx NUMERICAL ANALYSIS MEETS NUMBER THEORY: USINGROOTFINDINGMETHODSTOCALCULATE
More informationA Signal-Level Fusion Model for Image-Based Change Detection in DARPA's Dynamic Database System
SPIE Aerosense 001 Conference on Signl Processing, Sensor Fusion, nd Trget Recognition X, April 16-0, Orlndo FL. (Minor errors in published version corrected.) A Signl-Level Fusion Model for Imge-Bsed
More informationOperations with Polynomials
38 Chpter P Prerequisites P.4 Opertions with Polynomils Wht you should lern: How to identify the leding coefficients nd degrees of polynomils How to dd nd subtrct polynomils How to multiply polynomils
More informationCHAPTER 4a. ROOTS OF EQUATIONS
CHAPTER 4. ROOTS OF EQUATIONS A. J. Clrk School o Engineering Deprtment o Civil nd Environmentl Engineering by Dr. Ibrhim A. Asskk Spring 00 ENCE 03 - Computtion Methods in Civil Engineering II Deprtment
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 informationLecture 4: Piecewise Cubic Interpolation
Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 4: Piecewise Cubic Interpoltion Compiled 5 September In this lecture we consider piecewise cubic interpoltion
More informationLecture 17. Integration: Gauss Quadrature. David Semeraro. University of Illinois at Urbana-Champaign. March 20, 2014
Lecture 17 Integrtion: Guss Qudrture Dvid Semerro University of Illinois t Urbn-Chmpign Mrch 0, 014 Dvid Semerro (NCSA) CS 57 Mrch 0, 014 1 / 9 Tody: Objectives identify the most widely used qudrture method
More informationSOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı
Mthemticl nd Computtionl Applictions, Vol., No., pp. 37-5, 5 http://dx.doi.org/.99/mc-5- SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD Mehmet Pkdemirli nd Gözde
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 Bernstein polynomial approach for solution of nonlinear integral equations
Avilble online t wwwisr-publictionscom/jns J Nonliner Sci Appl, 10 (2017), 4638 4647 Reserch Article Journl Homepge: wwwtjnscom - wwwisr-publictionscom/jns A Bernstein polynomil pproch for solution of
More informationTHERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION
XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es
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 informationDuality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.
Dulity #. Second itertion for HW problem Recll our LP emple problem we hve been working on, in equlity form, is given below.,,,, 8 m F which, when written in slightly different form, is 8 F Recll tht we
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 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 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 information( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that
Arc Length of Curves in Three Dimensionl Spce If the vector function r(t) f(t) i + g(t) j + h(t) k trces out the curve C s t vries, we cn mesure distnces long C using formul nerly identicl to one tht we
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 informationStuff You Need to Know From Calculus
Stuff You Need to Know From Clculus For the first time in the semester, the stuff we re doing is finlly going to look like clculus (with vector slnt, of course). This mens tht in order to succeed, you
More informationAPPLICATIONS OF THE DEFINITE INTEGRAL
APPLICATIONS OF THE DEFINITE INTEGRAL. Volume: Slicing, disks nd wshers.. Volumes by Slicing. Suppose solid object hs boundries extending from x =, to x = b, nd tht its cross-section in plne pssing through
More informationJackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jckson 2.26 Homework Problem Solution Dr. Christopher S. Bird University of Msschusetts Lowell PROBLEM: The two-dimensionl region, ρ, φ β, is bounded by conducting surfces t φ =, ρ =, nd φ = β held t zero
More informationCOSC 3361 Numerical Analysis I Numerical Integration and Differentiation (III) - Gauss Quadrature and Adaptive Quadrature
COSC 336 Numericl Anlysis I Numericl Integrtion nd Dierentition III - Guss Qudrture nd Adptive Qudrture Edgr Griel Fll 5 COSC 336 Numericl Anlysis I Edgr Griel Summry o the lst lecture I For pproximting
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 information