# The Dimensional Fluxes of the Hypercylindrical Function

Authors: Dusko Letic; Branko Davidovic
Journal: Abstract and Applied Analysis
 (2011)

---

## Abstract

We present the results of the theoretical researches of some characteristics of the generalized hyperspherical function with two degrees as independent dimensions. Primarily are given the answers on the quantification of the dimensional potentials (fluxes) of this function in the domain of the integer natural numbers. Beside them, we have got the solutions for some continual fluxes of the contour hypercylindrical (HY) functions on the basis. The symbolical evaluation and numerical verification of the series values and integrals are realized through the program packages Mathcad Professional and Mathematica.


---

## Body

## 1. Introduction

Hypercylinder function is the hypothetical function connected to multidimensional space. The most significant value of this function is in the fact that it originated [1] on the properties of the cylindrical entities: point, diameter, square, surface, and volume of cylinder. Another property is generalizing of these functions from discretion to continuum. It belongs to the group of special functions, so its testing is performed on the basis of known functions of these types: gamma (Γ), psi (ψ0), error function (erf), and the like.Definition 1.1.
The generalized hypercylindrical function is defined by equality [2]  (1.1)HY(k,n,r)=2πk-1rk+n-3Γ(k+1)Γ(k+n-2)Γ((k+1)/2)(k,n∈R,r∈N),
where Γ(z) is the gamma function.Thanks to the interpolation properties of the gamma function, we can analytically pass from the field of the natural integer values on the set of real and noninteger values with which there is concurrence of events both for its graphic interpretation and more concise mathematical analysis. It is developed on the basis of the two freedom degreesk and n, as special (vectors) dimensions, besides r radius, as implicitly included freedom degree for every hypercylinder (Figures 1 and 2). The dominant theorem that is set is the one that refers to the recurrent property of this function (when the height is h=2r).It implicitly includes that the left vectors (n=2,1,0,-1,-2,…) of the M[HY]kxn matrix columns in Figure 1 we get on the basis of the reverent vector (n=3) deduction, and the right vectors (n=4,5,6,7,8,…) on the basis of integrals on r radius [2](1.2)∂∂rHY(k,n,r)=HY(k,n-1,r),HY(k,n+1,r)=∫0rHY(k,n,r)dr.Figure 1
The functional submatrixM[HY(k,n,r)] for the freedom degrees k∈-4,-3,…,4 and n∈-1,0,…,7 and with six characteristics hypercylindrical functions (undef. are nondefined, most often singular values of this function).Figure 2
The surface graphic of the hypercylindrical function with the constant radiusr=1 and the selected fields of freedom degrees (-4<n≤5) and (-2<k≤5) (using the software Mathematica).To the development of the theory of the multidimensional objects, especially have contributed: Conway and Sloane [3], Gwak et al. [4], Hinton [5], Hocking and Young [6], Manning [7], Maunder [8], Neville [9], Von R. Rucker [10], Sommerville [11], Sun and Bowman [12], and the others, and to its testing, Ramanujan and Hardy [16, 17]. Today the researches of the hypercylindrical function are represented both in Euclid's and Riemann's geometry (molecular dynamics, neural networks, hypercylindrical black holes and the like).


## 2. Dimensional Potentials: The Fluxes of the HY Function

### 2.1. Vertical Dimensional Flux of the Hypercylindrical Function

