Numerical Calculation of Area of Elliptical Segments
Total Page:16
File Type:pdf, Size:1020Kb
mathematics Article Numerical Calculation of Area of Elliptical Segments Dusan Vallo * , Jozef Fulier and Lucia Rumanova Department of Mathematics, Faculty of Natural Sciences, Constantine the Philosopher University in Nitra, Tr. A. Hlinku 1, 949 74 Nitra, Slovakia; [email protected] (J.F.); [email protected] (L.R.) * Correspondence: [email protected]; Tel.: +421-37-6408-690 Received: 9 August 2020; Accepted: 21 August 2020; Published: 23 August 2020 Abstract: In this paper, we introduce the notion of an elliptical segment as some analogy of the circular segment and we focus on the problem of calculation of its area. Based on the analytical method, we derive the formulas, which can be used for the numerical approximation of the area of the given segment. Keywords: ellipse; segment; area; approximation 1. Introduction In this paper, we present the problem of calculation of the area of elliptical segment. The notion elliptical segment is a generalization of the term circular segment. In accordance with this analogy, we consider the elliptical segment as a set of points surrounded by a given elliptical arc and its chord. The figures in Figure1 are the segments illustrated as the shaded regions. Figure 1. Circular and elliptical segment. There exist a few formulas about how to calculate the area A of the given circular segment. Their proving is based on the idea of subtraction of the area of the circular sector and the area of the subsistent triangle portion [1]. In this concept, the term central angle is crucial, but it is undefined for the ellipse [2]. The Gauss–Green formula is used to determine an elliptical sector area. It leads to an algorithm for calculating the area of the elliptical segment [3]. Our approach to the problem of the area of the elliptical segment is different. First, we prefer analytical method. Based on the idea of definite integral, we tessellate the given segment by polygons (trapezoids) and, consecutively, we derive a formula approximately evaluated the reflected area. Finally, we present an interesting corollaries of the results related to the area of the given ellipse. Mathematics 2020, 8, 1411; doi:10.3390/math8091411 www.mdpi.com/journal/mathematics Mathematics 2020, 8, 1411 2 of 8 2. Elliptical Segment First, we consider an ellipse in its central equation with center in the origin of Cartesian coordinate system Oxy. The straight line, determined by the chord of the conic, is defined in its slope-intercept form [4]. In Cartesian coordinate system Oxy, the ellipse is given by the equation x2 y2 + = 1 (1) a2 b2 for 0 < b < a, a, b 2 R and let p be the line of the form y = kx + q (2) for k, q 2 R. The line p must be a secant line. We show that h p p i q 2 − a2k2 + b2, a2k2 + b2 . Generally, let us consider an arbitrary parallel line in the equation y = kx + qi, (3) where qi 2 R. Suppose that the line in Equation (3) intersects the ellipse in points Ui, Vi. Their coordinates must be derived by solving an equation x2 (kx + q )2 + i = 1 a2 b2 . 2 2 2 2 2 2 2 2 x a k + b + 2a kqix + a qi − b = 0 We evaluate the discriminant D 2 2 2 2 2 2 2 2 D = 4a b a k + b − qi = 4a b Di (4) 2 2 2 2 where Di = a k + b − qi . The points of intersection exist, if and only if Di ≥ 0. It implies that p 2 2 2 p 2 2 2 − a k + b ≤ qi ≤ a k + b . (5) The line p also determines the chord UV of the given ellipse, if and only if p p − a2k2 + b2 < q < a2k2 + b2, (6) as it was necessary to prove. It also holds true that • if Di > 0, then there exists an one-to-one correspondence between the pair of points Ui, Vi and the intercept qi. These points of intersection have the coordinates p −ka2q ± ab D x = i i , (Ui,Vi) a2k2 + b2 Mathematics 2020, 8, 1411 3 of 8 p ±kab D + b2q y = kx + q = ··· = i i , (Ui,Vi) (Ui,Vi) i a2k2 + b2 specifically p p −kq a2 + ab D q b2 + kab D U i i , i i , (7) i a2k2 + b2 a2k2 + b2 p p −kq a2 − ab D q b2 − kab D V i i , i i . (8) i a2k2 + b2 a2k2 + b2 p 2 2 2 • If Di p= 0, then qi = ± a k + b and it means that Ui ≡ Vi. The lines in the equation y = kx ± a2k2 + b2 are tangent lines to the given ellipse. p 2 2 2 We label tP the tangentp line in the equation y = kx + a k + b and tQ will be the tangent line in the equation y = kx − a2k2 + b2. The tangent lines tP, tQ demarcate the strip in which the ellipse lies. The condition (5) also determines that each of the lines in Equation (3) belongs to this strip. Later, we use the points of intersection as the vertices of polygons which approximately tessellate the given elliptical segment. Without loss of generality, in the next section, we consider that the elliptical segment lies between the lines tQ and the line p. 3. Arithmetic Sequence of Intercepts To determine the vertices of the polygons, we introduce the intercepts qi as an arithmetic sequence n fqigi=0 such that p 2 2 2 p 2 2 2 q0 = − a k + b < q1 < ··· < qi < ··· < qn ≤ a k + b , (9) p 2 2 2 qi = − a k + b + id, (10) for i = 0, 1, 2, ... , n, where n 2 N is fixed and d > 0 is a difference of the sequence. Its value is bounded too. We determine it as follows: The inequalities (9) result in p 2 2 2 p 2 2 2 qi = − a k + b + id ≤ a k + b and we derive p a2k2 + b2 0 < d ≤ 2 i for i = 1, 2, . , n. It implies that p a2k2 + b2 0 < d ≤ 2 . (11) n This result allows us to put p q + a2k2 + b2 d = , (12) n This setting of the difference d has a geometric meaning. The line segment on the axis y, bounded p by the points 0, − a2k2 + b2 , (0, q), is fragmented to n particular segments of the equal length d. 4. Approximation of Elliptical Segment First, we construct the set of the lines, each in the equation y = kx + qi for the arithmetical n sequence fqigi=0 with its difference d. We determine the points of intersection with the ellipse and draw the triangle U0U1V1 and the trapezoids UiViVi+1Ui+1 for i = 1, 2, . , n − 1. They approximate the elliptical segment. Mathematics 2020, 8, 1411 4 of 8 It holds true that q 2 2 2ab p p jU V j = x − x + y − y = ··· = 1 + k2 D . i i Ui Vi Ui Vi a2k2 + b2 i Similarly, we derive 2ab p 2p jU + V + j = 1 + k D + . i 1 i 1 a2k2 + b2 i 1 In this setting of the intercept qi we derive that for Di it holds true 2 2 2 2 p 2 2 2 2 2 Di = a k + b − qi = ··· = 2id a k + b − i d , (13) where the difference d is determined by (12). The area Ae of the elliptical segment is approximately equal to the area A(n,d) n−1 jUiVij + jUi+1Vi+1j A(n,d) = ∑ hi, i=0 2 where hi is a height of the trapezoid. If we label k = tan j (Figure2), then we derive d hi = d cos j = p , (14) 1 + k2 − abd n 1 p p A = D + D . (15) (n,d) 2 2 + 2 ∑ i i+1 a k b i=0 Example 1. Approximately calculate the area of a semi-ellipse with axes a = 3, b = 2 for n ≤ 10. Solution. The semi-ellipse is a special case of the elliptical segment for q = 0. If we apply it in the Formulas (10), (12), and (13), then we derive p a2k2 + b2 i i d = , D = a2k2 + b2 2 − n i n n and we derive s s ! ab n−1 i i i + 1 i + 1 A(n,d) = ∑ 2 − + 2 − . n i=0 n n n n It is evident that k is a free parameter. It has also a geometric explanation-the line p passing through the center of the ellipse divides the ellipse in half. Due to the simplicity of the calculation, we put k = 0 and we organize data in Table1. Similarly, the reader can also verify some results listed in Table2. Mathematics 2020, 8, 1411 5 of 8 Figure 2. The secant lines of the ellipse and the points of intersection. Table 1. Practical calculation for n = 4. p p p p Parameters i i · d qi Di Di+1 Di Di+1 Di + Di+1 a = 3, b = 2 0 0.00 −2.00 0.00 1.75 0.00 1.32 1.32 k = 0, q = 0 1 0.50 −1.50 1.75 3.00 1.32 1.73 3.05 n = 4 2 1.00 −1.00 3.00 3.75 1.73 1.94 3.67 d = 0.50 3 1.50 −0.50 3.75 4.00 1.94 2.00 3.94 4 2.00 0.00 4.00 2.00 ∑ 11.98 A(n,d) 8.99 The visualizations with the different partition number n are shown in Figure3. Table 2.