Layer Potentials. Chapter Fundamental and Singular Solutions

Size: px
Start display at page:

Download "Layer Potentials. Chapter Fundamental and Singular Solutions"

Transcription

1 Chapter Layer Potentals We solve boundary value problems assocated wth system (1.1) (or, rather, ts homogeneous verson) by means of potental-type functons wth a sutably chosen kernel. In ths chapter we construct a matrx of fundamental solutons for the operator A ω ( x ), whch we can then use to defne generalzed sngle-layer and doublelayer plate potentals. The method used s analogous to the one employed n 14 to construct the correspondng matrx for the operator A( x ) defned by (1.11), whch occurs n the study of the equlbrum bendng of plates. Later, we wrte the matrx of fundamental solutons n a form that allows us to decompose t nto an nfnte seres of so-called wavefunctons. Ths form, constructed n Theorem.1, s smlar to the correspondng matrx n the theory of plane elastodynamcs 48 wth the added complcaton of certan computatonal constants. In Secton. we nvestgate the sngulartes of the matrx of fundamental solutons and of ts assocated matrx of sngular solutons. Thus, n Theorems.3 and.4 we fnd that these sngulartes concde wth those of the correspondng matrces from equlbrum plate theory. Therefore, the sngle-layer and double-layer potentals ntroduced n Secton.3 behave n the same way as the potentals consdered n 14. The mportant propertes of these functons, used extensvely n the subsequent analyss, are contaned n Theorems Fundamental and Sngular Solutons We construct a matrx of fundamental solutons for the operator A ω ( x ) usng the method descrbed n 13. If (ξ ) s the adjont of the matrx A ω (ξ ), then u(x) ( x )B(x), (.1) where B satsfes (deta ω )( x )B(x)H(x). (.) From (1.1) and (1.11) t follows that G.R. Thomson and C. Constanda, Statonary Oscllatons of Elastc Plates, A Boundary Integral Equaton Analyss, DOI 1.17/ _, Sprnger Scence+Busness Meda, LLC 11 7

2 8 Layer Potentals deta ω (ξ ) μξ 1 μ ξ 1 (h Δ 1)+ρω h μξ 1 + μξ μ ξ (h Δ 1)+ρω h μξ +(μδ + ρω ) h 4 μ(λ + μ)δ + h (λ + 3μ)(ρω h μ)δ +(ρω h μ) μ 3 (ξ1 + ξ )(h Δ 1)+ρω h μ (ξ1 + ξ ) + h 4 μ (λ + μ)δ 3 + h μ(λ + 3μ)(ρω h μ)+ρω h 4 μ(λ + μ) Δ + μ(ρω h μ) + ρω h (λ + 3μ)(ρω h μ) Δ + ρω (ρω h μ) h 4 μ (λ + μ)δ 3 + ρω h 4 μ(λ + 5μ) h μ (λ + μ) Δ + ρω h (ρω h μ)(λ + 4μ)Δ + ρω (ρω h μ). Factorng ths expresson, we obtan ( deta ω (ξ )h 4 μ (λ + μ) Δ + λ + 3μ λ + μ k Δ + μ ) λ + μ k k3 (Δ + k3), where and hence, where k ρω μ (.3) k 3 k 1 h ; (.4) deta ω (ξ )h 4 μ (λ + μ)(δ + k1 )(Δ + k )(Δ + k 3 ), (.5) k1 + k λ + 3μ λ + μ k, k1 k μ (.6) λ + μ k k3. Wthout loss of generalty, we assume that k1 k. Usng (.4) and (.6), we fnd that k1, k, and k 3 are connected by the equalty h (k1 k 3 )(k k 3 )h k1 k k 3 (k 1 + k )+k4 3 h ( μ λ + μ k k 3 λ + 3μ λ + μ k k 3 + k 4 3 h ( k k3 + k3) 4 ( h k3 k + k 1h ), ) from whch h (k 1 k 3)(k k 3) k 3. (.7)

3 .1 Fundamental and Sngular Solutons 9 We clam that, under assumptons (1.14) and (1.15), k 1, k, and k 3 are real, strctly postve, and dstnct. Frst, k1 and k are the two roots of the equaton x + λ + 3μ λ + μ k x + μ λ + μ k k3. The dscrmnant of ths quadratc s (λ + 3μ) (λ + μ) k4 4μ λ + μ k k 3 k (λ + μ) k (λ + μ) k (λ + μ) k (λ + 6λμ+ 9μ ) 4μ(λ + μ) k (λ + λμ+ μ )+ k (λ + μ) + 4μ(λ + μ) h 4μ(λ + μ) h >. (k 1h ) Consequently, k1 and k are real and dstnct. Also, by (.4), assumpton (1.15) mples that k3 >. By (.6) and (1.14), ths means that k 1 + k >, k 1 k >. Hence, k 1 and k are strctly postve. Fnally, from (.7) and the fact that k 3 t follows that k 1 k 3 and k k 3. Replacng, n turn, each component of H by δ( x y ), where δ s the Drac delta dstrbuton, and settng the other two equal to zero, from (.1) and (.) we obtan the matrx of fundamental solutons where, by (.) and (.5), t(x,y) s a soluton of D ω (x,y) ( x ) t(x,y)e 3, (.8) h 4 μ (λ + μ)(δ + k1 )(Δ + k )(Δ + k 3 )t(x,y) δ( x y ). (.9) We seek t(x,y) n the form t(x,y) 3 j1 b j (k j x y ), (.1) where s the Hankel functon of the frst knd of order zero and the b j are constants to be determned from (.9). Ths Hankel functon s a fundamental soluton of the Helmholtz operator and satsfes 75 (Δ + k j) (k j x y )4δ( x y ). (.11)

