# ORTHOGONAL POLYNOMIALS AND CUBATURE FORMULAE ON SPHERES AND ON SIMPLICES. Yuan Xu

Save this PDF as:

Size: px
Start display at page:

Download "ORTHOGONAL POLYNOMIALS AND CUBATURE FORMULAE ON SPHERES AND ON SIMPLICES. Yuan Xu"

## Transcription

1 ORTHOGONAL POLYNOMIALS AND CUBATURE FORMULAE ON SPHERES AND ON SIMPLICES Yuan Xu Abstract. Orthogonal polynomials on the standard simplex Σ d in R d are shown to be related to the spherical orthogonal polynomials on the unit sphere S d in R d+1 that are invariant under the group Z 2 Z 2. For a large class of measures on S d cubature formulae invariant under Z 2 Z 2 are shown to be characterized by cubature formulae on Σ d. Moreover, it is also shown that there is a correspondence between orthogonal polynomials and cubature formulae on Σ d and those invariant on the unit ball B d in R d. The results provide a new approach to study orthogonal polynomials and cubature formulae on spheres and on simplices. 1. Introduction The purpose of this paper is to study the connection between orthogonal polynomials and cubature formulae on spheres and on simplices. A special case of the results shows that the orthogonal polynomials with respect to the measure (u 1...u d (1 u 1... u d )) 1/2 du on the standard simplex Σ d on R d are related to the spherical harmonics on the unit sphere S d in R d+1 that are invariant under the group Z 2 Z 2, and there is a one-to-one correspondence between cubature formulae on Σ d and cubature formulae on S d that are invariant under Z 2 Z 2. The main results will be established for a large class of measures. Orthogonal polynomials on the sphere with respect to a measure other than the surface measure have been studied only recently (see [3, 4, 5, 8, 25] and the references there). The most important development has been a theory developed by Dunkl for measures invariant under a finite reflection group, in which the role of Laplacian operator is replaced by a differential-difference operator in the commutative algebra generated by a family of commuting first order differential-difference operators (Dunkl s operators). Motivated by the work of Dunkl, in [25] we used an elementary approach to study the connection between orthogonal polynomials and cubature formulae on the sphere S d and the unit ball B d of R d. It was shown that we can construct homogeneous orthogonal polynomials on S d from the corresponding orthogonal polynomials on B d for a large class of measures Mathematics Subject Classification. 33C50, 33C55, 65D32. Key words and phrases. Orthogonal polynomials in several variables, on spheres, on simplices, on balls, cubature formulae. Supported by the National Science Foundation under Grant DMS Typeset by AMS-TEX

2 2 YUAN XU The present work is a continuation of [25] in spirit. For a large class of measures we will show that orthogonal polynomials on Σ d and those on S d which are invariant under the group G = Z 2 Z 2 are connected by a simple transformation. The result offers a way to study the structure of orthogonal polynomials on the sphere; it may also lead to new understanding on orthogonal polynomials on the simplex. Together with the results in [25], we also obtain a correspondence between orthogonal polynomials on the simplex and those on the unit ball that are invariant under G. Our approach is elementary as in [25]; it is motivated by [3] in which the connection between orthogonal polynomials on S 2 and Σ 2 is used to study orthogonal polynomials with symmetry of order three; the approach does not involve differential or differential-difference operators. A closely related question is constructing cubature formulae on spheres and on simplices. According to a theorem of Sobolev, in order to establish a cubature formula invariant under a finite group, one only has to verify those polynomials that are invariant under the same group. In particular, to establish a Z 2 Z 2 invariant cubature formula on S d, we only have to consider the homogeneous polynomials that are invariant under the group G. This consideration leads us to establish a one-to-one correspondence between the cubature formulae on Σ d and the Z 2 Z 2 invariant cubature formula on S d. Moreover, using the result in [25], we also have a correspondence between the cubature formulae on Σ d and the Z 2 Z 2 invariant cubature formulae on B d. These results allow us to transform cubature formulae on one domain to another. Over the years, there has been a lot of efforts devoted to the construction of cubature formulae on the sphere, on the ball and on the simplex (mostly with respect to the Lebesgue measure); see [6, 12, 15, 17] for some of the references. It is remarkable that the simple correspondence between cubature formulae on these regions has not been noticed before. It yields, in particular, many new cubature formulae on spheres and on simplices; some of them can be derived from the known cubature formulae on a different domain by the correspondence. Because the main focus of this paper is on the relation between the orthogonal structure and cubature formulae on spheres and on simplices, we will concentrate on the theoretic side of the matter. The examples of new cubature formulae will be reported in a separate paper ([7]). The paper is organized as follows. In Section 2 we introduce notation and present the necessary preliminaries, where we also prove the basic lemma. In Section 3 we discuss the relation between orthogonal polynomials on S d and those on Σ d. In Section 4 we discuss the relation between cubature formulae on the unit sphere and those on the simplex. 2. Preliminary and Basic Lemma Basic notation. For x R d we denote by x = x x2 d the usual Euclidean norm. Let B d be the unit ball of R d and S d be the unit sphere on R d+1 ; that is, B d = {x R d : x 1} and S d = {y R d+1 : y = 1}. We also denote the l 1 norm of x R d by x 1 = x x d. Let Σ d be the standard simplex in R d ; that is, Σ d = {u R d : u 1 0,..., u d 0, 1 u 1 0}.

3 ORTHOGONAL POLYNOMIALS AND CUBATURE FORMULA 3 For d = 2, the simplex Σ 2 is the triangle with vertices at (0, 0), (1, 0) and (0, 1). Throughout this paper we fix the following notation: For y R d+1, we write (2.1) y = (y 1,..., y d, y d+1 ) = rx = r(x, x d+1 ), x S d, x B d, where r = y = y y2 d+1 and x = (x 1,..., x d ). Polynomial spaces. Let N 0 be the set of nonnegative integers. For α = (α 1,..., α d ) N d 0 and x = (x 1,..., x d ) R d we write x α = x α 1 1 xα d d. The number α 1 = α 1 + +α d is called the total degree of x α. We denote by Π d the set of polynomials in d variables on R d and by Π d n the subset of polynomials of total degree at most n. We also denote by Pn d the space of homogeneous polynomials of degree n on Rd and we let rn d = dim Pd n. It is well-known that ( ) ( ) n + d n + d 1 dim Π d n = and rn d n = dim Pd n =. d Weight functions. Throughout this paper a weight function on a compact region Ω is a nonnegative function on Ω whose integral over Ω is positive. We will take Ω = S d, B d or Σ d most of the time. On B d and Σ d we use dx (Lebesgue measure), on S d we use rotation invariant measure dω (surface measure) on S d. Orthogonal polynomials on B d and Σ d. Let Ω denote either B d or Σ d. Let W be a weight function on Ω. It is known that for each n N 0 the set of polynomials of degree n that are orthogonal to all polynomials of lower degree forms a vector space V n whose dimension is rn d. We denote by {P k n}, 1 k rd n and n N 0, one family of orthonormal polynomials with respect to W on Ω that forms a basis of Π d n, where the superscript n means that Pk n Πd n. The orthonormality means that Pk n (x)p j m (x)w(x)dx = δ j,kδ m,n. Ω For each n N 0, the polynomials Pk n, 1 k rd n, form an orthonormal basis of V n. We note that there are many bases of V n ; if Q is an invertible matrix of size rn d, then the components of QP n form another basis of V n which is orthonormal if Q is an orthogonal matrix. For general results on orthogonal polynomials in several variables, including some of the recent development, we refer to the survey [23] and the references there. Ordinary harmonics and h-harmonics. The ordinary harmonic polynomials on R d+1 are the homogeneous polynomials satisfying the Laplace equation P = 0, where = d+1 on Rd+1 and i is the ordinary partial derivative with respect to the i-th coordinate. The spherical harmonics are the restriction of harmonic polynomials on S d. They are orthogonal polynomials with respect to the surface measure dω on S d. The theory of h-harmonics has been established recently by Dunkl (see [3-5]). For a nonzero vector v R d+1 define the reflection σ v by xσ v := x 2(x v)v/ v 2, x R d+1.