The discrete dimensional potential orflux of the hypercylindrical function presents the total of all single functions in the (sub)matrix of this function that develops for the integer natural freedom degrees. Formally, flux can be quantification by twofold series that covers this area of HY function. The first step is to define the value of the infinite succession of functions ordered in columns (vectors) of the submatrix M[HY]kxn(k,n∈N). This is at the same time as well the definition of the vertical dimensional fluxes of HY function. The first value that is being calculated, refers to the fourth column (n=3) of the submatrix M[HY]kxn in Figure 1. In this case, it follows that the flux is equal
(2.1)∑k=0∞HY(k,3,r)=2π+2r+4r2+2πr3+83πr4+π2r4+⋯+ε(2πk-1rkΓ((k+1)/2)).
The analytical value of this series is
(2.2)∑k=0∞HY(k,3,r)=∑k=0∞2πk-1rkΓ((k+1)/2)=2(1π+reπr2erfc⁡(-rπ)).
In the paper are used three known error functions as follows: erf(z)-basic, erfc(z)-cumulative and erfi(z)-imaginary. When k values are even (0,2,4,…), in other words odd ones (1,3,5,…), the series can be divided as dichotomous, so we can now obtain two complementary series(2.3)∑k=0∞HY(k,3,1)=∑k=0,2,4,…∞HY(k,3,r)+∑k=1,3,5,…∞HY(k,3,r)=2(1π+reπr2erf⁡(rπ))+2reπr2.
Well, the result (2.3) can be presented in the from of series with even (k=2b) and odd members (k=2b+1) that complement (one another). In that sense follows(2.4)∑b=0∞2πb-1(2r)2bb!(2b)!=2(1π+reπr2erf⁡(rπ)),∑b=0∞2πbr2b+1b!=2reπr2.
On the basis of the solution (2.3), as the starting and reference (one), and applying the reference relations (1.2), we can obtain the series values for lower freedom degrees (n<3). So, we establish the connection of cylinder hypervolume (n=3) with its hypersurface (n=2). In that sense follows a new vector flux for n=2:(2.5)∂∂r∑k=0∞HY(k,3,r)=∑k=0∞HY(k,2,r),
so that
(2.6)∑k=0∞HY(k,2,r)=2[2r+eπr2erfc⁡(-rπ)(1+2πr2)].
Further, for the hypercylinder (n=1), we obtain a series on the basis of the derivative of the previous series, so that(2.7)∂∂r∑k=0∞HY(k,2,r)=∑k=0∞HY(k,1,r)=4[2(1+πr2)+πreπr2erfc⁡(-rπ)(3+2πr2)].
For the zero dimension n=0, the series value is as well located on the derivative basis, so it follows that
(2.8)∑k=0∞HY(k,0,r)=∑k=0∞2rk-3πk-1Γ(k+1)Γ(k-2)Γ((k+1)/2)=4π[r(10+4πr2)+eπr2erfc⁡(-rπ)[3+4πr2(3+πr2)]].
For greater freedom degree that n=3, series are found through an inverse operation, that is, by recurrent relation on the basis of integrating on the radius r. Well,
(2.9)∑k=0∞HY(k,n+1,r)=∫0r(∑k=0∞HY(k,n,r))dr,
and for the fourth dimension is valid the next integral form
(2.10)∑k=0∞HY(k,4,r)=∫0r(∑k=0∞HY(k,3,r))dr=∫0r2[1π+reπr2erf⁡(rπ)]dr=1π[eπr2erfc⁡(-rπ)-1].
For greater freedom degree than n=3, series are found through a recurrent relation (1.2), by integrating on radius(2.11)∑k=0∞HY(k,4,r)=∫0r(∑k=1,3,5,…∞HY(k,3,r))dr+∫0r(∑k=0,2,4,…∞HY(k,3,r))dr.
So, we come to the expression(2.12)∑k=0∞HY(k,4,r)=∫0r2reπr2dr+∫0r(∑b=0∞22b+1πb-1r2bb!(2b)!)dr.
Solving the integrals and using the relation on dichotomous series [2], we obtain the dimension flux of the fourth dimension
(2.13)∑k=0∞HY(k,4,r)=eπr2-1π+∑b=0∞22b+1πb-1r2b+1b!(2b+1)!.
For the freedom degree of n=5,series can be found through the integrating of the obtained solution (2.10)
(2.14)∑k=0∞HY(k,5,r)=∫0r1π[eπr2erfc⁡(-rπ)-1]dr.
The general solution is known, and it is
(2.15)∫ebz2erfc⁡(az)dz=πerf⁡(zb)2b-1bπ∑k=0ℵa2k+1Γ(k+1,-bz2)bkk!(2k+1)+C.
After settling the expression, we get the concrete solution
(2.16)∑k=0∞HY(k,5,r)=∑k=0∞2rk+2πk-1Γ(k+1)Γ(k+3)Γ((k+1)/2)=1π[erfi(rπ)2-r+1π∑k=0∞(-1)2k+12k+1(1-Γ(k+1,-πr2)k!)].
For the same freedom degree, the series can be found with the integrating of the dichotomous expression(2.17)∑k=0∞HY(k,5,r)=∫0reπr2-1πdr+∫0r(∑b=0∞22b+1πb-1r2b+1b!(2b+1)!)dr.
So that(2.18)∑k=0∞HY(k,5,r)=1π(erfi(rπ)2-r)+∑b=0∞22b+1πb-1r2(b+1)b!(2b+2)!.
The same result is obtained as well on the basis of the complementary dichotomous series
(2.19)∑k=0∞HY(k,5,r)=∑b=0∞2πbr2b+3(2b+2)(2b+3)b!+∑b=0∞22b+1πb-1r2(b+1)b!(2b+2)!.
Here, we use the imaginary error function that is equal to [14](2.20)erfi(z)=-i⋅erf⁡(i⋅z)=2π∑k=0∞z2k+1k!(2k+1).
For the last analysed vector flux, through similar procedures follow:
(2.21)∑k=0∞HY(k,6,r)=∫0r(∑b=0∞2πbr2b+3(2b+2)(2b+3)b!)dr+∫0r(∑b=0∞22b+1πb-1r2(b+1)b!(2b+2)!)dr.
It can be presented that the first subintegral member of the dichotomous series (2.19) is equal to
(2.22)∑b=0∞2πbr2b+3(2b+2)(2b+3)b!=1π(erfi(rπ)2-r).
So, its integral is
(2.23)∫0r1π(erfi(rπ)2-r)dr=12π[1-eπr2π+rerfi(rπ)-r2].
The dimensional flux for the sixth freedom degree is now
(2.24)∑k=0∞2rk+3πk-1Γ(k+1)Γ(k+4)Γ((k+1)/2)=12π[1-eπr2π+rerfi(rπ)-r2]+∑b=0∞22b+1πb-1r2b+3b!(2b+3)!.
For each freedom degree the recurrent relation would be formulated as
(2.25)∑k=0∞HY(k,n,r)=∫0r(∑k=0,2,…∞HY(k,n-1,r)+∑k=1,3,…∞HY(k,n-1,r))dr.
The similar formulation would be related as well to the recurrent relation of the type
(2.26)∑k=0∞HY(k,n,r)=∂∂r(∑k=0,2,…∞HY(k,n+1,r)+∑k=1,3,…∞HY(k,n+1,r)).
It can be supposed that the values of the vector fluxes are less and less with the increasing of the freedom degree n, so that the limit values are equal
(2.27)lim⁡n→∞∑k=0∞HY(k,n,r)=0,respectively,lim⁡n→∞∫0∞HY(k,n,r)dk=0.
The systematized numerical values of the discrete and continual fluxes, (for r=1), are given in Table 1.Table 1The freedom degreen∑k=0∞HY(k,n,1)∫0∞HY(k,n,1)dk046629.76368337472346629.88425089125415399.22768992989695399.16962705746822674.04323182871168674.09309651202785392.63527195113329292.399315367283533414.32373036585873914.03491997018803952.45719988167729042.299538951068193460.43798351716962270.382509391917891670.07531012467716440.060900692549048180.01191755293051970.008950280519384290.00170146817331110.0011952249074813100.00021823419520460.0001445005596724110.00002521961908790.0000158503684769120.00000264065257910.0000015847488146⋮⋮⋮306.610193523821e-292.6540266176391e-29⋮⋮⋮∞limn→∞∑k=0∞HY(k,n,r)=0limn→∞∫0∞HY(k,n,r)dk=0∑n∑k52812.977966365550—The dimensional fluxes can be studied as well for the complex part. So, for example with recurrence we get the series value for negative freedom degrees, and so forn=-1, follows [2](2.28)∑k=0∞HY(k,-1,r)=∑k=0∞2rk-4πk-1Γ(k+1)Γ(k-3)Γ((k+1)/2)=8π{2(1+2πr2)(4+πr2)+πreπr2erfc⁡(-rπ)[15+4πr2(5+πr2)]}.
For more lower freedom degree, that is n=-2, the flux is more complex and is defined by the next analytical value
(2.29)∑k=0∞HY(k,-2,r)=∂∂r∑k=0∞HY(k,-1,r)=∑k=0∞2rk-5πk-1Γ(k+1)Γ(k-4)Γ((k+1)/2)=8π2{2r(3+2πr2)(11+2πr2)+eπr2erfc⁡(-rπ)×[2πr2(2πr2(15+2πr2)+45)+15]}.
So, for example, the flux for radius r=1/2 is ∑k=0∞HY(k,-2,1/2)≈43770.014449969689. The other values, also can be established on the basis of the recurrent relations (2.25), in other words (2.26).


### 2.2. The Fluxes on the Basis of the Series of the Hypercylindrical Functional Matrix