4 1 Layer Potentals From (.1) and (.11) we fnd that (Δ + k3 )t(x,y)b 1 4δ( x y )+(k 3 k1 )H(1) (k 1 x y ) + b 4δ( x y )+(k 3 k )H(1) (k x y ) + 4b 3 δ( x y ). To elmnate the Drac dstrbuton n ths equaton, we requre that Now, f (.1) s satsfed, we see that b 1 + b + b 3. (.1) (Δ + k )(Δ + k3)t(x,y) b 1 (k3 k 1 ) 4δ( x y )+(k k 1 )H(1) (k 1 x y ) + 4b (k 3 k )δ( x y ), so we must have (k 3 k 1)b 1 +(k 3 k )b. (.13) Snce t(x, y) satsfes (.9), we obtan 1 h 4 μ (λ + μ) δ( x y )(Δ + k 1 )(Δ + k )(Δ + k 3 )t(x,y) 4b 1 (k1 k )(k 1 k 3 )δ( x y ), from whch we deduce that b 1 4h 4 μ (λ + μ)(k1 k )(k 1 (.14) k 3 ). Substtutng ths nto (.1) and (.13) yelds b b 3 4h 4 μ (λ + μ)(k k 1 )(k k 3 ), 4h 4 μ (λ + μ)(k 3 k 1 )(k 3 k ). (.15) The formula for b 3 can be smplfed by means of (.7): b 3 4h μ (λ + μ)k3. (.16) These constants are well defned snce the k j are dstnct. To calculate the matrx of fundamental solutons D ω (x,y) usng (.8), we frst need to compute ( x ). From (1.1) and (1.11) t follows that