4 4 YUAN XU The h-harmonics are homogeneous orthogonal polynomials on S d with respect to h 2 α dω, where the weight function h α is defined by (2.2) h α (x) := m x v i α i, α i 0, with α i = α j whenever σ i is conjugate to σ j in the reflection group G generated by the reflections {σ vi : 1 i m}. The function h α is a positively homogeneous G- invariant function. The key ingredient of the theory is a family of first-order differentialdifference operators, D i (Dunkl s operators), which generates a commutative algebra ([5]). The h-laplacian is defined by h = D D2 d+1, which plays the role of Laplacian in the theory of the ordinary harmonics. In particular, the h-harmonics are the homogeneous polynomials satisfying the equation h Y = 0. The h-spherical harmonics are the restriction of h-harmonics on the sphere. Orthogonal polynomials on spheres for S-symmetric functions. A weight function H defined on R d+1 is called S-symmetric if it is symmetric with respect to y d+1 and centrally symmetric with respect to the variables y = (y 1,..., y d ); i.e., H(y, y d+1 ) = H(y, y d+1 ) and H(y, y d+1 ) = H( y, y d+1 ). In particular, the weight functions of the form H(y) = W(y1 2,..., y2 d+1 ) are S-symmetric. It is proved in [25] that there is a unique decomposition [n/2] (2.3) Pn d+1 = y 2k H n 2k (H), k=0 where H n (H) is the subspace of homogeneous polynomials of degree n which are orthogonal to polynomials of lower degree. If H = h 2, where h is a reflection invariant weight function as in (2.2), then H n (H) is the same as the space of h-harmonics on R d+1 ; in particular, when H = 1, it is the same as the space of ordinary harmonics. It is also proved in [25] that an orthonormal basis of Hn d+1 can be constructed using orthogonal polynomials on B d as follows. Associated to the weight function H we define two weight functions W (1) (x) = H(x, 1 x 2 )/ 1 x 2 and W (2) (x) = H(x, 1 x 2 ) 1 x 2 on B d. We denote by {Pk n} and {Qn k } systems of orthonormal polynomials with respect to the weight functions W (1) and W (2), respectively, where we keep the convention that the superscript n means that Pk n and Qn k are polynomials in Πd n, and the subindex k has the range 1 k rn d. Keeping in mind the notation (2.1) we define (2.4) Y (1) k,n (y) = rn Pk n (x) and Y (2) j,n (y) = rn x d+1 Q n 1 j (x),

5 ORTHOGONAL POLYNOMIALS AND CUBATURE FORMULA 5 where 1 k rn d, 1 j rd (2) n 1 and we define Y j,0 (y) = 0. It is proved in [25] that these polynomials are homogeneous polynomials in y and they form an orthonormal basis for the space H n (H). The group Z 2 Z 2. In the rest of this paper we will reserve the letter G for the group G = Z 2 Z 2 = (Z 2 ) m. It is one of the simplest reflection groups. The elements of G take the form a = (ε 1,..., ε m ), where ε i = ±1. For a function f defined on R m, the action of a G on f is defined by R(a)f(x) = f(xa), x R m, a G; we have R(a)f(x) = f(ε 1 x 1,..., ε m x m ). If R(a)f = f for all a G, we say that f is invariant under G. We take the convention that the size of G agrees with the function it acts upon, so that we do not have to write G = G m. We will take m = d or m = d+1 in subsequent sections. Basic Lemma. We let dω = dω d denote the surface measure on S d, and the surface area ω(s d ) = ω d (S d ) = dω d = 2π (d+1)/2 /Γ((d + 1)/2). S d Let H be defined on R d+1. Assume that H is symmetric with respect to x d+1 ; i.e., H(x, x d+1 ) = H(x, x d+1 ), where x R d. It is proved in [25] that for any continuous function g defined on S d, we have (2.5) g(y)h(y)dω = S d B d [ g(x, 1 x 2 ) + g(x, ] 1 x 2 ) H(x, 1 x 2 )dx/ 1 x 2. This formula connects the integral on B d to S d and it plays an essential role in [25]. The following establishes a relation between the integrations over the unit sphere and over the standard simplex. Lemma 2.1. For any continuous function f defined on S d, (2.6) f(y1 2,..., y2 d+1 )dω = 2 du f(u 1,..., u d, 1 u 1 ) S d Σ d u1...u d (1 u 1 ). Proof. It follows from (2.5) with H = 1 and g(y) = f(y 2 1,..., y2 d+1 ) that f(y1 2,..., y2 d+1 )dω = 2 f(x 2 1,..., x2 d, 1 dx x 2 ). S d B d 1 x 2 Since the function g(x 2 1,..., x2 d ) is invariant under the action of G, we can write its integral over B d as 2 d times the integral over B+ d = {x Bd : x 1 0,..., x d 0}. Upon changing variables u 1 = x 2 1,..., u d = x 2 d in the integral over Bd +, we obtain that (2.7) g(x 2 1,..., x2 d )dx = du g(u 1,..., u d ), B d Σ d u1 u d