The discrete dimensional fluxes can be calculated as well on the “horizontal line”, that is, adding functions values on the submatrix seriesM[HY]k,n (Figure 1). For example, through the series development for k=3 the flux would contain the next members:
(2.30)∑n=0∞HY(3,n,r)=12π+12πr+6πr2+2πr3+12πr4+⋯+ε(12πrnΓ(n+1)).
Some values of discrete and continual fluxes (for r=1) are given in Table 2.Table 2Freedom degreek∑n=0∞HY(k,n,r)∑n=0∞HY(k,n,1)∫0+∞HY(k,n,1)dn02er/π1.7305119588645301.84608438006859712er5.436563656918095.24809906024688728er21.746254627672322.46216193622815312πer102.476810672082885.44633790197124464π(er-1)345.4810604213614237.53292390754055120π2(er-r-1)850.6988994458283508.72931334440466384π2[2er-(r+1)2-1]1654.544866434988892.20643559127927∑n=0∞(1680π3rn+4/Γ(n+5))2688.6618988873221338.1438558916638∑n=0∞(12288π3rn+5/Γ(n+6))3790.4276572623341768.5863448966209∑n=0∞(30240π4rn+6/Γ(n+7))4757.7028081498232104.22878227066410∑n=0∞(245760π4rn+7/Γ(n+8))5416.8054631352562289.34360590283511∑n=0∞(665280π5rn+8/Γ(n+9))5672.0250477771182304.87915812710812∑n=0∞(5898240π5rn+9/Γ(n+10))5520.7364051223062167.38351118232513∑n=0∞(17297280π6rn+10/Γ(n+11))5036.8154770260341917.751163982797⋮⋮⋮⋮50∑n=0∞HY(50,n,r)0.0000001178268380.000000029841827⋮⋮⋮⋮∞limk→∞∑n=0∞HY(k,n,r)=00limk→∞∫0∞HY(k,n,r)dn=0—∑k∑n52812.977966365550—


### 2.3. Some Continual Fluxes of the Hypercylindrical Function

The trend of the distribution of the vector fluxes is increasing, and then asymptotically falling, with the linear growth of the freedom degreen.From the standpoint of the functional analysis the most interesting series of the matrix M[HY]k,n are the ones referring to the freedom degrees k=2 and k=3. The first series (k=2) covers the known functions for the square size (8r) and surface (4r2).The members of the following series are, among the others, the cylinder functions of the surface (6πr2) and volume (2πr3) (Figure 1). The same series are interesting as well for calculatiing continual fluxes. So the continual natural flux for the hypercylinder surface is analysed in view of integrals, instead of series. This integral is specific, because its subintegral function is the reciprocal gamma function. Its value, as it is known, is equal to the value of Fransen-Robinson's constant [15]
(2.31)F=∫0∞1Γ(x)dx=e+∫0∞e-nπ2+ln⁡2ndn≈2.8077702420285.
The integral values of the wanted flux is now ∫0∞HY(2,n,1)dn=∑n=0∞HY(2,n,1)+∫0∞(8e-n/(π2+ln2n))dn, what concretely for the unit radius is
(2.32)∫0∞8Γ(n)dn=8(e+∫0∞e-nπ2+ln⁡2ndn)=8F≈22.462161936228.
In respect to the continual dimension n, the more generalized cylinder volume flux follows in view of Ramanujan-Hardy's integral [16, 17]:(2.33)∫0∞yxΓ(x+1)dx=ey-∫0∞e-xyx(π2+ln⁡2x)dx.
Ramanujan defined this integral, and it was analytically intensified by Hardy. In that sense, the previous expression can be applied in calculating the hypercylindrical function flux, when k=3, as(2.34)∫0∞HY(3,n,1)dn=∫0∞12πrnΓ(n+1)dn=12π(er-∫0∞e-nrn(π2+ln⁡2n)dn)(forx=n,y=r).
The integral can be calculated as well as the difference between the series and the integral with the value (for r=1)(2.35)∫0∞12πrnΓ(n+1)dn=∑n=0∞12πrnΓ(n+1)-∫0∞12πe-nrn(π2+ln⁡2n)dn|r=1≈278.6085228836.


### 2.4. The Progressions of the Vector Fluxes

The whole dimension flux in the freedom degree domain with the natural numbers, is obtained in the result of double amount with which are considered the integer values of the hypercylindrical functionHY(k,n,r), for all k,n,r≥0. This double series must be convergent, and this characteristic is in the function of hypercylinder radius. The flux can be watched as well for every column M[HY]k,n of the matrix individually. So, there is for the nth column (in the mark 〈n〉), the flux in the form of series
(2.36)ΦHY〈n〉(k,n,r)=∑k=0∞HY(k,n,r).


### 2.5. The Orthogonal Dimensional Fluxes

These fluxes are the fluxes of the all columns or of all series of the matrixM[HY]k,n. As the number of columns, that is of the series infinite, the total flux is as follows.Definition 2.1.
The dimensional flux (potentional) of the functional matrix with two freedom degreesk and n is defined as a double series(2.37)ΦHY(k,n,r)=∑n=0∞ΦHY〈n〉(k,n,r)=∑n=0∞∑k=0∞HY(k,n,r).
When the certain number of members is calculated, the flux has the next form:(2.38)ΦHY(k,n,r)=∑n=0∞∑k=0∞HY(k,n,r)=4π[r(10+4πr2)+eπr2erfc⁡(-rπ)[3+4πr2(3+πr2)]]+4[2(1+πr2)+πreπr2erfc⁡(-rπ)(3+2πr2)]+2[2r+eπr2erfc⁡(-rπ)(1+2πr2)]+2(1π+reπr2erfc⁡(-rπ))+1π[eπr2erfc⁡(-rπ)-1]+1π(erfi(rπ)2-r)+1π2∑k=0∞(-1)2k+12k+1(1-Γ(k+1,-πr2)k!)+12π[1-eπr2π+rerfi(rπ)-r2]+⋯.
The flux on the matrix series in the domain of the natural number is defined as a double series, but with the changed sequence of summing. This dimensional flux is, thus, defined as(2.39)ΩHY(k,n,r)=∑k=0∞∑n=0∞HY(k,n,r).
In view of previously established members, the matrix flux has the following form:(2.40)ΩHY(k,n,r)=∑k=0∞∑n=0∞HY(k,n,r)=2erπ+2er+8er+12πer+64π(er-1)+120π2(er-r-1)+768π2(er+r-r22+1)+∑n=0∞1680π3rn+4Γ(n+5)+∑n=0∞12288π3rn+5Γ(n+6)+⋯.
The equivalence of the orthogonal dimensional fluxes means the equality of the double series(2.41)ΦHY(k,n,r)=ΩHY(k,n,r).
So, for example, for r=1 the dimensional fluxes have the unambiguous numerical value that is:
(2.42)ΦHY(k,n,1)=ΩHY(k,n,1)=52812.977966365550.


### 2.6. The Application of the Recurrent Operators with Defining the Diagonal Dimensional Fluxes of theHY Function