5 .1 Fundamental and Sngular Solutons (ξ ) h μδ + h (λ + μ)ξ μ + ρω h μξ μξ μδ + ρω h μ ΔΔ + h μ(λ + μ)δξ +(ρω h μ)μδ + ρω h μδ In the same way, + ρω h (λ + μ)ξ + ρω (ρω h μ)+μ ξ h μ ΔΔ + h μ(λ + μ)δδ h μ(λ + μ)δξ 1 μ ξ 1 + ρω h μδ + ρω h (λ + μ)δ ρω h (λ + μ)ξ1 + ρω (ρω h μ) h μ(λ + μ)δδ h μ(λ + μ)δξ1 + ρω h (λ + 3μ)Δ ( μ + ρω h (λ + μ) ) ξ 1 + ρω (ρω h μ). (ξ )h μ(λ + μ)δδ h μ(λ + μ)δξ + ρω h (λ + 3μ)Δ Next, 1 (ξ ) h (λ + μ)ξ 1 ξ μξ 1 μξ μδ + ρω h (λ + μ)ξ 1 ξ μξ μξ 1 μδ + ρω Aω 1 (ξ ) ( μ + ρω h (λ + μ) ) ξ + ρω (ρω h μ). h μ(λ + μ)δξ 1 ξ ( μ + ρω h (λ + μ) ) ξ 1 ξ, 13 (ξ ) h (λ + μ)ξ 1 ξ μξ 1 h μδ + h (λ + μ)ξ μ + ρω h μξ h (λ + μ)ξ 1 ξ h μδ + h (λ + μ)ξ μ + ρω h μξ 1 μξ Aω 31 (ξ ) h μ(λ + μ)ξ 1 ξ + h μ Δξ 1 + h μ(λ + μ)ξ 1 ξ +(ρω h μ)μξ 1 h μ Δξ 1 + μ(ρω h μ)ξ 1, 3 (ξ ) h μδ + h (λ + μ)ξ1 μ + ρω h μξ 1 h (λ + 1 ξ μξ Aω 3 (ξ ) h μ Δξ + μ(ρω h μ)ξ. Fnally, 33 (ξ ) h μδ + h (λ + μ)ξ1 μ + ρω h h (λ + μ)ξ 1 ξ h (λ + 1 ξ h μδ + h (λ + μ)ξ μ + ρω h

6 1 Layer Potentals h 4 μ ΔΔ + h 4 μ(λ + μ)δξ + h μ(ρω h μ)δ + h 4 μ(λ + μ)δξ 1 + h 4 (λ + μ) ξ 1 ξ + h (λ + μ)(ρω h μ)ξ 1 + h μ(ρω h μ)δ + h (λ + μ)(ρω h μ)ξ +(ρω h μ) h 4 (λ + μ) ξ 1 ξ h 4 μ(λ + μ)δδ + h (λ + 3μ)(ρω h μ)δ +(ρω h μ). Therefore, the elements of the adjont of the matrx A ω (ξ ) are α (ξ )h μ(λ + μ)δ α ΔΔ h μ(λ + μ)δξ α ξ + ρω h (λ + 3μ)δ α Δ ( ρω h (λ + μ)+μ ) ξ α ξ + ρω (ρω h μ)δ α, (.17) α3 (ξ ) Aω 3α (ξ )h μ Δξ α + μ(ρω h μ)ξ α, (.18) 33 (ξ )h4 μ(λ + μ)δδ + h (λ + 3μ)(ρω h μ)δ +(ρω h μ). (.19) Theorem.1. The elements of the matrx of fundamental solutons D ω (x,y) are Dα ω (x,y) 4h μk3 α1 x α x (k 1 x y ) α x α x (k x y ) + x α x (k 3 x y )+δ α k3 H(1) (k 3 x y ), (.) Dα3 ω (x,y) Dω 3α (x,y) 4h μk3 γ x (k 1 x y )+γ α x (k x y ), (.1) α D33 ω (x,y) 4h μk3 where 1 (k 1 x y )+ H(1) (k x y ), (.) α 1 k μ k 3 k k 1 1 h k 3 (k 1 μ k 3 ) k 1 k, α k 1 μ k3 k1, (.3) k γ μ k3 k1, (.4) k, h k3 (k μ k3 ) k, (.5) k 1 wth μ μ λ + μ. (.6)

7 .1 Fundamental and Sngular Solutons 13 Proof. Frst, we rewrte the elements of the adjont matrx n a more convenent form. Usng (.3) and (.4), from (.17) we obtan α (ξ )δ α h μ(λ + μ)δδ + ρω h (λ + 3μ)Δ + ρω (ρω h μ) h μ(λ + μ)δ + ρω h (λ + μ)+μ ξ α ξ δ α h μ(λ + μ) ΔΔ + λ + 3μ λ + μ k Δ + μ λ + μ k k3 so, usng (.4) and (.6), we arrve at α (ξ )δ αh μ(λ + μ)(δ + k 1)(Δ + k ) Smlarly, from (.18) and (.19), h μ(λ + μ)(δ + k )ξ α ξ μ ξ α ξ, h μ(λ + μ)(δ + k 3)ξ α ξ μ(λ + μ)ξ α ξ. α3 (ξ ) Aω 3α (ξ ) h μ Δξ α + μ(ρω h μ)ξ α ( h μ Δ + ρω h ) μ h ξ α h μ (Δ + k μ 3)ξ α, 33 (ξ )h4 μ(λ + μ)δδ + h (λ + 3μ)(ρω h μ)δ +(ρω h μ) h 4 μ(λ + μ)(δ + μ k 3)(Δ + k 3). By the defnton of D ω (x,y), Dα ω (x,y)aω α ( x)t(x,y) δ α h μ(λ + μ)(k1 k3)(k k3)b 3 (k 3 x y ) h μ(λ + μ) (k3 k 1 )b 1 x α x (k 1 x y ) +(k3 k)b x α x (k x y ) μ(λ + μ) b 1 x α x (k 1 x y )+b x α x (k x y ) Takng (.14) (.16) nto account, we fnd that + b 3 x α x (k 3 x y ).

8 14 Layer Potentals Dα ω (x,y)δ α 4h μ H(1) (k 3 x y )+ 4h μk3 x α x (k 3 x y ) + h μ(λ + μ)(k1 k 3 ) μ(λ + μ) 4h 4 μ (λ + μ)(k1 k )(k 1 k 3 ) x α x (k 1 x y ) + h μ(λ + μ)(k k 3 ) μ(λ + μ) 4h 4 μ (λ + μ)(k k 1 )(k k 3 ) x α x (k x y ) 4h μk3 α1 x α x (k 1 x y ) α x α x (k x y ) + x α x (k 3 x y )+δ α k3h (1) (k 3 x y ), where α1 k 3 λ + μ k k 1 α k 3 k 1 k λ + μ 1 h (k1 k 3 ) λ + μ λ + μ 1 h (k k 3 ),. Now, usng (.7) and (.6), we deduce that α1 k 3 k 1 μ + k k 3 k 1 k3 1 k (k k 3 μ k3 + k k3) k μ k3 1 k, k 1 as requred. Smlarly, t can be shown that Next, α k 1 μ k3 k1. k D ω α3 (x,y) Dω 3α (x,y)aω α3 ( x)t(x,y) h μ (k3 k 1 ) 4h 4 μ (λ + μ)(k1 k )(k 1 k 3 ) x (k 1 x y ) α h μ (k3 + k ) 4h 4 μ (λ + μ)(k k 1 )(k k 3 ) x (k x y ) α 4h μk3 γ x (k 1 x y )+γ α x (k x y ), α where γ s defned by (.4).

9 .1 Fundamental and Sngular Solutons 15 Fnally, D ω 33 (x,y)aω 33 ( x)t(x,y) h4 μ(λ + μ)(μ k 3 k 1 )(k 3 k 1 ) 4h 4 μ (λ + μ)(k 1 k )(k 1 k 3 ) H(1) (k 1 x y ) + h4 μ(λ + μ)(μ k3 k )(k 3 k ) 4h 4 μ (λ + μ)(k k 1 )(k k 3 ) H(1) (k x y ) 4h μk3 1 (k 1 x y )+ (k x y ), where 1 and are defned by (.5). Ths completes the proof. The representaton of D ω (x,y) gven n Theorem.1 s suffcent for our purposes. There s no need to calculate explctly the frst-order and second-order dervatves of the Hankel functons n (.) (.). To smplfy the notaton, we wrte r x y. Snce and f (r) f (r) x α y α f (r) f (r), x α x y α y from (.) (.) t s easly seen that We ntroduce the matrx of sngular solutons D ω (x,y) D ω (y,x) T. (.7) P ω (x,y) T ( y )D ω (y,x) T. (.8) Theorem.. Each column of D ω (x,y) and P ω (x,y) satsfes the homogeneous system A ω ( x )u(x) at all ponts x R,x y. Proof. From (.8), (.5), and (.9) we see that for x y, A ω ( x )D ω (x,y)a ω ( x ) ( x ) t(x,y)e 3 Usng (.8) and (.7), we fnd that for x y, (deta ω )( x )t(x,y) E 3. Ak ω ( x)pkj ω (x,y)aω k ( x) T ( y )D ω (y,x) jk Ak ω ( x)t jm ( y )Dmk ω (y,x) T jm ( y )Ak ω ( x)dkm ω (x,y), whch proves the asserton.

10 16 Layer Potentals. Order of Sngularty The nvestgaton of the sngulartes of D ω and P ω plays an mportant role n the study of the behavor of the sngle-layer and double-layer plate potentals. It s known 1 that, as ξ, (ξ ) π so, as r, from (.1) we see that t(x, y)t(r) 3 j1 π π b j (k jr) 3 j1 3 j1 ( ξ ξ 4 ) lnξ ; b j (1 1 4 k j r + 1 ) 64 k4 j r4 ln(k j r) b j ( k j r k4 j r4 )lnr + d + O(r 6 lnr), where d s a constant. From (.14) and (.15) t follows that 3 j1 b j 3 j1 b j k j, 3 j1 b j k 4 j 4h 4 μ (λ + μ), whch leads to where t(r)c 1 r 4 lnr + d + O(r 6 lnr), (.9) 1 c 1 18πh 4 μ (λ + μ). (.3) Theorem.3. As r, ( D ω (x,y)lnr d 1 E γγ 1 ) πμ E (x α y α )(x y ) 33 + d r E α where C s a constant (3 3)-matrx, +C + D ω (x,y), (.31) λ + 3μ d 1 4πh μ(λ + μ), d λ + μ 4πh μ(λ + μ), (.3) and D ω (x,y) s a (3 3)-matrx whose elements are O(r lnr) as r.

11 . Order of Sngularty 17 Proof. By (.8), (.17), (.9), and (.3), we fnd that D ω α (x,y)aω α ( x)t(x,y) h μ(λ + μ)δ α ΔΔ(c 1 r 4 lnr) h μ(λ + μ) Δ(c 1 r 4 lnr) x α x hence, + ρω h (λ + 3μ)δ α Δ(c 1 r 4 lnr) ρω h (λ + μ)+μ (c 1 r 4 lnr) x α x + ρω (ρω h μ)δ α c 1 r 4 lnr +C α + O(r lnr) 64c 1 h μ(λ + μ)δ α lnr c 1 h μ(λ + μ) (16r lnr + 8r ) x α x +C α + O(r lnr); Dα ω (x,y)64c 1h μ(λ + μ)δ α lnr c 1 h μ(λ + μ) 3 (x α y α )(x y ) r + 3δ α lnr +C α + O(r lnr) λ + 3μ 4πh μ(λ + μ) δ λ + μ (x α y α )(x y ) α lnr + 4πh μ(λ + μ) r Also, usng (.18) and (.19), we see that and +C α + O(r lnr). D ω α3 (x,y) Dω 3α (x,y)aω α3 ( x)t(x,y) h μ x α Δ(c 1 r 4 lnr)+μ(ρω h μ) x α (c 1 r 4 lnr) C α3 + O(r lnr) +C α3 + O(r3 lnr) D33 ω (x,y)aω 33 ( x)t(x,y) h 4 μ(λ + μ)δδ(c 1 r 4 lnr)+h (λ + 3μ)(ρω h μ)δ(c 1 r 4 lnr) +(ρω h μ) c 1 r 4 lnr +C 33 + O(r lnr) 64c 1 h 4 μ(λ + μ)lnr +C 33 + O(r lnr) 1 πμ lnr +C 33 + O(r lnr). The asserton s proved.

12 18 Layer Potentals We ntroduce the alternatng symbol ε α α for α, {1,}. In what follows, ν(y) and denote the dervatves n the normal and tangental drectons, respectvely. s(y) Theorem.4. As r, P ω (x,y) 1 { μ ε π α s(y) lnr E α + ν(y) lnr E 3 (λ + μ (x α y α )(x y ) )ε αγ s(y) r E γ 1 (xα y α )lnr E 3α ν(y) + 1 ε α (xα y α )lnr ( λ E s(y) ) h E } 1 ν α(y) h E α3 + (λ + μ )E 3α where D s a constant (3 3)-matrx, λ λ (λ + μ), + D + P ω (x,y), (.33) μ s defned by (.6), and P ω (x,y) s a (3 3)-matrx whose elements are O(r lnr) as r. Proof. From (.8) (.3) we obtan so, P11 ω (x,y)t 11( y )D11 ω (y,x)+t 1( y )D1 ω (y,x)+t 13( y )D31 ω (y,x) h (λ + μ)ν 1 A ω 11 ( y)t(y,x), 1 +h μν A ω 11 ( y)t(y,x), + h μν A ω 1 ( y)t(y,x), 1 +h λν 1 A ω 1 ( y)t(y,x), h (λ + μ)ν 1 h μ(λ + μ)δδt, 1 h μ(λ + μ)δt, h μν h μ(λ + μ)δδt, h μ(λ + μ)δt, 11 + h μν h μ(λ + μ)δt, 11 + h λν 1 h μ(λ + μ)δt, 1 + D 11 + O(r lnr) h 4 μ(λ + μ)(λ + μ)ν 1 ΔΔt, 1 +h 4 μ (λ + μ)ν 1 ΔΔt, 1 + h 4 μ (λ + μ)ν ΔΔt, h 4 μ (λ + μ)ν 1 Δt, 111 h 4 μ (λ + μ)ν Δt, 11 h 4 λμ(λ + μ)ν 1 Δt, 111 h 4 λμ(λ + μ)ν 1 Δt, 1 +D 11 + O(r lnr);

13 . Order of Sngularty 19 P11 ω (x,y) h 4 μ (λ + μ) ν(y) ΔΔt + h4 μ (λ + μ)ν 1 Δt, 111 h 4 μ (λ + μ)ν 1 Δt, 111 +h 4 μ (λ + μ)ν 1 Δt, 1 h 4 μ (λ + μ)ν Δt, 11 + D 11 + O(r lnr) h 4 μ (λ + μ) ν(y) ΔΔt + h4 μ (λ + μ) s(y) Δt, 1 +D 11 + O(r lnr) 64h 4 μ (λ + μ)c 1 ν(y) lnr + 64h4 μ (λ + μ)c 1 s(y) whch means that P ω 11 (x,y) 1 π ν(y) lnr 1 π (λ + μ ) s(y) (x 1 y 1 )(x y ) r + D 11 + O(r lnr), (x 1 y 1 )(x y ) r + D 11 + O(r lnr). Analogous manpulaton yelds P1 ω (x,y) h4 λμ s(y) Δt, 11 +h 4 μ (λ + μ) s(y) Δt, +D 1 + O(r lnr); hence, P ω 1 (x,y) λ 4π(λ + μ) 1 Smlarly, we fnd that P ω π μ (x y ) s(y) r (x y ) 1 4π s(y) r + lnr + lnr + D 1 + O(r lnr) s(y) lnr 1 π (λ + μ (x y ) ) s(y) r + D 1 + O(r lnr). 13 (x,y)h μ (λ + μ) s(y) Δt, +D 13 + O(r lnr) 1 (x 4πh y )lnr 1 s(y) 4πh ν 1(y)+D 13 + O(r lnr), 1 (x,y) h4 μ (λ + μ) s(y) Δt, 11 +h 4 λμ s(y) Δt, +D 1 + O(r lnr) 1 π μ s(y) lnr + 1 π (λ + μ (x 1 y 1 ) ) s(y) r + D 1 + O(r lnr), P ω

14 Layer Potentals (x,y)h4 μ (λ + μ) ν(y) ΔΔt h4 μ (λ + μ) s(y) Δt, 1 P ω 1 π ν(y) lnr + 1 π (λ + μ ) s(y) + D + O(r lnr) (x 1 y 1 )(x y ) r + D + O(r lnr), P3 ω (x,y) h μ (λ + μ) s(y) Δt, 1 +D 3 + O(r lnr) 1 (x1 4πh y 1 )lnr 1 s(y) 4πh ν (y)+d 3 + O(r lnr), P31 ω (x,y)h4 μ (λ + μ) ν(y) Δt, 1 +h 4 λμ s(y) Δt, +D 31 + O(r lnr) 1 (x1 y 1 )lnr 1 4π ν(y) 4π λ (x y )lnr s(y) 1 π (λ + μ )ν 1 (y)+d 31 + O(r lnr), P3 ω (x,y) h4 λμ s(y) Δt, 1 +h 4 μ (λ + μ) ν(y) Δt, +D 3 + O(r lnr) 1 4π λ (x1 y 1 )lnr 1 (x y )lnr s(y) 4π ν(y) 1 π (λ + μ )ν (y)+d 3 + O(r lnr), P33 ω (x,y)h4 μ (λ + μ) ν(y) ΔΔt + D 33 + O(r lnr) 1 π ν(y) lnr + D 33 + O(r lnr), as requred. Remark.1. Expansons (.31) and (.33) of D ω (x,y) and P ω (x,y), respectvely, for y close to x concde wth those of the correspondng matrces arsng n the equlbrum bendng of plates, whch can be found n Propertes of the Potentals We defne the sngle-layer potental (V ω ϕ)(x) D ω (x,y)ϕ(y)ds(y) (.34) S

15 .3 Propertes of the Potentals 1 and the double-layer potental (W ω ϕ)(x) P ω (x,y)ϕ(y)ds(y), (.35) S where ϕ s a (3 1)-vector functon known as the densty. The propertes of these potentals play an mportant role n formulatng sutable ntegral equatons for varous boundary value problems. We are especally concerned wth the behavor of (V ω ϕ)(x) and (W ω ϕ)(x) as x approaches S. Ths requres a very detaled analyss based on the expansons (.31) and (.33). As was noted n Remark.1, the sngulartes of D ω and P ω are exactly the same as those occurrng n the correspondng matrces n the equlbrum bendng of plates. The behavor of the correspondng potentals s thoroughly nvestgated n 14, and many relevant propertes are smply quoted from there. We denote by C,α (X) the vector space of Hölder contnuous functons (wth ndex α) onx, and by C 1,α (X) the subspace of C 1 (X) of functons whose frstorder dervatves belong to C,α (X). In the sequel, when we menton that a functon ϕ C,α (X), α (,1), satsfes a certan equaton or condton, we mean that ϕ C,α (X) for all α n the nterval (,1). The next three assertons follow from Theorem. and the results n 14. Theorem.5. () If ϕ C(S), then V ω ϕ and W ω ϕ are analytc and satsfy A ω ( x )u(x) n S + S. () If ϕ C,α (S), α (,1), then the drect values V ωϕ and W ωϕ of V ω ϕ and W ω ϕ on S exst (the latter n the sense of prncpal value). Theorem.6. If ϕ C,α (S), α (,1), then the functons V ω+ (ϕ)(v ω ϕ) S +, V ω (ϕ)(v ω ϕ) S (.36) are of class C (S + ) C 1,α ( S + ) and C (S ) C 1,α ( S ), respectvely, and T V ω+ (ϕ) ( W ω + 1 I ) ϕ on S, (.37) T V ω (ϕ) ( W ω 1 I ) ϕ on S, (.38) where W ω s the adjont of W ω and I s the dentty operator. Theorem.7. If ϕ C 1,α (S), α (,1), then the functons { W ω+ (W ω ϕ) (ϕ) S + n S +, ( W ω 1 I ) ϕ on S, { W ω (W ω ϕ) (ϕ) S n S, ( W ω + 1 I ) ϕ on S, (.39) (.4) are of class C (S + ) C 1,α ( S + ) and C (S ) C 1,α ( S ), respectvely, and TW ω+ (ϕ)t W ω (ϕ) on S. (.41)

16 Layer Potentals Remark.. Theorems.5.7 are used n Chapters 6 11 to rgorously justfy the constructon of regular solutons of the fundamental boundary value problems n S + and S n the form of layer potentals.

17

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Appendix B. Criterion of Riemann-Stieltjes Integrability

Appendix B. Criterion of Riemann-Stieltjes Integrability Appendx B. Crteron of Remann-Steltes Integrablty Ths note s complementary to [R, Ch. 6] and [T, Sec. 3.5]. The man result of ths note s Theorem B.3, whch provdes the necessary and suffcent condtons for

More information

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS BOUNDEDNESS OF THE IESZ TANSFOM WITH MATIX A WEIGHTS Introducton Let L = L ( n, be the functon space wth norm (ˆ f L = f(x C dx d < For a d d matrx valued functon W : wth W (x postve sem-defnte for all

More information

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars

More information

Exercise Solutions to Real Analysis

Exercise Solutions to Real Analysis xercse Solutons to Real Analyss Note: References refer to H. L. Royden, Real Analyss xersze 1. Gven any set A any ɛ > 0, there s an open set O such that A O m O m A + ɛ. Soluton 1. If m A =, then there

More information

GELFAND-TSETLIN BASIS FOR THE REPRESENTATIONS OF gl n

GELFAND-TSETLIN BASIS FOR THE REPRESENTATIONS OF gl n GELFAND-TSETLIN BASIS FOR THE REPRESENTATIONS OF gl n KANG LU FINITE DIMENSIONAL REPRESENTATIONS OF gl n Let e j,, j =,, n denote the standard bass of the general lnear Le algebra gl n over the feld of

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

First day August 1, Problems and Solutions

First day August 1, Problems and Solutions FOURTH INTERNATIONAL COMPETITION FOR UNIVERSITY STUDENTS IN MATHEMATICS July 30 August 4, 997, Plovdv, BULGARIA Frst day August, 997 Problems and Solutons Problem. Let {ε n } n= be a sequence of postve

More information

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

COMPLEX NUMBERS AND QUADRATIC EQUATIONS COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not

More information

Lecture 5 Decoding Binary BCH Codes

Lecture 5 Decoding Binary BCH Codes Lecture 5 Decodng Bnary BCH Codes In ths class, we wll ntroduce dfferent methods for decodng BCH codes 51 Decodng the [15, 7, 5] 2 -BCH Code Consder the [15, 7, 5] 2 -code C we ntroduced n the last lecture

More information

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1. 7636S ADVANCED QUANTUM MECHANICS Soluton Set 1 Sprng 013 1 Warm-up Show that the egenvalues of a Hermtan operator  are real and that the egenkets correspondng to dfferent egenvalues are orthogonal (b)

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

Affine transformations and convexity

Affine transformations and convexity Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

REAL ANALYSIS I HOMEWORK 1

REAL ANALYSIS I HOMEWORK 1 REAL ANALYSIS I HOMEWORK CİHAN BAHRAN The questons are from Tao s text. Exercse 0.0.. If (x α ) α A s a collecton of numbers x α [0, + ] such that x α

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim Causal Damonds M. Aghl, L. Bombell, B. Plgrm Introducton The correcton to volume of a causal nterval due to curvature of spacetme has been done by Myrhem [] and recently by Gbbons & Solodukhn [] and later

More information

Solutions Homework 4 March 5, 2018

Solutions Homework 4 March 5, 2018 1 Solutons Homework 4 March 5, 018 Soluton to Exercse 5.1.8: Let a IR be a translaton and c > 0 be a re-scalng. ˆb1 (cx + a) cx n + a (cx 1 + a) c x n x 1 cˆb 1 (x), whch shows ˆb 1 s locaton nvarant and

More information

Formulas for the Determinant

Formulas for the Determinant page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use

More information

On Finite Rank Perturbation of Diagonalizable Operators

On Finite Rank Perturbation of Diagonalizable Operators Functonal Analyss, Approxmaton and Computaton 6 (1) (2014), 49 53 Publshed by Faculty of Scences and Mathematcs, Unversty of Nš, Serba Avalable at: http://wwwpmfnacrs/faac On Fnte Rank Perturbaton of Dagonalzable

More information

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space.

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space. Lnear, affne, and convex sets and hulls In the sequel, unless otherwse specfed, X wll denote a real vector space. Lnes and segments. Gven two ponts x, y X, we defne xy = {x + t(y x) : t R} = {(1 t)x +

More information

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

Georgia Tech PHYS 6124 Mathematical Methods of Physics I Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends

More information

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16 STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus

More information

9 Characteristic classes

9 Characteristic classes THEODORE VORONOV DIFFERENTIAL GEOMETRY. Sprng 2009 [under constructon] 9 Characterstc classes 9.1 The frst Chern class of a lne bundle Consder a complex vector bundle E B of rank p. We shall construct

More information

A summation on Bernoulli numbers

A summation on Bernoulli numbers Journal of Number Theory 111 (005 37 391 www.elsever.com/locate/jnt A summaton on Bernoull numbers Kwang-Wu Chen Department of Mathematcs and Computer Scence Educaton, Tape Muncpal Teachers College, No.

More information

Lecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2

Lecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2 P470 Lecture 6/7 (February 10/1, 014) DIRAC EQUATION The non-relatvstc Schrödnger equaton was obtaned by notng that the Hamltonan H = P (1) m can be transformed nto an operator form wth the substtutons

More information

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices Internatonal Mathematcal Forum, Vol 11, 2016, no 11, 513-520 HIKARI Ltd, wwwm-hkarcom http://dxdoorg/1012988/mf20166442 The Jacobsthal and Jacobsthal-Lucas Numbers va Square Roots of Matrces Saadet Arslan

More information

5 The Rational Canonical Form

5 The Rational Canonical Form 5 The Ratonal Canoncal Form Here p s a monc rreducble factor of the mnmum polynomal m T and s not necessarly of degree one Let F p denote the feld constructed earler n the course, consstng of all matrces

More information

A SURVEY OF PROPERTIES OF FINITE HORIZON DIFFERENTIAL GAMES UNDER ISAACS CONDITION. Contents

A SURVEY OF PROPERTIES OF FINITE HORIZON DIFFERENTIAL GAMES UNDER ISAACS CONDITION. Contents A SURVEY OF PROPERTIES OF FINITE HORIZON DIFFERENTIAL GAMES UNDER ISAACS CONDITION BOTAO WU Abstract. In ths paper, we attempt to answer the followng questons about dfferental games: 1) when does a two-player,

More information

5 The Laplace Equation in a convex polygon

5 The Laplace Equation in a convex polygon 5 Te Laplace Equaton n a convex polygon Te most mportant ellptc PDEs are te Laplace, te modfed Helmoltz and te Helmoltz equatons. Te Laplace equaton s u xx + u yy =. (5.) Te real and magnary parts of an

More information

On the symmetric character of the thermal conductivity tensor

On the symmetric character of the thermal conductivity tensor On the symmetrc character of the thermal conductvty tensor Al R. Hadjesfandar Department of Mechancal and Aerospace Engneerng Unversty at Buffalo, State Unversty of New York Buffalo, NY 146 USA ah@buffalo.edu

More information

TR/95 February Splines G. H. BEHFOROOZ* & N. PAPAMICHAEL

TR/95 February Splines G. H. BEHFOROOZ* & N. PAPAMICHAEL TR/9 February 980 End Condtons for Interpolatory Quntc Splnes by G. H. BEHFOROOZ* & N. PAPAMICHAEL *Present address: Dept of Matematcs Unversty of Tabrz Tabrz Iran. W9609 A B S T R A C T Accurate end condtons

More information

MULTIPLICATIVE FUNCTIONS: A REWRITE OF ANDREWS CHAPTER 6

MULTIPLICATIVE FUNCTIONS: A REWRITE OF ANDREWS CHAPTER 6 MULTIPLICATIVE FUNCTIONS: A REWRITE OF ANDREWS CHAPTER 6 In these notes we offer a rewrte of Andrews Chapter 6. Our am s to replace some of the messer arguments n Andrews. To acheve ths, we need to change

More information

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

General viscosity iterative method for a sequence of quasi-nonexpansive mappings Avalable onlne at www.tjnsa.com J. Nonlnear Sc. Appl. 9 (2016), 5672 5682 Research Artcle General vscosty teratve method for a sequence of quas-nonexpansve mappngs Cuje Zhang, Ynan Wang College of Scence,

More information

Geometry of Müntz Spaces

Geometry of Müntz Spaces WDS'12 Proceedngs of Contrbuted Papers, Part I, 31 35, 212. ISBN 978-8-7378-224-5 MATFYZPRESS Geometry of Müntz Spaces P. Petráček Charles Unversty, Faculty of Mathematcs and Physcs, Prague, Czech Republc.

More information

STEINHAUS PROPERTY IN BANACH LATTICES

STEINHAUS PROPERTY IN BANACH LATTICES DEPARTMENT OF MATHEMATICS TECHNICAL REPORT STEINHAUS PROPERTY IN BANACH LATTICES DAMIAN KUBIAK AND DAVID TIDWELL SPRING 2015 No. 2015-1 TENNESSEE TECHNOLOGICAL UNIVERSITY Cookevlle, TN 38505 STEINHAUS

More information

CHAPTER III Neural Networks as Associative Memory

CHAPTER III Neural Networks as Associative Memory CHAPTER III Neural Networs as Assocatve Memory Introducton One of the prmary functons of the bran s assocatve memory. We assocate the faces wth names, letters wth sounds, or we can recognze the people

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

The Feynman path integral

The Feynman path integral The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space

More information

Convexity preserving interpolation by splines of arbitrary degree

Convexity preserving interpolation by splines of arbitrary degree Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

Time-Varying Systems and Computations Lecture 6

Time-Varying Systems and Computations Lecture 6 Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS These are nformal notes whch cover some of the materal whch s not n the course book. The man purpose s to gve a number of nontrval examples

More information

Randić Energy and Randić Estrada Index of a Graph

Randić Energy and Randić Estrada Index of a Graph EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 5, No., 202, 88-96 ISSN 307-5543 www.ejpam.com SPECIAL ISSUE FOR THE INTERNATIONAL CONFERENCE ON APPLIED ANALYSIS AND ALGEBRA 29 JUNE -02JULY 20, ISTANBUL

More information

SL n (F ) Equals its Own Derived Group

SL n (F ) Equals its Own Derived Group Internatonal Journal of Algebra, Vol. 2, 2008, no. 12, 585-594 SL n (F ) Equals ts Own Derved Group Jorge Macel BMCC-The Cty Unversty of New York, CUNY 199 Chambers street, New York, NY 10007, USA macel@cms.nyu.edu

More information

Supplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso

Supplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso Supplement: Proofs and Techncal Detals for The Soluton Path of the Generalzed Lasso Ryan J. Tbshran Jonathan Taylor In ths document we gve supplementary detals to the paper The Soluton Path of the Generalzed

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

CSCE 790S Background Results

CSCE 790S Background Results CSCE 790S Background Results Stephen A. Fenner September 8, 011 Abstract These results are background to the course CSCE 790S/CSCE 790B, Quantum Computaton and Informaton (Sprng 007 and Fall 011). Each

More information

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning Journal of Machne Learnng Research 00-9 Submtted /0; Publshed 7/ Erratum: A Generalzed Path Integral Control Approach to Renforcement Learnng Evangelos ATheodorou Jonas Buchl Stefan Schaal Department of

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

Curvature and isoperimetric inequality

Curvature and isoperimetric inequality urvature and sopermetrc nequalty Julà ufí, Agustí Reventós, arlos J Rodríguez Abstract We prove an nequalty nvolvng the length of a plane curve and the ntegral of ts radus of curvature, that has as a consequence

More information

Binomial transforms of the modified k-fibonacci-like sequence

Binomial transforms of the modified k-fibonacci-like sequence Internatonal Journal of Mathematcs and Computer Scence, 14(2019, no. 1, 47 59 M CS Bnomal transforms of the modfed k-fbonacc-lke sequence Youngwoo Kwon Department of mathematcs Korea Unversty Seoul, Republc

More information

Bernoulli Numbers and Polynomials

Bernoulli Numbers and Polynomials Bernoull Numbers and Polynomals T. Muthukumar tmk@tk.ac.n 17 Jun 2014 The sum of frst n natural numbers 1, 2, 3,..., n s n n(n + 1 S 1 (n := m = = n2 2 2 + n 2. Ths formula can be derved by notng that

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

= z 20 z n. (k 20) + 4 z k = 4

= z 20 z n. (k 20) + 4 z k = 4 Problem Set #7 solutons 7.2.. (a Fnd the coeffcent of z k n (z + z 5 + z 6 + z 7 + 5, k 20. We use the known seres expanson ( n+l ( z l l z n below: (z + z 5 + z 6 + z 7 + 5 (z 5 ( + z + z 2 + z + 5 5

More information

Implicit Integration Henyey Method

Implicit Integration Henyey Method Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure

More information

Assortment Optimization under MNL

Assortment Optimization under MNL Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.

More information

3 Basic boundary value problems for analytic function in the upper half plane

3 Basic boundary value problems for analytic function in the upper half plane 3 Basc boundary value problems for analytc functon n the upper half plane 3. Posson representaton formulas for the half plane Let f be an analytc functon of z throughout the half plane Imz > 0, contnuous

More information

SELECTED SOLUTIONS, SECTION (Weak duality) Prove that the primal and dual values p and d defined by equations (4.3.2) and (4.3.3) satisfy p d.

SELECTED SOLUTIONS, SECTION (Weak duality) Prove that the primal and dual values p and d defined by equations (4.3.2) and (4.3.3) satisfy p d. SELECTED SOLUTIONS, SECTION 4.3 1. Weak dualty Prove that the prmal and dual values p and d defned by equatons 4.3. and 4.3.3 satsfy p d. We consder an optmzaton problem of the form The Lagrangan for ths

More information

ACTM State Calculus Competition Saturday April 30, 2011

ACTM State Calculus Competition Saturday April 30, 2011 ACTM State Calculus Competton Saturday Aprl 30, 2011 ACTM State Calculus Competton Sprng 2011 Page 1 Instructons: For questons 1 through 25, mark the best answer choce on the answer sheet provde Afterward

More information

Complete subgraphs in multipartite graphs

Complete subgraphs in multipartite graphs Complete subgraphs n multpartte graphs FLORIAN PFENDER Unverstät Rostock, Insttut für Mathematk D-18057 Rostock, Germany Floran.Pfender@un-rostock.de Abstract Turán s Theorem states that every graph G

More information

CHAPTER 5: Lie Differentiation and Angular Momentum

CHAPTER 5: Lie Differentiation and Angular Momentum CHAPTER 5: Le Dfferentaton and Angular Momentum Jose G. Vargas 1 Le dfferentaton Kähler s theory of angular momentum s a specalzaton of hs approach to Le dfferentaton. We could deal wth the former drectly,

More information

arxiv: v1 [math.ho] 18 May 2008

arxiv: v1 [math.ho] 18 May 2008 Recurrence Formulas for Fbonacc Sums Adlson J. V. Brandão, João L. Martns 2 arxv:0805.2707v [math.ho] 8 May 2008 Abstract. In ths artcle we present a new recurrence formula for a fnte sum nvolvng the Fbonacc

More information

Random Walks on Digraphs

Random Walks on Digraphs Random Walks on Dgraphs J. J. P. Veerman October 23, 27 Introducton Let V = {, n} be a vertex set and S a non-negatve row-stochastc matrx (.e. rows sum to ). V and S defne a dgraph G = G(V, S) and a drected

More information

REDUCTION MODULO p. We will prove the reduction modulo p theorem in the general form as given by exercise 4.12, p. 143, of [1].

REDUCTION MODULO p. We will prove the reduction modulo p theorem in the general form as given by exercise 4.12, p. 143, of [1]. REDUCTION MODULO p. IAN KIMING We wll prove the reducton modulo p theorem n the general form as gven by exercse 4.12, p. 143, of [1]. We consder an ellptc curve E defned over Q and gven by a Weerstraß

More information

SPECTRAL PROPERTIES OF IMAGE MEASURES UNDER THE INFINITE CONFLICT INTERACTION

SPECTRAL PROPERTIES OF IMAGE MEASURES UNDER THE INFINITE CONFLICT INTERACTION SPECTRAL PROPERTIES OF IMAGE MEASURES UNDER THE INFINITE CONFLICT INTERACTION SERGIO ALBEVERIO 1,2,3,4, VOLODYMYR KOSHMANENKO 5, MYKOLA PRATSIOVYTYI 6, GRYGORIY TORBIN 7 Abstract. We ntroduce the conflct

More information

Quantum Mechanics for Scientists and Engineers. David Miller

Quantum Mechanics for Scientists and Engineers. David Miller Quantum Mechancs for Scentsts and Engneers Davd Mller Types of lnear operators Types of lnear operators Blnear expanson of operators Blnear expanson of lnear operators We know that we can expand functons

More information

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Perron Vectors of an Irreducible Nonnegative Interval Matrix Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of

More information

Combined Wronskian solutions to the 2D Toda molecule equation

Combined Wronskian solutions to the 2D Toda molecule equation Combned Wronskan solutons to the 2D Toda molecule equaton Wen-Xu Ma Department of Mathematcs and Statstcs, Unversty of South Florda, Tampa, FL 33620-5700, USA Abstract By combnng two peces of b-drectonal

More information

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

Modelli Clamfim Equazioni differenziali 7 ottobre 2013 CLAMFIM Bologna Modell 1 @ Clamfm Equazon dfferenzal 7 ottobre 2013 professor Danele Rtell danele.rtell@unbo.t 1/18? Ordnary Dfferental Equatons A dfferental equaton s an equaton that defnes a relatonshp

More information

a b a In case b 0, a being divisible by b is the same as to say that

a b a In case b 0, a being divisible by b is the same as to say that Secton 6.2 Dvsblty among the ntegers An nteger a ε s dvsble by b ε f there s an nteger c ε such that a = bc. Note that s dvsble by any nteger b, snce = b. On the other hand, a s dvsble by only f a = :

More information

Solutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution.

Solutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution. Solutons HW #2 Dual of general LP. Fnd the dual functon of the LP mnmze subject to c T x Gx h Ax = b. Gve the dual problem, and make the mplct equalty constrants explct. Soluton. 1. The Lagrangan s L(x,

More information

Solutions to the 71st William Lowell Putnam Mathematical Competition Saturday, December 4, 2010

Solutions to the 71st William Lowell Putnam Mathematical Competition Saturday, December 4, 2010 Solutons to the 7st Wllam Lowell Putnam Mathematcal Competton Saturday, December 4, 2 Kran Kedlaya and Lenny Ng A The largest such k s n+ 2 n 2. For n even, ths value s acheved by the partton {,n},{2,n

More information

Lecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem.

Lecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem. prnceton u. sp 02 cos 598B: algorthms and complexty Lecture 20: Lft and Project, SDP Dualty Lecturer: Sanjeev Arora Scrbe:Yury Makarychev Today we wll study the Lft and Project method. Then we wll prove

More information

ON THE EXTENDED HAAGERUP TENSOR PRODUCT IN OPERATOR SPACES. 1. Introduction

ON THE EXTENDED HAAGERUP TENSOR PRODUCT IN OPERATOR SPACES. 1. Introduction ON THE EXTENDED HAAGERUP TENSOR PRODUCT IN OPERATOR SPACES TAKASHI ITOH AND MASARU NAGISA Abstract We descrbe the Haagerup tensor product l h l and the extended Haagerup tensor product l eh l n terms of

More information

where a is any ideal of R. Lemma 5.4. Let R be a ring. Then X = Spec R is a topological space Moreover the open sets

where a is any ideal of R. Lemma 5.4. Let R be a ring. Then X = Spec R is a topological space Moreover the open sets 5. Schemes To defne schemes, just as wth algebrac varetes, the dea s to frst defne what an affne scheme s, and then realse an arbtrary scheme, as somethng whch s locally an affne scheme. The defnton of

More information

FINITELY-GENERATED MODULES OVER A PRINCIPAL IDEAL DOMAIN

FINITELY-GENERATED MODULES OVER A PRINCIPAL IDEAL DOMAIN FINITELY-GENERTED MODULES OVER PRINCIPL IDEL DOMIN EMMNUEL KOWLSKI Throughout ths note, s a prncpal deal doman. We recall the classfcaton theorem: Theorem 1. Let M be a fntely-generated -module. (1) There

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

Determinants Containing Powers of Generalized Fibonacci Numbers

Determinants Containing Powers of Generalized Fibonacci Numbers 1 2 3 47 6 23 11 Journal of Integer Sequences, Vol 19 (2016), Artcle 1671 Determnants Contanng Powers of Generalzed Fbonacc Numbers Aram Tangboonduangjt and Thotsaporn Thanatpanonda Mahdol Unversty Internatonal

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

The exponential map of GL(N)

The exponential map of GL(N) The exponental map of GLN arxv:hep-th/9604049v 9 Apr 996 Alexander Laufer Department of physcs Unversty of Konstanz P.O. 5560 M 678 78434 KONSTANZ Aprl 9, 996 Abstract A fnte expanson of the exponental

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Homework & Solution. Contributors. Prof. Lee, Hyun Min. Particle Physics Winter School. Park, Ye

Homework & Solution. Contributors. Prof. Lee, Hyun Min. Particle Physics Winter School. Park, Ye Homework & Soluton Prof. Lee, Hyun Mn Contrbutors Park, Ye J(yej.park@yonse.ac.kr) Lee, Sung Mook(smlngsm0919@gmal.com) Cheong, Dhong Yeon(dhongyeoncheong@gmal.com) Ban, Ka Young(ban94gy@yonse.ac.kr) Ro,

More information

332600_08_1.qxp 4/17/08 11:29 AM Page 481

332600_08_1.qxp 4/17/08 11:29 AM Page 481 336_8_.qxp 4/7/8 :9 AM Page 48 8 Complex Vector Spaces 8. Complex Numbers 8. Conjugates and Dvson of Complex Numbers 8.3 Polar Form and DeMovre s Theorem 8.4 Complex Vector Spaces and Inner Products 8.5

More information

MAE140 - Linear Circuits - Winter 16 Midterm, February 5

MAE140 - Linear Circuits - Winter 16 Midterm, February 5 Instructons ME140 - Lnear Crcuts - Wnter 16 Mdterm, February 5 () Ths exam s open book. You may use whatever wrtten materals you choose, ncludng your class notes and textbook. You may use a hand calculator

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information