6 6 YUAN XU where we have changed the integral over B d + back to the integral over Bd. Putting these two formulae together with g(x) = f(x, 1 x 1 )/ 1 x 1 completes the proof. In the case d = 2 the formula (2.6) has appeared in [3, Lemma 1.6, p. 689]. We note, in particular, that the Lebesgue measure on S d is related to the weight function 1/ (u 1...u d (1 u 1 ) over Σ d. 3. Orthogonal polynomials on spheres and on simplices In this section we discuss the relation between orthogonal polynomials on spheres and on simplices. We need some notations. Definition 3.1. Let H(y) = W(y 2 1,..., y2 d+1 ) be a weight function defined on Rd+1. Associated to H we define a weight function W Σ on the simplex Σ d by (3.1) W Σ (u) = 2W(u 1,..., u d, 1 u 1 )/ u 1 u d (1 u 1 ), u Σ d, and we normalize the weight function W so that, recall the formula (2.6), (3.2) W(y1, 2..., yd+1)dω 2 = W Σ (u)du = 1. S d Σ d It is evident that H is S-symmetric on R d+1. Recall that the space of orthogonal homogeneous polynomials of degree n associated to the weight function H is denoted by H n (H). If a polynomial P in H n (H) is invariant under G, then P has to be of even degree, which means that n is even. We define a subspace H G 2n of H 2n(H) by H G 2n = {P H 2n(H) P(xa) = P(x), a G}; that is, H2n G contains polynomials in H 2n(H) that are invariant under the action of G. Similarly, we define P2n G to be the subspace of Pd+1 2n containing all invariant homogeneous polynomials on R d+1. Let Rk n denote a basis of orthonormal polynomials of degree n with respect to the weight function W Σ on Σ d, where we keep the convention that the superscript n means that Rk n are polynomials in Πd n, and the subindex k has the range 1 k rd n. Keep in mind the notation (2.1). We define (3.3) S 2n k (y) = r 2n R n k(x 2 1,..., x 2 d), 1 k r d n. These polynomials are in fact homogeneous polynomials in y R d+1. Theorem 3.2. Let W be a weight function on R d+1 defined as in Definition 3.1. Then the functions Sk 2n defined in (3.3) are homogeneous polynomials of degree 2n on R d+1, and the set {Sk 2n : 1 k rn d} forms an orthonormal basis for HG 2n. Moreover, n (3.4) P2n G = y 2k H2n 2k G ; k=0

7 ORTHOGONAL POLYNOMIALS AND CUBATURE FORMULA 7 that is, if P Pn G, then there is a unique decomposition P(y) = n y 2n 2k P 2k (y), P 2k H2k. G k=0 Proof. Recall the notation (2.1); we see that S 2n k S 2n k (y) = r2n R n k (x2 1,..., x2 d ) = r2n takes the form β 1 n a β x 2β 1 1 x 2β d d, where a β are constants. Since y i = rx i, it is not hard to see that Sk 2n is a homogeneous polynomial of degree 2n in the variables y R d+1. To show Sk 2n H n (H), we need to show that it is orthogonal to polynomials of degree less than 2n with respect to H(y)dω on S d. We consider P β (y) = y β for β N d+1 and β 1 2n 1. If one of the components of β = (β 1,..., β d+1 ) is odd, say β 1 is odd, then we conclude by symmetry of x 1 x 1 that Sk 2n (y)p β(y)h(y)dω = Rk n (y2 1,..., y2 d )yβ W(y1 2,..., y2 d+1 )dω = 0. S d S d Hence, we only need to consider the case when all components of β are even; in such a case we write β = 2γ for γ N d+1 and γ 1 n 1. Then the basic formula (2.6) shows that S d S 2n k (y)y2γ H(y)dω = Σ d R n k (x)xγ 1 1 xγ d d (1 x 1 x d ) γ d+1 W Σ (x)dx = 0 by the orthogonality of Rk n to polynomials of lower degrees. Since it is evident from (3.3 ) that Sk 2n are invariant under G, we have shown that Sk 2n are elements of H2n G. Moreover, if we replace y 2γ by Sj 2n (y) in the above formula, then it follows from the orthonormality of Rk n that {S2n k : 1 k rn} d forms an orthonormal set. Let P be a homogeneous polynomial on R d+1 that is invariant under G. Then P is necessarily of even degree, say degree 2n, which is of the form P(y) = p(y1 2,..., y2 d+1 ), where p is homogeneous of degree n. Using x 2 d+1 = 1 x2 1...x 2 d whenever possible, we can write P further as P(y) = r 2n q(x 2 1,..., x2 d ), q Πd n, where q is a polynomial of degree n, not homogeneous in general. Since {Rj k, 1 j rn d, 0 k n} forms a basis for all polynomials of degree n, we can write q in terms of. Therefore, we have the unique expansion R k j n r k d n r P(y) = r 2n b j,k Rj k (x2 1,..., x2 d ) = k d b j,k r 2n 2k Sj 2k (y). k=0 j=1 k=0 j=1

8 8 YUAN XU In particular, if P H2n G, then it follows from the orthogonality of P to polynomials of lower degree that P(y) = j a jsj 2n (y). This shows that {Sj n, 1 j rd n } forms an orthonormal basis of H2n, G and the unique decomposition (3.4) holds. As a special case of the theorem, a basis of orthonormal polynomials associated to the weight function W Σ (u) = (u 1 u d (1 u 1 )) 1/2 yields, by (3.3), an orthonormal basis for spherical harmonics that are invariant under G. Since {S n j } forms a basis of HG 2n, it follows that the dimension of HG 2n is rd n. We formulate this as a corollary. Corollary 3.3. Let W be a weight function on R d+1 defined as in Definition 3.1. Then ( ) n + d 1 dim H2n G = n ( ) n + d and dim P2n G =. n For the h-harmonics, the decomposition of P2n G in terms of HG 2n was established in [4] using the differential-difference operators. The importance of the above theorem lies in the relation between orthogonal polynomials on the sphere S d and on the simplex Σ d. In view of the result in [25], we can also establish a relation between orthogonal polynomials on the simplex Σ d and on the ball B d. Let U be a weight function defined on R d. We define a weight function U Σ on Σ d by (3.5) U Σ (u) = U(u 1,..., u d )/ u 1 u d, u Σ d. We denote the space of orthogonal polynomials of degree n with respect to the weight function U(x 2 1,..., x2 d ) on Bd by V n (B d ). If P V n (B d ) is invariant under G, then it must be of even degree. We define a subspace V G 2n (Bd ) of V 2n (B d ) by V G 2n(B d ) = {P V 2n (B d ) P(xa) = P(x), a G}. That is, V 2n (B d ) contains orthogonal polynomials in V 2n (B d ) that are invariant under G. Theorem 3.4. Let U be a weight function defined on R d. Let Rk n be orthonormal polynomials with respect to the weight function U Σ. Then the polynomials Tk 2n(x) = Rk n(x2 1,..., x2 d ) are orthogonal with respect to the weight function U(x2 1,..., x2 d ) on Bd and {Tk 2n, 1 k rd n } forms an orthonormal basis for VG 2n (Bd ). In particular, it follows that dim V2n G (Bd ) = rn d. The theorem follows as a corollary of Theorem 3.2 and the results in [25]. It can also be proved following a similar line as in the proof of Theorem 3.2, using the formula (2.7) instead of (2.6); we leave the details to the reader. These two theorems point out a rather basic relation between orthogonal polynomials on the simplex Σ d and those on the sphere S d and the ball B d. It is rather remarkable that, except for the case discussed in [3], they do not seem to have appeared in the literature. Let us illustrate them by an example, in which we look at the classical orthogonal polynomials.

9 ORTHOGONAL POLYNOMIALS AND CUBATURE FORMULA 9 Example 3.5. Let us consider the weight functions H(y) = H α y 1 2α 1 y 2 2α 2 y 3 2α 3 on S 2, where H α is a normalization constant so that H is normalized as in Definition 3.1. The corresponding weight function W Σ on the simplex Σ d is given by (3.6) W α Σ (x) = 2H αx α 1 1/2 1 x α 2 1/2 2 (1 x 1 x 2 ) α 3 1/2. In an effort to understand Dunkl s theory of h-harmonics, in [24] we study the orthogonal polynomials on S d associated to H(y) = H α y α 1 1 yα d+1 d+1, d 2, in detail. In particular, an orthonormal basis of h-harmonics is given in terms of the orthonormal polynomials D n (λ,µ) of one variable with respect to the measure (1 t 2 ) λ t 2µ on [ 1, 1], which in turn can be written in terms of Jacobi polynomials. For d = 2, using the spherical coordinates y 1 = r cos θ 1, y 2 = r sin θ 1 cos θ 2, y 3 = r sin θ 1 sin θ 2, r = y, some of the h-harmonics of degree n (Y n,(1) k in the notation of [24]) are given by Y n k (y) = An k rn D (k+α 2+α 3 +1/2,α 1 ) n k (cosθ 1 )(sinθ 1 ) k D (α 3,α 2 ) k (cos θ 2 ), where 0 k n and A n k are normalization constants. (It should be mentioned that the spherical coordinates adopted above are in the reverse order of the spherical coordinates used in [24].) On the other hand, an orthonormal basis for W Σ can be given in terms of the Jacobi polynomials P n (α,β) by (cf. [10]) Rk(x n 1, x 2 ) = Bk n P (2k+α 2+α 3,α 1 1/2) n k (2x 1 1)(1 x 1 ) k P (α 3 1/2,α 2 1/2) k (2 x 2 1), 1 x 1 where 0 k n, n 0, and Bk n 3.2, the polynomials are normalization constants. According to Theorem S 2n k (y) = r2n R n k (x2 1, x2 2 ), y = (y 1, y 2, y 3 ) = r(x 1, x 2, x 3 ), are h-harmonics with respect to H. In fact, using the spherical coordinates and the formula D (λ,µ) 2n (cos θ) = const P n (λ 1/2,µ 1/2) (cos2θ) ([24]), we see from the above formulae for Yk n and Rn 2n k that, up to a constant, the polynomials S2n k (y) and Y2k (y) are identical. In particular, this shows that Y2k 2n (y) forms a basis for h-harmonics invariant under G. Moreover, by Theorem 3.4, the polynomials Tk n(x 1, x 2 ) = Rk n(x2 1, x2 2 ) are orthogonal polynomials on B 2 with respect to the weight function U(x 1, x 2 ) = 2H α x 1 2α 1 x 2 2α 2 (1 x 2 1 x2 2 )2α 3 1/2. If α 1 = α 2 = 0, we can use the quadratic transform P (λ,λ) 2k (t) = const P (λ, 1/2) k (2t 2 1) [18, p. 59, (4.1.5)] to conclude that Tk n (x 1, x 2 ) = const P (2k+α 3,2k+α 3 ) 2n 2k (x 1 )(1 x 2 1 )k P (α 3 1/2,α 3 1/2) 2k (x 2 / 1 x 2 1 ).

10 10 YUAN XU Apart from the constant, these polynomials with α = α 3 1/2 are exactly the same as 2P2n,2k α given in [10, formula (3.8), p. 449], where the polynomials 2Pn,k α are given as an orthogonal basis for the weight function (1 x 2 1 x 2 2) α 1/2. The case α 1 = α 2 = α 3 = 0 in the above example is of particular interest, since h-harmonics reduce to the ordinary harmonics when α = 0. Another interesting case is α 1 = α 2 = α 3 ; Dunkl [3] studied the h-harmonics that are invariant under the symmetric group of order three in this case. Among other things, an orthogonal basis invariant under the permutations of x 1, x 2, and 1 x 1 x 2 is given, which leads to h-harmonics invariant under the octahedral group. 4. Cubature formula on spheres and on simplices In this section we discuss the connection between cubature formulae on spheres and on simplices. For a given integral L(f) := fdµ, where dµ is a nonnegative measure with support set on Σ d or B d, a cubature formula of degree M is a linear functional I M (f) = N λ k f(x k ), λ k R, x k R d, k=1 defined on Π d, such that L(f) = I M (f) whenever f Π d M, and L(f ) I M (f ) for at least one f Π d M+1. When the measure is supported on Sd, we need to replace Π d M by Π d+1 M in the above formulation and require x k S d. The points x 1,...,x N in the formula are called nodes and the numbers λ 1,..., λ N are called weights. Such a formula is called minimal, if N, the number of nodes, is minimal among all cubature formulae of degree M. We are interested in the cubature formulae that are invariant under the action of the group G. A linear functional I(f) is called invariant under G if I(R(a)f) = I(f) for all a G. In order to write down an invariant cubature formula explicitly, we need more notation. For u R m, m = d or d+1, we denote its G-orbit by G(u), which is defined by G(u) = {ua a G}; we also denote by G(u) the number of distinct elements in G(u). A cubature formula is invariant under G if the set of its nodes is a union of G-orbits and the nodes belonging to the same G-orbit have the same weight. In the case we are considering, the invariant cubature formula, denoted by IM G (f), takes the form N IM G (f) = λ k f(u i a) / G(u i ) = k=1 a G N k=1 λ k / f(ε 1 u i,1,...ε m u i,m ) G(ui ), ε { 1,1} m where the summation a G f(u ia) is understood as taken over all distinguished elements of G(u i ) as explained in the second equality. We note that G(u) = 2 k, where k is the number of non-zero elements of u in this case. The main result in this section establishes a correspondence between cubature formulae on the simplex Σ d and a Z 2 Z 2 invariant cubature formulae on the sphere S d. Because of their applications in numerical integration and in areas ranging from coding

11 ORTHOGONAL POLYNOMIALS AND CUBATURE FORMULA 11 theory to isometric embeddings between classical Banach spaces, cubature formulae on the sphere have been under intense study for years. One of the important approaches in the construction of cubature formulae is initiated by Sobolev [16], which deals with the invariant cubature formula under a finite group. The fundamental result of Sobolev states that a cubature formula invariant under a group is exact for all polynomials in a subspace P if, and only if, it is exact for all polynomials in P that are invariant under the same group. Since this result helps to reduce the number of polynomials that need to be evaluated by a great deal, it has been used to construct a number of cubature formulae on the unit sphere S 2 (cf. [6, 11, 15, 17] and the references there). In the case we are considering, the results in the previous section shows that the invariant polynomials on S d can be identified with the polynomials on the simplex Σ d under a proper change of variables. Therefore, Sobolev s result suggests that the Z 2 Z 2 invariant cubature formulae are related to the cubature formulae on the simplex. Indeed, we have the following result which does not seem to have been noticed before. Theorem 4.1. Let W Σ be a weight function on Σ d as defined in Definition 3.1. Suppose that there is a cubature formula of degree M on Σ d giving by (4.1) Σ d f(u)w Σ (u)du = N λ i f(u i ), f Π d M, whose N nodes lie on the simplex Σ d ; that is u i = (u i,1,..., u i,d ) and u i 0, 1 u 1 0. Then there is a cubature formula of degree 2M + 1 on the unit sphere S d, (4.2) S d g(y)w(y 2 1,..., y2 d+1 )dω = N λ i a G g(v i a) / G(v i ), g Π d+1 2M+1 where the nodes v i = (v i,1,..., v i,d+1 ) S d are defined in terms of u i by (4.3) v i = (v i,1,..., v i,d+1 ) = ( u i,1,..., u i,d, 1 u i 1 ). On the other hand, if there exists a cubature formula of degree 2M +1 on S d in the form of (4.2), then there is a cubature formula of degree M on the simplex Σ d in the form of (4.1) whose nodes u i Σ d are defined by (4.4) u i = (u i,1,..., u i,d ) = (v 2 i,1,..., v 2 i,d). Proof. Assuming (4.1), we prove (4.2) for the polynomials g β (y) = y β, β N d+1 and β 1 2M + 1. If one of the components of β = (β 1,..., β d ) is odd, say β 1 is odd, then the integral in the left hand side of (4.2) is zero by symmetry. Moreover, let a 0 = ( 1, 1,...1) G; we have that a Gg β (v i a) = a G g β (v i a 0 a 1 0 a) = g β (v i a 0 b) = g β (v i b), b G b G

12 12 YUAN XU which implies that the right hand side of (4.2) is also zero. Hence, we only need to verify (4.2) for g 2γ (y) = y 2γ, γ 1 M. Since g 2γ is clearly invariant under G, the right hand side of (4.2) with g = g 2γ becomes N λ i g 2γ (v i a) / G(v i ) = a G N λ i g 2γ (v i ). Hence, taking into consideration (4.3), it follows from the basic formula (2.6) that (4.2) with g = g 2γ is equivalent to (4.1) with f(u) = u γ 1 1 uγ d d (1 u 1) γ d+1, which holds true since γ 1 M. On the other hand, assuming (4.2), to prove (4.1) we only have to restrict (4.2) to invariant polynomials of the form g(y) = f(y1 2,..., y2 d ), y Rd+1, and use the basic formula (2.6). Although both the statement and the proof of this theorem are simple, its importance is apparent. It allows us to transform back and forth between cubature formulae on Σ d and Z 2 Z 2 invariant cubature formulae on S d for a large class of weight functions. Before we go any further, let us determine the number of nodes in the cubature formula (4.2). Evidently, the number of nodes of (4.2) depends on how many nodes of (4.1) are on the boundary of Σ d. Recall that Σ d is defined by d + 1 inequalities: u 1 0,..., u d 0 and 1 u 1 0 for u R d. A k-dimensional face of Σ d, 0 k d, contains elements of Σ d for which exactly d k inequalities become equalities. In particular, if k = d, then none of the inequalities becomes equality, so that the (unique) d-dimensional face of Σ d is the interior of Σ d. We also note that a 0-dimensional face is one of the vertices of the simplex Σ d. For 0 k d, we define (4.5) N k = number of nodes in (4.1) which lie on k-dimensional faces. Then the number of nodes of (4.1) is equal to N 0 + N N d. We have the following corollary. Corollary 4.2. Let N M (Σ d ) and N 2M+1 (S d ) denote the number of nodes of the formula (4.1) and the formula (4.2), respectively. Then the following relation holds N M (Σ d ) = N 0 + N N d N 2M+1 (S d ) = 2 d+1 N d N 1 + 2N 0. Proof. If a node u i of (4.1) lies on a k-dimensional face of Σ d, then it follows from the definition that the corresponding node v i in (4.3) has exactly k + 1 nonzero elements, which implies that G(u i ) = 2 k+1. Moreover, if v i is a node of a cubature formula of the form (4.2) which contains exactly k + 1 nonzero elements, then u i in (4.4) lies on a k-dimensional face of Σ d. There are lower bounds on the number of nodes of cubature formulae, which are used to test whether a given cubature formula is minimal. The lower bound for cubature formulae on Σ d is given by ( ) n + d (4.6) N M (Σ d ) := NM min n (Σd ), M = 2n or M = 2n + 1.

13 ORTHOGONAL POLYNOMIALS AND CUBATURE FORMULA 13 This lower bound is well-known (cf. [17]) and it holds for all cubature formulae with respect to weight functions on Σ d or on any other domain that has a positive measure. For cubature formulae of odd degree with respect to the weight function W α Σ on Σ2 in (3.6), there is an improved lower bound ([2]) ( ) n + 2 [ n + 1 ] (4.7) N 2n+1 (Σ 2 ) + n 2 := N min 2n+1(Σ 2 ), which is the same as the lower bound for cubature formulae of degree 2n+1 with respect to a centrally symmetric weight function, such as (1 x 2 ) α on B 2. This bound follows from a general result of Möller ([14]), but the verification ([2, 21]) uses explicit formulae for the orthogonal polynomials with respect to WΣ α. It is not known whether the lower bound (4.7) holds for other weight functions on Σ 2. For cubature formulae on S d we only give the lower bound for odd degree, since the cubature formula (4.2) is of degree 2M + 1. We have ( ) M + d (4.8) N 2M+1 (S d ) 2 := N2M+1 min M (Sd ), (cf. [14,15]) which holds for all centrally symmetric weight functions on S d. We note that for large M, the order of both NM min(σd ) and N2M+1 min (Sd ) is O(M d ). The following corollary states that if the cubature formula (4.1) is close to a minimal formula for M relatively large, then so is the Z 2 Z 2 invariant cubature formula (4.2). Corollary 4.3. Let the notation be the same as in Corollary 4.2. Then N M (Σ d ) = N min M (Σd ) + O(M d 1 ) N 2M+1 (S d ) = N min 2M+1 (Sd ) + O(M d 1 ). Proof. This follows easily from the following asymptotic formula ( ) n + λ = n Γ(n + λ + 1) Γ(λ + 1)Γ(n + 1) = and the formulae (4.6) and (4.8). n λ Γ(λ + 1) (1 + O(n 1 )), λ > 0, It is worth mentioning that an attractive feature of the correspondence in Theorem 4.1 is that the degree of the Z 2 Z 2 invariant cubature formula (4.2) is more than twice the degree of the cubature formula (4.1) on Σ d. Since it is easier to construct cubature formulae of lower degree, the theorem offers an effective method to construct cubature formulae on the sphere. However, most of the existing cubature formulae on Σ d are constructed for the Lebesgue measure dx (see the survey [12]), which corresponds to the measure y 1 y d+1 dω on S d ; and there are few cubature formulae known for the weight function (u 1 u d (1 u 1 u d )) 1/2 on Σ d, which corresponds to the surface measure dω. On the other hand, most of the existing cubature formulae on S d are constructed for the surface measure dω, from which we may obtain cubature formula on Σ d with respect

14 14 YUAN XU to the measure (u 1 u d (1 u 1 u d )) 1/2. Hence, Theorem 4.1 will also help us to obtain new cubature formulae on Σ d. We indicate the power of the method by the following example. Example 4.4. The following is a Z 2 Z 2 Z 2 invariant cubature formula of degree 5 on the unit sphere S 2 with 12 nodes ([1] and [17, p. 296, formula U 3 : 5-1]), (4.9) f(y 1, y 2, y 3 )dω = π [f(±r, ±s, 0) + f(0, ±r, ±s) + f(±s, 0, ±r)], S 3 2 where the sum is taken over all possible choices of (±1, ±1), and r, s are given by r 2 = (5 + 5)/10 and s 2 = (5 5)/10. It is a minimal cubature formula on S 2 by (4.8). We note that r 2 + s 2 = 1. Using Theorem 4.1, we can transform it to a cubature formula of degree 2 on Σ 2, dx 1 dx (4.10) f(x 1, x 2 ) 2 Σ 2 x1 x 2 (1 x 1 x 2 ) = 2π [ f(r 2, s 2 ) + f(0, r 2 ) + f(s 2, 0) ], 3 which is in fact a minimal cubature formula of degree 2 by (4.6). On the other hand, a minimal formula of degree 3 has four nodes; it follows from the general theory of minimal cubature formulae that these nodes are common zeros of two orthogonal polynomials of degree 2 ([14]). Since a basis of the orthogonal polynomials with respect to WΣ 0 is known explicitly (cf. Example 3.5), it is not hard to find such a formula. Let Rk n denote the orthonormal polynomials of degree n given in Example 3.5 with respect to the weight function WΣ 0. We have R 2 0 (x 1, x 2 ) =3(3 30x x 2 1 )/8, R 2 1 (x 1, x 2 ) =3 5( 1 + 7x 1 )( 1 + x 1 + 2x 2 )/4, R 2 2 (x 1, x 2 ) =3 35(1 2x 1 8x 2 + x x 1x 2 + 8x 2 2 )/8, where we take the subscript of Rk n as k = 0, 1, 2 instead of k = 1, 2, 3 as in Example 3.5. The two polynomials of degree 2 that have four common zeros are elements of the subspace spanned by these three polynomials. It turns out, however, that R0 2 and R2 2 have four common zeros. We denote these zeros by (s i, t i ), 1 i 4; they yield a minimal cubature formula of degree 3 in the form 4 (4.11) f(x 1, x 2 )WΣ(x 0 1, x 2 )dx 1 dx 2 = 2π λ i f(s i, t i ), Σ 2 where the nodes (s i, t i ) and the weights λ i are given by s 1 = s 2 = ( )/35, s 3 = s 4 = ( )/35, ( t 1 = ) ( 30 /35, t 2 = ) 30 /35, ( t 3 = ) ( 30 /35, t 4 = ) 30 /35 λ 1 = λ 2 = (18 30)/72, λ 3 = λ 4 = ( )/72,

15 ORTHOGONAL POLYNOMIALS AND CUBATURE FORMULA 15 and t 1, t 2, t 3, t 4 are exactly the zeros of the polynomial 1225t t t 2 40t + 1. We found the nodes and the weights with the help of the computer program Mathematica. For more on the minimal cubature formula, we refer to [14, 15, 21, 22] and the references there. Since all nodes of (4.11) are located inside the triangle Σ 2, it yields a Z 2 Z 2 Z 2 invariant cubature formula of degree 7 with 32 nodes, (4.12) S 2 f(x 1, x 2, x 3 )dω = π 2 4 λ i f(± si, ± t i, ± 1 s i t i ) where the inner sum is over all possible choices of signs. Both cubature formulae (4.11) and (4.12) appear to be new. Since N7 min (S 2 ) = 20 by (4.8), the formula (4.12) is far from a minimal formula. In [1] or [17, p. 299, formula U 3 : 7-2], a Z 2 Z 2 Z 2 invariant formula of degree 7 with 26 nodes is given, which yields via Theorem 4.1, however, a cubature formula of degree 3 on Σ 2 with 7 nodes. In order to construct a Z 2 Z 2 Z 2 invariant cubature formula on S 2 with respect to surface measure, we only need to work with a cubature formula with respect to W 0 Σ. Several methods used to construct cubature formulae for the Lebesgue measure on Σ 2 may be extended to construct cubature formulae for the weight function W 0 Σ, including the method of constructing minimal formulae. At the moment, there is a better understanding of the structure of the cubature formula on the simplex, especially the minimal cubature, and less of the structure of cubature formulae on the sphere. Using Theorem 4.1 in both directions, a number of new cubature formulae on S d as well as on Σ d can be constructed. We will present the bulk of the new formulae obtained through Theorem 4.1 in [7], where, among others, formulae on S 2 that are invariant under the octahedral group (called fully symmetric formulae in [9, 17]) will be constructed by considering symmetric formulae on the triangle Σ 2 (see [13]). As suggested by Theorem 3.4 and Sobolev s principle, there is also a correspondence between cubature formulae on the simplex Σ d and the cubature formulae on the unit ball B d, which we now describe. Let us denote by R d + the positive quadrant of Rd ; that is, if x = (x 1,..., x d ) R d +, then x i 0 for 1 i d. For x R d +, we define x1/2 R d + by x 1/2 = ( x 1,..., x d ). Theorem 4.5. Let U be a weight function defined on R d and U Σ be defined as in (3.5). If there is a cubature formula of degree M on Σ d giving by (4.13) Σ d f(u)u Σ (u)du = N λ i f(u i ), f Π d M, with all u i R d +, then there is a cubature formula of degree 2M + 1 on the unit ball B d, (4.14) B d g(x)u(x 2 1,..., x 2 d)dx = N λ i a G g(u 1/2 i a) / G(u i ), g Π d 2M+1.

16 16 YUAN XU Moreover, a cubature formula of degree 2M + 1 in the form of (4.14) implies a cubature formula of degree M in the form of (4.13). The proof of this theorem follows exactly as that of Theorem 4.1 if we use the basic formula (2.7) instead of (2.6). We notice that (4.13) is the same as (4.1), except that we no longer require its nodes to be on Σ d. If a node of (4.13) is not on Σ d, then the corresponding nodes of (4.14) will be outside of B d. Let N k be the number of nodes of (4.13) which lie on k-dimensional faces of Σ d as in (4.5). Then we also have Corollary 4.6. Let N M (Σ d ) and N 2M+1 (B d ) denote the number of nodes of the formula (4.13) and the formula (4.14), respectively. Then the following relation holds N M (Σ d ) = N 0 + N N d N 2M+1 (B d ) = 2 d N d N 1 + N 0. A corollary similar to that of Corollary 4.3 also holds. We note that in the correspondence of Theorem 4.5, the Lebesgue measure on Σ d corresponds to the measure x 1 x d dx on B d, and the Lebesgue measure on B d corresponds to the measure du/ u 1 u d on Σ d. In view of Theorem 4.1, another interesting case is the correspondence between dx/ 1 x 2 on B d and du/(u 1 u d (1 u 1 )) 1/2 on Σ d. Example 4.7. Let us consider the cubature formula of degree 2 in (4.10) of Example 4.4; by Theorem 4.5 we can transform it to dx 1 dx f(x 1, x 2 ) 2 = π [ ] B 2 1 x 2 f(±r, ±s) + f(0, ±r) + f(±s, 0), 1 x where the summations are over all possible choices of signs. This is a cubature formula of degree 5 with 8 nodes for W B (x) = 1/ 1 x 2. The number of nodes of this formula is one more than the lower bound (4.7); we have remarked that (4.7) holds for cubature formulae with respect to W B on B 2. Using Theorem 4.5, we also conclude that the cubature formula (4.11) yields a Z 2 Z 2 Z 2 invariant cubature formula of degree 7 with respect to W B on B 2 with 16 nodes. References 1. J. Albrecht and L. Collatz,, Zur numerischen auswertung mehrdimensionaler Integrale, Z. Angew. Math. Mech. 38 (1958), H. Berens and H. J. Schmid, On the number of nodes of odd degree cubature formulae for integrals with Jacobi weights on a simplex, Numerical Integration (T. O. Espelid and A. Genz, ed.), NATO ASI Series, vol. 367, Kluwer Acad. Publisher, 1992, pp C. Dunkl, Orthogonal polynomials with symmetry of order three, Can. J. Math. 36 (1984), C. Dunkl, Orthogonal polynomials on the sphere with octahedral symmetry, Trans. Amer. Math. Soc. 282 (1984), C. Dunkl, Differential-difference operators associated to reflection groups, Trans. Amer. Math. Soc. 311 (1989), H. Engles, Numerical quadrature and cubature, Academic Press, New York, S. Heo and Y. Xu, Constructing cubature formulae for the sphere and for the triangle, to appear. 8. E.G. Kalnins, W. Miller, Jr., and M. V. Tratnik, Families of orthogonal and biorthogonal polynomials on the n-sphere, SIAM J. Math. Anal. 22 (1991),

17 ORTHOGONAL POLYNOMIALS AND CUBATURE FORMULA P. Keast and J. C. Diaz, Fully symmetric integration formulas for the surface of the sphere in s dimensions, SIAM J. Numer. Anal. 20 (1983), T. Koornwinder, Two-variable analogues of the classical orthogonal polynomials, Theory and Application of Special Functions (R. A. Askey, ed.), Academic Press, New York, 1975, pp V.I. Lebedev, Quadrature on a sphere, USSR Comp. Math. and Math. Phys. 16 (1976), J.N. Lyness and R. Cools, A survey of numerical cubature over triangles, Mathematics of Computation : a half-century of computational mathematics (Vancouver, BC, 1993), Proc. Sympos. Appl. Math., 48, Amer. Math. Soc., Providence, RI, 1994, pp J.N. Lyness and D. Jespersen, Moderate degree symmetric quadrature rules for the triangle, J. Inst. Math. Anal. Appl. 15 (1975), H. M. Möller, Kubaturformeln mit minimaler Knotenzahl, Numer. Math. 35 (1976), I. P. Mysovskikh, Interpolatory cubature formulas, Nauka, Moscow, (Russian) 16. S. L. Sobolev, Cubature formulas on the sphere invariant under finite groups of rotations, Sov. Math. Dokl. 3 (1962), A. Stroud, Approximate calculation of multiple integrals, Prentice Hall, Englewood Cliffs, NJ, H. Szegő, Orthogonal polynomials, 4th ed., Amer. Math. Soc. Colloq. Publ. vol.23, Providence, RI, N. J. Vilenkin, Special functions and the theory of group representations, Amer. Math. Soc. Transl. of Math. Monographs, vol. 22, Providence, RI, Y. Xu, On zeros of multivariate quasi-orthogonal polynomials and Gaussian cubature formulae, SIAM J. Math. Anal. 25 (1994), Y. Xu, On a class of bivariate orthogonal polynomials and cubature formula, Numer. Math. 69 (1994), Y. Xu, Common zeros of polynomials in several variables and higher dimensional quadrature, Pitman Research Notes in Mathematics Series, Longman, Essex, Y. Xu, On orthogonal polynomials in several variables, Special functions, q-series, and related topics, Fields Institute Communications, vol. 14, 1997, pp Y. Xu, Orthogonal polynomials for a family of product weight functions on the spheres, Canadian J. Math. 49 (1997), Y. Xu, Orthogonal polynomials and cubature formulae on spheres and on balls, SIAM J. Math. Anal. (to appear). Yuan Xu, Department of Mathematics, University of Oregon, Eugene, Oregon address:

### BIVARIATE LAGRANGE INTERPOLATION AT THE PADUA POINTS: THE IDEAL THEORY APPROACH

BIVARIATE LAGRANGE INTERPOLATION AT THE PADUA POINTS: THE IDEAL THEORY APPROACH LEN BOS, STEFANO DE MARCHI, MARCO VIANELLO, AND YUAN XU Abstract The Padua points are a family of points on the square [,

### Numerical integration formulas of degree two

Applied Numerical Mathematics 58 (2008) 1515 1520 www.elsevier.com/locate/apnum Numerical integration formulas of degree two ongbin Xiu epartment of Mathematics, Purdue University, West Lafayette, IN 47907,

### Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product

Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )

### TOPOLOGICAL GROUPS MATH 519

TOPOLOGICAL GROUPS MATH 519 The purpose of these notes is to give a mostly self-contained topological background for the study of the representations of locally compact totally disconnected groups, as

### Orthogonal Polynomials, Quadratures & Sparse-Grid Methods for Probability Integrals

1/33 Orthogonal Polynomials, Quadratures & Sparse-Grid Methods for Probability Integrals Dr. Abebe Geletu May, 2010 Technische Universität Ilmenau, Institut für Automatisierungs- und Systemtechnik Fachgebiet

### Shayne Waldron ABSTRACT. It is shown that a linear functional on a space of functions can be described by G, a

UNIVERSITY OF WISCONSIN-MADISON CENTER FOR THE MATHEMATICAL SCIENCES Symmetries of linear functionals Shayne Waldron (waldron@math.wisc.edu) CMS Technical Summary Report #95-04 October 1994 ABSTRACT It

### Fully Symmetric Interpolatory Rules for Multiple Integrals over Infinite Regions with Gaussian Weight

Fully Symmetric Interpolatory ules for Multiple Integrals over Infinite egions with Gaussian Weight Alan Genz Department of Mathematics Washington State University ullman, WA 99164-3113 USA genz@gauss.math.wsu.edu

### Inner products. Theorem (basic properties): Given vectors u, v, w in an inner product space V, and a scalar k, the following properties hold:

Inner products Definition: An inner product on a real vector space V is an operation (function) that assigns to each pair of vectors ( u, v) in V a scalar u, v satisfying the following axioms: 1. u, v

### A RECONSTRUCTION FORMULA FOR BAND LIMITED FUNCTIONS IN L 2 (R d )

PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 127, Number 12, Pages 3593 3600 S 0002-9939(99)04938-2 Article electronically published on May 6, 1999 A RECONSTRUCTION FORMULA FOR AND LIMITED FUNCTIONS

### ON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT

ON A WEIGHTED INTERPOLATION OF FUNCTIONS WITH CIRCULAR MAJORANT Received: 31 July, 2008 Accepted: 06 February, 2009 Communicated by: SIMON J SMITH Department of Mathematics and Statistics La Trobe University,

### Polynomial approximation via de la Vallée Poussin means

Summer School on Applied Analysis 2 TU Chemnitz, September 26-3, 2 Polynomial approximation via de la Vallée Poussin means Lecture 2: Discrete operators Woula Themistoclakis CNR - National Research Council

### A linear algebra proof of the fundamental theorem of algebra

A linear algebra proof of the fundamental theorem of algebra Andrés E. Caicedo May 18, 2010 Abstract We present a recent proof due to Harm Derksen, that any linear operator in a complex finite dimensional

### A lower estimate for the Lebesgue constants of linear means of Laguerre expansions

A lower estimate for the Lebesgue constants of linear means of Laguerre expansions George Gasper 1 and Walter Trebels 2 Dedicated to P. L. Butzer on the occasion of his 70th birthday in gratitude (Oct.

### DIAGONAL TOEPLITZ OPERATORS ON WEIGHTED BERGMAN SPACES

DIAGONAL TOEPLITZ OPERATORS ON WEIGHTED BERGMAN SPACES TRIEU LE Abstract. In this paper we discuss some algebraic properties of diagonal Toeplitz operators on weighted Bergman spaces of the unit ball in

### Fredholmness of Some Toeplitz Operators

Fredholmness of Some Toeplitz Operators Beyaz Başak Koca Istanbul University, Istanbul, Turkey IWOTA, 2017 Preliminaries Definition (Fredholm operator) Let H be a Hilbert space and let T B(H). T is said

### Recurrence Relations and Fast Algorithms

Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations

### MULTIPLE SOLUTIONS FOR THE p-laplace EQUATION WITH NONLINEAR BOUNDARY CONDITIONS

Electronic Journal of Differential Equations, Vol. 2006(2006), No. 37, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) MULTIPLE

### Convexity Properties of The Cone of Nonnegative Polynomials

Convexity Properties of The Cone of Nonnegative Polynomials Grigoriy Blekherman November 00 Abstract We study metric properties of the cone of homogeneous non-negative multivariate polynomials and the

### Gaussian interval quadrature rule for exponential weights

Gaussian interval quadrature rule for exponential weights Aleksandar S. Cvetković, a, Gradimir V. Milovanović b a Department of Mathematics, Faculty of Mechanical Engineering, University of Belgrade, Kraljice

### An introduction to some aspects of functional analysis

An introduction to some aspects of functional analysis Stephen Semmes Rice University Abstract These informal notes deal with some very basic objects in functional analysis, including norms and seminorms

### Notes on Special Functions

Spring 25 1 Notes on Special Functions Francis J. Narcowich Department of Mathematics Texas A&M University College Station, TX 77843-3368 Introduction These notes are for our classes on special functions.

### b n x n + b n 1 x n b 1 x + b 0

Math Partial Fractions Stewart 7.4 Integrating basic rational functions. For a function f(x), we have examined several algebraic methods for finding its indefinite integral (antiderivative) F (x) = f(x)

### REAL AND COMPLEX ANALYSIS

REAL AND COMPLE ANALYSIS Third Edition Walter Rudin Professor of Mathematics University of Wisconsin, Madison Version 1.1 No rights reserved. Any part of this work can be reproduced or transmitted in any

### ADJOINTS, ABSOLUTE VALUES AND POLAR DECOMPOSITIONS

J. OPERATOR THEORY 44(2000), 243 254 c Copyright by Theta, 2000 ADJOINTS, ABSOLUTE VALUES AND POLAR DECOMPOSITIONS DOUGLAS BRIDGES, FRED RICHMAN and PETER SCHUSTER Communicated by William B. Arveson Abstract.

### Implicit Functions, Curves and Surfaces

Chapter 11 Implicit Functions, Curves and Surfaces 11.1 Implicit Function Theorem Motivation. In many problems, objects or quantities of interest can only be described indirectly or implicitly. It is then

### Journal of Inequalities in Pure and Applied Mathematics

Journal of Inequalities in Pure and Applied Mathematics THE EXTENSION OF MAJORIZATION INEQUALITIES WITHIN THE FRAMEWORK OF RELATIVE CONVEXITY CONSTANTIN P. NICULESCU AND FLORIN POPOVICI University of Craiova

### Algebra Exam Syllabus

Algebra Exam Syllabus The Algebra comprehensive exam covers four broad areas of algebra: (1) Groups; (2) Rings; (3) Modules; and (4) Linear Algebra. These topics are all covered in the first semester graduate

### Some notes on Coxeter groups

Some notes on Coxeter groups Brooks Roberts November 28, 2017 CONTENTS 1 Contents 1 Sources 2 2 Reflections 3 3 The orthogonal group 7 4 Finite subgroups in two dimensions 9 5 Finite subgroups in three

### Real Analysis Problems

Real Analysis Problems Cristian E. Gutiérrez September 14, 29 1 1 CONTINUITY 1 Continuity Problem 1.1 Let r n be the sequence of rational numbers and Prove that f(x) = 1. f is continuous on the irrationals.

### Bivariate Lagrange interpolation at the Padua points: The generating curve approach

Journal of Approximation Theory 43 (6) 5 5 www.elsevier.com/locate/jat Bivariate Lagrange interpolation at the Padua points: The generating curve approach Len Bos a, Marco Caliari b, Stefano De Marchi

### Infinite-dimensional Vector Spaces and Sequences

2 Infinite-dimensional Vector Spaces and Sequences After the introduction to frames in finite-dimensional vector spaces in Chapter 1, the rest of the book will deal with expansions in infinitedimensional

### ON WEIGHTED INEQUALITIES FOR FRACTIONAL INTEGRALS OF RADIAL FUNCTIONS. dy, 0 < γ < n. x y

ON WEIGHTED INEQUALITIES FOR FRACTIONAL INTEGRALS OF RADIAL FUNCTIONS PABLO L. DE NÁPOLI, IRENE DRELICHMAN, AND RICARDO G. DURÁN Abstract. We prove a weighted version of the Hardy-Littlewood-Sobolev inequality

### , > --1. [P(" )(x) (1 x) (1 x)dx. ( )(0) t( )P, (cos 0)(sin 0/2)+1/2(cos 0/2) +1/ A TRANSPLANTATION THEOREM FOR JACOBI SERIES

A TRANSPLANTATION THEOREM FOR JACOBI SERIES BY 1. Introduction Let P( )(x) be the Jacobi polynomial of degree n, order (a, ), defined by - [P(" )(x) (1 x) (1 x)dx (1 x)"(1 + x)p(" )(x) (-1) a [(1 x)+"(1

### Changing sign solutions for the CR-Yamabe equation

Changing sign solutions for the CR-Yamabe equation Ali Maalaoui (1) & Vittorio Martino (2) Abstract In this paper we prove that the CR-Yamabe equation on the Heisenberg group has infinitely many changing

### This ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0

Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x

Mathematical Olympiad Training Polynomials Definition A polynomial over a ring R(Z, Q, R, C) in x is an expression of the form p(x) = a n x n + a n 1 x n 1 + + a 1 x + a 0, a i R, for 0 i n. If a n 0,

### Convex Optimization & Parsimony of L p-balls representation

Convex Optimization & Parsimony of L p -balls representation LAAS-CNRS and Institute of Mathematics, Toulouse, France IMA, January 2016 Motivation Unit balls associated with nonnegative homogeneous polynomials

### Week 6 Notes, Math 865, Tanveer

Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,

### UPC condition in polynomially bounded o-minimal structures

Journal of Approximation Theory 132 (2005) 25 33 www.elsevier.com/locate/jat UPC condition in polynomially bounded o-minimal structures Rafał Pierzchała 1 Institute of Mathematics, Jagiellonian University,

### Locally Lipschitzian Guiding Function Method for ODEs.

Locally Lipschitzian Guiding Function Method for ODEs. Marta Lewicka International School for Advanced Studies, SISSA, via Beirut 2-4, 3414 Trieste, Italy. E-mail: lewicka@sissa.it 1 Introduction Let f

### Linear Algebra Massoud Malek

CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product

### Introduction to finite element exterior calculus

Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl

### The De Giorgi-Nash-Moser Estimates

The De Giorgi-Nash-Moser Estimates We are going to discuss the the equation Lu D i (a ij (x)d j u) = 0 in B 4 R n. (1) The a ij, with i, j {1,..., n}, are functions on the ball B 4. Here and in the following

### Aliprantis, Border: Infinite-dimensional Analysis A Hitchhiker s Guide

aliprantis.tex May 10, 2011 Aliprantis, Border: Infinite-dimensional Analysis A Hitchhiker s Guide Notes from [AB2]. 1 Odds and Ends 2 Topology 2.1 Topological spaces Example. (2.2) A semimetric = triangle

### Problems in Linear Algebra and Representation Theory

Problems in Linear Algebra and Representation Theory (Most of these were provided by Victor Ginzburg) The problems appearing below have varying level of difficulty. They are not listed in any specific

### A NOTE ON INVARIANT FINITELY ADDITIVE MEASURES

PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 93, Number 1, January 1985 A NOTE ON INVARIANT FINITELY ADDITIVE MEASURES S. G. DANI1 ABSTRACT. We show that under certain general conditions any

### YALE UNIVERSITY DEPARTMENT OF COMPUTER SCIENCE

fa it) On the Jacobi Polynomials PA ' Gregory Matviyenko Research Report YALETJ/DCS/RR-1076 June 13, 1995 DTI,* QUALITY DEFECTED 3 YALE UNIVERSITY DEPARTMENT OF COMPUTER SCIENCE & k.-. *,«.

### Math Camp Lecture 4: Linear Algebra. Xiao Yu Wang. Aug 2010 MIT. Xiao Yu Wang (MIT) Math Camp /10 1 / 88

Math Camp 2010 Lecture 4: Linear Algebra Xiao Yu Wang MIT Aug 2010 Xiao Yu Wang (MIT) Math Camp 2010 08/10 1 / 88 Linear Algebra Game Plan Vector Spaces Linear Transformations and Matrices Determinant

### Hankel Operators plus Orthogonal Polynomials. Yield Combinatorial Identities

Hanel Operators plus Orthogonal Polynomials Yield Combinatorial Identities E. A. Herman, Grinnell College Abstract: A Hanel operator H [h i+j ] can be factored as H MM, where M maps a space of L functions

### ANALYSIS QUALIFYING EXAM FALL 2017: SOLUTIONS. 1 cos(nx) lim. n 2 x 2. g n (x) = 1 cos(nx) n 2 x 2. x 2.

ANALYSIS QUALIFYING EXAM FALL 27: SOLUTIONS Problem. Determine, with justification, the it cos(nx) n 2 x 2 dx. Solution. For an integer n >, define g n : (, ) R by Also define g : (, ) R by g(x) = g n

### Inner product spaces. Layers of structure:

Inner product spaces Layers of structure: vector space normed linear space inner product space The abstract definition of an inner product, which we will see very shortly, is simple (and by itself is pretty

### MINIMAL REPRESENTING MEASURES ARISING FROM RANK-INCREASING MOMENT MATRIX EXTENSIONS

J. OPERATOR THEORY 42(1999), 425 436 c Copyright by Theta, 1999 MINIMAL REPRESENTING MEASURES ARISING FROM RANK-INCREASING MOMENT MATRIX EXTENSIONS LAWRENCE A. FIALKOW Communicated by Florian-Horia Vasilescu

### Linear Algebra. Workbook

Linear Algebra Workbook Paul Yiu Department of Mathematics Florida Atlantic University Last Update: November 21 Student: Fall 2011 Checklist Name: A B C D E F F G H I J 1 2 3 4 5 6 7 8 9 10 xxx xxx xxx

### A Simple Proof of the Generalized Cauchy s Theorem

A Simple Proof of the Generalized Cauchy s Theorem Mojtaba Mahzoon, Hamed Razavi Abstract The Cauchy s theorem for balance laws is proved in a general context using a simpler and more natural method in

### Introduction to Topology

Introduction to Topology Randall R. Holmes Auburn University Typeset by AMS-TEX Chapter 1. Metric Spaces 1. Definition and Examples. As the course progresses we will need to review some basic notions about

### The Symmetric Space for SL n (R)

The Symmetric Space for SL n (R) Rich Schwartz November 27, 2013 The purpose of these notes is to discuss the symmetric space X on which SL n (R) acts. Here, as usual, SL n (R) denotes the group of n n

### Problem 1A. Suppose that f is a continuous real function on [0, 1]. Prove that

Problem 1A. Suppose that f is a continuous real function on [, 1]. Prove that lim α α + x α 1 f(x)dx = f(). Solution: This is obvious for f a constant, so by subtracting f() from both sides we can assume

### PROBLEMS. (b) (Polarization Identity) Show that in any inner product space

1 Professor Carl Cowen Math 54600 Fall 09 PROBLEMS 1. (Geometry in Inner Product Spaces) (a) (Parallelogram Law) Show that in any inner product space x + y 2 + x y 2 = 2( x 2 + y 2 ). (b) (Polarization

### ON TYPES OF MATRICES AND CENTRALIZERS OF MATRICES AND PERMUTATIONS

ON TYPES OF MATRICES AND CENTRALIZERS OF MATRICES AND PERMUTATIONS JOHN R. BRITNELL AND MARK WILDON Abstract. It is known that that the centralizer of a matrix over a finite field depends, up to conjugacy,

### Numerische Mathematik

Numer. Math. (1998) 80: 39 59 Numerische Mathematik c Springer-Verlag 1998 Electronic Edition Numerical integration over a disc. A new Gaussian quadrature formula Borislav Bojanov 1,, Guergana Petrova,

### Math 421, Homework #9 Solutions

Math 41, Homework #9 Solutions (1) (a) A set E R n is said to be path connected if for any pair of points x E and y E there exists a continuous function γ : [0, 1] R n satisfying γ(0) = x, γ(1) = y, and

### On spectral analysis in varieties containing the solutions of inhomogeneous linear functional equations

On spectral analysis in varieties containing the solutions of inhomogeneous linear functional equations arxiv:1704.04753v1 [math.cv] 16 Apr 2017 Gergely Kiss gergely.kiss@uni.lu Csaba Vincze csvincze@science.unideb.hu

### 1. General Vector Spaces

1.1. Vector space axioms. 1. General Vector Spaces Definition 1.1. Let V be a nonempty set of objects on which the operations of addition and scalar multiplication are defined. By addition we mean a rule

### ON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES

Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 29, 2004, 167 176 ON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES Piotr Haj lasz and Jani Onninen Warsaw University, Institute of Mathematics

### HAIYUN ZHOU, RAVI P. AGARWAL, YEOL JE CHO, AND YONG SOO KIM

Georgian Mathematical Journal Volume 9 (2002), Number 3, 591 600 NONEXPANSIVE MAPPINGS AND ITERATIVE METHODS IN UNIFORMLY CONVEX BANACH SPACES HAIYUN ZHOU, RAVI P. AGARWAL, YEOL JE CHO, AND YONG SOO KIM

### disc f R 3 (X) in K[X] G f in K irreducible S 4 = in K irreducible A 4 in K reducible D 4 or Z/4Z = in K reducible V Table 1

GALOIS GROUPS OF CUBICS AND QUARTICS IN ALL CHARACTERISTICS KEITH CONRAD 1. Introduction Treatments of Galois groups of cubic and quartic polynomials usually avoid fields of characteristic 2. Here we will

### Wavelet analysis as a p adic spectral analysis

arxiv:math-ph/0012019v3 23 Feb 2001 Wavelet analysis as a p adic spectral analysis S.V.Kozyrev February 3, 2008 Institute of Chemical Physics, Russian Academy of Science Abstract New orthonormal basis

### Multidimensional Riemann derivatives

STUDIA MATHEMATICA 235 (1) (2016) Multidimensional Riemann derivatives by J. Marshall Ash and Stefan Catoiu (Chicago, IL) Abstract. The well-known concepts of nth Peano, Lipschitz, Riemann, Riemann Lipschitz,

### OPERATOR THEORY ON HILBERT SPACE. Class notes. John Petrovic

OPERATOR THEORY ON HILBERT SPACE Class notes John Petrovic Contents Chapter 1. Hilbert space 1 1.1. Definition and Properties 1 1.2. Orthogonality 3 1.3. Subspaces 7 1.4. Weak topology 9 Chapter 2. Operators

### GENERAL QUARTIC-CUBIC-QUADRATIC FUNCTIONAL EQUATION IN NON-ARCHIMEDEAN NORMED SPACES

U.P.B. Sci. Bull., Series A, Vol. 72, Iss. 3, 200 ISSN 223-7027 GENERAL QUARTIC-CUBIC-QUADRATIC FUNCTIONAL EQUATION IN NON-ARCHIMEDEAN NORMED SPACES M. Eshaghi Gordji, H. Khodaei 2, R. Khodabakhsh 3 The

### Mercer s Theorem, Feature Maps, and Smoothing

Mercer s Theorem, Feature Maps, and Smoothing Ha Quang Minh, Partha Niyogi, and Yuan Yao Department of Computer Science, University of Chicago 00 East 58th St, Chicago, IL 60637, USA Department of Mathematics,

### Representation theory and quantum mechanics tutorial Spin and the hydrogen atom

Representation theory and quantum mechanics tutorial Spin and the hydrogen atom Justin Campbell August 3, 2017 1 Representations of SU 2 and SO 3 (R) 1.1 The following observation is long overdue. Proposition

### Exercise Solutions to Functional Analysis

Exercise Solutions to Functional Analysis Note: References refer to M. Schechter, Principles of Functional Analysis Exersize that. Let φ,..., φ n be an orthonormal set in a Hilbert space H. Show n f n

### The coincidence Nielsen number for maps into real projective spaces

F U N D A M E N T A MATHEMATICAE 140 (1992) The coincidence Nielsen number for maps into real projective spaces by Jerzy J e z i e r s k i (Warszawa) Abstract. We give an algorithm to compute the coincidence

### PUTNAM TRAINING PROBLEMS

PUTNAM TRAINING PROBLEMS (Last updated: December 3, 2003) Remark This is a list of Math problems for the NU Putnam team to be discussed during the training sessions Miguel A Lerma 1 Bag of candies In a

### NEW BOUNDS FOR TRUNCATION-TYPE ERRORS ON REGULAR SAMPLING EXPANSIONS

NEW BOUNDS FOR TRUNCATION-TYPE ERRORS ON REGULAR SAMPLING EXPANSIONS Nikolaos D. Atreas Department of Mathematics, Aristotle University of Thessaloniki, 54006, Greece, e-mail:natreas@auth.gr Abstract We

### Diophantine equations for second order recursive sequences of polynomials

Diophantine equations for second order recursive sequences of polynomials Andrej Dujella (Zagreb) and Robert F. Tichy (Graz) Abstract Let B be a nonzero integer. Let define the sequence of polynomials

### EASY PUTNAM PROBLEMS

EASY PUTNAM PROBLEMS (Last updated: December 11, 2017) Remark. The problems in the Putnam Competition are usually very hard, but practically every session contains at least one problem very easy to solve

### 2 (Bonus). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?

MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due 9/5). Prove that every countable set A is measurable and µ(a) = 0. 2 (Bonus). Let A consist of points (x, y) such that either x or y is

### Real Analysis Notes. Thomas Goller

Real Analysis Notes Thomas Goller September 4, 2011 Contents 1 Abstract Measure Spaces 2 1.1 Basic Definitions........................... 2 1.2 Measurable Functions........................ 2 1.3 Integration..............................

### STRONG SINGULARITY OF SINGULAR MASAS IN II 1 FACTORS

Illinois Journal of Mathematics Volume 5, Number 4, Winter 2007, Pages 077 084 S 009-2082 STRONG SINGULARITY OF SINGULAR MASAS IN II FACTORS ALLAN M. SINCLAIR, ROGER R. SMITH, STUART A. WHITE, AND ALAN

### CHAPTER 0 PRELIMINARY MATERIAL. Paul Vojta. University of California, Berkeley. 18 February 1998

CHAPTER 0 PRELIMINARY MATERIAL Paul Vojta University of California, Berkeley 18 February 1998 This chapter gives some preliminary material on number theory and algebraic geometry. Section 1 gives basic

### ON BASIC EMBEDDINGS OF COMPACTA INTO THE PLANE

ON BASIC EMBEDDINGS OF COMPACTA INTO THE PLANE Abstract. A compactum K R 2 is said to be basically embedded in R 2 if for each continuous function f : K R there exist continuous functions g, h : R R such

### Error formulas for divided difference expansions and numerical differentiation

Error formulas for divided difference expansions and numerical differentiation Michael S. Floater Abstract: We derive an expression for the remainder in divided difference expansions and use it to give

### Journal of Inequalities in Pure and Applied Mathematics

Journal of Inequalities in Pure and Applied Mathematics ON SOME APPROXIMATE FUNCTIONAL RELATIONS STEMMING FROM ORTHOGONALITY PRESERVING PROPERTY JACEK CHMIELIŃSKI Instytut Matematyki, Akademia Pedagogiczna

### Auerbach bases and minimal volume sufficient enlargements

Auerbach bases and minimal volume sufficient enlargements M. I. Ostrovskii January, 2009 Abstract. Let B Y denote the unit ball of a normed linear space Y. A symmetric, bounded, closed, convex set A in

### Kernel Method: Data Analysis with Positive Definite Kernels

Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University

### A fixed point approach to orthogonal stability of an Additive - Cubic functional equation

Int. J. Adv. Appl. Math. and Mech. 3(4 (06 8 (ISSN: 347-59 Journal homepage: www.ijaamm.com IJAAMM International Journal of Advances in Applied Mathematics and Mechanics A fixed point approach to orthogonal

### arxiv: v1 [math.ca] 7 Aug 2015

THE WHITNEY EXTENSION THEOREM IN HIGH DIMENSIONS ALAN CHANG arxiv:1508.01779v1 [math.ca] 7 Aug 2015 Abstract. We prove a variant of the standard Whitney extension theorem for C m (R n ), in which the norm

### 1. Use the properties of exponents to simplify the following expression, writing your answer with only positive exponents.

Math120 - Precalculus. Final Review. Fall, 2011 Prepared by Dr. P. Babaali 1 Algebra 1. Use the properties of exponents to simplify the following expression, writing your answer with only positive exponents.

### 4.3 - Linear Combinations and Independence of Vectors

- Linear Combinations and Independence of Vectors De nitions, Theorems, and Examples De nition 1 A vector v in a vector space V is called a linear combination of the vectors u 1, u,,u k in V if v can be

### 2) Let X be a compact space. Prove that the space C(X) of continuous real-valued functions is a complete metric space.

University of Bergen General Functional Analysis Problems with solutions 6 ) Prove that is unique in any normed space. Solution of ) Let us suppose that there are 2 zeros and 2. Then = + 2 = 2 + = 2. 2)

### ETNA Kent State University

Electronic Transactions on Numerical Analysis. Volume 9, 1999, pp. 39-52. Copyright 1999,. ISSN 1068-9613. ETNA QUADRATURE FORMULAS FOR RATIONAL FUNCTIONS F. CALA RODRIGUEZ, P. GONZALEZ VERA, AND M. JIMENEZ

### GALOIS GROUPS OF CUBICS AND QUARTICS (NOT IN CHARACTERISTIC 2)

GALOIS GROUPS OF CUBICS AND QUARTICS (NOT IN CHARACTERISTIC 2) KEITH CONRAD We will describe a procedure for figuring out the Galois groups of separable irreducible polynomials in degrees 3 and 4 over

### RENORMINGS OF L p (L q )

RENORMINGS OF L p (L q ) by arxiv:math/9804002v1 [math.fa] 1 Apr 1998 R. Deville, R. Gonzalo and J. A. Jaramillo (*) Laboratoire de Mathématiques, Université Bordeaux I, 351, cours de la Libération, 33400

### On integral representations of q-gamma and q beta functions

On integral representations of -gamma and beta functions arxiv:math/3232v [math.qa] 4 Feb 23 Alberto De Sole, Victor G. Kac Department of Mathematics, MIT 77 Massachusetts Avenue, Cambridge, MA 239, USA

### THE DIRICHLET PROBLEM ON QUADRATIC SURFACES

THE DIRICHLET PROBLEM ON QUADRATIC SURFACES SHELDON AXLER, PAMELA GORKIN, AND KARL VOSS Abstract. We give a fast, exact algorithm for solving Dirichlet problems with polynomial boundary functions on quadratic

### Introduction to Signal Spaces

Introduction to Signal Spaces Selin Aviyente Department of Electrical and Computer Engineering Michigan State University January 12, 2010 Motivation Outline 1 Motivation 2 Vector Space 3 Inner Product