In the previous analyses the defining of the dimensional fluxes of the matrixM[HY]k,n is performed in view of the adding up of the HY function value for (on) the columns, in other words, the series of this matrix. The more detailed analysis would be very ample including the exponential function, then error functions erf(z), erfc(z), erfi(z), the incomplete gamma function Γ(a,z) and the like. When we use the ideas of transition operators from the reference function into the destined HY function, in the functional hypercylindrical matrix M[HY]k,n we can also establish the dimensional fluxes values on the diagonals (2.44), whose total would present the whole flux for the matrix in the field of freedom degrees for natural numbers, that is, k,n∈N. Such matrix contains infinite number of elements. For reference functions we take HY functions on the positions of the first series of the matrix, and they are so-called zero HY functions: HY(0,0,r),HY(0,1,r),…,HY(0,n,r),… The destined functions are arranged on the “gradual” growth law (+Δk) and decrease (-Δn) of increments.Definition 2.2.
The flux series operatorϑ(Δk,Δn,0) is defined by the relation
(2.43)ϑ(Δk,Δn,0)=HY(k+Δk,n+Δn,r)HY(k,n,r)=πΔk/2rΔk+ΔnΓ(k+n-2)Γ(k+Δk+1)Γ(k+1)Γ(k+n+Δk+Δn-2)Γ((k+Δk+1)/2)Γ(k+12).(2.44)Meanwhile, as increments for the absolute value are mutually equal and unit ones, that is,|Δk|=|-Δn|=1, to them is assigned a new joint argument u (Δk=Δn=u). In addition, the starting value of the k-og freedom degree states k=0, so the theta operator gets the form(2.45)θ(u,-u,0)=2uπuΓ(u2+1).
The destined function is now calculated as(2.46)HY(u,n-u,r)=θ(u,-u,0)⋅HY(0,n,r)=2u+1πu-2rn-3Γ(n-2)Γ(u2+1).
The dimensional flux on the diagonal presents the total of ist particular nembers. So, for the first diagonal (with the mark 〈0〉) flux is equal
(2.47)Π〈0〉(k,n,r)=HY(0,0,r)=0,
then for the second,
(2.48)Π〈1〉(k,n,r)=HY(0,1,r)+HY(1,0,r)=0.
For the third,
(2.49)Π〈2〉(k,n,r)=HY(0,2,r)+HY(1,1,r)+HY(2,0,r)=0.
Or for the fourth,
(2.50)Π〈3〉(k,n,r)=HY(0,3,r)+HY(1,2,r)+HY(2,1,r)+HY(3,0,r)=10+2π+12π.
The flux in the nth diagonal would be reckoned in the sum form(2.51)Π〈n〉(k,n,r)=∑u=0nHY(u,n-u,r)=rn-3Γ(n-2)∑u=0n2u+1πu-2Γ(u2+1)(n≠2).
The flux for the value n=2 is calculated on the basis of the highest function value. Considering that the equivalence of the double factorial and trigonometric functions [18] is known
(2.52)u!!=2(1/4)[1+2u-cos⁡(πu)]π(1/4)[cos⁡(πu)-1]Γ(u2+1).
The expression (2.19) can after settling be written as well in the equivalent form
(2.53)Π〈n〉(k,n,r)=rn-3(n-3)!∑u=0nu!!(2π)2u(2π)3+cos⁡(πu)4.
So, we get for the fifth diagonal (n=4), using the expressions (2.51) or (2.53)(2.54)Π〈4〉(k,n,r)=∑u=0nHY(u,n-u,r)=2rπ(5π+38π2+1).
For the sixth diagonal (n=5) follows:(2.55)Π〈5〉(k,n,r)=∑u=0n=5HY(u,n-u,r)=r2π(5π+38π2+60π3+1).
As the number of diagonals is infinite, the total flux is formed as the series of all diagonal fluxes
(2.56)ΠHY(k,n,r)=∑n=0∞Π〈n〉(k,n,r),orconcretely,ΠHY(k,n,r)=∑n=0∞∑u=0nHY(u,n-u,r)=∑n=0∞(rn-3Γ(n-2)∑u=0n2u+1πu-1Γ(u+12)).
For example, approximately, the flux for r=1 and n=11 it concretely is(2.57)ΠHY(k,11,1)=10960120160π+1096014032+143729910080π+4264371680π2+58497280π3+290128π4+332π5.
In the developed, form the total flux has the polynominal structure of members(2.58)ΠHY(k,n,r)=10+2π+12π+r(10+2π+76π)+r2(5+1π+38π+60π2)+r3(53+13π+383π+148π2)+r4(512+112π+37π2+70π3)+r5(112+160π+1930π+375π2+5825π3+112)+r6(172+1360π+19180π+3730π2+975π3+42π4)+⋯.
The diagonal flux of the hypercylindrical function can be expressed by the series of general form
(2.59)ΠHY(v,r)=∑v=0∞avrv.
Here, r is the summation index with which we take into consideration the order of elements from left to right and from top to bottom on (along) the diagonal (Figure 3). The polynomial coefficients contain the rational numbers and the graded constant π (Table 3). The first three coefficients are equal to zero, so they are not put in the summation sequence. Its other values are (v=0,1,…,12) given in Table 3.Table 3vav- polynomial coefficient ∑vavrv010+2π+12π12π+10+76π21π+5+38π+60π2313π+53+38π3+148π24112π+512+19π6+37π2+70π35160π+112+19π30+37π25+582π3561360π+172+19π180+37π230+97π35+42π4712520π+1504+19π1260+37π2210+97π335+1150π4218120160π+14032+19π10080+37π21680+97π3280+575π484+33π5591181440π+136288+19π90720+37π215120+97π32520+575π4756+2279π51261011814400π+1362880+19π907200+37π2151200+97π325200+115π41512+2279π51260+143π63011119958400π+13991680+19π9979200+37π21663200+97π3277200+115π416632+2279π513860+905π619812147900160+1239500800π+19π119750400+37π219958400+97π33326400+115π4199584+2279π5166320+905π62376+13π712Figure 3
The representative submatrix (above) and the addition principle of the matrix members on the diagonal.The series approximation with 13 coefficients in the decimal notation is in the form(2.60)∑v=012avrv≈48,33+249,4r+716,87r2+1502,27r3+2546,01r4+3684,24r5+4705,22r6+5422,03r7+5727,08r8+5610,38r9+5143,66r10+4445,23r11+3642,42r12.
Approximately, the double series leads us to the solution that is very near to the correct one. Namely, for the unit radius and reducing on ∞~n=30 the double series of the diagonal flux develops the following structure:(2.61)ΠHY(u,n,1)≈∑n=0n=30∑u=0nHY(u,n-u,r)=73997539898837593289987313727222173626045880401920000+739975398988375932899873137136110868130229402009600000π-1+5107360421374632558310731373581864950795510579200000π+23851430587700445181187313711342572344185783500800000π2+3950250824632679929196942891890428724030963916800000π3+103021449478899595399537986310638624850739200π4+345714978541889974399922394521436201548195840000π5+2633045028485350315237270061246290137702400π6+70440913248925148230273375765578626721280000π7+996435517964121615726791790306561280000π8+3543877234517627836251417738197632000π9+1285336197080807170493211041753600π10+1402624867171115786408429792000π11+5134803604675464525994928000π12+3516397543390323285352000π13+101418347130107761784000π14≈52810,5103716.
The diagonal dimension flux is characteristic with its coefficients that contain the πn constant in the graded series numbers, as distinguished from the vertical fluxes with the destination of the errors functions and π and e constants. The horizontal fluxes, as it is presented (2.40) contain the exponential functions. Meanwhile, the total flux for the unit radius is convergent and it can be calculated with considerably greater value. The value aberration (2.61) of the accurate value is only (or 0,00467%)
(2.62)ΔΠ=52812,97796636-52810,5103716≈2,467595.
The whole dimensional continual flux (k,n∈N) of the unit hypercylindrical function HY(k,n,1) is equal to the value of the twofold integral(2.63)∬0∞2πk-1Γ(k+1)Γ(k+n-2)Γ((k+1)/2)dkdn,
and we find its solution on the analytical and numerical bases.


