On the Miquel point of simplices
|
|
- Bridget Watson
- 5 years ago
- Views:
Transcription
1 Elem. Math. 63 (2008) /08 Swiss Mathematical Society, 2008 I Elemente er Mathematik Lothar Heinrich Lothar Heinrich promovierte im Jahr 1980 an er Technischen Universität Dresen. Von war er Professor für Wahrscheinlichkeitstheorie un Statistik an er Technischen Universität Bergakaemie in Freiberg. Seit 1997 ist er Professor für Mathematik an er Universität Augsburg. Seine Forschungsinteressen liegen in er angewanten Wahrscheinlichkeitstheorie (stochastische Geometrie, große Abweichungen, Grenzwertsätze, räumliche Statistik, probabilistische Zahlentheorie). 1 Introuction an preliminaries If points are marke on each sie of a planar triangle, one on each sie (or on a sie's extension), then the three circles (each passing through a vertex an the marke points on the ajacent sies) are concurrent at a point M. This interesting fact was first prove an publishe by Augueste Miquel [3] in 1838, see also Weisstein [4] for further etails an extensions. This result is well-known in planar geometry as Miquel's theorem, an M is calle the Miquel point. However, much less known (even amongst geometers) is the following multiimensional generalization of Miquel's theorem: If one point is marke on each of the ( + I) /2 eges of a -simplex S(xo,Xl,...,xa) : {*o+ipi (xi-xo),lu,. l, t"i> 0,i - L,..., (l) l a t i:i i:7 Der Satz von Miquel ist ein klassisches Resultat er Elementargeometrie: Wir auf jeer Seite eines gegebenen Dreiecks oer eren Verlängerung ein Funkt beliebig festgelegt, un wir urch jeweils eine Ecke un ie beien markierten Punk;te auf en Nachbarseiten ein Kreis gezeichnet, so schneien sich iese rei Kreise in einem Punk;t. In er vorliegenen Arbeit wir mit einfachen Hilfsmitteln er linearen Algebra erstmals ein vollstitniger Beweis er Verallgemeinerung es Miquelschen Satzes auf -imensionale Simplizes angegeben: Wir auf en ( + l) /ZKanten je ein Punkt beliebig markiert, ann haben ie * 1 Kugeln, wovon jee urch je einen Eckpunkt un ie markierten Punkte auf en in iesen Eckpunkt einlaufenen Kanten festgelegt
2 r27 with positive volume an a sphere S; is rawn through each vertex x; an the points marke on the eges which meet in x;, then these + l spheres S;, i - 0, 1,...,, all meet in a point M which will also be calle Miquel's point in this note. By the best of the author's knowlege the only known proof of this result seems to be that given by Konnully t2l. Konnully's proof is base on the fact that there exists a common orthogonal sphere with respect to the so-calle Miquel spheres So, Sr S7 an it is shown that the raius of this sphere equals zero. However, such a sphere oes not always exist even in the planar case, namely if the unique point which has the same circle power w.r.t. three pairwise non-concentric circles lies in the interior of each of these circles. The aim of the present note is to provie a rigorous analytical proof which requires only simple facts from analytic geometry an linear algebra. As a by-prouctwe obtain a family of upper bouns of Gram's eterminant (incluing Haamar's inequality) which seems to be of interest in its own rieht. This auxiliarv result in Section 3 is vali in all Eucliean vector spaces. It shoul be mentione that an analogous construction with points marke on the ( - I)- faces (instea of on the eges) of the simplex oes in general not yiel a common point that belongs to all spheres. Let the points of the Eucliean space R be represente by column vectors x _ (u,..., x)t having the Eucliean nonn llxll - /E;S, where the scalar prouct (.,.) is efine by (y,z) : y'z : yi zl+... * ya za for y: (yr,...,!)' an z- (2t,...,2)l Furthermore, we recall the well-known fact froq analytic geometry that the circumsphere of the -simplex (1) consists of all points x IR. satisfying the equation llxll2 llxoll2 llxrll2 x X g x l llxa ll2 x 1-0, (2) where the left-han sie is for a ( * 2) x ( t 2) eterminant. Without loss of generality, we shall assume the vertex xs of the -simplex (1) to coincie withtheorigino:(0,...,0),antheverticesxi:(x1; linearly inepenent vectors, i.e. A::etXl0 (3) where X - (xr x6) enotes the quaratic matrix with columns Xl,..., X. For 0 < i < j <, let xil be a fixe point on the ege joining the vertices xi an x; being istinct from its en-points, i.e. vij : xi * l;j (x: - x;) for some )"ii 4 {0, 1}. For notationaleaseput Li; -I-)"ij forl < j anlii:l/2fori:0,i,...,.
3 t28 L. Heinrich In accorance with (2), the unique -sphere Ss passing through the origin o an the points X01,..., xg consists of all points x - (xr x)t e IR satisfying As(x) :- llxll2 o rfr,llxrll2 xf;ollxall2 xr 0 lorxrr Loaxu - 0. x 0 Lofiat Loa)c 1 I 1... Since 117:o Loi * 0, the latter equation is equivalento llxll2 lor llxr ll2 xt xn xl loa llxall2-0 (4) *o xr ": x Likewise, a point x : (xr x)t e IR belongs to the unique -sphere $ containing thevertexx;anthemarkepointsxgi,...,xi-t,i,xi,i+l,...,xiontheajacentegesif an onlv if fori-i,...,. llxll2 llxi ll2,r,fr, llxi ll2 llxi + r;; (xj- xi) ll2 xr xri )"oixtt xu I j - xu) a; (x):- - 0 *O *,a, ),Oirxai xai * X;1@A1 - xai) ffi By appealing to the well-known transformation rules for eterminants we obtain A;(x)-(),e;-1) llxll2 - x?r,llxillz (1* r.o;)llxill2 Äo; llxill2 llxi + l;;(xj - xi)ll2 x1 - )"six1; xfi 0 xu I Lij (xt j - xu) x -)"oixar *a, ö xai i Lij('xaj - xai) : (Äo; - r) II r,; i*i llxll2-1,fr, llxi ll2 cii ro; llxi ll2 cij x1 - )"gix1i xfi 0 x1j x -'Loixai *,a, ö x'j \/ j e {1,...,}\{l} je {1,...,}\{,}
4 where llxi + l;; (xj- xi) ll2 - llxi ll2 + (t + ro;) ll*ill2 )uij - lo; ll*ill2 + llxill2 -Lii llxl -x;ll2 for i, i =1,...,. The latter equality follows by using the ientity llu - vll2 - llull2 + llvll2-2 {a, v). Since cii : (1* ),0;) llxi ll2 we may simplify the previous eterminant by multiplying the secon column by Äo; an aing it to the first column. Consequently, in view of l";; {0, 1} for i + j, the equation A;(x) - 0 can be expresse as,,,,),,,,) llxll" + ^o; llxi ll- cit ci xt xn xh *o xt ": x Obviously, the set of points x lr satisfying both equations (4) an (5) coincies with the ( - l)-imensional sphere So n S;. By subtracting equation (4) from equation (5) an applying the summation law for eterminants iffering in only one row we obtain the linear equation t ^o; llxill- ' ) r r t cir - ^0l llxtll- ^ r t t ) ci - ^0llxll' t t ) xi xtl xi D; (x) :- - 0 *o xt ": x that hols for all x - (xr x)t IR belonging to the (uniquely etermine) hyperplane H; containing the ( - l)-sphere 56 n S,. Combining the equations D; (x) : 0, i -, yiels a system of linear equations whose solution (if it exists!) coincies with the point of intersection of the hyperplanes Hr,..., H. To fin this point we introuce the matrix Al : (o4)!,1:1with entries aii::cij -ro;llx;ll2:äo; llxill2+(1 -ro;) llx;ll2 _ Lii llxi-";112. (6) Letusfirstnotetheremarkablefactthat,inviewof Lij+Lii - lforl { j anlit -I/2, aij laji: llxill2 + llxjll2 - (xii tl1) llxi -x;ll2-2(x1,x;) for i, i :L,...,, which can be expresse concisely by I (o" + A1\ ) - G(x) :- ((xi, xi)!,i:t, where G(X) an et G(X) are calle Gram's matrix an Gram's eterminant, respectively, of the vectors Xl,..., X. In IR we have G(X) - X/X an thus et G(X) : L2. On the (7)
5 130 L. Heinrich other han, the relation (7) is meaningful in any real vector space V on which a symmetric, positive ef,nite, bilinear form (.,.) - briefly calle scalar prouct - is efine. Such vector spaces are usually calle Euclieanvector spaces. Proposition. For arbitrary linearly inepenent vectors Xt,..., X R an real numbers)";i,i, j :0, 1,...,, satisfuing)"ij*)"ii - Ian).ij f {0,I}fori, j :0, 1,...,, there exists auniquepointof intersection H1n... iha: x* : ("f,...,*ä),,whichis given by x*: XAnlxl with XÄ: (Äor llxrll2,...,lo llxall2),. (g) The first step in proving this result is to show that x - (\,..., x)t obeys the linear equations Di(x) -, 0 for i : I x1. This is left to the reaer as an exercise. To see the invertibility of the matrix An we use (7) an ecompose An as follows: Ar : Bn * G(X) with "n: (ln - n'n ) ) (e) Here the matrix Bn is skew-symmetric, i.e. Bl : -Bn. This skew-symmetry an the positive efiniteness of the Gram matrix G(X) show the quaratic form x/a1x - x/ G(X) x to be strictly positive for any x # o. This in turn implies An x # o for any x t' o, which is equivalent to etan I 0 for ali ),i1 satisfying Lij + Lii - 1. This combine with a simple continuity argument tells us that etan < 0 contraicts et G(X) ) 0, leaving.as the only possibility etal cvg(x)) > 0for anya > 0,entailingetBzr > 0bylettingcv + 0.HoweveretBn may take positive values only for even > 2 ue to the very efinition of skew-symmetry. Moreover, we shall boun etan uniformly from below by Gram's eterminant et G(X) in Section 3. Finally, note that the proposition oes not answer the question whether x* & hols for some or even ail i - 0, 1,...,. This will be the subiect of the next section. 2 The main result an its proof Theorem. Uner the conitions of the above proposition the point x* given in (8) belongs toeachof thespheres S;, i :0, 1,...,,thatis, x* is theuniquepointof intersection So f-l Sr n... f) Sa. In other wors, x* coincies with Miquel's point M of the simplex S(o, x1,..., xo) w.rt. the marke points xii: Xi * I;; (xj - xi) on its eges. Proof. After transposing an expaning the eterminant on the left-han sie of (4) along the first row we recognize that Ä9(x*) : 0 is equivalent to ll**ll2^-i "i j:1 xn lor llxr ll2 xr xr2 Lozllxzllz xz xr loa llxa ll2 x - 0. (10)
6 Miquel point M = (3.41,0.89, 1.54) xr = (3,6, 0) Fig. 1: Tetraheron,S(xg,x1,x2,x3)anthefourMiquelsphereswithcommonMiquel point M for ),61 : l/2, ),g2 - I/4, ),ß : 2/3, Lt2 : I/3, ),p : I/4, )'Y :2/3 Diviing by a t'. 0 the latter equation takes the form (X*, x* - zo) - components of zs - (zt za)' ut" given by 0, where the 1 zj:t xtl xt2 r.or llxr ll2 Lozllxzllz xl x2 for j:1,..., xt roa llxa ll2 x Applying Cramer's rule we see that z6 satisfies the linear equation X'zo _ xtr, i.e. (xi/llxill,zo) - Äo; llxill for i : I, sothatzs - (X,)-1 x,,". The geometric interpretation of the above relations reveals that the orthogonal projection of z6 onto the ege joining o an xi equals Xgi : Äo; xi for i : 1,...,. Furtherrnore, zg lies on the sphere Ss an, provie that x* So, the point zs/2 coincies with the centre of Ss by appealing to the converse of Thales' theorem. Hence, Ao(x*) : 0 hols if an only if (x*, x* - (x')-1xr,) : xi 1An1)'*' (*Ant - (x')-r ) *^ - xi crrx). : 0, where Cn :: (A;t)'X'XA;1 - (Ant),.
7 r32 L. Heinrich Fig. 2: Intersection of the tetraheron an the Miquel spheres in Fig. 1 with the plane z: parallel to the xy-plane Using (7) we fin that cn : I t,t;tl' (nn +n'n) lnt - (Ant)' :i (ant - (Ant)')- -air, z thus proving the skew-symmetry of the matrix Cn which is necessary an sufficient for the quaiatic f*- ti Crr x1 to isappear for any real.1"s1 )"sa anflr1l X1,..., x R, as we wishe to prove. ItremainstoshowthatA;(x*) -0fori - l,...,.forthiswetransposetheeterminant on the left-han sie of (5) an expan it along the first row leaing to the system of equations (llx* ll2 * Äo; llxi ll2) A - I j:i ti xtt cit xi x't ci x v j-th column In the next step we subtract equation (9) from the latter equation an ivie the ifference by A I 0. Finally, using the abbreviation (6) we arrive at the equation lo;llxillt : I j:l,j rti with wi1 l A xn ait xl ;I a'i xaa \./ j-th column, j:i,...,,
8 r33 which is equivalentto A;(x*) - 0 for eachi - 1,...,. Again, accoring to Cramer's rule the vector wi : (wn,..., tt)i)t sasfies the equation X/wi - (an,..., ai)',that is, we may write wi : (Xr)-t A'rr.t with ei - (0,..., 1 0)/ enoting the i-th column vector of the x ientity matrix. Hence, the above equations can be rewritten as Äo;llxillz: wi'x* -ei'anx-l x* for i-1,...,. However, these equations follow irectly from (8) an vice versa. Thus the proof of the Theorem is complete. I Note that the point Zi - wi* zl-xiwhich the orthogonality relation can be shown to belong to the sphere S; satisfies (x*-zi,x*-xi) -0 for i-1,..., This is easily verifie by straightforwar computations using the expressions of x*, wi, an zs given in the above proof. As a result, the sphere S; has the centre (zi + x) /2 - (wi + zü/2. 3 Bouns for eterminants In the subsequent lemma we establish a lower boun for the eterminant of An that is uniform in all the varying parameters )";y, 0 ( i < j <, an even positive provie the vectors Xl,..., x are linearly inepenent. Note that the below inequality (11) is vali for any (> 2) elements of an arbitrary Eucliean vector space. Lemma. Let V be a Eucliean vector space equippe with scalar prouct (.,.). For any Xl,...,X V ananyrealnumberslij,i, j :0,I,...,, satisfiing),i1 *L1i - lfor i, j :0, 1,...,, the inequality etan Z etbr * etc(x) > etc(x) (11) hols, where Al, B71, an G(X) are efine by (6), (9), an (7), respectively. Equality is attaine in (lt) if,1.;; llxi -*;ll2 - (1-ro;) llxillz +xo,;1x;ll2 - (x1,x;) fo, <i < j <. Proof. As alreay pointe out at the en of Section 1 the eterminant of any skewsymmetric x matrixbisnon-negativeanequals zerolf iso. Thus,theseconpart of (11) is trivial an instea of the first part we prove the slightly more general inequality et(b + G(X)) > etb * et G(X) (r2) By applying the well-known principal axis theorem to the non-negative efinite Gram matrix G(X) we fin an orthogonal x maftix O (with eto - 1) such that D - O'G(X) O is a iagonal matrix with non-negative iagonal elements. The multiplication rule for eterminants enables us to replace G(X) by D an B by the skew-symmetric matrix
9 134 L. Heinrich O'B O without changing the inequality (12). Therefore it suffices to verify et(b + D) > etb*etd only for iagonal matrices D with non-negative iagonal elements )l1,... z la an any skew-symmetric x matrix B. The Taylor expansion of the function f (yt,..., ya) - et(b + D) - lr bn bu -bn lz bza -bta -bza io leas to -2 f (Yr,...,Y) * etb + I lh'''!i* etbi1,...,ik k:l l<it<...<ip< * Yr ' ''Y, wherethe(-k)x(-k)matrixbir,...,iuemergesfrombbyeletingtherowsil,...,ik an columns i 1,..., ik.obviously, all these matrices are skew-symmetric which, together with y1,..., la > 0, implies f(yr,...,ya) > etb *yt"'l: etb*etd. Thus the inequality (11) is prove. The proof of the lemma is complete by noting that Al : A/n an (7) imply equality in (11). n Remark. The above proof turns out that the inequality (I2) remains vali if G(X) is replace by any other non-negative efinite x matrix. Corollary. If xt,..., X e V are linearly inepenent, then etal is positive which in turn implies the existence of the inverse Alr for all real )";1 satisfying ),ii * ).1i : I for i, j :0, 1,...,. As a special case, (II) inclues thewell-knownhaamarinequality etg(x) < llxr P..... llxall2 for any xr,..., x4 e V, whichreas lals llxrll... llxall in R. The first part of the corollary follows from the well-known fact that et G(X) > 0 characterizes the linear inepenence of xr,..., x V. The secon part follows by setting _ L1 :0an )rij : llxill2 /llxi-xjll2 sothatby(6) aji :0 for < i < j <. In other wors, Azr is an upper triangular matrix entailing etan - afi '... ' a. 4 Concluing remarks 1. Formula (8) an the above theorem reveal that the Miquel point M of the simplex S(o,x1,...,x)w.r.t.themarkepointsxij:xi*l;j(x:-xi)onthesimplex'egescan be represente as a linear combination of the eges Xl,..., X, **:I i:l uixi - Xu with weight vector a - (ut,..., ü)t - All xt (13)
10 135 From (6) an (8) it is seen that the weights ur,...,u4 epen only on the parameters )"ii anthesquareegelengthsllxi-x;ll2for0 < i < j < (withxo: o). For this reason the Miquel point can be efine in a meaningful way for any finite family of linearly inepenent elements of a Eucliean vector space V. By (13) the Miquel points *i,i - 0, I,...,, of those ( - l)-simplices having the vertices x:, i f i, can be easily etermine. Geometrically spoken,xi is the point of intersectionof themiquelspheres Sj, j t' i,withthehyperplanecontainingthe ( - l)-tace of the original simplex which x; oes not belong to. In this way a new - simplex S(xö,x{,...,*ä) arises an relations between it an the original one for given Ä;7's coul be of interest. 2. Since the matrix An is invertible for any )"ii the point x* is also well-efine for \ij e {0, 1}. For special choices of the,1.;y's we have interesting geometric intepretations, e.g. for:2,lettinglor + l,loz ) 0, an),p + l entailsthateachof thelimitingmiquel circles touches one sie of the triangle A x6x1x2 at a vertex an passes through the opposite vertex; the corresponing Miquel point turns into the Brocar point, see [1]. By means of the Theorem in Section 2 several generalizations of this point to higher imensions are possible. In particular, for a tetraheron S(o, xr,x2,x3) in IR3 the choice l0r : Loz : 103 : Ä13 : 1 an Ltz : \23 :0 (in the above setting) yiels x* as common point of the circumsphere 56 an 51 fl Sz fl 53, where e.g. Sr is the unique sphere through x3 touching the face triangle A ox1x2 atx1, an 52, 53 are eflne analogously. 3. The vertices of a further -simplex associate with,s(o, x1,..., xa) an the given,l.;;'s coinciewiththe mipointsxf of themiquel spheres S;, i :0, 1,...,.In Section 2 we have erive the followins formulas: *3: jtx')-t *, anxf -)r*'r-tan., **3, i : r,..., In the planar case simple geometric arguments show that the triangles A xsx1x2 an A xfixixi are similar, cf. e.g. [1]. It is natural to ask whether the simplex S(xfi, x!,..., *;) an the original simplex are similar for any > 2. This fact can be expresse analytically by an orthogonal matrix O an some scaling factor y > 0 by requiring *i -* lt*'l-tan., yox; for i - r,...,, where' 2 TetAn 1t7a \- ^r-l ln view of O/ - O-1 an (7) these relations are equivalento G-l(X) - 2y2 (n;t + (A'n)-t ). fnis equality hols actually only for. - 2 without aitional restrictions on the.1";;'s. 4. Onthe other han, for any > 2, the particular choice )"ij l/2, 0 < i < j <, yiels the Miquel point x* _ \.(x')-t (llxrll2,...,llxa;12;', which coincies with the circumcentre of the simplex,s(o, Xl,..., xa).
11 136 L. Heinrich Acknowlegements I wish to thank Prof. Dr. H. Martini (TU Chemnitz) who brought Konnully's paper l2l to my attention. I am grateful to Dr. L. Muche for computer simulations concerning Miquel's theorem. Many thanks also to Dr. T. Klein for his helpful comments on an earlier version of this paper an to stu. math. oec. G. Freuenreich who assiste me in the 3D visualization of Miquel's theorem. References [1] Donath, E.: Die merhuürigen Punkte un Linien es ebenen Dreiecks. VEB Deutscher Verlag er Wissenschaften, Berlin [2] Konnully, A.O.: Pivot theorems in n-space. Beitröge Algebra Geom. 10 (1980), [3] Miquel, A.: M6moire e G6om6trie, J. Math. Pures Appl. e Liouville 1 (1838), [4] Weisstein, E.W.: Miquel's Theorem. Lothar Heinrich University of Augsburg Institute of Mathematics D Augsburg, Germany he inri ch@math. uni - augsburg. e
On the Miquel point of simplices
Elem Math 63 (2008) 126 136 0013-6018/08/030126-11 c Swiss Mathematical Society, 2008 Elemente der Mathematik On the Miquel point of simplices Lothar Heinrich Lothar Heinrich promovierte im Jahr 1980 an
More informationLecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012
CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration
More informationProblem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs
Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable
More informationSturm-Liouville Theory
LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationExam 2 Review Solutions
Exam Review Solutions 1. True or False, an explain: (a) There exists a function f with continuous secon partial erivatives such that f x (x, y) = x + y f y = x y False. If the function has continuous secon
More informationAcute sets in Euclidean spaces
Acute sets in Eucliean spaces Viktor Harangi April, 011 Abstract A finite set H in R is calle an acute set if any angle etermine by three points of H is acute. We examine the maximal carinality α() of
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationRobust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k
A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationSYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is
SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. Uniqueness for solutions of ifferential equations. We consier the system of ifferential equations given by x = v( x), () t with a given initial conition
More informationTwo formulas for the Euler ϕ-function
Two formulas for the Euler ϕ-function Robert Frieman A multiplication formula for ϕ(n) The first formula we want to prove is the following: Theorem 1. If n 1 an n 2 are relatively prime positive integers,
More informationMath 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions
Math 3 Winter 2 Avance Bounary Value Problems I Bessel s Equation an Bessel Functions Department of Mathematical an Statistical Sciences University of Alberta Bessel s Equation an Bessel Functions We use
More informationInternational Journal of Pure and Applied Mathematics Volume 35 No , ON PYTHAGOREAN QUADRUPLES Edray Goins 1, Alain Togbé 2
International Journal of Pure an Applie Mathematics Volume 35 No. 3 007, 365-374 ON PYTHAGOREAN QUADRUPLES Eray Goins 1, Alain Togbé 1 Department of Mathematics Purue University 150 North University Street,
More informationSome vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10
Some vector algebra an the generalize chain rule Ross Bannister Data Assimilation Research Centre University of Reaing UK Last upate 10/06/10 1. Introuction an notation As we shall see in these notes the
More informationREAL ANALYSIS I HOMEWORK 5
REAL ANALYSIS I HOMEWORK 5 CİHAN BAHRAN The questions are from Stein an Shakarchi s text, Chapter 3. 1. Suppose ϕ is an integrable function on R with R ϕ(x)x = 1. Let K δ(x) = δ ϕ(x/δ), δ > 0. (a) Prove
More informationOn the enumeration of partitions with summands in arithmetic progression
AUSTRALASIAN JOURNAL OF COMBINATORICS Volume 8 (003), Pages 149 159 On the enumeration of partitions with summans in arithmetic progression M. A. Nyblom C. Evans Department of Mathematics an Statistics
More informationIterated Point-Line Configurations Grow Doubly-Exponentially
Iterate Point-Line Configurations Grow Doubly-Exponentially Joshua Cooper an Mark Walters July 9, 008 Abstract Begin with a set of four points in the real plane in general position. A to this collection
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationu!i = a T u = 0. Then S satisfies
Deterministic Conitions for Subspace Ientifiability from Incomplete Sampling Daniel L Pimentel-Alarcón, Nigel Boston, Robert D Nowak University of Wisconsin-Maison Abstract Consier an r-imensional subspace
More informationALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS
ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an
More informationApplications of the Wronskian to ordinary linear differential equations
Physics 116C Fall 2011 Applications of the Wronskian to orinary linear ifferential equations Consier a of n continuous functions y i (x) [i = 1,2,3,...,n], each of which is ifferentiable at least n times.
More informationCHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold
CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the
More informationPermanent vs. Determinant
Permanent vs. Determinant Frank Ban Introuction A major problem in theoretical computer science is the Permanent vs. Determinant problem. It asks: given an n by n matrix of ineterminates A = (a i,j ) an
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationSome functions and their derivatives
Chapter Some functions an their erivatives. Derivative of x n for integer n Recall, from eqn (.6), for y = f (x), Also recall that, for integer n, Hence, if y = x n then y x = lim δx 0 (a + b) n = a n
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationLinear Algebra- Review And Beyond. Lecture 3
Linear Algebra- Review An Beyon Lecture 3 This lecture gives a wie range of materials relate to matrix. Matrix is the core of linear algebra, an it s useful in many other fiels. 1 Matrix Matrix is the
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationNotes on Lie Groups, Lie algebras, and the Exponentiation Map Mitchell Faulk
Notes on Lie Groups, Lie algebras, an the Exponentiation Map Mitchell Faulk 1. Preliminaries. In these notes, we concern ourselves with special objects calle matrix Lie groups an their corresponing Lie
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationINRADII OF SIMPLICES
INRADII OF SIMPLICES ULRICH BETKE, MARTIN HENK, AND LYDIA TSINTSIFA Abstract. We stuy the following generalization of the inraius: For a convex boy K in the -imensional Eucliean space an a linear k-plane
More informationUnit vectors with non-negative inner products
Unit vectors with non-negative inner proucts Bos, A.; Seiel, J.J. Publishe: 01/01/1980 Document Version Publisher s PDF, also known as Version of Recor (inclues final page, issue an volume numbers) Please
More informationPure Further Mathematics 1. Revision Notes
Pure Further Mathematics Revision Notes June 20 2 FP JUNE 20 SDB Further Pure Complex Numbers... 3 Definitions an arithmetical operations... 3 Complex conjugate... 3 Properties... 3 Complex number plane,
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationarxiv: v2 [math.cv] 2 Mar 2018
The quaternionic Gauss-Lucas Theorem Riccaro Ghiloni Alessanro Perotti Department of Mathematics, University of Trento Via Sommarive 14, I-38123 Povo Trento, Italy riccaro.ghiloni@unitn.it, alessanro.perotti@unitn.it
More informationSolutions to Math 41 Second Exam November 4, 2010
Solutions to Math 41 Secon Exam November 4, 2010 1. (13 points) Differentiate, using the metho of your choice. (a) p(t) = ln(sec t + tan t) + log 2 (2 + t) (4 points) Using the rule for the erivative of
More informationThe Sokhotski-Plemelj Formula
hysics 25 Winter 208 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),
More informationLower bounds on Locality Sensitive Hashing
Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,
More informationMath 115 Section 018 Course Note
Course Note 1 General Functions Definition 1.1. A function is a rule that takes certain numbers as inputs an assigns to each a efinite output number. The set of all input numbers is calle the omain of
More informationTOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH
English NUMERICAL MATHEMATICS Vol14, No1 Series A Journal of Chinese Universities Feb 2005 TOEPLITZ AND POSITIVE SEMIDEFINITE COMPLETION PROBLEM FOR CYCLE GRAPH He Ming( Λ) Michael K Ng(Ξ ) Abstract We
More informationWitt#5: Around the integrality criterion 9.93 [version 1.1 (21 April 2013), not completed, not proofread]
Witt vectors. Part 1 Michiel Hazewinkel Sienotes by Darij Grinberg Witt#5: Aroun the integrality criterion 9.93 [version 1.1 21 April 2013, not complete, not proofrea In [1, section 9.93, Hazewinkel states
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More informationA new proof of the sharpness of the phase transition for Bernoulli percolation on Z d
A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition
More informationCharacterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations
Characterizing Real-Value Multivariate Complex Polynomials an Their Symmetric Tensor Representations Bo JIANG Zhening LI Shuzhong ZHANG December 31, 2014 Abstract In this paper we stuy multivariate polynomial
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationPhysics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1
Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems
More informationDiscrete Mathematics
Discrete Mathematics 309 (009) 86 869 Contents lists available at ScienceDirect Discrete Mathematics journal homepage: wwwelseviercom/locate/isc Profile vectors in the lattice of subspaces Dániel Gerbner
More informationOptimal Variable-Structure Control Tracking of Spacecraft Maneuvers
Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign
More information1 Lecture 20: Implicit differentiation
Lecture 20: Implicit ifferentiation. Outline The technique of implicit ifferentiation Tangent lines to a circle Derivatives of inverse functions by implicit ifferentiation Examples.2 Implicit ifferentiation
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationLarge Triangles in the d-dimensional Unit-Cube (Extended Abstract)
Large Triangles in the -Dimensional Unit-Cube Extene Abstract) Hanno Lefmann Fakultät für Informatik, TU Chemnitz, D-0907 Chemnitz, Germany lefmann@informatik.tu-chemnitz.e Abstract. We consier a variant
More informationSimilar Operators and a Functional Calculus for the First-Order Linear Differential Operator
Avances in Applie Mathematics, 9 47 999 Article ID aama.998.067, available online at http: www.iealibrary.com on Similar Operators an a Functional Calculus for the First-Orer Linear Differential Operator
More informationDECOMPOSITION OF POLYNOMIALS AND APPROXIMATE ROOTS
DECOMPOSITION OF POLYNOMIALS AND APPROXIMATE ROOTS ARNAUD BODIN Abstract. We state a kin of Eucliian ivision theorem: given a polynomial P (x) an a ivisor of the egree of P, there exist polynomials h(x),
More informationA note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz
A note on asymptotic formulae for one-imensional network flow problems Carlos F. Daganzo an Karen R. Smilowitz (to appear in Annals of Operations Research) Abstract This note evelops asymptotic formulae
More informationDiophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations
Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything
More informationLower Bounds for the Smoothed Number of Pareto optimal Solutions
Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.
More informationPhysics 251 Results for Matrix Exponentials Spring 2017
Physics 25 Results for Matrix Exponentials Spring 27. Properties of the Matrix Exponential Let A be a real or complex n n matrix. The exponential of A is efine via its Taylor series, e A A n = I + n!,
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationDIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10
DIFFERENTIAL GEOMETRY, LECTURE 15, JULY 10 5. Levi-Civita connection From now on we are intereste in connections on the tangent bunle T X of a Riemanninam manifol (X, g). Out main result will be a construction
More informationYear 11 Matrices Semester 2. Yuk
Year 11 Matrices Semester 2 Chapter 5A input/output Yuk 1 Chapter 5B Gaussian Elimination an Systems of Linear Equations This is an extension of solving simultaneous equations. What oes a System of Linear
More informationarxiv: v4 [cs.ds] 7 Mar 2014
Analysis of Agglomerative Clustering Marcel R. Ackermann Johannes Blömer Daniel Kuntze Christian Sohler arxiv:101.697v [cs.ds] 7 Mar 01 Abstract The iameter k-clustering problem is the problem of partitioning
More information1 Math 285 Homework Problem List for S2016
1 Math 85 Homework Problem List for S016 Note: solutions to Lawler Problems will appear after all of the Lecture Note Solutions. 1.1 Homework 1. Due Friay, April 8, 016 Look at from lecture note exercises:
More informationA Sketch of Menshikov s Theorem
A Sketch of Menshikov s Theorem Thomas Bao March 14, 2010 Abstract Let Λ be an infinite, locally finite oriente multi-graph with C Λ finite an strongly connecte, an let p
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More information4.2 First Differentiation Rules; Leibniz Notation
.. FIRST DIFFERENTIATION RULES; LEIBNIZ NOTATION 307. First Differentiation Rules; Leibniz Notation In this section we erive rules which let us quickly compute the erivative function f (x) for any polynomial
More informationQubit channels that achieve capacity with two states
Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March
More informationDeterminant and Trace
Determinant an Trace Area an mappings from the plane to itself: Recall that in the last set of notes we foun a linear mapping to take the unit square S = {, y } to any parallelogram P with one corner at
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationDifferentiation ( , 9.5)
Chapter 2 Differentiation (8.1 8.3, 9.5) 2.1 Rate of Change (8.2.1 5) Recall that the equation of a straight line can be written as y = mx + c, where m is the slope or graient of the line, an c is the
More informationProblem set 2: Solutions Math 207B, Winter 2016
Problem set : Solutions Math 07B, Winter 016 1. A particle of mass m with position x(t) at time t has potential energy V ( x) an kinetic energy T = 1 m x t. The action of the particle over times t t 1
More informationTractability results for weighted Banach spaces of smooth functions
Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March
More information1 Definition of the derivative
Math 20A - Calculus by Jon Rogawski Chapter 3 - Differentiation Prepare by Jason Gais Definition of the erivative Remark.. Recall our iscussion of tangent lines from way back. We now rephrase this in terms
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationA. Exclusive KL View of the MLE
A. Exclusive KL View of the MLE Lets assume a change-of-variable moel p Z z on the ranom variable Z R m, such as the one use in Dinh et al. 2017: z 0 p 0 z 0 an z = ψz 0, where ψ is an invertible function
More informationThe effect of dissipation on solutions of the complex KdV equation
Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,
More informationInterconnected Systems of Fliess Operators
Interconnecte Systems of Fliess Operators W. Steven Gray Yaqin Li Department of Electrical an Computer Engineering Ol Dominion University Norfolk, Virginia 23529 USA Abstract Given two analytic nonlinear
More informationMarkov Chains in Continuous Time
Chapter 23 Markov Chains in Continuous Time Previously we looke at Markov chains, where the transitions betweenstatesoccurreatspecifietime- steps. That it, we mae time (a continuous variable) avance in
More informationMath Chapter 2 Essentials of Calculus by James Stewart Prepared by Jason Gaddis
Math 231 - Chapter 2 Essentials of Calculus by James Stewart Prepare by Jason Gais Chapter 2 - Derivatives 21 - Derivatives an Rates of Change Definition A tangent to a curve is a line that intersects
More information3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes
Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationGeneralized Tractability for Multivariate Problems
Generalize Tractability for Multivariate Problems Part II: Linear Tensor Prouct Problems, Linear Information, an Unrestricte Tractability Michael Gnewuch Department of Computer Science, University of Kiel,
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationMARTIN HENK AND MAKOTO TAGAMI
LOWER BOUNDS ON THE COEFFICIENTS OF EHRHART POLYNOMIALS MARTIN HENK AND MAKOTO TAGAMI Abstract. We present lower bouns for the coefficients of Ehrhart polynomials of convex lattice polytopes in terms of
More informationLINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form
LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1 We consier ifferential equations of the form y + a()y = b(), (1) y( 0 ) = y 0, where a() an b() are functions. Observe that this class of equations inclues equations
More informationFebruary 21 Math 1190 sec. 63 Spring 2017
February 21 Math 1190 sec. 63 Spring 2017 Chapter 2: Derivatives Let s recall the efinitions an erivative rules we have so far: Let s assume that y = f (x) is a function with c in it s omain. The erivative
More informationd dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1
Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of
More informationMATH Topics in Applied Mathematics Lecture 12: Evaluation of determinants. Cross product.
MATH 311-504 Topics in Applied Mathematics Lecture 12: Evaluation of determinants. Cross product. Determinant is a scalar assigned to each square matrix. Notation. The determinant of a matrix A = (a ij
More information4. Important theorems in quantum mechanics
TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this
More informationLinear and quadratic approximation
Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function
More informationCalculus in the AP Physics C Course The Derivative
Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.
More informationThe Chain Rule. d dx x(t) dx. dt (t)
The Chain Rule The Problem You alreay routinely use the one imensional chain rule t f xt = f x xt x t t in oing computations like t sint2 = cost 2 2t In this example, fx = sinx an xt = t 2. We now generalize
More information