A Technique for the Numerical Solution of Certain Integral Equations of the First Kind*
|
|
- David Lucas
- 5 years ago
- Views:
Transcription
1 A Technique for the Numerical Solution of Certain Integral Equations of the First Kind* DAVID L. PHH~LIPSt Argonne National Laboratory, Argonne, Illinois Introduction The general linear equation may be written as h(x)f(x) --~ I p Ja b K(x, y)f(y) dy = g(x) (a<=x<=b) where the known functions h(x), K(x, y) and g(x) are assumed to be bounded and usually to be continuous. If h(x) ~- the equation is of first kind; if h(x) ~ for a -<_ x ~ b, the equation is of second kind; if h(x) vanishes somewhere but not identically, the equation is of third kind. If the range of integration is infinite or if the kernel K(x, y) is not bounded, the equation is singular. Here we will consider only nonsingular linear integral equations of the first kind: b ~ K(x,y)f(y) dy = g(x) (a ~ x ~ b) (l) There is extensive literature on equations of the second kind, but literature on linear equations of the first kind is sparse. However, several, methods for solving equations of the first kind numerically have been proposed [1-1]. No method has been very successful for arbitrary kernels when the function g(x) is known with only modest accuracy. The reason for this is inherent in the equation itself. Think of the equation as a linear operator, operating on f(y) to produce g(x). This operator does not have a bounded inverse (it may not even have an inverse, but we will assume here that it does) which can be seen as follows. Let f(y) be the solution to (1) and add to it the function fm = sin my. /, b For any integrable kernel it is known that gm= J. K(x, y) sin (my) dy -~ a as m -~ ~. Hence only an infinitesimal change g,~ in g causes s finite change f~ inf (i.e. the equation is unstable).,also, one would expect that g,~ -~ as m --~ faster for flat smooth kernels than for sharply peaked kernels (indeed if K(x, y) were the a-function, K(x, y) = a(x -- y), then g,~ = fm would not approach zero). Hence we conclude that the success in solving equation (1) by any method depends to a large extent on the accuracy of g(x) and the shape of K(x, y). * Received June, t Based on work performed under the auspices of the U. S. Atomic Energy Commission. 84
2 CERTAIN INTEGRAL EQUATIONS OF THE FIRST KIND 85 Description of Technique If the straightforward matrix approximation to equation (1) is used, it is found that as the mesh width decreases, the solutions at first become more accurate, but eventually begin to get worse. How soon the solutions begin to get worse depends on the accuracy of g(x). Larger errors in g cause the solutions to get worse sooner. Also, the error in each of the approximate solutions tends to be an oscillatory function of x. Since the function g(x) is not known accurately, we should state the problem as f v K(x, y)f(y) dy = g(x) + e(x) (a <- x <- b) (2) a where e(x) is an arbitrary function except for some condition on the size of ~(x), such as I ~(x) I =< M or p(x)e2(x) dx < 2~, p(x) >. Instead of a F a unique solution of (2) we get a family ~ of solutions. The real problem then is to pick out of the family of functions ~ the true solution f. This cannot be done without more information about the problem than is given in equation (2). We will assume here that the functional form of f is not known. If it were known we could use a least square fit to find a best fit to f. However, we will assume here that f is a reasonably smooth function. With this assumption the best approximation tof we can choose is probably the function f8 C ~ which is the smoothest in some sense. Of the various smoothness conditions, we choose the following (assuming the f have pieeewise continuous second derivatives): f b (f~t)2dx = min f b (ff)2dx. (3) In order to solve numerically (2) and (3) we make a matrix approximation to them. We subdivide the interval into n parts by the points a = x < xl < x2 < < x r, = b, and replace the integral equation (2) by the following linear system, i= wikiifi = gi + ei (j = O, 1,...,n) (4) where f~ = f(xi), gi = g(xi), ei = e(xj), lcii =- K(xi, x~), and the w~ are weight factors whose values depend on the quadrature formula used. However, for simplicity, we will assume that the x~ are uniformly spaced. For the condition on the magnitude of e(x) we take i~ ~ e ~ (5) 2 where the where d is a constant (more generally one can take ~= p~e~ = e 2 p~ = are weights). The anmogous problem to (3) is to look for the vector
3 86 DAVID L. PHILLIPS f~ = (H, fl ~, "'", f~) such that 2 " ~ ~ ~ (f,+l - 2 f,~-l) (6) i= f~3* i= where ~* is the set of vectors satisfying (4) and (5) (here we assume f is zero Outside the interval (a, b) and set f-1 = f~+l = ). LeE us introduce the following matrix notation: A =- (w~lcj,), A -~ =- (a,j). Equation (4) can be written in the form f = A-lg + A-I~. (7) This shows that the f~ are linear functions of the ~i and that of, _ ~ (i,j = o, 1,..., n). (s) Oei From (6), (8), and the constraint (5), we see that the conditions f, must satisfy ean be written as follows: i= El = e i= (f~+l -- 2H + ]~-1) (ai+l.~ -- 2aii + ai-l,s) + 7-1ej = (j = O, 1,..., n) (9) where ~-1 is the Lagrangian multiplier. It will be much simpler to solve these equations for f, if we take ~, as known and e 2 as an unknown instead of the other way around. This is equivalent to replacing the conditions (5) and (6) by the single condition that f, be the vector (satisfying (4)) that minimizes the expression ~' ~-~ (S,+~ -- 2f, +f,-~)2 + ~ ~,2 (1) i~ i~ where ~/is a given constant. From (1) we see tha~ ~/should be non-negative and that for ~ = the whole problem reduces to just the straightforward matrix approximation to (1), Af = g. Equations (9) can be written in the matrix form where B = (~,) and ~,BL+ ~ = (11) fhk = ak-2,z -- 4ak--l.l + 6akz -- 4ak+l,, + ak+~,l (k,l =, 1, "'., n). (12) In (12) we define a_2,~ = --aoa, a-la =, a~+t,z =, and a~+%z = -a~t. Using (11) and the fact that f, satisfies (7) we solve for f, and eand gee f. = (A + ~,B)-~g, r = -~,Bf,. (13)
4 CERTAIN INTEGRAL EQUATIONS OF THE FIRST KIND 87 From equations (13), the matrix method described here merely replaces A by A + ~B where B is a certain matrix whose elements depend only on A and ~/ is an arbitrary non-negative parameter which controls the amount of smoothing. Increasing ~ produces greater smoothing. It follows from the second equation in (13) and (5) that e is approximately proportional to ~/so that only a very few v~flues of ~ need to be used in order to find one giving e the desired magnitude. However, one additional matrix inversion is needed for each new value of - used. The value of e is determined from the accuracy of the g~-.,vzamples Example 1. Let the problem be the following: where ff K(x - 3 X)f(z) dx = g(x) 7rz K(z) = 1+ cosy, l z] =< 3, =, lzl>3; ( 1 3) 9 ~rk g(z) = (6+X) 1--~cos - ~sin-~-, Izl =< 6, =, Izl>6. The solution to this problem is f(x) = K(x). Hence we can easily check the numerical solution against the true solution. Let us first take n = 12 (13 points) and use Simpson's rule for the quadrature formula. The truncation error is about.4. The values of g are rounded off so that the maximum error in g(kj) is.5. Table I and Figure 1 give the comparison between the true solution and the numerical solution for several choices of % Since the solution is symmetric TABLE I k true vaiues of/,y ~, =.11 ~'=.11 "t' ~.3 ~, =.1 T ~.5 ~'= ~ 6~.1.(.{].G [ [ [ ) ) ) ave. 1~5 I max. ]e~- [
5 ~ DAVID L. t'h[liaps about zero, outy uonmegative vaues of ;', are listed in the ta, blc, ;\1~(), t t.~ average ej i and maximum %. i are listed for each % Here the trtmc::~tion error, which acf, s much like an error i~.q, is bound(.,d by.4. He~ee we expect to get, the best results for some value of 7 for which avg. ei i <.4. Here v... 3 is the best choice even t, hough a, ve.. I ei] =.26, which is eov, siderabty sma, ller Wan.4. Now consider the same problem with -n (25 points). Here the truncation error is bmmded by.25. Table II "TRUE CURVE, ~.y =... Y:.11 j,~j " +7 :.@3,'=.5 I j '\ -a -4 J 2 I _ I I t 1 J I ]" I I FIG, 1 TABLE II [Pde values off ~ = v=.25 ~ =.2 7 =.1 v =.5 ~ = I 2, (~}.5 1.8(~ ~g fN 2.5O % i.;~ ( ~) (K) O , , ,312 2, (X), , , ,8531 1,48(51 1.9, ,18 --,531 --, E , t!.t7.3r P.155 i.28 max. '~ ~i i....(x)27.59,122 i k :... :...
6 CERTAIN INTEGRAL EQUATIONS OF THE FIRST KIND 89 =, = '1 I ] t- + + ~ ' t I I i I! I t I I I I I FIG. 2 TABLE III X true values off ~, = ~'=.1 ~ =.1 3' ~.1 '~ =.1 ~' = ave. [ ej" [ max. [ ej I and Figure 2 give the comparison between the true solution and the numerical solution for several values of ~. Here the smoothing effect of ~, on the oscillation due to the ill-condition of the matrix is much greater than the smoothing effect on the solution. Thus here it is possible to practically eliminate the oscillation due to the ill-conditioned matrix without appreciably affecting the solution. Now consider the same problem where the g's have been further rounded off.
7 9 DAVID L. PHILLIPS In this case the maximum error in g is.4 and the average error in g is.14. Table IIt and Figure 3 give the comparison between the true solution and the numerical solution for severm values of 7. 7: + 7 "=.1 TRUE I I I FIG. 3 TABLE IV X true values of f 7 = 7 =.5 7 =.5 7 =.5 7 = , i , ave. I e~ I... } I max. I ~J!... ] ,84
8 t5 --!4 -- I3 -- CERTAIN INTEGRAL EQUATIONS OF THE FIRST KIND ~ Z =o +~'= II -- I I I i FIG. 4 Notice that for 1' = there is considerably more oscillation than in the previous case. However, the oscillation is smoothed out just as before by using a suitable value of % Now consider the same problem where the g's have been even further rounded off. In this case the average error in I g I is.2 and the maximum error is.41. Table IV and Figure 4 give the comparison between the true solution and the numerical solution for several values of % Notice that for 7 = the solution is completely dominated by oscillation. The errors in g have been magnified in f so much that the solution is completely useless. However, for ~, =.5 the solution is reasonably good, especially near the peak which is the most important region. Also notice that the average ] ej] for =.5 is of the same order of magnitude as the average error of the g's.
9 92 DAVID L. PHILLIPS Example 2, The problem is the following: 3 dz f t<(x - x)f(z) = g(x) 3 where K(z), g(z) and f(z) are given in tabular form in "'t' ta)m V. K(z) and g(z) are also plotted in ]~,gure -a. o 5. The values of g were obtained from tile tabulated values of f(z) and K(z) by numerical integration. The values of g have errors which average about.t in magnitude. The maximum error in g is about,.2 Let us first take n - 14 (15 points) and use Simpson's rule for the quadrature formula and solve for f(z) assuming K(z), g(z) given. The truncation error is less than.66. Table VI and Figure 6 give the comparison between the true solution and the numerical solution for several choices of 1/. g (z) K (z) I TABLE V f(z) g(z) K(z) l(z) K(z) = for l z i > ~K (z) 3 - I IO 2 5 FIG. 5
10 CERTAIN INTEGRAL EQUATIONS OF THE FIRST KIND 93 TABLE VI z true values off ~= --I i 7 = ~/ = = "/= = 1? = ave. [ ej : max. ] ej s: /I~.~ = O, n = I --.I TRUE +x:l, n:3 _j-y= O, n=14 14, O FzG. 6
11 - -.] 94 DAVID L. PHILLIPS The poor results here are due to a large truncation error. The parameter 3' has to be chosen so large that the true shape of the function is noticeably af- fected. Hence t.he oscillation and the true shape cannot be separated. Now consider the same problem where we take n = 3 (31 points). Here the upper bound for the truncation error is.2, which is the same order of mag- nitude as the errors in g. Table VII and Figure 6 give the comparison between the true solution and the numerical solution for several choices of 7. TABLE VII z true values off O ] = 7 =.1 7 =.1 i =2 7~i ave. l~s I... i max. t osl... i --O oo i.4.2 i
12 CERTAIN INTEGRAL EQUATIONS OF TIIE FIRST KIND 95 TABLE VIII I trueo~llles 3,= 7=.1 7=I 3,~2 3,=5 T=IO O OO28.5 l I , : O ave. I ej { max. l ej! Notice that the best values of 3' are the ones for which the max I e~] is approximately the same as the maximum error in g. Also notice that the peak is rounded off a small amount. This is unavoidable with thepresent method. Finally, consider the same problem where the g's are given less accurately. Let the errors be about 5 per cent of the peak value (This simulates a certain problem from experimental physics ) The maximum error in g is about 15. Table VIII and Figure 7 give the comparison between the true solution (for accurate g) and the numerical solution for several values of 7.
13 96 DAVID L. PtlI:LLIP~ = +y: TRUE I , l I -.t8 }-- -zo ~ I -28 I I I 1. I I I t t I 1 _ -2 < Fro. 7 Notice that for T = the solution is dominated by oscillation. The best values of v are 1 and 2. For these values the error in the peak value is about 7per cent, a reasonable value since the maximum error in g is about 5 per cent. For the best values of T the max l eji is about.1, somewhat lower than the maximum error in g. This is probably because of the unusually large errors in ~?. Sumrr~ary The method described here should work reasonably well on problems where the unknown function f is assumed to be a relatively smooth function. The truncation error should be as small a.s, or smaller than the errors in g. The most difficult task is to choose v. When the errors in g are relatively small, ~/should probably be chosen such that. e is approximately the same magtfitude as t:hese errors. When the errors in g are relatively large, as they are in the last part of example 2, ~/should probably be chosen such that e is somewhat smaller l~han
14 CERTAIN INTEGRAL EQUATIONS OF THE FIRST KIND 97 these errors. In any case several values of 7 should be tried and the best value should be the one that appears to take out the oscillation without appreciably smoothing the function f. REFERENCES 1. C~ou% P. D. J. Math. Phys. 19 (194), CROUT, P. D.; AND HILDEBRAND, F.B. J. Math. Phys. 2 (1941), DIXON, W. R.; AND AITKEN, J. t-i. Canad. J. Phys. 36 (1958), FOX, T.; AND GOODWIN, E.J. Philos. Trans. Roy. Soc. London. A 245 (1953), GOLD, R.; AND SCOFIELD, N. E. Bull. Amer. Phys. Soc. 2 (196), t~reisel, G. Proc. Roy. Soc. London. A 197 (1949), N~(STnO~, E.J. Acta Math. 5~ (193), R~z, A. Ark. Mat. Astr. Fys. 29A No. 29 (1943). 9. VAN DE t-iulst, H. C. Bull. Astr. Inst. Netherlands 1 (1946), YouNo, A. Proc. Roy. Soc. London. A224 (1954), 561.
ESTIMATES FOR MAXIMAL SINGULAR INTEGRALS
ESTIMATES FOR MAXIMAL SINGULAR INTEGRALS LOUKAS GRAFAKOS Abstract. It is shown that maximal truncations of nonconvolution L -bounded singular integral operators with kernels satisfying Hörmander s condition
More information18.02 Multivariable Calculus Fall 2007
MIT OpenCourseWare http://ocw.mit.edu 18.02 Multivariable Calculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. V9. Surface Integrals Surface
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory
More informationNumerical integration and differentiation. Unit IV. Numerical Integration and Differentiation. Plan of attack. Numerical integration.
Unit IV Numerical Integration and Differentiation Numerical integration and differentiation quadrature classical formulas for equally spaced nodes improper integrals Gaussian quadrature and orthogonal
More informationLECTURE 5, FRIDAY
LECTURE 5, FRIDAY 20.02.04 FRANZ LEMMERMEYER Before we start with the arithmetic of elliptic curves, let us talk a little bit about multiplicities, tangents, and singular points. 1. Tangents How do we
More informationHere, logx denotes the natural logarithm of x. Mrsky [7], and later Cheo and Yien [2], proved that iz;.(»)=^io 8 "0(i).
Curtis Cooper f Robert E* Kennedy, and Milo Renberg Dept. of Mathematics, Central Missouri State University, Warrensburg, MO 64093-5045 (Submitted January 1997-Final Revision June 1997) 0. INTRODUCTION
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationTHE ADJOINT OF A DIFFERENTIAL OPERATOR WITH INTEGRAL BOUNDARY CONDITIONS
738 a. m. krall [August References 1. L. V. Ahlfors and Leo Sario, Riemann surfaces, Princeton Univ. Press, Princeton, N. J., 1960. 2. H. B. Curtis, Jr., Some properties of functions which uniformize a
More informationPointwise convergence rates and central limit theorems for kernel density estimators in linear processes
Pointwise convergence rates and central limit theorems for kernel density estimators in linear processes Anton Schick Binghamton University Wolfgang Wefelmeyer Universität zu Köln Abstract Convergence
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More informationLimits at Infinity. Horizontal Asymptotes. Definition (Limits at Infinity) Horizontal Asymptotes
Limits at Infinity If a function f has a domain that is unbounded, that is, one of the endpoints of its domain is ±, we can determine the long term behavior of the function using a it at infinity. Definition
More information: Does there exist a t-space which is not bornologic?
294 [Vol. 30, 63. On Locally Convex Vector Spaces of Continuous Functions By Taira SmOTA Department of Mathematics, Osaka University (Comm. by K. KUNUGI, M.J.A., April 12, 1954) 1. In the present note
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More informationLAMINAR FLOW BEHIND A TWO-DIMENSIONAL GRID
LAMINAR FLOW BEHIND A TWO-DIMENSIONAL GRID BY L. I. G. KOVASZNAY Communicated by Sir GEOFFREY TAYLOR Received 20 May 1947 INTRODUCTION The production of a 'wake' behind solid bodies has been treated by
More informationMATH 104: INTRODUCTORY ANALYSIS SPRING 2008/09 PROBLEM SET 8 SOLUTIONS
MATH 04: INTRODUCTORY ANALYSIS SPRING 008/09 PROBLEM SET 8 SOLUTIONS. Let f : R R be continuous periodic with period, i.e. f(x + ) = f(x) for all x R. Prove the following: (a) f is bounded above below
More informationScientific Computing
2301678 Scientific Computing Chapter 2 Interpolation and Approximation Paisan Nakmahachalasint Paisan.N@chula.ac.th Chapter 2 Interpolation and Approximation p. 1/66 Contents 1. Polynomial interpolation
More informationON THE CONVERGENCE OF APPROXIMATE SOLUTIONS OF
ON THE CONVERGENCE OF APPROXIATE SOLUTIONS OF THE HEAT EQUATION TO THE EXACT SOLUTION WERNER LEUTERT Introduction. In many cases it is practically impossible to solve an initial value problem for a partial
More informationBanach Journal of Mathematical Analysis ISSN: (electronic)
Banach J. Math. Anal. 7 (2013), no. 1, 59 72 Banach Journal of Mathematical Analysis ISSN: 1735-8787 (electronic) www.emis.de/journals/bjma/ WEIGHTED COMPOSITION OPERATORS BETWEEN VECTOR-VALUED LIPSCHITZ
More informationQualifying Examination
Summer 24 Day. Monday, September 5, 24 You have three hours to complete this exam. Work all problems. Start each problem on a All problems are 2 points. Please email any electronic files in support of
More information8 Singular Integral Operators and L p -Regularity Theory
8 Singular Integral Operators and L p -Regularity Theory 8. Motivation See hand-written notes! 8.2 Mikhlin Multiplier Theorem Recall that the Fourier transformation F and the inverse Fourier transformation
More informationNUMERICAL METHODS. x n+1 = 2x n x 2 n. In particular: which of them gives faster convergence, and why? [Work to four decimal places.
NUMERICAL METHODS 1. Rearranging the equation x 3 =.5 gives the iterative formula x n+1 = g(x n ), where g(x) = (2x 2 ) 1. (a) Starting with x = 1, compute the x n up to n = 6, and describe what is happening.
More informationSMALL ZEROS OF QUADRATIC FORMS WITH LINEAR CONDITIONS. Lenny Fukshansky. f ij X i Y j
SALL ZEROS OF QUADRATIC FORS WITH LINEAR CONDITIONS Lenny Fukshansky Abstract. Given a quadratic form and linear forms in N + 1 variables with coefficients in a number field K, suppose that there exists
More informationMonte Carlo Methods for Statistical Inference: Variance Reduction Techniques
Monte Carlo Methods for Statistical Inference: Variance Reduction Techniques Hung Chen hchen@math.ntu.edu.tw Department of Mathematics National Taiwan University 3rd March 2004 Meet at NS 104 On Wednesday
More informationMATH 220 solution to homework 4
MATH 22 solution to homework 4 Problem. Define v(t, x) = u(t, x + bt), then v t (t, x) a(x + u bt) 2 (t, x) =, t >, x R, x2 v(, x) = f(x). It suffices to show that v(t, x) F = max y R f(y). We consider
More informationDiffraction by a Half-Plane LL. G. CHAMBERS. {Received 29th March Read 5th May 1950.)
Diffraction by a Half-Plane By LL. G. CHAMBERS {Received 29th March 1950. Read 5th May 1950.) 1. Introduction. The diffraction of a simple harmonic wave train by a straightedged semi-infinite screen was
More informationDavid A. Stephens Department of Mathematics and Statistics McGill University. October 28, 2006
556: MATHEMATICAL STATISTICS I COMPUTING THE HYPEBOLIC SECANT DISTIBUTION CHAACTEISTIC FUNCTION David A. Stephens Department of Mathematics and Statistics McGill University October 8, 6 Abstract We give
More informationAn idea how to solve some of the problems. diverges the same must hold for the original series. T 1 p T 1 p + 1 p 1 = 1. dt = lim
An idea how to solve some of the problems 5.2-2. (a) Does not converge: By multiplying across we get Hence 2k 2k 2 /2 k 2k2 k 2 /2 k 2 /2 2k 2k 2 /2 k. As the series diverges the same must hold for the
More information(2) ^max^lpm g M/?"[l - ^-(l - /T1)2}
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 56, April 1976 SOME INEQUALITIES FOR POLYNOMIALS Q. I. RAHMAN Abstract. Let pn(z) be a polynomial of degree n. Given that pn(z) has a zero on the
More informationA Nonstandard Fourier Inequality
A Nonstandard Fourier Inequality Jonas Azzam, Bobbie Buchholz, Ian Grooms Gretchen Hagge, Kyle Hays, Greg Norgard Abstract We consider a class of functions given by series of the form h(x) := f k e a kx
More informationSMOOTHNESS OF FUNCTIONS GENERATED BY RIESZ PRODUCTS
SMOOTHNESS OF FUNCTIONS GENERATED BY RIESZ PRODUCTS PETER L. DUREN Riesz products are a useful apparatus for constructing singular functions with special properties. They have been an important source
More informationTHE POISSON TRANSFORM^)
THE POISSON TRANSFORM^) BY HARRY POLLARD The Poisson transform is defined by the equation (1) /(*)=- /" / MO- T J _M 1 + (X t)2 It is assumed that a is of bounded variation in each finite interval, and
More informationarxiv: v1 [math.ca] 29 Dec 2018
A QUANTITATIVE WEIGHTED WEAK-TYPE ESTIMATE FOR CALDERÓN-ZYGMUND OPERATORS CODY B. STOCKDALE arxiv:82.392v [math.ca] 29 Dec 208 Abstract. The purpose of this article is to provide an alternative proof of
More informationWeighted norm inequalities for singular integral operators
Weighted norm inequalities for singular integral operators C. Pérez Journal of the London mathematical society 49 (994), 296 308. Departmento de Matemáticas Universidad Autónoma de Madrid 28049 Madrid,
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 informationTHE MAXIMUM OF SUMS OF STABLE RANDOM VARIABLES
THE MAXIMUM OF SUMS OF STABLE RANDOM VARIABLES BY D. A. DARLING(') 1. Introduction. Let Xu X2, be identically distributed independent random variables and set Sn=Xi+ +X. In this paper we obtain the limiting
More informationDepartment of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016
Department of Mathematics, University of California, Berkeley YOUR 1 OR 2 DIGIT EXAM NUMBER GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016 1. Please write your 1- or 2-digit exam number on
More informationOdds find E n d s ENTRIES. and A L T O SOLO. iiik places for many species of liik JMtHSIDKVrS r ITCH I
LG Bx L G L L FY-F Y LLL G Y 593 JBL F GBB 3 P P P P P B P P P g J g -g ( J z x 000 P gg - F g L g g L YB J 933 g 95 Y LL g g g g- PY L FG LLL F LLX F q F PP06 BG P P P F g 000 g F L 93 PL - L g g G gg
More informationAPSIDAL ANGLES FOR SYMMETRICAL DYNAMICAL SYSTEMS WITH TWO DEGREES OF FREEDOM
APSIDAL ANGLES FOR SYMMETRICAL DYNAMICAL SYSTEMS WITH TWO DEGREES OF FREEDOM J. L. SYNGE The bob of a spherical pendulum (or particle on a smooth sphere under gravity) oscillates between two levels, and
More informationBOUNDARY-VALUE TECHNIQUES FOR THE NUMERICAL SOLUTION OF INITIAL-VALUE PROBLEMS IN ORDINARY DIFFERENTIAL EQUATIONS
BOUNDARY-VALUE TECHNIQUES FOR THE NUMERICAL SOLUTION OF INITIAL-VALUE PROBLEMS IN ORDINARY DIFFERENTIAL EQUATIONS By L. FOX (The National Physical Laboratory, Teddington, Middlesex) and A. R. MITCHELL
More informationHence a root lies between 1 and 2. Since f a is negative and f(x 0 ) is positive The root lies between a and x 0 i.e. 1 and 1.
The Bisection method or BOLZANO s method or Interval halving method: Find the positive root of x 3 x = 1 correct to four decimal places by bisection method Let f x = x 3 x 1 Here f 0 = 1 = ve, f 1 = ve,
More informationVolterra integral equations solved in Fredholm form using Walsh functions
ANZIAM J. 45 (E) ppc269 C282, 24 C269 Volterra integral equations solved in Fredholm form using Walsh functions W. F. Blyth R. L. May P. Widyaningsih (Received 8 August 23; revised 6 Jan 24) Abstract Recently
More informationEXAMPLES OF PROOFS BY INDUCTION
EXAMPLES OF PROOFS BY INDUCTION KEITH CONRAD 1. Introduction In this handout we illustrate proofs by induction from several areas of mathematics: linear algebra, polynomial algebra, and calculus. Becoming
More informationSolutions Final Exam May. 14, 2014
Solutions Final Exam May. 14, 2014 1. Determine whether the following statements are true or false. Justify your answer (i.e., prove the claim, derive a contradiction or give a counter-example). (a) (10
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Simplifying integration for logarithmic singularities R.N.L. Smith Department ofapplied Mathematics & OR, Cranfield University, RMCS, Shrivenham, Swindon, Wiltshire SN6 SLA, UK Introduction Any implementation
More informationLinear DifferentiaL Equation
Linear DifferentiaL Equation Massoud Malek The set F of all complex-valued functions is known to be a vector space of infinite dimension. Solutions to any linear differential equations, form a subspace
More informationSection 5.5 More Integration Formula (The Substitution Method) 2 Lectures. Dr. Abdulla Eid. College of Science. MATHS 101: Calculus I
Section 5.5 More Integration Formula (The Substitution Method) 2 Lectures College of Science MATHS : Calculus I (University of Bahrain) Integrals / 7 The Substitution Method Idea: To replace a relatively
More informationMath 106 Calculus 1 Topics for first exam
Chapter 2: Limits and Continuity Rates of change and its: Math 06 Calculus Topics for first exam Limit of a function f at a point a = the value the function should take at the point = the value that the
More informationMath 140A - Fall Final Exam
Math 140A - Fall 2014 - Final Exam Problem 1. Let {a n } n 1 be an increasing sequence of real numbers. (i) If {a n } has a bounded subsequence, show that {a n } is itself bounded. (ii) If {a n } has a
More informationParts Manual. EPIC II Critical Care Bed REF 2031
EPIC II Critical Care Bed REF 2031 Parts Manual For parts or technical assistance call: USA: 1-800-327-0770 2013/05 B.0 2031-109-006 REV B www.stryker.com Table of Contents English Product Labels... 4
More informationTHE EVALUATION OF LINEAR FORMS P. S. DWYER
PSYCHOMETRIKA~VOL. 6, NO. 6 DECEMBER. 1941 THE EVALUATION OF LINEAR FORMS P. S. DWYER UNIVERSITY OF MICHIGAN This presentation deals with the evaluation and transformation of linear forms. Especial emphasis
More informationLecture 8: Boundary Integral Equations
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 8: Boundary Integral Equations Gunnar Martinsson The University of Colorado at Boulder Research support by: Consider
More informationOSCILLATORY SINGULAR INTEGRALS ON L p AND HARDY SPACES
POCEEDINGS OF THE AMEICAN MATHEMATICAL SOCIETY Volume 24, Number 9, September 996 OSCILLATOY SINGULA INTEGALS ON L p AND HADY SPACES YIBIAO PAN (Communicated by J. Marshall Ash) Abstract. We consider boundedness
More informationMath Stochastic Processes & Simulation. Davar Khoshnevisan University of Utah
Math 5040 1 Stochastic Processes & Simulation Davar Khoshnevisan University of Utah Module 1 Generation of Discrete Random Variables Just about every programming language and environment has a randomnumber
More information1. Introduction. Let P be the orthogonal projection from L2( T ) onto H2( T), where T = {z e C z = 1} and H2(T) is the Hardy space on
Tohoku Math. Journ. 30 (1978), 163-171. ON THE COMPACTNESS OF OPERATORS OF HANKEL TYPE AKIHITO UCHIYAMA (Received December 22, 1976) 1. Introduction. Let P be the orthogonal projection from L2( T ) onto
More informationEquipartition of energy in wave motion
Carnegie Mellon University Research Showcase @ CMU Department of Mathematical Sciences Mellon College of Science 1969 Equipartition of energy in wave motion Richard James Duffin Carnegie Mellon University
More informationThe backward continued fraction map and geodesic flow
Ergod. Th. & Dynam. Sys. (1984), 4, 487-492 Printed in Great Britain The backward continued fraction map and geodesic flow ROY L. ADLER AND LEOPOLD FLATTO IBM Thomas J. Watson Research Center, Yorktown
More informationON THE AVERAGE NUMBER OF REAL ROOTS OF A RANDOM ALGEBRAIC EQUATION
ON THE AVERAGE NUMBER OF REAL ROOTS OF A RANDOM ALGEBRAIC EQUATION M. KAC 1. Introduction. Consider the algebraic equation (1) Xo + X x x + X 2 x 2 + + In-i^" 1 = 0, where the X's are independent random
More informationElliptic Curves and Public Key Cryptography
Elliptic Curves and Public Key Cryptography Jeff Achter January 7, 2011 1 Introduction to Elliptic Curves 1.1 Diophantine equations Many classical problems in number theory have the following form: Let
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... find the approximate solutions of derivative (first- and second-order) and antiderivative (definite integral only). Numerical Differentiation
More information2.12: Derivatives of Exp/Log (cont d) and 2.15: Antiderivatives and Initial Value Problems
2.12: Derivatives of Exp/Log (cont d) and 2.15: Antiderivatives and Initial Value Problems Mathematics 3 Lecture 14 Dartmouth College February 03, 2010 Derivatives of the Exponential and Logarithmic Functions
More informationReading. w Foley, Section 11.2 Optional
Parametric Curves w Foley, Section.2 Optional Reading w Bartels, Beatty, and Barsky. An Introduction to Splines for use in Computer Graphics and Geometric Modeling, 987. w Farin. Curves and Surfaces for
More informationA NOTE ON TRIGONOMETRIC MATRICES
A NOTE ON TRIGONOMETRIC MATRICES garret j. etgen Introduction. Let Q(x) he an nxn symmetric matrix of continuous functions on X: 0^x
More informationMulti-Linear Mappings, SVD, HOSVD, and the Numerical Solution of Ill-Conditioned Tensor Least Squares Problems
Multi-Linear Mappings, SVD, HOSVD, and the Numerical Solution of Ill-Conditioned Tensor Least Squares Problems Lars Eldén Department of Mathematics, Linköping University 1 April 2005 ERCIM April 2005 Multi-Linear
More informationOn Additive Polynomials
On Additive Polynomials Zachary Scherr 1 Preliminaries The purpose of this article is to explore the following question: Question 1. Let K be a field. For which f(z) K[z] does f(a+b) = f(a)+f(b) for all
More informationMath 480 The Vector Space of Differentiable Functions
Math 480 The Vector Space of Differentiable Functions The vector space of differentiable functions. Let C (R) denote the set of all infinitely differentiable functions f : R R. Then C (R) is a vector space,
More informationMath 671: Tensor Train decomposition methods II
Math 671: Tensor Train decomposition methods II Eduardo Corona 1 1 University of Michigan at Ann Arbor December 13, 2016 Table of Contents 1 What we ve talked about so far: 2 The Tensor Train decomposition
More informationPrinciple of Mathematical Induction
Advanced Calculus I. Math 451, Fall 2016, Prof. Vershynin Principle of Mathematical Induction 1. Prove that 1 + 2 + + n = 1 n(n + 1) for all n N. 2 2. Prove that 1 2 + 2 2 + + n 2 = 1 n(n + 1)(2n + 1)
More informationSolution to Set 7, Math 2568
Solution to Set 7, Math 568 S 5.: No. 18: Let Q be the set of all nonsingular matrices with the usual definition of addition and scalar multiplication. Show that Q is not a vector space. In particular,
More informationMotion in Three Dimensions
Motion in Three Dimensions We ve learned about the relationship between position, velocity and acceleration in one dimension Now we need to extend those ideas to the three-dimensional world In the 1-D
More informationNumerical Quadrature Over a Rectangular Domain in Two or More Dimensions
Numerical Quadrature Over a Rectangular Domain in Two or More Dimensions Part 1. Quadrature over a Square, Using up to Sixteen Equally Spaced Points By J. C. P. Miller 1. Introduction. Except for a section
More informationNumerical Analysis: Solutions of System of. Linear Equation. Natasha S. Sharma, PhD
Mathematical Question we are interested in answering numerically How to solve the following linear system for x Ax = b? where A is an n n invertible matrix and b is vector of length n. Notation: x denote
More informationMathematics Course 111: Algebra I Part I: Algebraic Structures, Sets and Permutations
Mathematics Course 111: Algebra I Part I: Algebraic Structures, Sets and Permutations D. R. Wilkins Academic Year 1996-7 1 Number Systems and Matrix Algebra Integers The whole numbers 0, ±1, ±2, ±3, ±4,...
More informationLEBESGUE INTEGRATION. Introduction
LEBESGUE INTEGATION EYE SJAMAA Supplementary notes Math 414, Spring 25 Introduction The following heuristic argument is at the basis of the denition of the Lebesgue integral. This argument will be imprecise,
More informationNOWHERE LOCALLY UNIFORMLY CONTINUOUS FUNCTIONS
NOWHERE LOCALLY UNIFORMLY CONTINUOUS FUNCTIONS K. Jarosz Southern Illinois University at Edwardsville, IL 606, and Bowling Green State University, OH 43403 kjarosz@siue.edu September, 995 Abstract. Suppose
More information18 BESSEL FUNCTIONS FOR LARGE ARGUMENTS. By M. Goldstein and R. M. Thaler
18 BESSEL FUNCTIONS FOR LARGE ARGUMENTS Bessel Functions for Large Arguments By M. Goldstein R. M. Thaler Calculations of Bessel Functions of real order argument for large values of the argument can be
More informationTHE JOHN-NIRENBERG INEQUALITY WITH SHARP CONSTANTS
THE JOHN-NIRENBERG INEQUALITY WITH SHARP CONSTANTS ANDREI K. LERNER Abstract. We consider the one-dimensional John-Nirenberg inequality: {x I 0 : fx f I0 > α} C 1 I 0 exp C 2 α. A. Korenovskii found that
More informationChapter 4. Continuous Random Variables
Chapter 4. Continuous Random Variables Review Continuous random variable: A random variable that can take any value on an interval of R. Distribution: A density function f : R R + such that 1. non-negative,
More informationOn a combinatorial method for counting smooth numbers in sets of integers
On a combinatorial method for counting smooth numbers in sets of integers Ernie Croot February 2, 2007 Abstract In this paper we develop a method for determining the number of integers without large prime
More informationSPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS
SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly
More informationSystems of Equations
Prerequisites: Solving simultaneous equations in 2 variables; equation and graph of a straight line. Maths Applications: Finding circle equations; finding matrix inverses; intersections of lines and planes.
More informationOn Systems of Diagonal Forms II
On Systems of Diagonal Forms II Michael P Knapp 1 Introduction In a recent paper [8], we considered the system F of homogeneous additive forms F 1 (x) = a 11 x k 1 1 + + a 1s x k 1 s F R (x) = a R1 x k
More informationh(x) lim H(x) = lim Since h is nondecreasing then h(x) 0 for all x, and if h is discontinuous at a point x then H(x) > 0. Denote
Real Variables, Fall 4 Problem set 4 Solution suggestions Exercise. Let f be of bounded variation on [a, b]. Show that for each c (a, b), lim x c f(x) and lim x c f(x) exist. Prove that a monotone function
More informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Taylor s Theorem Can often approximate a function by a polynomial The error in the approximation
More informationON A CHARACTERIZATION OF THE KERNEL OF THE DIRICHLET-TO-NEUMANN MAP FOR A PLANAR REGION
ON A CHARACTERIZATION OF THE KERNEL OF THE DIRICHLET-TO-NEUMANN MAP FOR A PLANAR REGION DAVID INGERMAN AND JAMES A. MORROW Abstract. We will show the Dirichlet-to-Neumann map Λ for the electrical conductivity
More information8.3 Partial Fraction Decomposition
8.3 partial fraction decomposition 575 8.3 Partial Fraction Decomposition Rational functions (polynomials divided by polynomials) and their integrals play important roles in mathematics and applications,
More informationAnalysis III. Exam 1
Analysis III Math 414 Spring 27 Professor Ben Richert Exam 1 Solutions Problem 1 Let X be the set of all continuous real valued functions on [, 1], and let ρ : X X R be the function ρ(f, g) = sup f g (1)
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 informationbriefly consider in connection with the differential equations of the type G catenaries, (5) the co I
VOL. 28, 1942 MA THEMA TICS: E. KA SNER 333 Kasner, "Transformation Theory of Isothermal Families and Certain Related Trajectories," Revista de Matematica, 2, 17-24 (1941). De Cicco, "The Two Conformal
More informationUses of duality. Geoff Gordon & Ryan Tibshirani Optimization /
Uses of duality Geoff Gordon & Ryan Tibshirani Optimization 10-725 / 36-725 1 Remember conjugate functions Given f : R n R, the function is called its conjugate f (y) = max x R n yt x f(x) Conjugates appear
More informationOPERATORS WITH SINGULAR CONTINUOUS SPECTRUM, V. SPARSE POTENTIALS. B. Simon 1 and G. Stolz 2
OPERATORS WITH SINGULAR CONTINUOUS SPECTRUM, V. SPARSE POTENTIALS B. Simon 1 and G. Stolz 2 Abstract. By presenting simple theorems for the absence of positive eigenvalues for certain one-dimensional Schrödinger
More informationAn Introduction to Differential Algebra
An Introduction to Differential Algebra Alexander Wittig1, P. Di Lizia, R. Armellin, et al. 1 ESA Advanced Concepts Team (TEC-SF) SRL, Milan Dinamica Outline 1 Overview Five Views of Differential Algebra
More informationNumerical Quadrature over a Rectangular Domain in Two or More Dimensions
Numerical Quadrature over a Rectangular Domain in Two or More Dimensions Part 2. Quadrature in several dimensions, using special points By J. C. P. Miller 1. Introduction. In Part 1 [1], several approximate
More information,,,,..,,., {. (, ),, {,.,.,..,,.,.,,....... {.. : N {, Z {, Q {, Q p { p{ {. 3, R {, C {. : ord p {. 8, (k) {.42,!() { {. 24, () { {. 24, () { {. 25,., () { {. 26,. 9, () { {. 27,. 23, '() { ( ) {. 28,
More informationNumerical Methods of Approximation
Contents 31 Numerical Methods of Approximation 31.1 Polynomial Approximations 2 31.2 Numerical Integration 28 31.3 Numerical Differentiation 58 31.4 Nonlinear Equations 67 Learning outcomes In this Workbook
More information10-701/ Recitation : Kernels
10-701/15-781 Recitation : Kernels Manojit Nandi February 27, 2014 Outline Mathematical Theory Banach Space and Hilbert Spaces Kernels Commonly Used Kernels Kernel Theory One Weird Kernel Trick Representer
More informationQuantum Mechanics for Scientists and Engineers. David Miller
Quantum Mechanics for Scientists and Engineers David Miller The particle in a box The particle in a box Linearity and normalization Linearity and Schrödinger s equation We see that Schrödinger s equation
More informationAlmost Sure Convergence in Extreme Value Theory
Math. Nachr. 190 (1998), 43-50 Almost Sure Convergence in Extreme Value Theory By SHIHONG CHENG of Beijing, LIANG PENG of Rotterdam, and YONGCHENG &I of Beijing (Received April 18, 1995) (Revised Version
More informationMY PUTNAM PROBLEMS. log(1 + x) dx = π2
MY PUTNAM PROBLEMS These are the problems I proposed when I was on the Putnam Problem Committee for the 984 86 Putnam Exams. Problems intended to be A or B (and therefore relatively easy) are marked accordingly.
More informationYears. Marketing without a plan is like navigating a maze; the solution is unclear.
F Q 2018 E Mk l lk z; l l Mk El M C C 1995 O Y O S P R j lk q D C Dl Off P W H S P W Sl M Y Pl Cl El M Cl FIRST QUARTER 2018 E El M & D I C/O Jff P RGD S C D M Sl 57 G S Alx ON K0C 1A0 C Tl: 6134821159
More information