## 2.1. Vertical Dimensional Flux of the Hypercylindrical Function

The discrete dimensional potential orflux of the hypercylindrical function presents the total of all single functions in the (sub)matrix of this function that develops for the integer natural freedom degrees. Formally, flux can be quantification by twofold series that covers this area of HY function. The first step is to define the value of the infinite succession of functions ordered in columns (vectors) of the submatrix M[HY]kxn(k,n∈N). This is at the same time as well the definition of the vertical dimensional fluxes of HY function. The first value that is being calculated, refers to the fourth column (n=3) of the submatrix M[HY]kxn in Figure 1. In this case, it follows that the flux is equal
(2.1)∑k=0∞HY(k,3,r)=2π+2r+4r2+2πr3+83πr4+π2r4+⋯+ε(2πk-1rkΓ((k+1)/2)).
The analytical value of this series is
(2.2)∑k=0∞HY(k,3,r)=∑k=0∞2πk-1rkΓ((k+1)/2)=2(1π+reπr2erfc⁡(-rπ)).
In the paper are used three known error functions as follows: erf(z)-basic, erfc(z)-cumulative and erfi(z)-imaginary. When k values are even (0,2,4,…), in other words odd ones (1,3,5,…), the series can be divided as dichotomous, so we can now obtain two complementary series(2.3)∑k=0∞HY(k,3,1)=∑k=0,2,4,…∞HY(k,3,r)+∑k=1,3,5,…∞HY(k,3,r)=2(1π+reπr2erf⁡(rπ))+2reπr2.
Well, the result (2.3) can be presented in the from of series with even (k=2b) and odd members (k=2b+1) that complement (one another). In that sense follows(2.4)∑b=0∞2πb-1(2r)2bb!(2b)!=2(1π+reπr2erf⁡(rπ)),∑b=0∞2πbr2b+1b!=2reπr2.
On the basis of the solution (2.3), as the starting and reference (one), and applying the reference relations (1.2), we can obtain the series values for lower freedom degrees (n<3). So, we establish the connection of cylinder hypervolume (n=3) with its hypersurface (n=2). In that sense follows a new vector flux for n=2:(2.5)∂∂r∑k=0∞HY(k,3,r)=∑k=0∞HY(k,2,r),
so that
(2.6)∑k=0∞HY(k,2,r)=2[2r+eπr2erfc⁡(-rπ)(1+2πr2)].
Further, for the hypercylinder (n=1), we obtain a series on the basis of the derivative of the previous series, so that(2.7)∂∂r∑k=0∞HY(k,2,r)=∑k=0∞HY(k,1,r)=4[2(1+πr2)+πreπr2erfc⁡(-rπ)(3+2πr2)].
For the zero dimension n=0, the series value is as well located on the derivative basis, so it follows that
(2.8)∑k=0∞HY(k,0,r)=∑k=0∞2rk-3πk-1Γ(k+1)Γ(k-2)Γ((k+1)/2)=4π[r(10+4πr2)+eπr2erfc⁡(-rπ)[3+4πr2(3+πr2)]].
For greater freedom degree that n=3, series are found through an inverse operation, that is, by recurrent relation on the basis of integrating on the radius r. Well,
(2.9)∑k=0∞HY(k,n+1,r)=∫0r(∑k=0∞HY(k,n,r))dr,
and for the fourth dimension is valid the next integral form
(2.10)∑k=0∞HY(k,4,r)=∫0r(∑k=0∞HY(k,3,r))dr=∫0r2[1π+reπr2erf⁡(rπ)]dr=1π[eπr2erfc⁡(-rπ)-1].
For greater freedom degree than n=3, series are found through a recurrent relation (1.2), by integrating on radius(2.11)∑k=0∞HY(k,4,r)=∫0r(∑k=1,3,5,…∞HY(k,3,r))dr+∫0r(∑k=0,2,4,…∞HY(k,3,r))dr.
So, we come to the expression(2.12)∑k=0∞HY(k,4,r)=∫0r2reπr2dr+∫0r(∑b=0∞22b+1πb-1r2bb!(2b)!)dr.
Solving the integrals and using the relation on dichotomous series [2], we obtain the dimension flux of the fourth dimension
(2.13)∑k=0∞HY(k,4,r)=eπr2-1π+∑b=0∞22b+1πb-1r2b+1b!(2b+1)!.
For the freedom degree of n=5,series can be found through the integrating of the obtained solution (2.10)
(2.14)∑k=0∞HY(k,5,r)=∫0r1π[eπr2erfc⁡(-rπ)-1]dr.
The general solution is known, and it is
(2.15)∫ebz2erfc⁡(az)dz=πerf⁡(zb)2b-1bπ∑k=0ℵa2k+1Γ(k+1,-bz2)bkk!(2k+1)+C.
After settling the expression, we get the concrete solution
(2.16)∑k=0∞HY(k,5,r)=∑k=0∞2rk+2πk-1Γ(k+1)Γ(k+3)Γ((k+1)/2)=1π[erfi(rπ)2-r+1π∑k=0∞(-1)2k+12k+1(1-Γ(k+1,-πr2)k!)].
For the same freedom degree, the series can be found with the integrating of the dichotomous expression(2.17)∑k=0∞HY(k,5,r)=∫0reπr2-1πdr+∫0r(∑b=0∞22b+1πb-1r2b+1b!(2b+1)!)dr.
So that(2.18)∑k=0∞HY(k,5,r)=1π(erfi(rπ)2-r)+∑b=0∞22b+1πb-1r2(b+1)b!(2b+2)!.
The same result is obtained as well on the basis of the complementary dichotomous series
(2.19)∑k=0∞HY(k,5,r)=∑b=0∞2πbr2b+3(2b+2)(2b+3)b!+∑b=0∞22b+1πb-1r2(b+1)b!(2b+2)!.
Here, we use the imaginary error function that is equal to [14](2.20)erfi(z)=-i⋅erf⁡(i⋅z)=2π∑k=0∞z2k+1k!(2k+1).
For the last analysed vector flux, through similar procedures follow:
(2.21)∑k=0∞HY(k,6,r)=∫0r(∑b=0∞2πbr2b+3(2b+2)(2b+3)b!)dr+∫0r(∑b=0∞22b+1πb-1r2(b+1)b!(2b+2)!)dr.
It can be presented that the first subintegral member of the dichotomous series (2.19) is equal to
(2.22)∑b=0∞2πbr2b+3(2b+2)(2b+3)b!=1π(erfi(rπ)2-r).
So, its integral is
(2.23)∫0r1π(erfi(rπ)2-r)dr=12π[1-eπr2π+rerfi(rπ)-r2].
The dimensional flux for the sixth freedom degree is now
(2.24)∑k=0∞2rk+3πk-1Γ(k+1)Γ(k+4)Γ((k+1)/2)=12π[1-eπr2π+rerfi(rπ)-r2]+∑b=0∞22b+1πb-1r2b+3b!(2b+3)!.
For each freedom degree the recurrent relation would be formulated as
(2.25)∑k=0∞HY(k,n,r)=∫0r(∑k=0,2,…∞HY(k,n-1,r)+∑k=1,3,…∞HY(k,n-1,r))dr.
The similar formulation would be related as well to the recurrent relation of the type
(2.26)∑k=0∞HY(k,n,r)=∂∂r(∑k=0,2,…∞HY(k,n+1,r)+∑k=1,3,…∞HY(k,n+1,r)).
It can be supposed that the values of the vector fluxes are less and less with the increasing of the freedom degree n, so that the limit values are equal
(2.27)lim⁡n→∞∑k=0∞HY(k,n,r)=0,respectively,lim⁡n→∞∫0∞HY(k,n,r)dk=0.
The systematized numerical values of the discrete and continual fluxes, (for r=1), are given in Table 1.Table 1The freedom degreen∑k=0∞HY(k,n,1)∫0∞HY(k,n,1)dk046629.76368337472346629.88425089125415399.22768992989695399.16962705746822674.04323182871168674.09309651202785392.63527195113329292.399315367283533414.32373036585873914.03491997018803952.45719988167729042.299538951068193460.43798351716962270.382509391917891670.07531012467716440.060900692549048180.01191755293051970.008950280519384290.00170146817331110.0011952249074813100.00021823419520460.0001445005596724110.00002521961908790.0000158503684769120.00000264065257910.0000015847488146⋮⋮⋮306.610193523821e-292.6540266176391e-29⋮⋮⋮∞limn→∞∑k=0∞HY(k,n,r)=0limn→∞∫0∞HY(k,n,r)dk=0∑n∑k52812.977966365550—The dimensional fluxes can be studied as well for the complex part. So, for example with recurrence we get the series value for negative freedom degrees, and so forn=-1, follows [2](2.28)∑k=0∞HY(k,-1,r)=∑k=0∞2rk-4πk-1Γ(k+1)Γ(k-3)Γ((k+1)/2)=8π{2(1+2πr2)(4+πr2)+πreπr2erfc⁡(-rπ)[15+4πr2(5+πr2)]}.
For more lower freedom degree, that is n=-2, the flux is more complex and is defined by the next analytical value
(2.29)∑k=0∞HY(k,-2,r)=∂∂r∑k=0∞HY(k,-1,r)=∑k=0∞2rk-5πk-1Γ(k+1)Γ(k-4)Γ((k+1)/2)=8π2{2r(3+2πr2)(11+2πr2)+eπr2erfc⁡(-rπ)×[2πr2(2πr2(15+2πr2)+45)+15]}.
So, for example, the flux for radius r=1/2 is ∑k=0∞HY(k,-2,1/2)≈43770.014449969689. The other values, also can be established on the basis of the recurrent relations (2.25), in other words (2.26).


