A collocation method for solving some integral equations in distributions
|
|
- Nickolas Butler
- 5 years ago
- Views:
Transcription
1 A collocation method for solving some integral equations in distributions Sapto W. Indratno Department of Mathematics Kansas State University, Manhattan, KS , USA A G Ramm Department of Mathematics Kansas State University, Manhattan, KS , USA ramm@math.ksu.edu Abstract A collocation method is presented for numerical solution of a typical integral equation Rh := R(x, y)h(y)dy = f(x), x D of the class D R, whose kernels are of positive rational functions of arbitrary selfadjoint elliptic operators defined in the whole space R r, D R r is a bounded domain. Several numerical examples are given to demonstrate the efficiency stability of the proposed method. MSC: 45A05, 45P05, 46F05, 62M40, 65R20, 74H15 Key words: integral equations in distributions, signal estimation, collocation method. 1 Introduction In [4] a general theory of integral equations of the class R was developed. The integral equations of the class R are written in the following form: Rh := R(x, y)h(y)dy = f(x), x D := D Γ, (1) D where D R r is a (bounded) domain with a (smooth) boundary Γ. Here the kernel R(x, y) has the following form [4, 5, 6, 7]: R(x, y) = P (λ)q (λ)φ(x, y, λ)dρ(λ), (2) Λ where P (λ), Q(λ) > 0 are polynomials, degp = p, degq = q, q > p, Φ, ρ, Λ are the spectral kernel, spectral measure, spectrum of a selfadjoint 1
2 elliptic operator L on L 2 (R r ) of order s. It was also proved in [4] that R : Ḣ α (D) H α (D) is an isomorphism, where H α (D) is the Sobolev space Ḣ α (D) its dual space with respect to the L 2 (D) inner product, α = s(q p) 2. In this paper we consider a particular type of integral equations of the class R with D = (, 1), r = 1, L = i, :=, Λ (, ), dρ(λ) = dλ, Φ(x, y, λ) = eiλ(x y) 2π, P (λ) = 1, Q(λ) = λ2 +1 2, s = 1, p = 0,q = 2, α = 1, i.e., Rh(x) := d dx e x y h(y)dy = f(x), (3) where h Ḣ ([, 1]) f H 1 ([, 1]). We denote the inner product norm in H 1 ([, 1]) by u, v 1 := ( ) u(x)v(x) + u (x)v (x)dx dx u, v H 1 ([, 1]), (4) u 2 1 := ( u(x) 2 + u (x) 2) dx, (5) respectively, where the primes denote derivatives the bar sts for complex conjugate. If u v are real-valued functions in H 1 ([, 1]) then the bar notations given in (4) can be dropped. Note that if f is a complex valued function then solving equation (3) is equivalent to solving the equations: e x y h k (y)dy = f k (x), k = 1, 2, (6) where h 1 (x) :=Reh(x), h 2 (x) :=Imh(x), f 1 (x) :=Ref(x), f 2 (x) :=Imf(x) h(x) = h 1 (x)+ih 2 (x), i =. Therefore, without loss of generality we assume throughout that f(x) is real-valued. It was proved in [5] that the operator R defined in (3) is an isomorphism between Ḣ ([, 1]) H 1 ([, 1]). Therefore, problem (3) is well posed in the sense that small changes in the data f(x) in the H 1 ([, 1]) norm will result in small in Ḣ ([, 1]) norm changes to the solution h(y). Moreover, the solution to (3) can be written in the following form: where h(x) = a δ(x + 1) + a 0 δ(x 1) + g(x), (7) a := f() f () 2, a 0 := f (1) + f(1), (8) 2 g(x) := f (x) + f(x), (9) 2 δ(x) is the delta function. Here throughout this paper we assume that f C α ([, 1]), α 2. It follows from (8) that h(x) = g(x) if only if f() = f () f(1) = f (1). 2
3 In [6, 7] the problem of solving equation (3) numerically have been posed solved. The least squares method was used in these papers. The goal of this paper is to develop a version of the collocation method which can be applied easily numerically efficiently. In [8] some basic ideas for using collocation method are proposed. In this paper some of these ideas are used new ideas, related to the choice of the basis functions, are introduced. In this paper the emphasis is on the development of methodology for solving basic equation (1) of the estimation theory by a version of the collocation method. The novelty of this version consists in minimization of a discrepancy functional (26), see below. This methodology is illustrated by a detailed analysis applied to solving equation (3), but is applicable to general equations of the class R. One of the goals of this paper is to demonstrate that collocation method can be successfully applied to numerical solution of some integral equations whose solutions are distributions, provided that the theoretical analysis gives sufficient information about the singular part of the solutions. Since f C α ([, 1]), α 2, it follows from (9) that g C([, 1]). Therefore, there exist basis functions ϕ j (x) C([, 1]), j = 1, 2,..., m, such that where max g(x) g m(x) 0 as m, (10) x [,1] g m (x) := m j ϕ j (x), (11) j, j = 1, 2,..., m, are constants. Hence the approximate solution of equation (3) can be represented by h m (x) = δ(x + 1) + c(m) 0 δ(x 1) + g m (x), (12) where j, j =, 0, are constants g m (x) is defined in (11). The basis functions ϕ j play an important role in our method. It is proved in Section 3 that the basis functions ϕ j in (12) can be chosen from the linear B-splines. The usage of the linear B-splines reduces the computation time, because computing (12) at a particular point x requires at most two out of the m basis functions ϕ j. For a more detailed discussion of the family of B-splines we refer to [10]. In Section 2 we derive a method for obtaining the coefficients j, j =, 0, 1,..., m, given in (12). This method is based on solving a finite-dimensional least squares problem ( see equation (33) below ) differs from the usual collocation method discussed in [2] [3]. We approximate f Rh m 2 1 by a quadrature formula. The resulting finite-dimensional linear algebraic system depends on the choice of the basis functions. Using linear B-splines as the basis functions, we prove the existence uniqueness of the solution to this linear algebraic system for all m = m(n) depending on the number n of collocation points used in the left rectangle quadrature rule. The convergence of our collocation method is proved in this Section. An example of the choice of the basis functions which yields the convergence of our version of the collocation method is given in Section 3. In 3
4 Section 4 we give numerical results of applying our method to several problems that discussed in [7]. 2 The collocation method In this Section we derive a collocation method for solving equation (3). From equation (3) we get Rh(x) = a e (x+1) + a 0 e (1 x) x + (e x e y g(y)dy + e x x ) e y g(y)dy = f(x). (13) Assuming that f C 2 ([, 1]) differentiating the above equation, one obtains (Rh) (x) = a e (x+1) + a 0 e (1 x) x + (e x e y g(y)dy e x x ) e y g(y)dy = f (x). (14) Thus, f(x) f (x) are continuous in the interval [, 1]. Let us use the approximate solution given in (12). From (13), (14) (12) we obtain Rh m (x) = e (x+1) + 0 e (1 x) n x + (e x e y ϕ j (y)dy + e x j x ) e y ϕ j (y)dy := f m (x), (15) (Rh m ) (x) = e (x+1) + 0 e (1 x) m x + [e x e y ϕ j (y)dy e x j x ] e y ϕ j (y)dy := (f m ) (x). (16) Thus, Rh m (x) (Rh m ) (x) are continuous in the interval [, 1]. Since f(x) Rh m (x) f (x) (Rh m ) (x) are continuous in the interval [, 1], we may assume throughout that the functions are Riemann-integrable over the interval [, 1]. J 1,m := [f(x) Rh m (x)] 2 (17) J 2,m := [f (x) (Rh m ) (x)] 2 (18) 4
5 Let us define q (x) := e (x+1), q 0 (x) := e (1 x), q j (x) := e x y ϕ j (y)dy, j = 1, 2,..., m, (19) define a mapping H n : C 2 ([, 1]) R n 1 by the formula: φ(x 1 ) φ(x 2 ) H n φ =., φ(x) C2 ([, 1]), (20) φ(x n ) where R n 1 := z 1 z 2. z n Rn z j := z(x j ), z C 2 ([, 1]) (21) x j are some collocation points which will be chosen later. We equip the space R n 1 with the following inner product norm u, v w (n),1 := u 2 w (n),1 := n n w (n) j (u j v j + u jv j), u, v R n 1, (22) w (n) j [u 2 j + (u j) 2 ], u R n 1, (23) respectively, where u j := u(x j ), u j := u (x j ), v j := v(x j ), v j := v (x j ), w j > 0 are some quadrature weights corresponding to the collocation points x j, j = 1, 2,..., n. Applying H n to Rh m, one gets (H n Rh m ) i = e (1+xi) + 0 e (1 xi) + m m j e (y xi) ϕ j (y)dy, i = 1, 2,..., n, x i where m = m(n) is an integer depending on n such that xi j e (xi y) ϕ j (y)dy+ (24) Let m(n) + 2 n, lim m(n) =. (25) n G n ( ) := H n (f Rh m ) 2 w (n),1, (26) 5
6 where f Rh m are defined in (3) (15), respectively, H n defined in (20), c c 0 w (n),1 defined in (23) = c 1 R m+2. Let us choose. c m w (n) j = 2, j = 1, 2,..., n, (27) n x j = + (j 1)s, s := 2, j = 1, 2,..., n, (28) n so that H n (f Rh m ) 2 w (n),1 is the left Riemannian sum of f Rh m 2 1, i.e., f Rh m 2 1 H n (f Rh m ) 2 w (n),1 := δ n 0 as n. (29) Remark 2.1. If J 1,m J 2,m are in C 2 ([, 1]), where J 1,m J 2,m are defined in (17) (18), respectively, then one may replace the weights w (n) j with the weights of the compound trapezoidal rule, get the estimate n δ n = (J 1,m (x) + J 2,m (x))dx w (n) j (J 1,m (x j ) + J 2,m (x j )) (30) 1 3n 2 D J, where δ n is defined in (29) D J := J 1,m(1) + J 2,m(1) (J 1,m() + J 2,m()). (31) Here we have used the following estimate of the compound trapezoidal rule [1, 9]: b n η(x)dx w (n) j η(x j ) (b b a)2 a 12n 2 η (b a)2 (x)dx = a 12n 2 η (b) η (a), (32) where η C 2 ([a, b]). Therefore, if D J C for all m, where C > 0 is a constant, then δ n = O ( ) 1 n. 2 The constants j in the approximate solution h m, see (12), are obtained by solving the following least squares problem: where G n is defined in (26). min G n ( ), (33) 6
7 A necessary condition for the minimum in (33) is n ( 0 = w (n) E m,l l E m,l + E E ) m,l m,l k =, 0, 1,..., m, (34) c k c k where l=1 E m,l := (f Rh m )(x l ), E m,l := (f Rh m ) (x l ), l = 1, 2,..., n. (35) Necessary condition (34) yields the following linear algebraic system (LAS): A m+2 = F m+2, (36) where R m+2, A m+2 is a square, symmetric matrix with the following entries: n (A m+2 ) 1,1 := 2 w (n) l e 2(xl+1), (A m+2 ) 1,2 = 0, (A m+2 ) 1,j := 2 (A m+2 ) 2,2 := 2 (A m+2 ) 2,j := 2 l=1 n l=1 n l=1 n l=1 w (n) l C l,j 2 e (x l+1), j = 3,..., m + 2, w (n) l e 2(1 x l), w (n) l B l,j 2 e (1 x l), j = 3,..., m + 2, (37) (38) (A m+2 ) i,j := n l=1 2w (n) l (B l,i 2 B l,j 2 + C l,i 2 C l,j 2 ), i = 3,..., m + 2, j = i,..., m + 2, (A m+2 ) j,i = (A m+2 ) i,j, i, j = 1, 2,..., m + 2, F m+2 is a vector in R m with the following elements: n (F m+2 ) 1 := w (n) l (f(x l ) f (x l ))e (xl+1) = H n Rh, Hq w (n),1 (F m+2 ) 2 := (F m+2 ) i := l=1 n l=1 n l=1 w (n) l (f(x l ) + f (x l ))e (1 x l) = H n Rh, Hq 0 w (n),1, w (n) l [f(x l )(C l,i 2 + B l,i 2 ) + f (x l )(B l,i 2 C l,i 2 )] = H n Rh, Hq i w (n),1, i = 3,..., m + 2, (39) (40) B l,j := x l e (y x l) ϕ j (y)dy, C l,j := xl e (x l y) ϕ j (y)dy. (41) 7
8 Theorem 2.2. Assume that the vectors H n q j, j =, 0, 1,..., m are linearly independent. Then linear algebraic system (36) is uniquely solvable for all m, where m is an integer depending on n such that (25) holds. Proof. Consider q j H 1 ([, 1]) defined in (19). Using the inner product in R n 1, one gets (A m+2 ) i,j = H n q i 2, H n q j 2 w (n),1, i, j = 1, 2,..., m + 2, (42) i.e., A m+2 is a Gram matrix. We have assumed that the vectors H n q j R n 1, j =, 0, 1,..., m, are linearly independent. Therefore, the determinant of the matrix A m+2 is nonzero. This implies linear algebraic system (36) has a unique solution. Theorem 2.2 is proved. It is possible to choose basis functions ϕ j such that the vectors H n q j, j =, 0, 1,..., m, are linearly independent. An example of such choice of the basis functions is given in Section 3. Lemma 2.3. Let y m := min be the unique minimizer for problem (33). Then G n (y m ) 0 as n, (43) where G n is defined in (26) m is an integer depending on n such that (25) holds. Proof. Let h(x) = a δ(x + 1) + a 0 δ(x 1) + g(x) be the exact solution to (3), Rh = f, where g(x) C([, 1]), define h m (x) = a δ(x + 1) + a 0 δ(x 1) + g m (x), (44) where Choose g m (x) so that m g m (x) := a j ϕ j (x). (45) Then max g(x) g m(x) 0 as m. (46) x [,1] G n (y m ) H n (f R h m ) 2 w (n),1, (47) because y m is the unique minimizer of G n. Let us prove that H(f R h m ) 2 0 as n. Let w (n),1 W 1,m (x) := f(x) R h m (x), W 2,m := f (x) (R h m ) (x). 8
9 Then x W 1,m (x) = e x e y (g(y) g m (y))dy + e x e y (g(y) g m (y))dy (48) x W 2,m (x) = e x e y (g(y) g m (y))dy e x e y (g(y) g m (y))dy. (49) x Thus, the functions [W 1,m (x)] 2 [W 2,m (x)] 2 are Riemann-integrable. Therefore, δ n := f R h m 2 1 H n (f R h m ) 2 w (n),1 0 as n. (50) Formula (50) the triangle inequality yield H n (f R h m ) 2 w (n),1 δ n + f R h m 2 1. (51) Let us derive an estimate for f R h m 2 1. From (48) (49) we obtain the estimates: x ) W 1,m (x) max g(y) g m(y) (e x e y dy + e x e y dy y [,1] x = max g(y) g m(y) [ e x (e x e ) + e x (e x e ) ] (52) y [,1] = max y [,1] g(y) g m(y) [ (2 e x e +x ) ] δ m,1 W 2,m (x) max g(y) g m(y) (e x y [,1] where Therefore, it follows from (52) (53) that x x e y dy + e x x ) e y dy δ m,1, (53) δ m,1 := 2 max g(y) g m(y). (54) y [,1] f R h m 2 1 = W 1,m (x) 2 dx + W 2,m (x) 2 dx 4δm,1, 2 (55) where δ m,1 is defined in (54). Using relation (46), we obtain lim m δ m,1 = 0. Since m = m(n) lim n m(n) =, it follows from (51) (55) that H(f R h m ) 2 w,1 0 as n. This together with (47) imply G n (y m ) 0 as n. Lemma 2.3 is proved. 9
10 Theorem 2.4. Let the vector min := system (36) 0 1. m h m (x) = δ(x + 1) + c(m) 0 δ(x 1) + R m+2 solve linear algebraic m j ϕ j (x). Then h h m H 0 as n. (56) Proof. We have h h m 2 H ([,1]) = R (f Rh m ) 2 H ([,1]) R 2 H 1 ([,1]) Ḣ ([,1]) f Rh m 2 1 ( ) C G n ( f min ) + Rhm 2 1 G n ( min ) C[G n ( min ) + δ n] 0 as n, (57) where C > 0 is a constant Lemma 2.3 was used. Theorem 2.4 is proved. 3 The choice of collocation points basis functions In this section we give an example of the collocation points x i, i = 1, 2,..., n, basis functions ϕ j, j = 1, 2,..., m, such that the vectors H n q j, j =, 0, 1,..., m, are linearly independent, where m + 2 n, H n q j = H n q = e x1 e x2. e xn, H nq 0 = e x1 ey ϕ j (y)dy + e x1 e y ϕ j (y)dy e x2 ey ϕ j (y)dy + e x2 e y ϕ j (y)dy. e xn ey ϕ j (y)dy + e xn e y ϕ j (y)dy e +x1 e +x2. e +xn (58), j = 1, 2,..., m. (59) 10
11 Figure 1: The structure of the basis functions ϕ j Let us choose w (n) j x j as in (27) (28), respectively, with an even number n 6. As the basis functions in C([, 1]) we choose the following linear B-splines: { ψ2 (x) x ϕ 1 (x) = 1 x x 3, 0, otherwise, ψ 1 (x (j 1)), x 2j 3 x x 2j, ϕ j (x) = ψ 2 (x (j 1)), x 2j x x 2j+1, 0, otherwise, j = 2,..., m 1, { ψ1 (x (m 1)), x ϕ m (x) = n x 1, 0, otherwise, (60) where m = n 2 + 1, s := 2 n, ψ 1 (x) := x x 1 +, ψ 2 (x) := (x x 1 ). Here we have chosen x 2j, j = 1, 2,..., n 2 + 1, as the knots of the linear B- splines. From Figure 1 we can see that at each j = 2,..., m, ϕ j (x) is a hat function. The advantage of using these basis functions is the following: at most two basis functions are needed for computing the solution h m (x), because { h m (x) = l, x = x 2l, l = 2,..., n 2, l ϕ l (x) + l+1 ϕ (62) l+1(x), x 2l < x < x 2l+1. (61) 11
12 From the structure of the basis functions ϕ j we have ϕ 1 (x) = 0, x 3 x 1, ϕ j (x) = 0, x x 2j 3 x 2j+1 x 1, j = 2, 3,..., m 1, ϕ m (x) = 0, x x n. (63) Let (H n q j ) i be the i-th element of the vector H n q j, j =, 0, 1,..., m. Then Using (63) in (59), we obtain where (H n q ) i = e (1+xi), i = 1, 2,..., n, (64) (H n q 0 ) i = e (1 xi), i = 1, 2,..., n. (65) +e, i = 1, (H n q 1 ) i = 1 e s, i = 2, e 1 (i)s C 1, 3 i n, (H n q j ) i = e +(i)s D j, 1 i 2j 3, j = 2, 3,..., m 1, (H n q j ) 2j 2 = e 3s e s +, j = 2, 3,..., m 1, (H n q j ) 2j = + e +, j = 2, 3,..., m 1, s (H n q j ) 2j = (H n q j ) 2j 2, j = 2, 3,..., m 1, (H n q i ) j = e 1 (i)s C j, 2j + 1 i n, j = 2, 3,..., m 1, { e (H n q m ) i = +(i)s D m, 1 i n 1, 1 e s, i = n, (66) C 1 := C j := x3 x2j+1 e y ϕ 1 (y)dy = + e, 2es e y ϕ j (y)dy = e+2( 2+j)s ( + e ) 2 x 2j 3, j = 2, 3,..., m 1, (67) D 1 := D j := D m := x3 x2j+1 e y ϕ 1 (y)dy = e( 1 + e ), e y ϕ j (y)dy = e1 2js ( + e ) 2, j = 2, 3,..., m 1, x 2j 3 e y ϕ m (y)dy = + e x n 2es. (68) 12
13 Theorem 3.1. Consider q j defined in (19) with ϕ j defined in (60). Let the collocation points x j, j = 1, 2,..., n, be defined in (28) with an even number n 6. Then the vectors H n q j, j =, 0, 1, 2,..., m, m = s + 1, s = n 2, are linearly independent, where H n is defined in (20). Proof. Let V 0 := {H n q, H n q 0 }, V j := V j {H n q j }, j = 1, 2,..., m. (69) We prove that the elements of the sets V j, j = 0, 1,..., m, are linearly independent. The elements of the set V 0 are linearly independent. Indeed, H n q j 0 j, assuming that there exists a constant α such that H n q = αh n q 0, (70) one gets a contradiction: consider the first the n-th equations of (70), i.e., respectively. It follows from (64), (65) (71) that From (72), (64), (65) (73) it follows that (H n q ) 1 = α(h n q 0 ) 1 (71) (H n q ) n = α(h n q 0 ) n, (72) α = e 2. (73) e 2+s = e 2 s. (74) This is a contradiction, which proves that H n q H n q 0 are linearly independent. Let us prove that the element of the set V j are linearly independent, j = 1, 2, 3,..., m 2. Assume that there exist constants α k, k = 1, 2,..., j + 1, such that H n q j = j k= α k+2 H n q k. (75) Using relations (64)-(66) one can write the (2j 1)-th equation of linear system (75) as follows: (H n q j ) 2j = j k= α k+2 (H n q k ) 2j j = α 1 e (2j 2)s + α 2 e 2+(2j 2)s + e 1 (2j 2)s α k+2 C k. k=1 (76) 13
14 Similarly, by relations (64)-(66) the (n 1)-th n-th equations of linear system (75) can be written in the following expressions: j (H n q j ) n = e + C j = α 1 e 2+ + α 2 e + e + α k+2 C k (77) j (H n q j ) n = e +s C j = α 1 e 2+s + α 2 e s + e +s α k+2 C k, (78) respectively. Multiply (78) by e s compare with (77) to conclude that α 2 = 0. From (78) with α 2 = 0 one obtains k=1 k=1 j α 1 = ec j e α k+2 C k. (79) k=1 Substitute α 1 from (79) α 2 = 0 into (76) get (H n q j ) 2j = e 1 (2j 2)s C j. (80) From (67) (66) one obtains for 0 < s < 1, j = 1, 2, 3,..., m 2, the following relation e 1 (2j 2)s C j (H n q j ) 2j = e ( + e ) 2 + e + s = e e 4s sinh() = > 0, s (81) which contradicts relation (80). This contradiction proves that the elements of the set V j are linearly independent, j = 1, 2, 3,..., m 2, for 0 < s < 1. Let us prove that the elements of the set V m, are linearly independent. Assume that there exist constants α k, k = 1, 2,..., m, such that H n q m = m 2 k= α k+2 H n q k. (82) Using (64)-(66), the (n 3)-th equation of (82) can be written as follows: (H n q m ) n 3 =e 1 4s D m = m 3 k= α k+2 (H n q k ) n 3 m 3 =α 1 e 2+4s + α 2 e 4s + e +4s k=1 α k+2 C k + α m (H n q m 2 ) n 3. (83) 14
15 Similarly we obtain the (n 2)-th, (n 1)-th n-th equations, corresponding to vector equation (82): m 3 (H n q m ) n 2 = α 1 e 2+3s + α 2 e 3s + e +3s k=1 α k+2 C k + α m (H n q m 2 ) n 2, (84) m 2 (H n q m ) n = α 1 e 2+ + α 2 e + e + k=1 α k+2 C k (85) m 2 (H n q m ) n = α 1 e 2+s + α 2 e s + e +s k=1 respectively. Multiply (86) by e s compare with (85) to get α k+2 C k, (86) α 2 = (H nq m ) n e s (H n q m ) n e 1 = 1 e + 4se e 3s (1 e, (87) ) where formula (66) was used. Multiplying (86) by e 3s, comparing with equation (83), using (87), we obtain α m = e1 4s D m e 3s (H n q m ) n α 2 (e 4s e ) (H n q m 2 ) n 3 e +4s C m 2 = 2 + 4s 2es s + 4e s 2e 3s s + e 4s ( 2 + 4s) + e 4s 4e. s (88) Another expression for α m is obtained by multiplying (86) by e comparing with (84): α m = (H nq m ) n 2 e (H n q m ) n α 2 (e 3s e s ) (H n q m 2 ) n 2 e +3s C m 2 = 2 + 4s 4es s + e ( 2 + 4s) + e 2e s, s (89) where α 2 is given in (87). In deriving formulas (88) (89) we have used the relation m = 1 s + 1 equation (66). Let us prove that equations (88) (89) lead to a contradiction. Define r 1 := 2 + 4s 2e s s + 4e s 2e 3s s + e 4s ( 2 + 4s), r 2 := + e 4s 4e s, r 3 := 2 + 4s 4e s s + e ( 2 + 4s), r 4 := + e 2e s s. (90) 15
16 Then from (88) (89) we get We have r 3 r 2 r 1 r 4 = 0. (91) r 3 r 2 r 1 r 4 = 2e s ( + e s ) 2 s(3 + 4s + (4s 3)e e s ) > 0 for s (0, 1). (92) The sign of the right side of equality (92) is the same as the sign of 3+4s+(4s 3)e e s := β(s). Let us check that β(s) > 0 for s (0, 1). One has β(0) = 0, β (0) = 0, β (s) = 4 2e +8se 2e s e s, β = 4e +16se 4e s e s > 0. If β (s) > 0 for s (0, 1) β(0) = 0, β (0) = 0, then β(s) > 0 for s (0, 1). Inequality (92) contradicts relation (91) which proves that H n q j, j =, 0, 1, 2,..., m 1, are linearly independent. Similarly, to prove that H n q j, j =, 0, 1, 2,..., m, are linearly independent, we assume that there exist constants α k, k = 1, 2,..., m + 1, such that H n q m = m k= α k+2 H n q k. (93) Using formulas (64)-(66), one can write the (n 5)-th equation: (H n q m ) n 5 = e 1 6s D m = m k= α k+2 (H n q k ) n 5 m 4 = α 1 e 2+6s + α 2 e 6s + e +6s α j+2 C j + α m (H n q m 3 ) n 5 + α m e 1 6s D m 2 + α m+1 e 1 6s D m, (94) Similarly one obtains the (n 4)-th, (n 3)-th, (n 2)-th, (n 1)-th n-th equations corresponding to the vector equation (93): m 4 (H n q m ) n 4 = e 1 5s D m = α 1 e 2+5s + α 2 e 5s + e +5s α j+2 C j + α m (H n q m 3 ) n 4 + α m (H n q m 2 ) n 4 + α m+1 e 1 5s D m, (95) m 3 (H n q m ) n 3 = e 1 4s D m = α 1 e 2+4s + α 2 e 4s + e +4s + α m (H n q m 2 ) n 3 + α m+1 e 1 4s D m, m 3 (H n q m ) n 2 = e 1 3s D m = α 1 e 2+3s + α 2 e 3s + e +3s + α m (H n q m 2 ) n 2 + α m+1 (H n q m ) n 2, α j+2 C j α j+2 C j (96) (97) 16
17 m 2 (H n q m ) n = e 1 D m = α 1 e 2+ + α 2 e + e + α j+2 C j (98) + α m+1 (H n q m ) n m 2 (H n q m ) n = α 1 e 2+s + α 2 e s + e +s α j+2 C j + α m+1 (H n q m ) n, (99) respectively. Here we have used the assumption n 6. From (99) one gets m 2 α 1 = (H n q m ) n e 2 s α 2 e 2 e α k+2 C k α m+1 (H n q m ) n e 2 s. (100) k=1 If one substitutes (100) into equations (98), (97) (96), then one obtains the following relations: α 2 = p 1 p 2 α m+1, (101) respectively, where p 1 := e1 D m (H n q m ) n e s e 1 α 2 = p 3 p 4 α m p 5 α m+1 (102) α 2 = p 6 p 7 α m p 8 α m+1, (103), p 2 := (H nq m ) n (H n q m ) n e s e, 1 p 3 := e1 3s D m (H n q m ) n e e 3s e s, p 4 := (H nq m 2 ) n 2 e +3s C m 2 e 3s e s, p 5 := (H nq m ) n 2 (H n q m ) n e e 3s e s, p 6 := e1 4s D m (H n q m ) n e 3s e 4s e, p 7 := (H nq m 2 ) n 3 e +4s C m 2 e 4s e, p 8 := e1 4s D m (H n q m ) n e 3s e 4s e. Another formula for α 1 one gets from equation (96): (104) m 3 α 1 = e 3 8s D m α 2 e 2 8s e α j+2 C j α m (H n q m 2 ) n 3 e 2 4s α m+1 e 3 8s D m. (105) Substituting (105) into equations (95) (94), yields α 2 = p 9 p 10 α m p 11 α m p 12 α m+1 (106) 17
18 respectively, where p 9 := ed m, α 2 = p 9 p 13 α m p 14 α m p 12 α m+1, (107) p 10 := (H nq m 3 ) n 4 e +5s C m 3 e 5s e 3s, p 11 := (H nq m 2 ) n 4 e s (H n q m 2 ) n 3 e 5s e 3s, p 12 := ed m, p 13 := (H nq m 3 ) n 5 e +6s C m 3 e 6s e, p 14 := e1 6s D m 2 e (H n q m 2 ) n 3 e 6s e. (108) Let us prove that the equations (101) (106) lead to a contradiction. From equations (101) (102) we obtain This together with (103) (101) yield α m = p 3 p 1 p 4 + p 2 p 5 p 4 α m+1. (109) α m+1 = p 3 p 1 p 4 p 2 p 8 p 7 Equations (106), (107) (109) yield α m = p 14 p 11 p 10 p 13 α m = p 14 p 11 p 10 p 13 p6 p1 p 7 p2 p5 p 4. (110) ( p3 p 1 p 4 + p 2 p 5 p 4 α m+1 ). (111) This together with (106) imply ( ) ( p 14 p 11 p3 p 1 α 2 = p 9 p 10 + p 11 + p ) 2 p 5 α m+1 p 12 α m+1, (112) p 10 p 13 p 4 p 4 where α m+1 is given in (110). Let ( ) ( p 14 p 11 p3 p 1 L 1 := p 9 p 10 + p 11 + p ) 2 p 5 α m+1 p 12 α m+1 (113) p 10 p 13 p 4 p 4 Then, from (101) (106) one gets L 2 := p 1 p 2 α m+1. (114) L 1 L 2 = 0. (115) Applying formulas (64)-(66) in (113) (114) using the relation e s =, we obtain j=0 sj j! L 1 L 2 = 2e4s [1 + e ( + s) + s]( + e 2e s s) ( + e )s[3 + 4s 2e s s + e ( 3 + 4s)] [ ] [ ] 2e 4s 2 j ( j ) j=3 (j)! s j 2( 2j j+1 ) j=2 j! s j+1 = ( + e )s[3 + 4s 2e s s + e > 0, ( 3 + 4s)] 18 (116)
19 because e > 1 for all 0 < s < 1, 2 j ( j j ) > 0 for all j 3, 2( j+1 1) > 0 for all j 2 3+4s 2e s s+e ( 3+4s) > 0 which was proved in (92). Inequality (116) contradicts relation (115) which proves that H n q j, j =, 0, 1, 2,..., m, are linearly independent. Theorem 3.1 is proved. 4 Numerical experiments Note that for all w, u, v R n we have n w l u l v l = v t W u, (117) l=1 where t sts for transpose w w W := w n 0 (118) w n Then DP := H n (f Rh m ) 2 w (n),1 n = w (n) l [(f(x l ) Rh m (x l )) 2 + (f (x l ) (Rh m ) (x l )) 2 ] l=1 = [H n (f Rh m )] t W H n (f Rh m ) + [H n (f (Rh m ) )] t W H n (f (Rh m ) ), (119) where W is defined in (118) with w j = w (n) j, j = 1, 2,..., n, defined in (27). The vectors H n Rh m H n (Rh m ) are computed as follows. Using (64)-(66), the vector H n Rh m can be represented by where = entries: 0 1. m H n Rh m = S m, (120) S m is an n (m + 2) matrix with the following (S m ) i,1 = (H n q ) i, (S m ) i,2 = (H n q 0 ) i, i = 1, 2,..., n, (S m ) i,j = (H n q j 2 ) i, i = 1, 2,..., n, j = 3, 4,..., m + 2. (121) 19
20 The vector H n (Rh m ) is computed as follows. Let J i,j := e xi x i e y ϕ j (y)dy e xi xi This together with (60) yield e e(+e ), i = 1, e J i,1 = s (1 e s +s) s, i = 2, e x i (+e ) 2es, i 3, e xi e1 2js (+e ) 2, i 2j 3; 2+e 3s 3e s, i = 2j 2; J i,j = 0, i = 2j 1; 2+e 3s 3e s, i = 2j; e xi e+2( 2+j)s (+e ) 2, i 2j + 1, { e x i (+e ) J i,m = 2es, i n 1, e s (+e s s) s, i = n. e y ϕ j (y)dy, i = 1, 2,..., n, j = 1, 2,..., m. (122) 1 i n, 1 < j < m, Then, using (123), the vector H n (Rh m ) can be rewritten as follows: where = entries: 0 1. m (123) H n (Rh m ) = T m, (124) T m is an n (m + 2) matrix with the following (T m ) i,1 = e (1+xi), (T m ) i,2 = e (1 xi), i = 1, 2,..., n, (T m ) i,j = J i,j 2, i = 1, 2,..., n, j = 3, 4,..., m + 2. (125) We consider the following examples discussed in [7]: (1) f(x) = cos(π(x + 1)) with the exact solution h(x) = + (1 + π 2 ) cos(π(x + 1)). (2) f(x) = 2e x + 2 π sin(π(x + 1)) + 2 cos(π(x + 1)) with the exact solution h(x) = 1 π sin(π(x + 1)) + (1 + π2 ) cos(π(x + 1)). (3) f(x) = cos( π(x+1) 2 ) + 4 cos(2π(x + 1)) 1.5 cos( 7π(x+1) 2 ) with the exact solution h(x) = 1 π2 2 (1 + 4 ) cos( π(x+1) 2 ) + (2 + 8π 2 ) cos(2π(x + 1)) 0.75( π 2 ) cos( 7π(x+1) 2 ) δ(x + 1) δ(x 1). 20
21 (4) f(x) = e x + 2 sin(2π(x + 1)) with the exact solution h(x) = (1 + 4π 2 ) sin(2π(x + 1)) + (e 2π)δ(x + 1) + 2πδ(x 1). In all the above examples we have f C 2 ([, 1]). Therefore, one may use the basis functions given in (60). In each example we compute the relative pointwise errors: RP E(t i ) := g m(t i ) g(t i ) max 1 i M g(t i ), (126) where g(x) g m (x) are defined in (9) (11), respectively, 2 t i := + (i 1), i = 1, 2,..., M. (127) M 1 The algorithm can be written as follows. Step 0. Set k = 3, n = 2k, m = n 2 defined in (119). + 1, ɛ (0, 1) DP 10, where DP is Step 1. Construct the weights w (n) j, j = 1, 2,..., n, defined in (27). Step 2. Construct the matrix A m+2 the vector F m+2 which are defined in (36). 0 Step 3. Solve for c := the linear algebraic system A m+2 c = F m m Step 4. Compute DP = H n (f Rh m ) 2 w (n),1. (128) Step 5. If DP > ɛ then set k = k + 1, n = 2k m = n 2 + 1, go to Step 1. Otherwise, stop the iteration use h m (x) = m j= c(m) j ϕ j (x) as the approximate solution, where ϕ (x) := δ(1 + x), ϕ 0 (x) := δ(x 1) ϕ j (x), j = 1, 2,..., m, are defined in (60) j, j =, 0, 1,..., m, are obtained in Step 3. In all the experiments the following parameters are used: M = 200 ɛ = 10 4, 10 6, We also compute the relative error RE := max 1 i M RP E(t i), (129) where RP E is defined in (120). Let us discuss the results of our experiments. Example 1. In this example the coefficients a a 0, given in (8) (9), respectively, are zeros. Our experiments show, see Table 1, that the approximate 21
22 Table 1: Example 1 n m ɛ a a n m ɛ DP RE coefficients c(m) 0 converge to a a 0, respectively, as ɛ 0. Here to get DP 10 6, we need n = 32 collocation points distributed uniformly in the interval [, 1]. Moreover, the matrix A m+2 is of the size 19 by 19 which is small. For ɛ = 10 6 the relative error RE is of order The RP E at the points t j are distributed in the interval [0, 0.018) as shown in Figure 2. In computing the approximate solution h m at the points t i, i = 1, 2,..., M, one needs at most two out of m = 17 basis functions ϕ j (x). The reconstruction of the continuous part of the exact solution can be seen in Figure 1. One can see from this Figure that for ɛ = 10 6 the continuous part g(x) of the exact solution h(x) can be recovered very well by the approximate function g m (x) at the points t j, j = 1, 2,..., M. Figure 2: A reconstruction of the continuous part g(x) (above) of Example 1 with ɛ = 10 6, the corresponding Relative Pointwise Errors (RPE) (below) 22
23 Table 2: Example 2 n m ɛ a a n m ɛ DP RE Example 2. This example is a modification of Example 1, where the constant 2 is replaced with the function 2e x + 2 π sin(π(x + 1)). In this case the coefficients a a 0 are also zeros. The results can be seen in Table 2. As in Example 1, both approximate coefficients c(m) 0 converge to 0 as ɛ 0. The number of collocation points at each case is equal to the number of collocation points obtained in Example 1. Also the RP E at each observed point is in the interval [0, 0.02). One can see from Figure 3 that the continuous part g(x) of the exact solution h(x) can be well approximated by the approximate function g m (x) with ɛ = 10 6 RE = O(10 2 ). Figure 3: A reconstruction of the continuous part g(x) (above) of Example 2 with ɛ = 10 6, the corresponding Relative Pointwise Errors (RPE) (below) Example 3. In this example the coefficients of the distributional parts a 23
24 Table 3: Example 3 n m ɛ a a n m ɛ DP RE a 0 are not zeros. The function f is oscillating more than the functions f given in Examples 1 2, the number of collocation points is larger than in the previous two examples, as shown in Table 3. In this table one can see that the approximate coefficients c(m) 0 converge to a a 0, respectively. The continuous part of the exact solution can be approximated by the approximate function g m (x) very well with ɛ = 10 6 RE = O(10 3 ) as shown in Figure 4. In the same Figure one can see that the RP E at each observed point is distributed in the interval [0, ). Figure 4: A reconstruction of the continuous part g(x) (above) of Example 3 with ɛ = 10 6, the corresponding Relative Pointwise Errors (RPE) (below) Example 4. Here we give another example of the exact solution h having non-zero coefficients a a 0. In this example the function f is oscillating less than the f in Example 3, but more than the f in examples 1 2. As 24
25 shown in Table 4 the number of collocation points n is smaller than the the number of collocation points given in Example 3. In this example the exact coefficients a a 0 are obtained at the error level ɛ = 10 8 which is shown in Table 4. Figure 5 shows that at the level ɛ = 10 6 we have obtained a good approximation of the continuous part g(x) of the exact solution h(x). Here the relative error RE is of order O(10 2 ). Table 4: Example 4 n m ɛ a a n m ɛ DP RE Figure 5: A reconstruction of the continuous part g(x) (above) of Example 4 with ɛ = 10 6, the corresponding Relative Pointwise Errors (RPE) (below) 25
26 References [1] P.J. Davis P. Rabinowitz, Methods of numerical integration, Academic Press, London, [2] L.Kantorovich G.Akilov, Functional Analysis, Pergamon Press, New York, [3] S. Mikhlin, S. Prössdorf, Singular integral operators, Springer-Verlag, Berlin, [4] A.G. Ramm, Theory Applications of Some New Classes of Integral Equations, Springer-Verlag, New York, [5] A.G. Ramm, Rom Fields Estimation, World Sci. Publishers, Singapore, [6] A.G. Ramm, Numerical solution of integral equations in a space of distributions, J. Math. Anal. Appl., 110, (1985), [7] A.G. Ramm Peiqing Li, Numerical solution of some integral equations in distributions, Computers Math. Applic., 22, (1991), [8] A.G. Ramm, Collocation method for solving some integral equations of estimation theory, Internat. Journ. of Pure Appl. Math., 62, N1, (2010). [9] E. Rozema, Estimating the error in the trapezoidal rule, The American Mathematical Monthly, 87, 2, (1980), [10] L.L. Schumaker, Spline Functions: Basic Theory, Cambridge University Press, Cambridge,
Electromagnetic wave scattering by many small bodies and creating materials with a desired refraction coefficient
168 Int. J. Computing Science and Mathematics, Vol. Electromagnetic wave scattering by many small bodies and creating materials with a desired refraction coefficient A.G. Ramm Mathematics epartment, Kansas
More informationInverse scattering problem with underdetermined data.
Math. Methods in Natur. Phenom. (MMNP), 9, N5, (2014), 244-253. Inverse scattering problem with underdetermined data. A. G. Ramm Mathematics epartment, Kansas State University, Manhattan, KS 66506-2602,
More informationComm. Nonlin. Sci. and Numer. Simul., 12, (2007),
Comm. Nonlin. Sci. and Numer. Simul., 12, (2007), 1390-1394. 1 A Schrödinger singular perturbation problem A.G. Ramm Mathematics Department, Kansas State University, Manhattan, KS 66506-2602, USA ramm@math.ksu.edu
More informationCreating materials with a desired refraction coefficient: numerical experiments
Creating materials with a desired refraction coefficient: numerical experiments Sapto W. Indratno and Alexander G. Ramm Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA
More informationA G Ramm, Implicit Function Theorem via the DSM, Nonlinear Analysis: Theory, Methods and Appl., 72, N3-4, (2010),
A G Ramm, Implicit Function Theorem via the DSM, Nonlinear Analysis: Theory, Methods and Appl., 72, N3-4, (21), 1916-1921. 1 Implicit Function Theorem via the DSM A G Ramm Department of Mathematics Kansas
More informationHeat Transfer in a Medium in Which Many Small Particles Are Embedded
Math. Model. Nat. Phenom. Vol. 8, No., 23, pp. 93 99 DOI:.5/mmnp/2384 Heat Transfer in a Medium in Which Many Small Particles Are Embedded A. G. Ramm Department of Mathematics Kansas State University,
More informationHow large is the class of operator equations solvable by a DSM Newton-type method?
This is the author s final, peer-reviewed manuscript as accepted for publication. The publisher-formatted version may be available through the publisher s web site or your institution s library. How large
More informationA method for creating materials with a desired refraction coefficient
This is the author s final, peer-reviewed manuscript as accepted for publication. The publisher-formatted version may be available through the publisher s web site or your institution s library. A method
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationA numerical algorithm for solving 3D inverse scattering problem with non-over-determined data
A numerical algorithm for solving 3D inverse scattering problem with non-over-determined data Alexander Ramm, Cong Van Department of Mathematics, Kansas State University, Manhattan, KS 66506, USA ramm@math.ksu.edu;
More informationQ. Zou and A.G. Ramm. Department of Mathematics, Kansas State University, Manhattan, KS 66506, USA.
Computers and Math. with Applic., 21 (1991), 75-80 NUMERICAL SOLUTION OF SOME INVERSE SCATTERING PROBLEMS OF GEOPHYSICS Q. Zou and A.G. Ramm Department of Mathematics, Kansas State University, Manhattan,
More informationPhys.Let A. 360, N1, (2006),
Phys.Let A. 360, N1, (006), -5. 1 Completeness of the set of scattering amplitudes A.G. Ramm Mathematics epartment, Kansas State University, Manhattan, KS 66506-60, USA ramm@math.ksu.edu Abstract Let f
More informationThus f is continuous at x 0. Matthew Straughn Math 402 Homework 6
Matthew Straughn Math 402 Homework 6 Homework 6 (p. 452) 14.3.3, 14.3.4, 14.3.5, 14.3.8 (p. 455) 14.4.3* (p. 458) 14.5.3 (p. 460) 14.6.1 (p. 472) 14.7.2* Lemma 1. If (f (n) ) converges uniformly to some
More informationON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX- ON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS N. S. HOANG AND A. G. RAMM (Communicated
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationOutline of Fourier Series: Math 201B
Outline of Fourier Series: Math 201B February 24, 2011 1 Functions and convolutions 1.1 Periodic functions Periodic functions. Let = R/(2πZ) denote the circle, or onedimensional torus. A function f : C
More informationFinite-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 (, )
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationFebruary 12, :52 WSPC/Book Trim Size for 9in x 6in book. Publishers page
Publishers page Publishers page Publishers page Publishers page To the memory of my parents vi Preface This book presents analytic theory of random fields estimation optimal by the criterion of minimum
More information7: FOURIER SERIES STEVEN HEILMAN
7: FOURIER SERIES STEVE HEILMA Contents 1. Review 1 2. Introduction 1 3. Periodic Functions 2 4. Inner Products on Periodic Functions 3 5. Trigonometric Polynomials 5 6. Periodic Convolutions 7 7. Fourier
More informationAn 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
More informationFinding discontinuities of piecewise-smooth functions
Finding discontinuities of piecewise-smooth functions A.G. Ramm Mathematics Department, Kansas State University, Manhattan, KS 66506-2602, USA ramm@math.ksu.edu Abstract Formulas for stable differentiation
More informationOrthonormal series expansion and finite spherical Hankel transform of generalized functions
Malaya Journal of Matematik 2(1)(2013) 77 82 Orthonormal series expansion and finite spherical Hankel transform of generalized functions S.K. Panchal a, a Department of Mathematics, Dr. Babasaheb Ambedkar
More informationCHAPTER 10 Shape Preserving Properties of B-splines
CHAPTER 10 Shape Preserving Properties of B-splines In earlier chapters we have seen a number of examples of the close relationship between a spline function and its B-spline coefficients This is especially
More informationA CHARACTERIZATION OF STRICT LOCAL MINIMIZERS OF ORDER ONE FOR STATIC MINMAX PROBLEMS IN THE PARAMETRIC CONSTRAINT CASE
Journal of Applied Analysis Vol. 6, No. 1 (2000), pp. 139 148 A CHARACTERIZATION OF STRICT LOCAL MINIMIZERS OF ORDER ONE FOR STATIC MINMAX PROBLEMS IN THE PARAMETRIC CONSTRAINT CASE A. W. A. TAHA Received
More informationMath The Laplacian. 1 Green s Identities, Fundamental Solution
Math. 209 The Laplacian Green s Identities, Fundamental Solution Let be a bounded open set in R n, n 2, with smooth boundary. The fact that the boundary is smooth means that at each point x the external
More informationGENERAL 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
More informationPolynomial Approximation: The Fourier System
Polynomial Approximation: The Fourier System Charles B. I. Chilaka CASA Seminar 17th October, 2007 Outline 1 Introduction and problem formulation 2 The continuous Fourier expansion 3 The discrete Fourier
More informationSpectra of Multiplication Operators as a Numerical Tool. B. Vioreanu and V. Rokhlin Technical Report YALEU/DCS/TR-1443 March 3, 2011
We introduce a numerical procedure for the construction of interpolation and quadrature formulae on bounded convex regions in the plane. The construction is based on the behavior of spectra of certain
More informationDynamical Systems Method for Solving Ill-conditioned Linear Algebraic Systems
Dynamical Systems Method for Solving Ill-conditioned Linear Algebraic Systems Sapto W. Indratno Department of Mathematics Kansas State University, Manhattan, KS 6656-6, USA sapto@math.ksu.edu A G Ramm
More informationStarting from Heat Equation
Department of Applied Mathematics National Chiao Tung University Hsin-Chu 30010, TAIWAN 20th August 2009 Analytical Theory of Heat The differential equations of the propagation of heat express the most
More informationSolving Analytically Singular Sixth-Order Boundary Value Problems
Int. Journal of Math. Analysis, Vol. 3, 009, no. 39, 1945-195 Solving Analytically Singular Sixth-Order Boundary Value Problems Fazhan Geng 1 Department of Mathematics, Changshu Institute of Technology
More informationNonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm
More information1 Continuity Classes C m (Ω)
0.1 Norms 0.1 Norms A norm on a linear space X is a function : X R with the properties: Positive Definite: x 0 x X (nonnegative) x = 0 x = 0 (strictly positive) λx = λ x x X, λ C(homogeneous) x + y x +
More informationCreating materials with a desired refraction coefficient: numerical experiments. Sapto W. Indratno and Alexander G. Ramm*
76 Int. J. Computing Science and Mathematics, Vol. 3, Nos. /2, 200 Creating materials with a desired refraction coefficient: numerical experiments Sapto W. Indratno and Alexander G. Ramm* Department of
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More informationwhere c R and the content of f is one. 1
9. Gauss Lemma Obviously it would be nice to have some more general methods of proving that a given polynomial is irreducible. The first is rather beautiful and due to Gauss. The basic idea is as follows.
More informationMATH34032 Mid-term Test 10.00am 10.50am, 26th March 2010 Answer all six question [20% of the total mark for this course]
MATH3432: Green s Functions, Integral Equations and the Calculus of Variations 1 MATH3432 Mid-term Test 1.am 1.5am, 26th March 21 Answer all six question [2% of the total mark for this course] Qu.1 (a)
More informationProblem 1A. Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1. (Hint: 1. Solution: The volume is 1. Problem 2A.
Problem 1A Find the volume of the solid given by x 2 + z 2 1, y 2 + z 2 1 (Hint: 1 1 (something)dz) Solution: The volume is 1 1 4xydz where x = y = 1 z 2 This integral has value 16/3 Problem 2A Let f(x)
More informationTitle: Localized self-adjointness of Schrödinger-type operators on Riemannian manifolds. Proposed running head: Schrödinger-type operators on
Title: Localized self-adjointness of Schrödinger-type operators on Riemannian manifolds. Proposed running head: Schrödinger-type operators on manifolds. Author: Ognjen Milatovic Department Address: Department
More informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More informationTEST CODE: PMB SYLLABUS
TEST CODE: PMB SYLLABUS Convergence and divergence of sequence and series; Cauchy sequence and completeness; Bolzano-Weierstrass theorem; continuity, uniform continuity, differentiability; directional
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationCourse 212: Academic Year Section 1: Metric Spaces
Course 212: Academic Year 1991-2 Section 1: Metric Spaces D. R. Wilkins Contents 1 Metric Spaces 3 1.1 Distance Functions and Metric Spaces............. 3 1.2 Convergence and Continuity in Metric Spaces.........
More informationLinear 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
More informationSOLUTIONS OF SELECTED PROBLEMS
SOLUTIONS OF SELECTED PROBLEMS Problem 36, p. 63 If µ(e n < and χ En f in L, then f is a.e. equal to a characteristic function of a measurable set. Solution: By Corollary.3, there esists a subsequence
More informationEncyclopedia of Mathematics, Supplemental Vol. 3, Kluwer Acad. Publishers, Dordrecht,
Encyclopedia of Mathematics, Supplemental Vol. 3, Kluwer Acad. Publishers, Dordrecht, 2001, 328-329 1 Reproducing kernel Consider an abstract set E and a linear set F of functions f : E C. Assume that
More informationMathematics Department Stanford University Math 61CM/DM Inner products
Mathematics Department Stanford University Math 61CM/DM Inner products Recall the definition of an inner product space; see Appendix A.8 of the textbook. Definition 1 An inner product space V is a vector
More informationbe any ring homomorphism and let s S be any element of S. Then there is a unique ring homomorphism
21. Polynomial rings Let us now turn out attention to determining the prime elements of a polynomial ring, where the coefficient ring is a field. We already know that such a polynomial ring is a UFD. Therefore
More informationOverview of normed linear spaces
20 Chapter 2 Overview of normed linear spaces Starting from this chapter, we begin examining linear spaces with at least one extra structure (topology or geometry). We assume linearity; this is a natural
More informationA Simple Proof of the Fredholm Alternative and a Characterization of the Fredholm Operators
thus a n+1 = (2n + 1)a n /2(n + 1). We know that a 0 = π, and the remaining part follows by induction. Thus g(x, y) dx dy = 1 2 tanh 2n v cosh v dv Equations (4) and (5) give the desired result. Remarks.
More informationNOTES ON SCHAUDER ESTIMATES. r 2 x y 2
NOTES ON SCHAUDER ESTIMATES CRISTIAN E GUTIÉRREZ JULY 26, 2005 Lemma 1 If u f in B r y), then ux) u + r2 x y 2 B r y) B r y) f, x B r y) Proof Let gx) = ux) Br y) u r2 x y 2 Br y) f We have g = u + Br
More informationMath 413/513 Chapter 6 (from Friedberg, Insel, & Spence)
Math 413/513 Chapter 6 (from Friedberg, Insel, & Spence) David Glickenstein December 7, 2015 1 Inner product spaces In this chapter, we will only consider the elds R and C. De nition 1 Let V be a vector
More informationIr O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v )
Section 3.2 Theorem 3.6. Let A be an m n matrix of rank r. Then r m, r n, and, by means of a finite number of elementary row and column operations, A can be transformed into the matrix ( ) Ir O D = 1 O
More informationbe the set of complex valued 2π-periodic functions f on R such that
. Fourier series. Definition.. Given a real number P, we say a complex valued function f on R is P -periodic if f(x + P ) f(x) for all x R. We let be the set of complex valued -periodic functions f on
More informationMTH 503: Functional Analysis
MTH 53: Functional Analysis Semester 1, 215-216 Dr. Prahlad Vaidyanathan Contents I. Normed Linear Spaces 4 1. Review of Linear Algebra........................... 4 2. Definition and Examples...........................
More informationHow to create materials with a desired refraction coefficient?
How to create materials with a desired refraction coefficient? A. G. Ramm Mathematics Department, Kansas State University, Manhattan, KS66506, USA ramm@math.ksu.edu www.math.ksu.edu/ ramm Abstract The
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
More informationWEYL S LEMMA, ONE OF MANY. Daniel W. Stroock
WEYL S LEMMA, ONE OF MANY Daniel W Stroock Abstract This note is a brief, and somewhat biased, account of the evolution of what people working in PDE s call Weyl s Lemma about the regularity of solutions
More informationMATH 101A: ALGEBRA I PART C: TENSOR PRODUCT AND MULTILINEAR ALGEBRA. This is the title page for the notes on tensor products and multilinear algebra.
MATH 101A: ALGEBRA I PART C: TENSOR PRODUCT AND MULTILINEAR ALGEBRA This is the title page for the notes on tensor products and multilinear algebra. Contents 1. Bilinear forms and quadratic forms 1 1.1.
More informationThe Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods
The Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods by Hae-Soo Oh Department of Mathematics, University of North Carolina at Charlotte, Charlotte, NC 28223 June
More informationTHE INVERSE FUNCTION THEOREM
THE INVERSE FUNCTION THEOREM W. PATRICK HOOPER The implicit function theorem is the following result: Theorem 1. Let f be a C 1 function from a neighborhood of a point a R n into R n. Suppose A = Df(a)
More informationIterative Solution of a Matrix Riccati Equation Arising in Stochastic Control
Iterative Solution of a Matrix Riccati Equation Arising in Stochastic Control Chun-Hua Guo Dedicated to Peter Lancaster on the occasion of his 70th birthday We consider iterative methods for finding the
More informationLecture 4: Fourier Transforms.
1 Definition. Lecture 4: Fourier Transforms. We now come to Fourier transforms, which we give in the form of a definition. First we define the spaces L 1 () and L 2 (). Definition 1.1 The space L 1 ()
More information(3) Let Y be a totally bounded subset of a metric space X. Then the closure Y of Y
() Consider A = { q Q : q 2 2} as a subset of the metric space (Q, d), where d(x, y) = x y. Then A is A) closed but not open in Q B) open but not closed in Q C) neither open nor closed in Q D) both open
More information4 Riesz Kernels. Since the functions k i (ξ) = ξ i. are bounded functions it is clear that R
4 Riesz Kernels. A natural generalization of the Hilbert transform to higher dimension is mutiplication of the Fourier Transform by homogeneous functions of degree 0, the simplest ones being R i f(ξ) =
More informationAN INVERSE PROBLEM FOR THE WAVE EQUATION WITH A TIME DEPENDENT COEFFICIENT
AN INVERSE PROBLEM FOR THE WAVE EQUATION WITH A TIME DEPENDENT COEFFICIENT Rakesh Department of Mathematics University of Delaware Newark, DE 19716 A.G.Ramm Department of Mathematics Kansas State University
More informationKey Words: cardinal B-spline, coefficients, moments, rectangular rule, interpolating quadratic spline, hat function, cubic B-spline.
M a t h e m a t i c a B a l a n i c a New Series Vol. 24, 2, Fasc.3-4 Splines in Numerical Integration Zlato Udovičić We gave a short review of several results which are related to the role of splines
More informationExistence and Multiplicity of Solutions for a Class of Semilinear Elliptic Equations 1
Journal of Mathematical Analysis and Applications 257, 321 331 (2001) doi:10.1006/jmaa.2000.7347, available online at http://www.idealibrary.com on Existence and Multiplicity of Solutions for a Class of
More informationImmerse Metric Space Homework
Immerse Metric Space Homework (Exercises -2). In R n, define d(x, y) = x y +... + x n y n. Show that d is a metric that induces the usual topology. Sketch the basis elements when n = 2. Solution: Steps
More informationReal 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.
More informationProblem Set 5: Solutions Math 201A: Fall 2016
Problem Set 5: s Math 21A: Fall 216 Problem 1. Define f : [1, ) [1, ) by f(x) = x + 1/x. Show that f(x) f(y) < x y for all x, y [1, ) with x y, but f has no fixed point. Why doesn t this example contradict
More informationARCS IN FINITE PROJECTIVE SPACES. Basic objects and definitions
ARCS IN FINITE PROJECTIVE SPACES SIMEON BALL Abstract. These notes are an outline of a course on arcs given at the Finite Geometry Summer School, University of Sussex, June 26-30, 2017. Let K denote an
More information1. Introduction Boundary estimates for the second derivatives of the solution to the Dirichlet problem for the Monge-Ampere equation
POINTWISE C 2,α ESTIMATES AT THE BOUNDARY FOR THE MONGE-AMPERE EQUATION O. SAVIN Abstract. We prove a localization property of boundary sections for solutions to the Monge-Ampere equation. As a consequence
More informationSOLUTION OF POISSON S EQUATION. Contents
SOLUTION OF POISSON S EQUATION CRISTIAN E. GUTIÉRREZ OCTOBER 5, 2013 Contents 1. Differentiation under the integral sign 1 2. The Newtonian potential is C 1 2 3. The Newtonian potential from the 3rd Green
More informationLecture 7.3: Ring homomorphisms
Lecture 7.3: Ring homomorphisms Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4120, Modern Algebra M. Macauley (Clemson) Lecture 7.3:
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationCHAPTER VIII HILBERT SPACES
CHAPTER VIII HILBERT SPACES DEFINITION Let X and Y be two complex vector spaces. A map T : X Y is called a conjugate-linear transformation if it is a reallinear transformation from X into Y, and if T (λx)
More informationFunctional Analysis HW #3
Functional Analysis HW #3 Sangchul Lee October 26, 2015 1 Solutions Exercise 2.1. Let D = { f C([0, 1]) : f C([0, 1])} and define f d = f + f. Show that D is a Banach algebra and that the Gelfand transform
More informationFourier Series. 1. Review of Linear Algebra
Fourier Series In this section we give a short introduction to Fourier Analysis. If you are interested in Fourier analysis and would like to know more detail, I highly recommend the following book: Fourier
More informationINVERSE FUNCTION THEOREM and SURFACES IN R n
INVERSE FUNCTION THEOREM and SURFACES IN R n Let f C k (U; R n ), with U R n open. Assume df(a) GL(R n ), where a U. The Inverse Function Theorem says there is an open neighborhood V U of a in R n so that
More informationAlgebra 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
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More informationWhere is matrix multiplication locally open?
Linear Algebra and its Applications 517 (2017) 167 176 Contents lists available at ScienceDirect Linear Algebra and its Applications www.elsevier.com/locate/laa Where is matrix multiplication locally open?
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationMetric Spaces. Exercises Fall 2017 Lecturer: Viveka Erlandsson. Written by M.van den Berg
Metric Spaces Exercises Fall 2017 Lecturer: Viveka Erlandsson Written by M.van den Berg School of Mathematics University of Bristol BS8 1TW Bristol, UK 1 Exercises. 1. Let X be a non-empty set, and suppose
More informationContinuous Functions on Metric Spaces
Continuous Functions on Metric Spaces Math 201A, Fall 2016 1 Continuous functions Definition 1. Let (X, d X ) and (Y, d Y ) be metric spaces. A function f : X Y is continuous at a X if for every ɛ > 0
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationThere are two things that are particularly nice about the first basis
Orthogonality and the Gram-Schmidt Process In Chapter 4, we spent a great deal of time studying the problem of finding a basis for a vector space We know that a basis for a vector space can potentially
More information1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer 11(2) (1989),
Real Analysis 2, Math 651, Spring 2005 April 26, 2005 1 Real Analysis 2, Math 651, Spring 2005 Krzysztof Chris Ciesielski 1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer
More informationarxiv:math/ v1 [math.dg] 6 Jul 1998
arxiv:math/9807024v [math.dg] 6 Jul 998 The Fundamental Theorem of Geometric Calculus via a Generalized Riemann Integral Alan Macdonald Department of Mathematics Luther College, Decorah, IA 520, U.S.A.
More informationQUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1)
QUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday September 21, 2004 (Day 1) Each of the six questions is worth 10 points. 1) Let H be a (real or complex) Hilbert space. We say
More informationEXTENSIONS OF ABSOLUTE VALUES
CHAPTER III EXTENSIONS OF ABSOLUTE VALUES 1. Norm and Trace Let k be a field and E a vector space of dimension N over k. We write End k (E) for the ring of k linear endomorphisms of E and Aut k (E) = End
More informationMath 209B Homework 2
Math 29B Homework 2 Edward Burkard Note: All vector spaces are over the field F = R or C 4.6. Two Compactness Theorems. 4. Point Set Topology Exercise 6 The product of countably many sequentally compact
More informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationThis article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author s institution, sharing
More informationHomework 8 Solutions to Selected Problems
Homework 8 Solutions to Selected Problems June 7, 01 1 Chapter 17, Problem Let f(x D[x] and suppose f(x is reducible in D[x]. That is, there exist polynomials g(x and h(x in D[x] such that g(x and h(x
More information