Immediate Calculation of Some Poisson Type Integrals Using Supermathematics Circular Ex-Centric Functions

IMMEDIATE CALCULATION OF SOME POISSON TYPE INTEGRALS USING SUPERMATHEMATICS CIRCULAR EX-CENTRIC FUNCTIONS 1 Florentin SMARANDACHE1 & Mircea Eugen SELARIU2 Chair of Department of Math & Sciences, University of New Mexico-Gallup, USA 2 Polytechnic University of Timişoara, Romania 0. ABSTRACT This article presents two methods, in parallel, of solving more complex integrals, among which is the Poisson’s integral, in order to emphasize the obvious advantages of a new method of integration, which uses the supermathematics circular ex-centric functions. We will specially analyze the possibilities of easy passing/changing of the supermathematics circular ex-centric functions of a centric variable α to the same functions of ex-centric variable θ. The angle α is the angle at the center point O(0,0), which represents the centric variable and θ is the angle at the ex-center E(k,ε), representing the ex-centric variable. These are the angles from which the points W1 and W2 are visible on the unity circle – resulted from the intersection of the unity/trigonometric circle with the revolving straight line d around the ex-centric E(k,ε) – from O and from E, respectively. KEYWORDS AND ABBREVIATIONS C - Centric, Circular, CC - Circular Centric, E - Ex-centric, EC - Ex-centric Circular, F - Function, H - Hyperbolic, PI - Poisson Integral, M - Mathematics, CM - Centric Mathematics, EM - Ex-centric Mathematics, SM – Supermathematics, FSM - F & SM, FSM_EC- FSM & EC, FSM_EH-FSM & EH. 1. INTRODUCTION The discovery of the ex-centric mathematics (EM), as a vast extension of the centric common/mathematics (CM), which together form the SuperMathematics (SM), allows new simpler approaches, for resolving more complex integrals, among which we present (11) the Poisson integral (PI) [1]. To emphasize the new integration method, we will present, in parallel, the classic method of solving, only for PI, presented in [1] and the new method which utilizes SM’s ex-centric circular functions (EC) [2], [3], [4] . The SM-EC functions, which will be in the center of our attention, are the radial ex-centric functions rex θ and Rex α and the derivatives ex-centric dex θ and Dex α, functions which are independent of the reference system selected. 1 The functions rex θ, of ex-centric variable θ, of the principal determination 1 and secondary 2, defined on the whole real axis for numeric ex-centricity k2 <1, and for k 2 > 1 exist only in the interval ℑ ∈ (θi, θf), in which θf,i = π + ε ± arcsin(1/k), αf,i = θf,i + βf,i are (1) rex1,2 θ = rex1,2 (θ, Ε(k,ε)) = − k . cos(θ − ε ) ± 1 − k 2 sin 2 (θ − η ) , where E(k, ε) is a pole, called ex-center, which divides the straight line d (d = d+ U d - ), revolving around this point, in the positive semi straight line d+, on which is situated the first principal determination rex1 θ , as function of ex-centric variable θ and, respectively, Rex α1 , of centric variable α of the function and in the negative semi straight line d -, on which is situated, along it, the second determination, secondary, of the function rex2 θ and Rex α2. The expressions of the same entities (1), as functions of centric variable α, which exist on the whole real axis, no matter which is the numeric ex-centricity k, are (2) Rex α1,2 = ± 1 + k 2 − 2k . cos(α 1, 2 − ε ) These functions represent, as Prof. Dr. Math. Octav Gheorghiu observed, the distance in plane, as oriented segments, in polar coordinates, between two points: the ex-center E(k, ε) and the intersection points W1,2 (1, α1,2) – between the straight line d and the unity circle CT[1,O(0,0)] with the center in the origin O of the system of coordinates axis, right Cartesian or polar reference point. For an E which is interior to the unity disc, the segment EW1 is situated on the positive direction of the semi straight line d +, being, in this case, positive, that is Rex α1 = rex1θ > 0, while the oriented segment EW2, positive oriented on the negative semi straight line is negative, that is Rex α2 = rex2 θ < 0, as it can be seen in the Figure 1. For k = ± 1 at α1 ∈ (0, 2π) ⇒ θ ∈ (0, π) and at α2 ∈ (0, 2π) ⇒ θ ∈ (π, 2π ). In other words, if, the straight line d rotates around E(k,ε) ⊂ C (1,O) with an angular speed Ω (θ = Ωt ), the points W1,2 rotate on the unity circle C(1,O) with a double angular speed ( α1,2 = 2Ωt) in half of the period and is stationary in the second half of the period (α1,2 = 0) taking turns in E (k, ε). If E is exterior to the unity disk, that is |k | > 1, then both determinations will be on the same semi straight line, being, successive, both positive and then, after the rotation of d of π, both negative, therefore are of the same sign, and this will make their product, in this case, positive, and, while in the precedent case, the product of the two determinations of the function was always negative (see Fig.1). We must observe also that at k > 1 and for α1,2 ∈ (0, 2π) ⇒ θ ∈ (θi , θf); the ex-centric variable θ diminishes the interval of existence of FSM-CE, between an initial value θi and a final one θf, with as much as the numeric growth of ex-centricity k. For k → ∞ the interval is reduced to a single point on the real axis R, for each determination. The results presented so far, will also be obtained from the relations that will follow. 2 W<1 W>1+ W<1+ ≡ W>2+ θ> k >+ θ< k <- E<- θ> α k <+ θ ε E<+ O(0,0 E>+ x A(R,0) k >- E>- W<2- ≡ W>2W<2 - W>1- Fig. 1 The representation of the functions rex1,2 θ and Rexα1,2 (EW1,2 ) on the unity circle. The dependency between the two variables is: (3) α1,2 (θ) = θ - β1,2 (θ) = θ m arcsin[e.sin (θ - ε)] and, respectively (4) θ ( α1,2 ) = α1,2 + β (α1,2) = α1,2 + arcsin( = α1,2 + arcsin( or (4’) where (5) k . sin(α 1, 2 − ε ) Re xα 1, 2 θ (α1,2 ) = α1,2 + arctan ( cosβ (α1,2 ) = k . sin(α 1, 2 − ε ) ± 1 + k 2 − 2k . cos(α 1, 2 − ε ) )= ) k . sin(α 1, 2 − ε ) 1 − k . cos(α 1, 2 − ε ) 1 − k . cos(α 1, 2 − ε ) Re xα 1, 2 and 3 ) = α1,2 + arctan ( sin β (α 1, 2 ) cos β (α 1, 2 ) , (6) sin β (α1,2) = k . sin(α 1, 2 − ε ) Re xα 1, 2 , and the derivative of d[β(α)]/dα is dβ (α ) k[cos(α − ε ) − k ] k[cos(α − ε ) − k ] = (7) = 2 dα Re x 2α 1, 2 1 + k − 2.k . cos(α − ε ) From (1), it results, without difficulty, that the sum, the difference, the product, and the ratio of the two determinations of the functions rex are: (8) + ⎧ ⎪ ∑ = rex1θ + rex2θ = −2k . cos(θ − ε ) ⎪ − 2. 2 ⎪∑ = rex1θ − rex 2θ = 2 1 − k sin (θ − ε ) ⎪ ⎧k 2 − 1 < 0, → k < 1 ⎪ ⎨ ∏ = rex 1θ .rex2θ = ⎪⎨ 0, ⇒ k = 1 ⎪ ⎪ ± (1 − k 2 ), → k > 1 ⎪ ⎩ ⎪ 2 ⎪ = rex2θ = dα 2 = 1 − k = d (rex2θ ) ⎪⎩U rex1θ dα 1 Re x 2α 1 d (rex1θ ) A function similarly useful, in this article, is the ex-centric derivative function of a centric variable α, for which the form of expression is invariable at the position of the ex-center E is: 1 − k . cos(α 1, 2 − ε ) d (α 1, 2 + β1, 2 ) dθ 1 =: , = = (9) Dex α1,2 = 2 dα 1, 2 dα dex1, 2θ 1 + k − 2k . cos(α 1, 2 − ε ) and k . cos(θ − ε ) (10) dex1,2 = 1− 1 − k 2 sin 2 (θ − ε ) and the nucleus of Poisson integral 1− k 2 1− k 2 dγ d (θ + β ) dβ = 1+ 2 (11) Nip α1,2 = = = = = dα 1, 2 dα 1, 2 dα 1, 2 1 + k 2 − 2k (α 1, 2 − ε ) Re x 2α 1, 2 = − Re xα 2 Re xα 1 2. THE INTEGRATION USING THE CLASSIC METHOD [1] The Poisson’s integral, with modified notations, in accordance with the supermathematics ex-centric circular functions (SM - EC), is π dα (12) PI (k, ε) = ∫ , 2 1 k 2 . k . cos( α ε ) + − − −π in which k ∈ ℜ and ε ∈ [−π , π ] are the parameters and, in the same time, the polar coordinates of the ex-center E. This is resolved in [1] as a simple integral which is dependent of a real parameter λ ≡ k, which will be further denoted as k, and representing, in EM, numeric excentricity k = e/R, the ratio between the real ex-centricity e and the circle radius R on which are placed the intersection points W1 and W2. The integral is simple, but the integration is quite laborious, as we will see later, and it will become indeed simple, only when passed from CM to EM with the utilization of the new supermathematics functions. 4 Classical Solution: The periodic real function 1 1 + k − 2k cos(α − ε ) is, as it is easily observed, the square of the radial ex-centric function of α (13) (14) f(α) = 2 f(α) = 1 / (Rex2 α), defined for any k ∈ ℜ − {±1} and ϕ ∈ [-π, π ] . Remark: Only one from the two determinations of the function Rex α1,2 is null (!) when E belongs to unity circle, that is /k/ = 1; the second determination having the expression which will be presented bellow. Based on the new knowledge from EM, now we can assert that the radial ex-centric function is defined also for k = ± 1. If k = +1, then (15) rex1,2 θ = - cos(θ-ε) ± 1 − sin 2 (θ − ε ) → rex1 θ = Rex α1 = 0 and rex2θ = Rex α2 = − 2cos(θ− ε) and, for k = −1, it results (16) rex1,2 θ = cos(θ-ε) ± 1 − sin 2 (θ − ε ) such that, now, rex1θ = Rex α1 = 2 cos (θ-ε) and rex 2θ = Rex α2 = 0, which it results and it can be seen, equally easy, also from the graphic. Because (17) Rex2 α = [k− ei(α - ε)].[k – e -i.(α - ε)] = [k−rad(α-ε)].[k – rad-(α- ε)] , in which, the radial centric functions [5], or in short, radial (denoted rad), equivalent to the exponential functions are unitary vectors, of symmetric directions, in relation to the straight line which contains the points O and E, therefore: (18) rad (α- ε) − rad[- (α - ε)] = 2 cos(α - ε) and (19) rad(α-ε) . rad-(α-ε) = rad (α − ε ) = 1, rad (α − ε ) in which (20) rad α = ei α is equivalent, in centric (for k = 0, when α1 = θ and α2 = θ + π) of functions rex θ and Rexα [ 5]. The function Rex2α (16) has the roots: (21) e ±( α - ε) = rad[ ±(α−ε )] which, for α = ε and also for α = ε −π, din (14) it results (22) k = ±1. 5 By introducing in PI the variable α’ = α + π, the change will lead to the integral: (23) PI (− k) = ∫ 2π 0 dα ' , 1 + k 2 + 2k . cos(α ' + ε ) in which the numeric ex-centricity changes the sign, that is k → − k , which is equivalent to the rotation of the ex-center E(k, ε) around the origin O (0,0) with π, on the circle with the radius k, that is ε → ± ( ε ± π ), or, yet, because of the inter-conversion properties of α with ε in the function cosine from ( 12 ), α → ± ( α ± π ). Suppose that k ≠ ± 1, the change of the variable α’ = α + π (24) z = ei(α’ + ε) , for which (25) dz / z = der (α ' + ε )dα ' i.rad (α ' + ε )dα ' = i dα’ = rad (α ' + ε ) rad (α ' + ε ) it will transform the segment [−π, + π] in the unity circumference, going in trigonometric positive sense ( sinistrorum / levogin). Then: (26) PI(k, α) = i ∫ dz C (1 + k ) z + (1 + z 2 ) k 2 in which m = k +1/k. The poly-functions f(z) from under the sign of ∫ = i k ∫ are z’ = − k and z” = − 1 / k with the C C dz , z + mz + 1 2 residues a’-1 = Rez[f(z), -k] = k / (k2 - 1) and a”-1 = Rez[f(z), -1 / k] = k / (1 – k2), such that a’-1 + a”-1 = 0. By applying the residues and semi-residues theorems, it results that for any angle ε ∈ [ −π, +π ] ⎧ 2π ⎪1 − k 2 , for k < 1; ⎪ (27) IP( k, ε) = ⎨ 0, for k = ±1 ⎪ 2π ⎪ , for k > 1 ⎩ k 2 −1 The zero value for lim PI (k ,ε ) k →±1 can be found choosing the contour Γ made of the circumferences C and γ (Fig. 2) , the last having the center in z” = 1 / k and the radius r < 1, from which we suppressed the interior portions of the reunion of the two circles. In this dz is null even when k → 1 (or –1), appearing as a principal conditions, the integral ∫ 2 Γ z − kz + 1 value in the Cauchy sense. We can then write: 6 ⎧ 2π , for ⇒ k = ℜ − {±1} ⎪ 2 PI(k, ε) = ⎨ 1 − k ⎪ 0, for ⇒ k = ±1 ⎩ (28) Fig. 3 The functions Rex1,2 and rex1,2 of k = ± 1 and ε = 0 Fig. 2 Explanatory Drawing The result (28), presented in [1], can also be established directly, knowing that from (14), for k = 1, rex1θ = 0, and for k = −1, rex2 θ = 2cos(α-ε) such that: 1 1 dα dα π ∫−π Re x 2α 2 = −∫π 4 cos 2 (α − ε ) = 4 tan(α − ε ) −π = 4 [tan(π − ε ) − tan(−π − ε )] = 0 π (28) π For k ≠ ± 1, will present the integral below. 3. THE INTEGRATION WITH THE HELP OF CIRCULAR EX-CENTRIC SUPERMATHEMATICS FUNCTIONS. Multiplying PI(k, ε) with (1−k2) / (1−k2) it results (29) 1 PI (k, ε) = 1− k 2 = for k < 1 and (30) 1 1− k 2 1− k 2 1 ∫−π Re x 2α1,2 dα = 1 − k 2 π 1 d (θ + β ) ∫−π dα dα = 1 − k 2 π 1 PI (k, ε) = 1− k 2 Re xα 2 ∫π Re xα π − 1 = ∫π d (θ + β ) π − − 2π , 1− k 2 7 Re x 2α ∫π − Re x α dα π − = = 1 1 θ +β 1− k 2 π −π = 2π , 1− k 2 for k > 1, in which we took in consideration the relation (9) and the sign of the functions Rex α1,2, for k < 1 and for k > 1, that is an ex-center interior or exterior to unity disk and of the relation, for k < 1. The relations between the integration limits, taking in account of the dependencies [2] (31) ⎧ α1 = θ − β ⎨ ⎩α 2 = θ + β + π knowing that β1 + β2 = π [3] are : (32) If α1 ∈ [−π, π] ⇒ then θ ∈ [ π − β1, π + β1 ] and their difference is +2π, and if (33) α2 ∈ [ −π, π ], then ⇒ θ ∈ [−2π −β, −β], as it can be seen also in Fig. 1, and their difference is −2π. (34) 1+2 dβ dα d (2β ) d (α + β + β ) d (θ + β ) 1− k 2 e[(cos(α − ε ) − e] = =1+2 = + = = 2 2 dα dα dα dα Re x α 1, 2 Re x α 1, 2 dα because θ = α + β, and for ⎧γ 1 = −π → θ = −π + 2 β 1, 2 (35) α= ⎨ , = → = + 2 γ π θ π β 2 1 , 2 ⎩ as it results also from the figure, therefore (36) γ2 - γ1 = 2 π CONCLUSIONS Because of the labor volume in the two variants, the conclusion is, evidently, in the favor of the new method of integration, taking in account, firstly, the degree of complexity of the integration. By utilizing the existing relations in EM, as, for example the relation (28), which can be written by denoting γ = θ + β , from which dγ = d(θ +β), but α = (θ - β ) and dα = d(γ - 2β) or dα = d(θ - β) , such that Re xα 2 k . cos(θ − ε ) 1− k 2 dγ / dα = 1 +2.dβ / dα = − =− = rex1θ Re xα 1 Re x 2α 1, 2 and PI is an immediate integral, k being a constant parameter, as we saw before. Furthermore, from the relation (29) it results the Poisson’s integral value undefined as being: k sin(αε ) α + 2 arcsin dα 1 1 Re xα = [θ (α ) + β (α )] = [α + 2β (α )] = (37) IPN = ∫ 2 2 Re x 2α 1 − k 2 1− k 1− k The integrals calculated in [1] with the help of the residues theorem 8 ∫R R − r cos(α − ε ) , + r 2 − 2r cos(α − ε ) 2π (38) I1 = 2 0 in which with r = k.R we denoted the real ex-centricity and with R the radius of a certain circle and 2π r sin(α − ε ) (39) I 2 = ∫ 2 2 0 R + r − 2r cos(α − ε ) which, by the classical method presented in [1, pp. 186-187] are equally laborious and, unfortunately, wrong; by the new method, from EM, these integrals are immediate. Reducing R from (38) and (39) it results the functions to be integrated: (40) F1 = dθ R − r cos(α − ε ) 1 − k cos(α − ε ) = = Dex α1,2 = 2 2 dα R + r − 2r cos(α − ε ) Re x α 2 such that the undefined integral is (41) I1N = dθ ∫ dα dα = ∫ dθ = θ (α 1, 2 ) = α + arcsin[ k sin(α 1, 2 − ε ) ± 1 + k 2 − 2k cos(α 1, 2 − ε ) ], such that, the defined integral (38) will be: ¾ For k = +1 ⇒ I1 = π, because in the first determination 1 (principal) θ(α = 0) = π/2 and θ(α=2π) = 3π/2 and the difference is π. If k = − 1 for the first determination θ(α=0) = π and θ(α = 2π) = 2π, such that the difference is the same π. It results that for |k| = 1 ⇒ I1 = π. ¾ For k > 1, the integral value I1 is 0. ¾ For k < 1, the integral value is 2π. The undefined integral I2N is: (42) I2N = ∫1+ k k . sin(α − ε ) dα = − 2k cos(α − ε ) 2 ∫ k . sin(α − ε ) dα = Re x 2α ∫ Re xα 1 d (Re xα ) dα = ln|Rexα dα Therefore, the defined integral I2 is: (43) I2 = ln Re xα 2π 0 = 0, for any k and ε, knowing that Rex0 = rex0 = Rex2π = rex2π. More integrals can be resolved immediately in this way without difficulties, if one knows the expressions of some supermathematics functions. More integrals are presented in [6]. 9 REFERENCES [1] [2] [3] [4] [5] [6] [7] Bădescu, Radu; Maican, C-tin - Integrale utilizate in mecanică, fizică, tehnică şi calculul lor - Ed. Tehnică, Bucharest, 1968, 696 pages. Selariu, Mircea Eugen – Funcţii Circulare Excentrice - Conf. Nat. Vibr. In Constr. de Maşini, Timişoara, pp 101–108. Selariu, Mircea Eugen – Funcţii Circulare Excentrice şi Extensia lor - Bul. Şt. şi Tehn. al IPTVT, Timişoara, Seria Mec., Tom 25 (39), Fasc. 1-1980, pp. 95–100. Selariu, Mircea Eugen – Supermatematica - Com. VII Conf. Internat. de Ing. Man. Şi Tehnologică, TEHNO’95, Timişoara, Vol.9, Matematica Aplicată, pp. 41–64. Selariu, Mircea Eugen – Forma Trigonometrică a Sumei şi a Diferenţei Numerelor Complexe-Com. VII Conf. Internaţ. de Ing. Man. şi Tehn. TEHNO’95, Vol. 9: Matematică Aplicată, pp. 65–72. Selariu, Mircea Eugen; Ajiduah, Crist.; Bozantan, Emil; Filipescu, Avram – Integralele Unor Funcţii Supermatematice - Com. VII Conf. Internaţ. de Ing. Man. şi Tehn. TEHNO’95, Timişoara, 1995, Vol. IX: Matem. Aplic. pp. 73–82. Smarandache, Florentin - A Triple Inequality with Series and Improper Integrals Bulletin of Pure and Applied Sciences, Vol. 25E, No. 1, pp. 215-217, 2006. 10