## 2.2. The Fluxes on the Basis of the Series of the Hypercylindrical Functional Matrix

The discrete dimensional fluxes can be calculated as well on the “horizontal line”, that is, adding functions values on the submatrix seriesM[HY]k,n (Figure 1). For example, through the series development for k=3 the flux would contain the next members:
(2.30)∑n=0∞HY(3,n,r)=12π+12πr+6πr2+2πr3+12πr4+⋯+ε(12πrnΓ(n+1)).
Some values of discrete and continual fluxes (for r=1) are given in Table 2.Table 2Freedom degreek∑n=0∞HY(k,n,r)∑n=0∞HY(k,n,1)∫0+∞HY(k,n,1)dn02er/π1.7305119588645301.84608438006859712er5.436563656918095.24809906024688728er21.746254627672322.46216193622815312πer102.476810672082885.44633790197124464π(er-1)345.4810604213614237.53292390754055120π2(er-r-1)850.6988994458283508.72931334440466384π2[2er-(r+1)2-1]1654.544866434988892.20643559127927∑n=0∞(1680π3rn+4/Γ(n+5))2688.6618988873221338.1438558916638∑n=0∞(12288π3rn+5/Γ(n+6))3790.4276572623341768.5863448966209∑n=0∞(30240π4rn+6/Γ(n+7))4757.7028081498232104.22878227066410∑n=0∞(245760π4rn+7/Γ(n+8))5416.8054631352562289.34360590283511∑n=0∞(665280π5rn+8/Γ(n+9))5672.0250477771182304.87915812710812∑n=0∞(5898240π5rn+9/Γ(n+10))5520.7364051223062167.38351118232513∑n=0∞(17297280π6rn+10/Γ(n+11))5036.8154770260341917.751163982797⋮⋮⋮⋮50∑n=0∞HY(50,n,r)0.0000001178268380.000000029841827⋮⋮⋮⋮∞limk→∞∑n=0∞HY(k,n,r)=00limk→∞∫0∞HY(k,n,r)dn=0—∑k∑n52812.977966365550—


## 2.3. Some Continual Fluxes of the Hypercylindrical Function

