u = 0 in Ω, u = f on B,
|
|
- Quentin Wilkins
- 6 years ago
- Views:
Transcription
1 Preprint April 1991 Colordo Stte University Deprtment of Mthemtics COMPUTATION OF WEAKLY AND NEARLY SINGULAR INTEGRALS OVER TRIANGLES IN R 3 EUGENE L. ALLGOWER 1,2,4, KURT GEORG 1,3,4 nd KARL KALIK 3,5 Abstrct. We study the pproximtion of wekly singulr integrls over tringles in generl position in R 3, giving explicit formule where convenient nd numericl qudrture in more generl cses. Prticulr models considered concern the colloction nd Glerkin methods in the boundry integrl pproch to the Dirichlet problem for Lplce s eqution. Key words. singulr integrls, numericl qudrture, boundry element method, symptotic estimtes AMS(MOS) subject clssifictions. 65D30, 65R20, 65N35, 45E05 1. Introduction. In the boundry integrl pproch to solving prtil differentil equtions over domin Ω R 3 one converts to n integrl eqution over B = Ω nd pproximtes the domin Ω by polytope whose boundry pproximtes B. The fces of the polytope re polygons, so there is no loss of generlity in ssuming tht consists of tringles. See, e.g., [2, sec. 15.4] for n introduction to piecewise liner methods for pproximting surfces. An efficient implementtion including smoothing steps hs been developed in [11]. To motivte our generl discussion, let us consider the problem of finding solution u to the first kind integrl eqution (1.1) u(x) x ξ d x = f(ξ), ξ, where denotes the Eucliden norm in R 3. Eqution (1.1) stems from the boundry integrl pproch for the Dirichlet problem u = 0 in Ω, u = f on B, where B is pproximted by, see, e.g., the books [3,4,5,7,8]. The numericl solution of (1.1) is usully crried out vi colloction or Glerkin method. Our discussion below pertins to both methods. Let Q 1,...,Q N be the vertices of nd let 1,..., M be the fces (tringles) of. We will ssume tht the domin Ω nd hence the pproximte boundry re bounded. 1 Prtilly supported by the Ntionl Science Foundtion under grnt number DMS Prtilly supported by the Deutscher Akdemischer Austuschdienst under study visit. 3 Prtilly supported by the Deutsche Forschungsgemeinschft. 4 Deprtment of Mthemtics, Colordo Stte University, Ft. Collins, Colordo Mthemtisches Institut A, Universität Stuttgrt, D(W)-7000 Stuttgrt, Germny. 1
2 2 e. l. llgower, k. georg nd k. klik Let us seek n pproximte solution to (1.1) which is of the form u A i ϕ i, i=1 where the ϕ i re chosen bsis functions nd the A i R re unknown coefficients which must be clculted. Typiclly, the ϕ i re piecewise polynomils. Tht is, ϕ i (x) = B i,q,α x α for x q, where we dopt the convention α n x =(x 1,x 2,x 3 ) q, Z + : = set of non-negtive integers, α =(α 1,α 2,α 3 ) Z 3 +, α :=α 1 + α 2 + α 3, x α :=x α 1 1 xα 2 2 xα 3 3. The non-negtive integer n denotes the degree of the polynomil ϕ i. In the colloction method for (1.1) the unknown coefficients A 1,...,A n re determined by solving the system of liner equtions (1.2) A i i=1 ϕ(x) x Q p d x = f(q p ), p =1,...,N. In (1.2) it is necessry to compute the entries of the coefficient mtrix: i,p := ϕ(x) x Q p d x = M q=1 α n B i,q,α q x α x Q p d x. Consequently, it becomes necessry to clculte integrls of the form x (1.3) I(p, q, α) := α x Q p d x. In the Glerkin method for (1.1) system of equtions of the form q (1.4) A i i=1 ϕ j (ξ) ϕ i (x) x ξ d x d ξ = ϕ j (ξ) f(ξ) d ξ must be stisfied for j =1,...,N. Here the ϕ i re the sme generl bsis functions s those described in the colloction method.
3 Now writing computtion of singulr integrls over tringles 3 I i (ξ) := ϕ i (x) x ξ d x, the coefficients of the system (1.4) ssume the form i,j = ϕ j (ξ) I i (ξ) d ξ. Typiclly, the coefficients i,j re computed vi qudrture formul: i,j c p ϕ(q p ) I i (Q p ) p=1 where the c p, p =1,...,N re weights. Since I i (ξ) = = q=1 ϕ i (x) x ξ d x α n B i,q,α q x α x Q p d x, we re gin led to the tsk of computing integrls of the form (1.3) lso in the Glerkin method. Our im therefore is to develop methods for evluting or pproximting integrls of the form (1.3). 2. Geometric Trnsformtions. Recent relted works concerning pproximtion of surfce integrls in the context of boundry element methods re [1,6,9,10]. The book of Hckbusch [5] discusses the pnel method in detil nd gives some cses in which certin integrls over tringles cn be integrted exctly. Since the tringles q re generlly in n rbitrry position in spce, we begin the tretment of the computtion of (1.3) by describing briefly trnsformtions which bring the generl configurtion of Q p nd q into more stndrd configurtion. Let us suppose tht Q p nd q re described in stndrd 3-dimensionl Crtesin co-ordintes (x 1,x 2,x 3 ), nd the vertices of q re denoted by v 1,v 2,v 3. An ffine trnsformtion x y =(y 1,y 2,y 3 ) brings the generl configurtion into stndrd configurtion hving the following properties: (i) q lies in the plne y 3 = 0. Tht is, y q implies y =(y 1,y 2, 0). (ii) The point Q p lies on the y 3 -xis. Tht is, Q p =(0, 0,c) for some c R. By introducing the orthonorml bsis (2.1) ρ 1 := v2 v 1 v 2 v 1,ρ3 := (v2 v 1 ) (v 3 v 1 ) (v 2 v 1 ) (v 3 v 1 ),ρ2 := ρ 3 ρ 1,
4 4 e. l. llgower, k. georg nd k. klik we obtin y = Ax b or x = A t (y + b), where A :=(ρ 1,ρ 2,ρ 3 ) t,b:= ρ1 v 1 ρ 2 Q p. ρ 3 Q p It is immeditely seen tht x α = d γ y γ 1 1 yγ 2 2 for x q, γ α where γ =(γ 1,γ 2 ) Z 2 + nd d γ R re uniquely determined nd esily computed from the bove trnsformtion. Hence we cn concentrte upon pproximting integrls of the form (2.2) Î(p, q, γ) := q γ y 1 1 yγ 2 2 y Q p d y, since (1.3) is merely liner combintion of the integrls in (2.2). Let us denote by v 0 := (0, 0, 0) the projection of Q p into the y 3 = 0 plne, nd the vertices v i of q re understood to be expessed in the y-co-ordintes. The indexing of the vertices v 1, v 2, v 3 will be chosen to correspond to the orienttion of, i.e., ρ 1 defined in (2.1) is the outer norml of on q. The tringle q oriented in this wy will be written s q =[v 1,v 2,v 3 ]. The prmeter c ppering in the y-co-ordintes of Q p =(0, 0,c) will ply n importnt rôle below. By mking use of the orienttion, it follows tht for ny integrnd F which is defined on ll the domins involved, the following formul holds: F (y) d y = F (y) d y + F (y) d y + F (y) d y. [v 1,v 2,v 3 ] In prticulr, we hve where (2.3) Î i,j (p, q, γ) := [v 0,v 1,v 2 ] [v 0,v 2,v 3 ] Î(p, q, γ) =Î1,2(p, q, γ)+î2,3(p, q, γ)+î3,1(p, q, γ) [v 0,v i,v j ] [v 0,v 3,v 1 ] y γ 1 1 yγ 2 2 y Q p d y for i, j {1, 2, 3}, i j. For convenience, let us denote ny integrl of the form (2.3) by J(γ,c), where c is the y 3 -co-ordinte of Q p nd the indices of the vertices of the tringles hve been supressed. Note tht we hve reduced our problem to the evlution of such integrls. 3. Numericl Approximtion of J(γ,c). We begin the clcultion of J(γ,c) by introducing the polr co-ordintes (3.1) y 1 = ρ cos(θ), y 2 = ρ sin(θ).
5 computtion of singulr integrls over tringles 5 Then (3.2) J(γ,c)= v i b R(θ) (θ) sin γ 2 (θ) 0 ρ γ +1 dρ dθ, ρ2 + c2 where the limits of the integrtion R(θ),, b re derived corresponding to the following figure : R(θ) d µ v j b v 0 β θ α Figure (3.3) Hence we hve π 2 = ( ) v tn 1 i 2 v1 i π 2 ( ) b = tn 1 v j 2 R(θ) = v j 1 d cos(θ β) = if v i 1 =0, if v i 1 0, if v j 1 =0, if v j 1 0, d cos(µ), β = π 2 α, µ = θ β, π 2 ( ) α = tn 1 v j 2 vi 2 v j 1 vi 1 if v j 1 = vi 1, if v j 1 vi 1.
6 6 e. l. llgower, k. georg nd k. klik We will consider the integrl J(γ,c) in the form (3.4) J(γ,c)= where b (3.5) G p (r, c) := (θ) sin γ 2 (θ) G γ +1 (R(θ),c) dθ, r Integrtion by prts yields the recursion formul G 0 (r, c) =ln ( r + r 2 + c 2) ln( c ), 0 ρ p ρ2 + c 2 dρ. (3.6) G 1 (r, c) = r 2 + c 2 c, G p (r, c) = 1 p rp 1 r 2 + c 2 p 1 p c2 G p 2 (r, c) for p =2, 3,... Note tht the function G p (r, c) cn be efficiently clculted vi this recursion. The formule hold for r>0 nd c 0, nd re supplemented by G p (r, 0) = 1 p rp for r>0 nd p>0. It is cler tht for fixed c R nd p =1, 2,..., the function r G p (r, c) isc, nd consequently the integrnd in (3.4) is C. Hence, the integrl J(γ,c) cn be pproximted by highly ccurte nd efficient numericl qudrture methods. We do not elborte this point but refer to the stndrd literture on numericl qudrture. 4. Explicit Integrtion for the Cse c =0. For the cse c = 0, the integrl J(γ,c) in (3.4) cn be explicitly clculted s follows. We hve G p (r, 0) = 1 p rp nd consequently, see (3.3), = J(γ,0) = 1 b (θ) sin γ 2 (θ) R γ +1 γ +1 (θ) dθ = dγ +1 γ +1 b (θ) sin γ 2 (θ) cos γ +1 (θ β) dθ b β = dγ +1 (µ + β) sin γ 2 (µ + β) dµ. γ +1 β cos γ +1 (µ) For convenience, let us set A := sin(β) nd B := cos(β). Then (µ + β) sin γ 2 (µ + β) γ 1 p 1 =0 p 2 =0 cos γ +1 (µ) γ 2 ( γ1 ( 1) γ 2 p 2 p 1 )( γ2 p 2 ) sin p 1+p 2 (µ) cos γ p 1 p 2 (µ) A γ 2+p 1 p 2B γ 1 +p 2 p 1. cos γ +1 (µ)
7 From this we hve computtion of singulr integrls over tringles 7 (4.1) J(γ,0) = dγ +1 γ +1 γ 1 γ 2 p 1 =0 p 2 =0 where we denoted by F k the integrl F k := ( 1) γ 2 p 2 ( γ1 b β β Now integrtion by prts yields the recursion p 1 )( γ2 p 2 sin k (µ) cos k+1 (µ) dµ. ) A γ 2+p 1 p 2B γ 1 +p 2 p 1 Fp1 +p 2, ( ) 1 + tn(µ/2) µ=b β F 0 =ln, 1 tn(µ/2) µ= β µ=b β 1 F 1 = cos(µ), µ= β F k = 1 sin k 1 µ=b β (µ) k cos k k 1 (µ) k F k 2 for k =2, 3,..., which mkes (4.1) explicit. µ= β 5. Asymptotic Estimtes for Smll c. In colloction or Glerkin methods, the entries of the system mtrix i,j re typiclly pproximted only to such precision tht the order of the discretiztion error is not worsened. Hence, it is of interest to investigte the error which is introduced when replcing the integrl J(γ,c) byj(γ,0) for smll c. If this error is smll, then the (more efficient) explicit clcultions described in section 4 cn be used for pproximting J(γ,c). We differentite the recursion (3.6) twice with respect to c nd let c tend to zero. Thus the following limits re obtined: 2 G 1 (r, 0 ± 0) = 1, 2 G p (r, 0) = 0 for p =2, 3,..., 2 2 G p(r, 0) = rp 2 p 2 for p =3, 4,..., Here 2 indictes the prtil derivtive with respect to the c-vrible, nd 0 ± 0 indictes whether the limit c 0 is tken from the right or the left. As consequence, we obtin the following symptotic estimtes for the integrl in (3.4): O ( c ) for γ =0, J(γ,c) =J(γ,0) + o ( c ) for γ =1, O ( c 2) for γ =2, 3,...
8 8 e. l. llgower, k. georg nd k. klik Let us recll tht the prmeter c describes the distnce of the colloction point Q p from the plne generted by the tringle q. On the other hnd, the discretiztion error of the colloction or Glerkin method is usully known nd expressed s O(h τ ) where h indictes the mximl dimeter of the tringles q. Hence, if for exmple γ 2, then the pproximtion J(γ,c) J(γ,0) is permissible for colloction points Q p which re ner to q in the sense tht the distnce c does not exceed the order O ( h τ/2). REFERENCES [1] M. H. Alibdi, W. S. Hll nd T. T. Hibbs, Exct singulrity cncelling for the potentil kernel in the boundry element method, Communictions in Applied Numericl Methods, 3 (1987), [2] E.L. Allgower nd K. Georg, Numericl continution methods: An introduction. Springer Verlg, Berlin, Heidelberg, New York, 1990 [3] P. K. Bnerjee nd R. Butterfield, Boundry Element Methods in Engineering Science, [4] R. B Guenther nd J. Lee, Prtil differentil equtions of mthemticl physics nd integrl equtions, Prentice Hll, Englewood Cliffs, 1987 [5] W. Hckbusch, Integrlgleichungen. Theorie und Numerik, B. G. Teubner, Stuttgrt, 1989 [6] K. Hymi nd C. A. Brebbi, A new coordinte trnsformtion method for singulr nd nerly singulr integrls over generl curved boundry elements, Boundry Elements IX, vol 1, editors: C. A. Brebbi, W. L. Wendlnd, G. Kutiu, Springer Verlg, Berlin, Heidelberg, New York [7] M. A. Jswon nd G. T. Symm, Integrl Eqution Methods in Potentil Theory nd Electrosttics, Acdemic Press, New York, 1977 [8] E. Mrtensen, Potentiltheorie, B. G. Teubner, Stuttgrt, 1968 [9] Ch. Schwb nd W. L. Wendlnd, 3-D BEM nd numericl integrtion, in Proc. Seventh. Conf. on Boundry Element Methods in Engineering, C. Brebbi, ed., Springer, 1985 [10] W. L. Wendlnd, Asymptotic ccurcy nd convergence progress in boundry element methods, editor: C. A. Brebbi, Pentech-Press, London, Plymouth, (1981), [11] R. Widmnn, An efficient lgorithm for the tringultion of surfces in R 3, preprint, Colordo Stte University, October 1990.
New 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 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 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 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 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 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 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 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 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 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 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 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 informationTheoretische Physik 2: Elektrodynamik (Prof. A.-S. Smith) Home assignment 4
WiSe 1 8.1.1 Prof. Dr. A.-S. Smith Dipl.-Phys. Ellen Fischermeier Dipl.-Phys. Mtthis Sb m Lehrstuhl für Theoretische Physik I Deprtment für Physik Friedrich-Alexnder-Universität Erlngen-Nürnberg Theoretische
More informationNotes on length and conformal metrics
Notes on length nd conforml metrics We recll how to mesure the Eucliden distnce of n rc in the plne. Let α : [, b] R 2 be smooth (C ) rc. Tht is α(t) (x(t), y(t)) where x(t) nd y(t) re smooth rel vlued
More informationMath 100 Review Sheet
Mth 100 Review Sheet Joseph H. Silvermn December 2010 This outline of Mth 100 is summry of the mteril covered in the course. It is designed to be study id, but it is only n outline nd should be used s
More informationNumerical integration
2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter
More informationSummary: Method of Separation of Variables
Physics 246 Electricity nd Mgnetism I, Fll 26, Lecture 22 1 Summry: Method of Seprtion of Vribles 1. Seprtion of Vribles in Crtesin Coordintes 2. Fourier Series Suggested Reding: Griffiths: Chpter 3, Section
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationNumerical Integration
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 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 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 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 information1.9 C 2 inner variations
46 CHAPTER 1. INDIRECT METHODS 1.9 C 2 inner vritions So fr, we hve restricted ttention to liner vritions. These re vritions of the form vx; ǫ = ux + ǫφx where φ is in some liner perturbtion clss P, for
More informationSpace Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space.
Clculus 3 Li Vs Spce Curves Recll the prmetric equtions of curve in xy-plne nd compre them with prmetric equtions of curve in spce. Prmetric curve in plne x = x(t) y = y(t) Prmetric curve in spce x = x(t)
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 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 information1.3 The Lemma of DuBois-Reymond
28 CHAPTER 1. INDIRECT METHODS 1.3 The Lemm of DuBois-Reymond We needed extr regulrity to integrte by prts nd obtin the Euler- Lgrnge eqution. The following result shows tht, t lest sometimes, the extr
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 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 informationk ) and directrix x = h p is A focal chord is a line segment which passes through the focus of a parabola and has endpoints on the parabola.
Stndrd Eqution of Prol with vertex ( h, k ) nd directrix y = k p is ( x h) p ( y k ) = 4. Verticl xis of symmetry Stndrd Eqution of Prol with vertex ( h, k ) nd directrix x = h p is ( y k ) p( x h) = 4.
More informationConvex Sets and Functions
B Convex Sets nd Functions Definition B1 Let L, +, ) be rel liner spce nd let C be subset of L The set C is convex if, for ll x,y C nd ll [, 1], we hve 1 )x+y C In other words, every point on the line
More informationPractice final exam solutions
University of Pennsylvni Deprtment of Mthemtics Mth 26 Honors Clculus II Spring Semester 29 Prof. Grssi, T.A. Asher Auel Prctice finl exm solutions 1. Let F : 2 2 be defined by F (x, y (x + y, x y. If
More informationREGULARITY OF NONLOCAL MINIMAL CONES IN DIMENSION 2
EGULAITY OF NONLOCAL MINIMAL CONES IN DIMENSION 2 OVIDIU SAVIN AND ENICO VALDINOCI Abstrct. We show tht the only nonlocl s-miniml cones in 2 re the trivil ones for ll s 0, 1). As consequence we obtin tht
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 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 informationPlates on elastic foundation
Pltes on elstic foundtion Circulr elstic plte, xil-symmetric lod, Winkler soil (fter Timoshenko & Woinowsky-Krieger (1959) - Chpter 8) Prepred by Enzo Mrtinelli Drft version ( April 016) Introduction Winkler
More informationPartial Differential Equations
Prtil Differentil Equtions Notes by Robert Piché, Tmpere University of Technology reen s Functions. reen s Function for One-Dimensionl Eqution The reen s function provides complete solution to boundry
More informationLine Integrals. Partitioning the Curve. Estimating the Mass
Line Integrls Suppose we hve curve in the xy plne nd ssocite density δ(p ) = δ(x, y) t ech point on the curve. urves, of course, do not hve density or mss, but it my sometimes be convenient or useful to
More informationThe Wave Equation I. MA 436 Kurt Bryan
1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string
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 information. Double-angle formulas. Your answer should involve trig functions of θ, and not of 2θ. sin 2 (θ) =
Review of some needed Trig. Identities for Integrtion. Your nswers should be n ngle in RADIANS. rccos( 1 ) = π rccos( - 1 ) = 2π 2 3 2 3 rcsin( 1 ) = π rcsin( - 1 ) = -π 2 6 2 6 Cn you do similr problems?
More informationMathematics. Area under Curve.
Mthemtics Are under Curve www.testprepkrt.com Tle of Content 1. Introduction.. Procedure of Curve Sketching. 3. Sketching of Some common Curves. 4. Are of Bounded Regions. 5. Sign convention for finding
More informationAnalytical Approximate Solution of Carleman s Equation by Using Maclaurin Series
Interntionl Mthemticl Forum, 5, 2010, no. 60, 2985-2993 Anlyticl Approximte Solution of Crlemn s Eqution by Using Mclurin Series M. Yghobifr 1 Institute for Mthemticl Reserch University Putr Mlysi Serdng
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 informationMath Advanced Calculus II
Mth 452 - Advnced Clculus II Line Integrls nd Green s Theorem The min gol of this chpter is to prove Stoke s theorem, which is the multivrible version of the fundmentl theorem of clculus. We will be focused
More informationAN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION
Applied Mthemtics E-Notes, 5(005), 53-60 c ISSN 1607-510 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ AN INTEGRAL INEQUALITY FOR CONVEX FUNCTIONS AND APPLICATIONS IN NUMERICAL INTEGRATION
More information. Double-angle formulas. Your answer should involve trig functions of θ, and not of 2θ. cos(2θ) = sin(2θ) =.
Review of some needed Trig Identities for Integrtion Your nswers should be n ngle in RADIANS rccos( 1 2 ) = rccos( - 1 2 ) = rcsin( 1 2 ) = rcsin( - 1 2 ) = Cn you do similr problems? Review of Bsic Concepts
More informationVariational Techniques for Sturm-Liouville Eigenvalue Problems
Vritionl Techniques for Sturm-Liouville Eigenvlue Problems Vlerie Cormni Deprtment of Mthemtics nd Sttistics University of Nebrsk, Lincoln Lincoln, NE 68588 Emil: vcormni@mth.unl.edu Rolf Ryhm Deprtment
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 informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More information1 The 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 informationThe Basic Functional 2 1
2 The Bsic Functionl 2 1 Chpter 2: THE BASIC FUNCTIONAL TABLE OF CONTENTS Pge 2.1 Introduction..................... 2 3 2.2 The First Vrition.................. 2 3 2.3 The Euler Eqution..................
More informationDisclaimer: This Final Exam Study Guide is meant to help you start studying. It is not necessarily a complete list of everything you need to know.
Disclimer: This is ment to help you strt studying. It is not necessrily complete list of everything you need to know. The MTH 33 finl exm minly consists of stndrd response questions where students must
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 informationCandidates must show on each answer book the type of calculator used.
UNIVERSITY OF EAST ANGLIA School of Mthemtics My/June UG Exmintion 2007 2008 ELECTRICITY AND MAGNETISM Time llowed: 3 hours Attempt FIVE questions. Cndidtes must show on ech nswer book the type of clcultor
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 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 informationINTRODUCTION TO INTEGRATION
INTRODUCTION TO INTEGRATION 5.1 Ares nd Distnces Assume f(x) 0 on the intervl [, b]. Let A be the re under the grph of f(x). b We will obtin n pproximtion of A in the following three steps. STEP 1: Divide
More informationDiscrete Least-squares Approximations
Discrete Lest-squres Approximtions Given set of dt points (x, y ), (x, y ),, (x m, y m ), norml nd useful prctice in mny pplictions in sttistics, engineering nd other pplied sciences is to construct curve
More informationIntegral points on the rational curve
Integrl points on the rtionl curve y x bx c x ;, b, c integers. Konstntine Zeltor Mthemtics University of Wisconsin - Mrinette 750 W. Byshore Street Mrinette, WI 5443-453 Also: Konstntine Zeltor P.O. Box
More informationMath 32B Discussion Session Session 7 Notes August 28, 2018
Mth 32B iscussion ession ession 7 Notes August 28, 28 In tody s discussion we ll tlk bout surfce integrls both of sclr functions nd of vector fields nd we ll try to relte these to the mny other integrls
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 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 informationOrthogonal Polynomials and Least-Squares Approximations to Functions
Chpter Orthogonl Polynomils nd Lest-Squres Approximtions to Functions **4/5/3 ET. Discrete Lest-Squres Approximtions Given set of dt points (x,y ), (x,y ),..., (x m,y m ), norml nd useful prctice in mny
More informationSTEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.
STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA STEPHEN SCHECTER. The unit step function nd piecewise continuous functions The Heviside unit step function u(t) is given by if t
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 informationThe First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).
The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples
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 Fall 2006 Sample problems for the final exam: Solutions
Mth 42-5 Fll 26 Smple problems for the finl exm: Solutions Any problem my be ltered or replced by different one! Some possibly useful informtion Prsevl s equlity for the complex form of the Fourier series
More information(4.1) D r v(t) ω(t, v(t))
1.4. Differentil inequlities. Let D r denote the right hnd derivtive of function. If ω(t, u) is sclr function of the sclrs t, u in some open connected set Ω, we sy tht function v(t), t < b, is solution
More informationIntegration Techniques
Integrtion Techniques. Integrtion of Trigonometric Functions Exmple. Evlute cos x. Recll tht cos x = cos x. Hence, cos x Exmple. Evlute = ( + cos x) = (x + sin x) + C = x + 4 sin x + C. cos 3 x. Let u
More 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 informationn f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1
The Fundmentl Theorem of Clculus As we continue to study the re problem, let s think bck to wht we know bout computing res of regions enclosed by curves. If we wnt to find the re of the region below the
More information4. Calculus of Variations
4. Clculus of Vritions Introduction - Typicl Problems The clculus of vritions generlises the theory of mxim nd minim. Exmple (): Shortest distnce between two points. On given surfce (e.g. plne), nd the
More information12 TRANSFORMING BIVARIATE DENSITY FUNCTIONS
1 TRANSFORMING BIVARIATE DENSITY FUNCTIONS Hving seen how to trnsform the probbility density functions ssocited with single rndom vrible, the next logicl step is to see how to trnsform bivrite probbility
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 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 informationMATH 174A: PROBLEM SET 5. Suggested Solution
MATH 174A: PROBLEM SET 5 Suggested Solution Problem 1. Suppose tht I [, b] is n intervl. Let f 1 b f() d for f C(I; R) (i.e. f is continuous rel-vlued function on I), nd let L 1 (I) denote the completion
More informationKinematic Waves. These are waves which result from the conservation equation. t + I = 0. (2)
Introduction Kinemtic Wves These re wves which result from the conservtion eqution E t + I = 0 (1) where E represents sclr density field nd I, its outer flux. The one-dimensionl form of (1) is E t + I
More informationNumerical Analysis. 10th ed. R L Burden, J D Faires, and A M Burden
Numericl Anlysis 10th ed R L Burden, J D Fires, nd A M Burden Bemer Presenttion Slides Prepred by Dr. Annette M. Burden Youngstown Stte University July 9, 2015 Chpter 4.1: Numericl Differentition 1 Three-Point
More 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 informationIndefinite Integral. Chapter Integration - reverse of differentiation
Chpter Indefinite Integrl Most of the mthemticl opertions hve inverse opertions. The inverse opertion of differentition is clled integrtion. For exmple, describing process t the given moment knowing the
More informationIntegrals along Curves.
Integrls long Curves. 1. Pth integrls. Let : [, b] R n be continuous function nd let be the imge ([, b]) of. We refer to both nd s curve. If we need to distinguish between the two we cll the function the
More informationProblem Set 3 Solutions
Msschusetts Institute of Technology Deprtment of Physics Physics 8.07 Fll 2005 Problem Set 3 Solutions Problem 1: Cylindricl Cpcitor Griffiths Problems 2.39: Let the totl chrge per unit length on the inner
More informationMath 6455 Oct 10, Differential Geometry I Fall 2006, Georgia Tech
Mth 6455 Oct 10, 2006 1 Differentil Geometry I Fll 2006, Georgi Tech Lecture Notes 12 Riemnnin Metrics 0.1 Definition If M is smooth mnifold then by Riemnnin metric g on M we men smooth ssignment of n
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 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 informationg i fφdx dx = x i i=1 is a Hilbert space. We shall, henceforth, abuse notation and write g i f(x) = f
1. Appliction of functionl nlysis to PEs 1.1. Introduction. In this section we give little introduction to prtil differentil equtions. In prticulr we consider the problem u(x) = f(x) x, u(x) = x (1) where
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 informationMethod of stationary phase
Physics 4 Spring 16 Method of sttionry phse Lecture notes by M. G. Rozmn Lst modified: April 13, 16 There is n immedite generliztion of the Lplce integrls f t)e φt) dt 1) which we obtin by llowing the
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 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 informationNote 16. Stokes theorem Differential Geometry, 2005
Note 16. Stokes theorem ifferentil Geometry, 2005 Stokes theorem is the centrl result in the theory of integrtion on mnifolds. It gives the reltion between exterior differentition (see Note 14) nd integrtion
More informationUnit 5. Integration techniques
18.01 EXERCISES Unit 5. Integrtion techniques 5A. Inverse trigonometric functions; Hyperbolic functions 5A-1 Evlute ) tn 1 3 b) sin 1 ( 3/) c) If θ = tn 1 5, then evlute sin θ, cos θ, cot θ, csc θ, nd
More informationReview of basic calculus
Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below
More informationComposite Mendeleev s Quadratures for Solving a Linear Fredholm Integral Equation of The Second Kind
Globl Journl of Pure nd Applied Mthemtics. ISSN 0973-1768 Volume 12, Number (2016), pp. 393 398 Reserch Indi Publictions http://www.ripubliction.com/gjpm.htm Composite Mendeleev s Qudrtures for Solving
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 informationThe Dirac distribution
A DIRAC DISTRIBUTION A The Dirc distribution A Definition of the Dirc distribution The Dirc distribution δx cn be introduced by three equivlent wys Dirc [] defined it by reltions δx dx, δx if x The distribution
More informationS. S. Dragomir. 2, we have the inequality. b a
Bull Koren Mth Soc 005 No pp 3 30 SOME COMPANIONS OF OSTROWSKI S INEQUALITY FOR ABSOLUTELY CONTINUOUS FUNCTIONS AND APPLICATIONS S S Drgomir Abstrct Compnions of Ostrowski s integrl ineulity for bsolutely
More informationTRAPEZOIDAL TYPE INEQUALITIES FOR n TIME DIFFERENTIABLE FUNCTIONS
TRAPEZOIDAL TYPE INEQUALITIES FOR n TIME DIFFERENTIABLE FUNCTIONS S.S. DRAGOMIR AND A. SOFO Abstrct. In this pper by utilising result given by Fink we obtin some new results relting to the trpezoidl inequlity
More information