The trend of the distribution of the vector fluxes is increasing, and then asymptotically falling, with the linear growth of the freedom degreen.From the standpoint of the functional analysis the most interesting series of the matrix M[HY]k,n are the ones referring to the freedom degrees k=2 and k=3. The first series (k=2) covers the known functions for the square size (8r) and surface (4r2).The members of the following series are, among the others, the cylinder functions of the surface (6πr2) and volume (2πr3) (Figure 1). The same series are interesting as well for calculatiing continual fluxes. So the continual natural flux for the hypercylinder surface is analysed in view of integrals, instead of series. This integral is specific, because its subintegral function is the reciprocal gamma function. Its value, as it is known, is equal to the value of Fransen-Robinson's constant [15]
(2.31)F=∫0∞1Γ(x)dx=e+∫0∞e-nπ2+ln⁡2ndn≈2.8077702420285.
The integral values of the wanted flux is now ∫0∞HY(2,n,1)dn=∑n=0∞HY(2,n,1)+∫0∞(8e-n/(π2+ln2n))dn, what concretely for the unit radius is
(2.32)∫0∞8Γ(n)dn=8(e+∫0∞e-nπ2+ln⁡2ndn)=8F≈22.462161936228.
In respect to the continual dimension n, the more generalized cylinder volume flux follows in view of Ramanujan-Hardy's integral [16, 17]:(2.33)∫0∞yxΓ(x+1)dx=ey-∫0∞e-xyx(π2+ln⁡2x)dx.
Ramanujan defined this integral, and it was analytically intensified by Hardy. In that sense, the previous expression can be applied in calculating the hypercylindrical function flux, when k=3, as(2.34)∫0∞HY(3,n,1)dn=∫0∞12πrnΓ(n+1)dn=12π(er-∫0∞e-nrn(π2+ln⁡2n)dn)(forx=n,y=r).
The integral can be calculated as well as the difference between the series and the integral with the value (for r=1)(2.35)∫0∞12πrnΓ(n+1)dn=∑n=0∞12πrnΓ(n+1)-∫0∞12πe-nrn(π2+ln⁡2n)dn|r=1≈278.6085228836.


## 2.4. The Progressions of the Vector Fluxes

The whole dimension flux in the freedom degree domain with the natural numbers, is obtained in the result of double amount with which are considered the integer values of the hypercylindrical functionHY(k,n,r), for all k,n,r≥0. This double series must be convergent, and this characteristic is in the function of hypercylinder radius. The flux can be watched as well for every column M[HY]k,n of the matrix individually. So, there is for the nth column (in the mark 〈n〉), the flux in the form of series
(2.36)ΦHY〈n〉(k,n,r)=∑k=0∞HY(k,n,r).


## 2.5. The Orthogonal Dimensional Fluxes

These fluxes are the fluxes of the all columns or of all series of the matrixM[HY]k,n. As the number of columns, that is of the series infinite, the total flux is as follows.Definition 2.1.
The dimensional flux (potentional) of the functional matrix with two freedom degreesk and n is defined as a double series(2.37)ΦHY(k,n,r)=∑n=0∞ΦHY〈n〉(k,n,r)=∑n=0∞∑k=0∞HY(k,n,r).
When the certain number of members is calculated, the flux has the next form:(2.38)ΦHY(k,n,r)=∑n=0∞∑k=0∞HY(k,n,r)=4π[r(10+4πr2)+eπr2erfc⁡(-rπ)[3+4πr2(3+πr2)]]+4[2(1+πr2)+πreπr2erfc⁡(-rπ)(3+2πr2)]+2[2r+eπr2erfc⁡(-rπ)(1+2πr2)]+2(1π+reπr2erfc⁡(-rπ))+1π[eπr2erfc⁡(-rπ)-1]+1π(erfi(rπ)2-r)+1π2∑k=0∞(-1)2k+12k+1(1-Γ(k+1,-πr2)k!)+12π[1-eπr2π+rerfi(rπ)-r2]+⋯.
The flux on the matrix series in the domain of the natural number is defined as a double series, but with the changed sequence of summing. This dimensional flux is, thus, defined as(2.39)ΩHY(k,n,r)=∑k=0∞∑n=0∞HY(k,n,r).
In view of previously established members, the matrix flux has the following form:(2.40)ΩHY(k,n,r)=∑k=0∞∑n=0∞HY(k,n,r)=2erπ+2er+8er+12πer+64π(er-1)+120π2(er-r-1)+768π2(er+r-r22+1)+∑n=0∞1680π3rn+4Γ(n+5)+∑n=0∞12288π3rn+5Γ(n+6)+⋯.
The equivalence of the orthogonal dimensional fluxes means the equality of the double series(2.41)ΦHY(k,n,r)=ΩHY(k,n,r).
So, for example, for r=1 the dimensional fluxes have the unambiguous numerical value that is:
(2.42)ΦHY(k,n,1)=ΩHY(k,n,1)=52812.977966365550.


## 2.6. The Application of the Recurrent Operators with Defining the Diagonal Dimensional Fluxes of theHY Function

In the previous analyses the defining of the dimensional fluxes of the matrixM[HY]k,n is performed in view of the adding up of the HY function value for (on) the columns, in other words, the series of this matrix. The more detailed analysis would be very ample including the exponential function, then error functions erf(z), erfc(z), erfi(z), the incomplete gamma function Γ(a,z) and the like. When we use the ideas of transition operators from the reference function into the destined HY function, in the functional hypercylindrical matrix M[HY]k,n we can also establish the dimensional fluxes values on the diagonals (2.44), whose total would present the whole flux for the matrix in the field of freedom degrees for natural numbers, that is, k,n∈N. Such matrix contains infinite number of elements. For reference functions we take HY functions on the positions of the first series of the matrix, and they are so-called zero HY functions: HY(0,0,r),HY(0,1,r),…,HY(0,n,r),… The destined functions are arranged on the “gradual” growth law (+Δk) and decrease (-Δn) of increments.Definition 2.2.
The flux series operatorϑ(Δk,Δn,0) is defined by the relation
(2.43)ϑ(Δk,Δn,0)=HY(k+Δk,n+Δn,r)HY(k,n,r)=πΔk/2rΔk+ΔnΓ(k+n-2)Γ(k+Δk+1)Γ(k+1)Γ(k+n+Δk+Δn-2)Γ((k+Δk+1)/2)Γ(k+12).(2.44)Meanwhile, as increments for the absolute value are mutually equal and unit ones, that is,|Δk|=|-Δn|=1, to them is assigned a new joint argument u (Δk=Δn=u). In addition, the starting value of the k-og freedom degree states k=0, so the theta operator gets the form(2.45)θ(u,-u,0)=2uπuΓ(u2+1).
The destined function is now calculated as(2.46)HY(u,n-u,r)=θ(u,-u,0)⋅HY(0,n,r)=2u+1πu-2rn-3Γ(n-2)Γ(u2+1).
The dimensional flux on the diagonal presents the total of ist particular nembers. So, for the first diagonal (with the mark 〈0〉) flux is equal
(2.47)Π〈0〉(k,n,r)=HY(0,0,r)=0,
then for the second,
(2.48)Π〈1〉(k,n,r)=HY(0,1,r)+HY(1,0,r)=0.
For the third,
(2.49)Π〈2〉(k,n,r)=HY(0,2,r)+HY(1,1,r)+HY(2,0,r)=0.
Or for the fourth,
(2.50)Π〈3〉(k,n,r)=HY(0,3,r)+HY(1,2,r)+HY(2,1,r)+HY(3,0,r)=10+2π+12π.
The flux in the nth diagonal would be reckoned in the sum form(2.51)Π〈n〉(k,n,r)=∑u=0nHY(u,n-u,r)=rn-3Γ(n-2)∑u=0n2u+1πu-2Γ(u2+1)(n≠2).
The flux for the value n=2 is calculated on the basis of the highest function value. Considering that the equivalence of the double factorial and trigonometric functions [18] is known
(2.52)u!!=2(1/4)[1+2u-cos⁡(πu)]π(1/4)[cos⁡(πu)-1]Γ(u2+1).
The expression (2.19) can after settling be written as well in the equivalent form
(2.53)Π〈n〉(k,n,r)=rn-3(n-3)!∑u=0nu!!(2π)2u(2π)3+cos⁡(πu)4.
So, we get for the fifth diagonal (n=4), using the expressions (2.51) or (2.53)(2.54)Π〈4〉(k,n,r)=∑u=0nHY(u,n-u,r)=2rπ(5π+38π2+1).
For the sixth diagonal (n=5) follows:(2.55)Π〈5〉(k,n,r)=∑u=0n=5HY(u,n-u,r)=r2π(5π+38π2+60π3+1).
As the number of diagonals is infinite, the total flux is formed as the series of all diagonal fluxes
(2.56)ΠHY(k,n,r)=∑n=0∞Π〈n〉(k,n,r),orconcretely,ΠHY(k,n,r)=∑n=0∞∑u=0nHY(u,n-u,r)=∑n=0∞(rn-3Γ(n-2)∑u=0n2u+1πu-1Γ(u+12)).
For example, approximately, the flux for r=1 and n=11 it concretely is(2.57)ΠHY(k,11,1)=10960120160π+1096014032+143729910080π+4264371680π2+58497280π3+290128π4+332π5.
In the developed, form the total flux has the polynominal structure of members(2.58)ΠHY(k,n,r)=10+2π+12π+r(10+2π+76π)+r2(5+1π+38π+60π2)+r3(53+13π+383π+148π2)+r4(512+112π+37π2+70π3)+r5(112+160π+1930π+375π2+5825π3+112)+r6(172+1360π+19180π+3730π2+975π3+42π4)+⋯.
The diagonal flux of the hypercylindrical function can be expressed by the series of general form
(2.59)ΠHY(v,r)=∑v=0∞avrv.
Here, r is the summation index with which we take into consideration the order of elements from left to right and from top to bottom on (along) the diagonal (Figure 3). The polynomial coefficients contain the rational numbers and the graded constant π (Table 3). The first three coefficients are equal to zero, so they are not put in the summation sequence. Its other values are (v=0,1,…,12) given in Table 3.Table 3vav- polynomial coefficient ∑vavrv010+2π+12π12π+10+76π21π+5+38π+60π2313π+53+38π3+148π24112π+512+19π6+37π2+70π35160π+112+19π30+37π25+582π3561360π+172+19π180+37π230+97π35+42π4712520π+1504+19π1260+37π2210+97π335+1150π4218120160π+14032+19π10080+37π21680+97π3280+575π484+33π5591181440π+136288+19π90720+37π215120+97π32520+575π4756+2279π51261011814400π+1362880+19π907200+37π2151200+97π325200+115π41512+2279π51260+143π63011119958400π+13991680+19π9979200+37π21663200+97π3277200+115π416632+2279π513860+905π619812147900160+1239500800π+19π119750400+37π219958400+97π33326400+115π4199584+2279π5166320+905π62376+13π712Figure 3
The representative submatrix (above) and the addition principle of the matrix members on the diagonal.The series approximation with 13 coefficients in the decimal notation is in the form(2.60)∑v=012avrv≈48,33+249,4r+716,87r2+1502,27r3+2546,01r4+3684,24r5+4705,22r6+5422,03r7+5727,08r8+5610,38r9+5143,66r10+4445,23r11+3642,42r12.
Approximately, the double series leads us to the solution that is very near to the correct one. Namely, for the unit radius and reducing on ∞~n=30 the double series of the diagonal flux develops the following structure:(2.61)ΠHY(u,n,1)≈∑n=0n=30∑u=0nHY(u,n-u,r)=73997539898837593289987313727222173626045880401920000+739975398988375932899873137136110868130229402009600000π-1+5107360421374632558310731373581864950795510579200000π+23851430587700445181187313711342572344185783500800000π2+3950250824632679929196942891890428724030963916800000π3+103021449478899595399537986310638624850739200π4+345714978541889974399922394521436201548195840000π5+2633045028485350315237270061246290137702400π6+70440913248925148230273375765578626721280000π7+996435517964121615726791790306561280000π8+3543877234517627836251417738197632000π9+1285336197080807170493211041753600π10+1402624867171115786408429792000π11+5134803604675464525994928000π12+3516397543390323285352000π13+101418347130107761784000π14≈52810,5103716.
The diagonal dimension flux is characteristic with its coefficients that contain the πn constant in the graded series numbers, as distinguished from the vertical fluxes with the destination of the errors functions and π and e constants. The horizontal fluxes, as it is presented (2.40) contain the exponential functions. Meanwhile, the total flux for the unit radius is convergent and it can be calculated with considerably greater value. The value aberration (2.61) of the accurate value is only (or 0,00467%)
(2.62)ΔΠ=52812,97796636-52810,5103716≈2,467595.
The whole dimensional continual flux (k,n∈N) of the unit hypercylindrical function HY(k,n,1) is equal to the value of the twofold integral(2.63)∬0∞2πk-1Γ(k+1)Γ(k+n-2)Γ((k+1)/2)dkdn,
and we find its solution on the analytical and numerical bases.


## 3. Conclusion

On view of the supposition of recurrent relations (1.2), namely, (2.25) and (2.26) that exist in the scope of the hypercylindrical function, we can calculate the discrete dimension flux of this function in the domain of integer freedom degrees. The quantitative flux value at the most depends of the formulated value of the hypercylinder radius. Meanwhile, as the function HY(k,n,r) is the three variables function, its dependence is by all means restricted by the values of k variables, namely n. In the paper are calculated several continual fluxes for the contour hypercylindrical functions, in view of Ramanujan-Hardy's integrals. With the continual flux in the domain k,n∈0,∞¯ the problem is more complex, because we must, for its defining, to perform the twofold integration (2.63). Although it is not yet calculated, we suppose that its value is very close to the discrete flux, obtained in view of the double series.Calculating the dimensional flux by thediagonal algorithmis much simpler and faster by computer, because the total flux is now defined as the convergent graded series and does not contain as components the special functions. In any case its value is identical with the fluxes that are calculated on the basis of the series, relatively to HY-matrix columns, so there is valid the numerically verified statement that(3.1)ΦHY(k,n,r)=ΩHY(k,n,r)=ΠHY(k,n,r)|r=1=52812.977966365550.
In any case, this solution is initial for solving the other dimensional fluxes, both for hypercylindrical and hyperspherical [3, 19], in other words hypercube function [20].



---

*Source: 245326-2011-06-18.xml*