The present paper tries to find out the common ground for the zeros of the Riemann zeta function and of the modified Bessel functions (or Bessel functions of imaginary argument z) for imaginary argument z and, furthermore, for the absence of zeros of the Gaussian Bell function . For the function now called Riemann zeta function which was known already to Euler but was extended by Riemann to the complex plane Riemann expressed the hypothesis that all nontrivial zeros of this function lie on the axis that means on the axis through and parallel to the imaginary axis y (Riemann hypothesis)    (both with republication of Riemann’s paper) and many others, e.g.      . Riemann never proved his hypothesis. He introduced in  also a Xi function which excludes the only singularity of the function at and its trivial zeros at and possesses more symmetry than the zeta function . Concerning their zeros it is equivalent to the nontrivial zeros of the zeta function. In present paper we will mainly have to do only with this Xi function which we displaced in a way that its zeros lie on the imaginary axis provided; the Riemann Hypothesis is correct and we denote as Xi function . With respect to the position of the zeros the function is fully equivalent to the nontrivial zeros of the Riemann zeta function only with displacement of the imaginary axis to these zeros.
The content of this article was not intended as a proof of the Riemann hypothesis but during the work we found a further, as it seems, essential building stone for its proof by the second mean-value approach which is represented in Appendix. The article is merely intended as an illustration to the zeros of a function with a possible representation in an integral form given in Section 2 (Equation (2.8)) with monotonically decreasing functions satisfied by the Riemann Xi function and by the modified Bessel functions and explains why the Gauss bell function although it can be represented in such form does not possess zeros. Other kinds of interesting illustrations of the Riemann zeta function (and of other functions) by the Newton flow are given in   . A main purpose was to understand how the zeros in the Taylor series approximations of such functions behave when we go from one order of the approximation to the next higher one. To get the possibility of a comparison with the pictures of zeros for functions without an integral representation of the mentioned form we made an analogous picture for an unorthodox entire function in Section 9.
2. Basic Equations for the Considered Functions
The Xi function to the Riemann zeta function is defined by
where is the Riemann Xi function  which is related to the Riemann zeta function by1
The Riemann zeta function is basically defined by the following Euler product
where is the sequence of prime numbers beginning with . The definition of the Riemann Xi function (2.3) is equivalent to the definition by the following (Dirichlet) series for complex variable
which is convergent for and arbitrary y and can be analytically continued to the whole complex z-plane. The function is an entire function which excludes the only singularity of at and its “trivial” zeros at .
Next we consider the whole class of modified Bessel functions of imaginary argument which is connected with the basic class of Bessel functions in the following slightly modified form by ( )
The functions are entire functions which satisfy the differential equation
In comparison to the functions exclude the zeros or infinities of the first ones at but the other zeros remain the same for bothfunctions.
Finally, we consider the Gaussian functions with parameter which can be represented by the following integral representation (continued from the imaginary axis y to the whole complex z-plane)
which become Gaussian bell functions for imaginary argument . Clearly, the functions do not possess zeros on the imaginary axis and zeros at all.
The three mentioned types of functions written as have in common that they are symmetrical functions in z and that they possess a representation by an integral of the type
with monotonically decreasing functions for that means
The Taylor series of can be written in the form
where are defined as the moments of the symmetrical function with respect to the reference point . A consequence of the definitions in (2.8) and (2.10) is
The odd moments of the function in the definition (2.10) vanish. The function at is equal to the zeroth moment of the function and is independent of the chosen reference point. In the following the moments of the function play an important role.
That is a symmetrical function in u is, in principle, not necessary since the integration over u in (2.8) is restricted by but the symmetry permits to extend the integral over negative values of u using it in the form
which for imaginary is a Fourier transformation of with the inversion
if the integral exists in some sense (e.g. weak convergence). We consider this now more explicitly.
The explicit representation of the Xi function to the Riemann Xi function in this form (2.8) together with (2.9) is
with the special values
This was derived in detail in  . The function together with its first derivative is represented in Figure 1.
The function is monotonically decreasing from on and decreases rapidly for . The function increases rapidly for real and decreases rapidly for imaginary . The symmetry
Figure 1. Representation of and of its first derivative . The function is monotonically decreasing for and together with its derivatives vanishes rapidly for .
of is not easily to see from (2.14) and it was a genuine surprise for us to meet such a kind of a symmetrical function (see discussion in  ). The function , for example, is already not a symmetrical function for .
In case of the modified Bessel functions of imaginary argument the following basic integral representations are known which for the functions proportional to may be written as follows (taking into account )
or expressed by the Confluent Hypergeometric function
The Taylor series expansion is
The functions possess the principal form
with the following explicit expressions for
where is Heaviside’s step function defined by
It restricts the upper limit of integration in (2.19) according to the choice in (2.20) to . The functions are equivalent to functions of only and are in this sense symmetrical functions of u and for real we have to choose the real value of in case of non-integer .
The first six cases of the function with integer or semi-integer index are
The function does not exist since the corresponding integral in (2.16) is divergent or is not finite for but exists. Furthermore, we have introduced in (2.16)
amplitudes and we will soon see that it is favorable to choose them for our purpose as constant for the whole class of functions . We illustrate the functions and the corresponding functions with a certain important modification in next Section 3.
It is well known that the Bessel functions possess zeros only on the real axis x and, correspondingly, only on the imaginary axis y. For this there exists a direct proof using their differential equations and which is similar to the derivation of duality and orthogonality relations (e.g.   ).
The considered Omega functions for the Riemann Xi function and for the modified Bessel functions and also for the Gaussian bell function have in common that they are monotonically decreasing functions for up to .
3. Modified Bessel Functions with Stretched Argument of the Kernel Function and Limiting Transition to Gaussian Bell Function
We now calculate the even moments of the functions in (2.20) which lead to well-known integrals
and in the special case of zeroth moments
These are the areas under the curves on the positive u-axes. The following considerations show that it is favorable to make equal independently on parameter and to keep constant and we choose
To keep in addition also constant we have now the only possibility to introduce a stretch factor to the variable u and we make the transformation (now already with choice (3.3))
with the definition of
from which follows for the zeroth moments of
The functions are graphically represented in Figure 2 for some low
Figure 2. Functions with stretched variable in comparison to functions . G rid lines are set at the values where the finite functions end. In the chosen scale the curve for is visibly already hardly to distinguish from a Gaussian bell function and, therefore, is not drawn here.
values of up to .
The corresponding transformation of the functions using (2.16) and (3.3) are
Their Taylor series are
or with substitution (see also (2.5))
The first three functions with integer and with semi-integer are explicitly
The functions on the imaginary axes for integer and semi-integer index up to are illustrated in Figure 3.
Figure 2 and Figure 3 admit the conjecture that and for real variable y become Gaussian bell functions in the limit . For in Figure 2 the function is in visible way already hardly to distinguish from the Gaussian bell function (below we see to ).
That the mentioned approach for to a Gaussian bell function is really true we establish exactly and determine these limits.
We now show that the new functions in (3.7) for go to a Gaussian function that becomes a Gaussian bell function for imaginary . As auxiliary formulae for in case of follows
Figure 3. Functions for values of up to . W ith increasing from a certain on the first zeros of increase (see also text) and the amplitudes of the oscillations decrease and in the limiting case they vanish and the resulting Gaussian function (on the imaginary y-axis ) does not possess zeros. This Gaussian curve is not well to distinguish in the bulk of the other curves and is not drawn here.
and analogously for
Applying these approximations we find from (3.8)
This means that in the limiting transition the functions in (3.7) approach to a Gaussian function according to
or by substitution with, in general, complex variable y
where we use the identity . From the corresponding limiting transition using the definition in (3.4) we find
We may check the transition from to via (2.8) using the auxiliary formula (2.7) with .
Furthermore, we find using the approximation (3.11) that the factor between and in (3.5) approaches for according to
with a high precision and monotonically increasing without a finite limit.
We now show that for high increasing the first positive roots of and thus also the higher roots increase. For this is clear since it is proportional to and it is well known that their roots are situated approximately at
However in the argument in the Bessel functions is stretched by the factor and therefore on the imaginary axis we have the functions that diminishes the values for the roots (3.18) by the factor and with the approximation (3.17) the zeros of are now situated approximately at
Thus the first roots of go with proportional to also to infinity. This explains more in detail why the first zeros of the curves for in Figure 3 go to the limiting case for .
We remind that the starting point for the derivation of the approximations was making equal the zeroth moments by stretching the argument of the functions to new functions leaving constant the amplitudes . The zeros of for move in this process also to infinity but “very slowly”. This is a good illustration why the Gaussian function which on the imaginary axis y becomes a Gaussian bell function does not possess zeros at all although its Omega functio is monotonically decreasing. The first zero and in this way all other zeros are moved by the limiting transition to infinity although very slowly. The same is the case with the discontinuities of derivatives in the functions .
4. Graphical Representations to Zeros of the Xi Function to Riemann Zeta Function in Approximations by Truncated Taylor Series
In this Section we consider the Taylor series expansions of the Xi function to the Riemann Xi function in powers of z defined in (2.14) according to
The coefficients in this Taylor series are the even moments of the function as defined.
We truncate the Taylor series of at upper summation numbers
and calculate all zeros of these approximations up to a certain maximal M and make graphical representations of their zeros. As mentioned the coefficients of the series are determined by the even moments of the function where the odd moments vanish automatically if we determine as symmetrical functions according to (2.8) and (2.10). In the following Sections we discuss the same also for modified Bessel functions (in particular for ) and for the limiting transition to a Gaussian function and compare this with the zeros for the Xi function (2.14) to the Riemann zeta function. In each approximation of the Taylor series we find 2M complex zeros. It was very interesting to see which difference appears between the zeros on the imaginary axis and the bulk of the other zeros with increasing number M. This is best seen in graphical representations of the zeros.
For the function related to the Riemann hypothesis which is explicitly given in (2.14) we calculated this series numerically up to with sufficiently high precision and obtained the following numerical coefficients (we write here explicitly down only a few of the obtained coefficients)2
It is interesting to mention that the coefficients possess an absolute minimum for with and, apparently, are (“slowly”) monotonically increasing for .
The first term in (4.3) can also be written (see  )
corresponding to the splitting of terms in the representation (2.14) of after the substitution .
In the truncation of the series in powers of z with the highest term proportional to we found by numerical solution of the corresponding algebraic equations of degree 2M the 2M complex solutions from which the following pairs of zeros lie on the imaginary axis 3:
Table 1. Zeros on the imaginary axis for the first 2M approximations up to 2M = 80.
In the following we give graphical illustrations of all zeros of the Taylor approximations in the complex z-plane where all zeros up to a certain order 2M are taken into account and where one may see how the zeros change from an order to a higher order.
We explain first how the following figures are made. We take a certain order 2M of the Xi function given by the truncated Taylor series (4.2) and determine numerically all of its zeros and represents them by points in the complex -plane where we choose the same scale on the x- and y-axis. Two variants are made, first the representation by isolated points and second the representation by connected neighbored points. The obtained partial pictures are a little different for odd and even M that is represented in Figure 4 for and .
Then we calculate and represent all zeros of the Taylor approximations of (i.e. of ) up to a certain maximal M and represent the zeros in described way by isolated points and by connected neighbored points in each of the approximations up to the maximal M. This is made in Figure 5 and Figure 6 for maximal . These approximations capture already approximations of the first two nontrivial zeros of the Riemann zeta function on the positive y-axis at and at seen in Figure 5 by some accumulation of points at these values. Figure 6 shows the same picture with all neighbored points in each approximation joined as described. Since not all details are well recognizable in Figure 6 the same is made in Figure 7 but
Figure 4. Zeros of for according to (2.14) (i.e. for Riemann hypothesis) in approximation of the Taylor series with and . The obtained 58, respectively, 60 zeros are shown in the complex plane as points without mutual distortion of the axis lengths. In the pictures to the right-hand sides we have joined neighbored numbers. On the imaginary axes where it is not clear which zeros are neighbored to zeros off the axis we went in clockwise sense on the positive part first to the highest zero and then to the next lower zeros and so on and then from the lowest zero on the imaginary axis clockwise to the next complex zero. This shows also the way we went in the next picture where the details on the imaginary axis are not so clearly visible. The two zeros at and at on the positive and negative imaginary axes are already decoupled from the main bulk of zeros and the zeros at begin only to decouple.
Figure 5. Zeros of Xi function to Riemann zeta function in the first 30 approximation of its Taylor series with . The neighbored zeros are not joined in each approximation. We see already the beginning accumulation of points at the first two genuine zeros of the Riemann zeta function at and .
Figure 6. Zeros of Xi function to Riemann zeta function in the first 30 approximation of its Taylor series with . The neighbored zeros are joined in each approximation separately. In the immediate neighborhood of the axis the picture becomes a little confusing since then happens a big step to the smallest genuine zero of the Xi function and then rises up to higher zeros on the imaginary axis before it goes to the main bulk of zeros outside the imaginary axis.
Figure 7. Zeros of Xi function to Riemann zeta function in the first 20 approximation of its Taylor series with . The neighbored zeros are joined in each approximation separately. We see that the zeros of the main bulk in the approximations go slowly but with great regularity (though not proved) to infinity for and vanish in this way as zeros of the whole function. Only the first zero at as genuine zeros of the Riemann zeta function are seen here as decoupled and stabilized.
only for the first 40 Taylor approximations where this is clearer to see. These pictures show that the zeros on the imaginary axis stabilize from order to higher orders at the genuine zeros of the Riemann zeta function on the y-axis and separate themselves from the main bulk of zeros in a considered order which do not lie on the imaginary axis. That this remains in this way for is, clearly, only a conjecture but in Section 8 we try to understand this by some analytic approximations.
In Figure 7 for the case we see only one accumulation point of the zeros on the positive and negative y-axis corresponding to the two zeros . The four higher zeros on the positive y-axis (correspondingly negative y-axis) belong already to not yet stabilized approximations to the second zeros at .
It is interesting to compare the functions and to the Riemann zeta function with a corresponding Gaussian functions and where the two parameters of the Gaussian function which are the amplitude and the stretching of the parameter u are chosen in the way that the first two terms of the Taylor series approximation are equal. For the function
this means for the Gaussian function
One may check that generally
For the moments of these functions result the inequalities ( )
The resulting function is a Gaussian function of the complex variable
It possesses another amplitude in comparison to (i.e. to 0.497120778). Clearly, as a Gaussian function it does not possess zeros on the y-axis and zeros at all. The function to the Riemann Xi function and the considered modified Bessel functions possess the common property that they vanish in infinity more rapidly (or are even finite) than the Omega function to a Gauss function. In principle, this does not exclude Omega functions which vanish less rapidly. For example, a function provides without zeros of at all in finite regions of the complex plane but with poles that regarding the zeros is the same as for a Gauss function (see Section 8 where this is explained by motion to infinity from finite approximations). However, by comparison with Gauss functions we may get inequalities for the moments of the considered Omega functions.
5. Zeros of the Taylor Series Approximations of the Function
In this Section we now come back to the modified Bessel functions in (3.7) to the Omega functions in (3.4). From this class of functions which in the limit go to a Gaussian function we choose the function
and give graphical representations of the zeros for its Taylor series approximations
They are represented in Figure 8 and in Figure 9 up to the approximation for with the difference that in the first all zeros are presented together and in the second we have joined the neighbored zeros in each approximations. The genuine zeros of lie at according to
Figure 8. Zeros of the function in the first 30 approximations of its Taylor series with . The neighbors within an approximations are not joined and it is not fully easy to see which point belongs to a certain approximation.
Figure 9. Zeros of the function in the first 30 approximations of its Taylor series with . The neighbors within each approximations are here joined. Since in each approximation a point on the positive imaginary goes approximately to the first zeros at the connections on the imaginary axis are not easily recognizable.
We may expect that for other low values of index in the functions we get similar illustrations with small topological distortions of the pictures for .
The function possesses a peculiar importance since for monotonically decreasing functions (for with ) we may apply the second mean-value theorem (Gauss-Bonnet theorem; see, e.g. Courant  (chap. IV), Widder  ) to bring the integral (2.14) to the form
Herein, is a mean value function from which we may assume that it is an analytic function of z since the integral on the left-hand side depends analytically on z. In next Section we consider shortly the problem of zeros of the function of the form (5.4).
6. General Conditions for Zeros of Functions
The problem of zeros of functions of the form (5.4) leads essentially to the problem of zeros of with exclusion of the zero at . If we separate the real and imaginary part of according to
and if we then separate the real and imaginary part of the function (5.4) can be written in the form (the separation of real and imaginary part in the first factor is uninteresting since a zero at is excluded)
From both forms of the right-hand side in (6.2) we find that for zeros of the following two conditions 
have to be satisfied at the same time and this is necessary and sufficient.
We now consider the special case of on the imaginary axis y and find from (6.2)
Since due to symmetries (2.8) the function has to be a real-valued function that for in (6.4) is only possible if vanishes we find on the imaginary axis
This follows also immediately by application of the second mean-value theorem to the integral for where the mean value can only take on real values here in dependence on y as parameter. Thus for zeros on the imaginary axis the two conditions (6.3) reduce to only one condition
The condition in the second line of (6.3) is then identically satisfied. If one knows the function one may determine the zeros on the imaginary axis analytically or, at least, numerically but for many problems including the considered one it is not necessary to know these zeros explicitly.
We return to the general case of general values x on the real axis with the conditions (6.3) for zeros. From (2.8) follows that should be an analytic function for all z for which the integral exists that means to satisfy the Cauchy-Riemann equations. For a general analytic function one can derive by integration of the Cauchy-Riemann equations the following relations in case of 
Now come into play the following operator identities (operator identities are such identities which can be applied to arbitrary functions to provide function identities) 
More general identities of such kind can be derived by representing the Cosine and Sine functions by Exponential functions and using that applied to analytic functions displace the argument of these functions to that is discussed in  . Using (6.7) and (6.8) we may write the conditions for zeros (6.3) in the following way
If we now apply the operator to the first of these conditions and the operator to the second one we obtain
By addition of both equalities using the operator identity
follows from (6.10)
This consequence for zeros results from both conditions (6.3) and provides a necessary condition for the imaginary values y of all zeros for which their imaginary value has to agree with one of the solutions for zeros on the imaginary axis. Thus we have to take all solutions y for zeros on the imaginary axis which we denote now by and which satisfy the conditions
and have to look for solutions x at the same time of both equations
If one finds such solutions of one of the conditions then it is almost in all cases not a solution of the other condition and, therefore, is then not a zero. Only in very few cases of depending on the function in (2.8) one would expect to find values as solution of both conditions (6.14). This is a very strong restriction to the satisfaction of both conditions (6.14). It seems that the set of functions with such a property can be only a set of measure zero within the set of all possible functions in (2.8). One such case is when is stepwise constant and monotonically decreasing with equal lengths of the steps and this is the only case (Appendix A). The corresponding functions are then superpositions of functions with different amplitudes. This was not correctly discussed in  and Katsnelson  showed an error in a short Email4 and we recognized it  but it was already seen in  that the possible zeros off the imaginary axis, i.e. of with , must possess imaginary parts which agree with one of the zeros on the imaginary axis. Reactions regarding concern about the applicability of the Bonnet theorem for present case were expressed by others, in particular, in a nice Email by Gélinas  with appended file but I could not find it published now. This also means that the Riemann hypothesis which is the absence of zeros of the Riemann Xi function off the imaginary axis through was not correctly solved by the second mean-value approach to this time although it was very improbable that a nontrivial zero of the Riemann zeta function off the imaginary axis has exactly the same imaginary value as that of one on the imaginary axis through . We add now in Appendix A a further important building stone to a full proof of the Riemann hypothesis by the second mean-value theorem which seems to be “deciding” (?).
7. Zeros of Taylor Series Approximations of Gauss Function and Absence of Genuine Zeros Understood in Uncommon Way
We now consider the Gauss function and calculate all zeros of its low-order Taylor series approximations
In Figure 10 and Figure 11 we illustrate then their zeros up to a certain maximal M in the same way as we did this for the Xi function to the Riemann zeta function in Section 4 and for the Xi functions to the modified Bessel function in Section 5. This provides good analogies and shows the
Figure 10. Zeros of Gaussian Function in the first 30 approximations of its Taylor series with . The neighbors within an approximations are here not joined and it is not easily to see which point belongs to a certain approximation.
Figure 11. Zeros of Gaussian Function in the first 30 approximations of its Taylor series with . The neighbors within an approximations are joined in each approximation. The whole bulk of zeros in each approximation goes here with increasing higher approximation to infinity.
essential differences in the pictures for the Gauss function to the other discussed functions with genuine zeros when we increase the order of approximation. It becomes obvious that in the limiting transition to the Gaussian function all zeros of the finite-order Taylor approximations go to infinity with great uniformity from order to higher order and the Gaussian function does not possess genuine zeros in finite regions of the complex z-plane. In Section 3 we established in (3.13) a limiting transition from modified Bessel functions with stretched arguments to a Gauss function and it was interesting to see how the zeros of finite-order Taylor series approximations move from approximation to the next higher approximation in the limiting case to infinity. The pictures for the Gaussian functions are very homogenous in contrast to the considered functions with genuine zeros. This provides an uncommon view onto the absence of zeros of Gaussian functions despite a representation in the form (2.7).
8. Approximations of the Zeros from One to the Next Higher Orders
The zeros in each approximation for the considered functions are either pairs on the imaginary axis or in majority quadruples and there was no doubt which are pairs and which are quadruples even in case that their real part is small compared with maximal modulus of the zeros and since the whole number of zeros has to be 2M. It is noticeable that new zeros on the imaginary axis when they first appear in the 2M-th approximation may disappear in the next higher -th approximation from the imaginary axis and reappear then in the -th approximation as can be also seen from Table 1. In every case when there appeared a new zero on the imaginary axis the next lower zeros began to stabilize and to decouple from the main bulk of zeros in the complex domain and stabilize there in each new -th approximation as we may see in Figures 4-7 for the Xi function to the Riemann hypothesis and in Figure 8 and Figure 9 for a modified Bessel function . In Figure 10 and Figure 11 for the function the zeros on the imaginary axis show a similar picture with alternatingly generating and not generating zeros on the imaginary axis from one approximation 2M to the next higher approximation . The main bulk of zeros in the complex domain in all these pictures drifts with their modulus to infinity although very slowly that we can see in Figure 11. In the other pictures this is the main bulk of zeros which does not correspond to genuine zeros of the considered functions whereas the lower zeros on the imaginary axis stabilized more and more to the genuine zeros.
To understand the discussed behavior of the zeros from order to next higher order we try to discuss this now in some approximation. We suppose that we have the Taylor series approximation of a Xi function in the order 2M in the form
We assume that is an exact zero in 2M-th approximation that means a solution of the equation
Then we try to calculate next higher solutions which are near to from the next higher approximation which satisfies the equation
In full generality this equation would provide the solutions of zeros from only one arbitrary solution in the considered approximation but in such generality we cannot and do not want to solve it. Well soluble is the equation for additions to leading to a quadratic equation for . This provides two solutions in a neighborhood for the considered . Thus we now make an expansion of the left-hand side in (8.3) up to quadratic terms in
If we neglect from the additional terms of the -th approximation the terms proportional to and as small terms we obtain the following quadratic equation for
The two solutions of this equation are
or more compactly written
where are the n-th derivatives of taken at . One may express this also by the logarithmic derivative of that, however, is inconvenient since it goes into the formula as a denominator.
For the most interesting case of points on the imaginary axis we find from (8.6)
The sum terms in braces are real ones but likely change their signs at the zeros (we think that it can be proved). The first sum term in braces changes a little the imaginary value of the root but the second term with the root with the two possible signs must be real or imaginary. In case of negative values of the content of the root it gives two imaginary values and together with the whole expression it provides two corrections off the imaginary axis. This can be seen in many of the picture for the roots. However a full discussion of the behavior from approximation to next higher approximation by (8.8) is complicated and has to describe how the roots stabilize on the imaginary axis in dependence on the moments. Such a discussion we cannot give to this time.
9. About the Zeros of an Unorthodox Function in Their Taylor Series Approximations
We consider here shortly for comparison with the pictures for the up to now discussed functions with a representation of the principal form (2.8) the following unorthodox entire function of an essential other kind
with respect to the zeros in its finite Taylor series approximations. This function plays a role for the calculation of the properties of coherent phase states  .
The function in (9.1) possesses even and odd powers of variable z and, therefore it is not symmetrical with mirror symmetry to the real and imaginary axis but only symmetrical to the real axis. Taking separately the even and odd powers of z and applying the duplication formula for the factorials one may represent (9.1) in the form
with the possible approximations in the coefficients for (see (3.11))
This shows that for real positive grows moderately faster than the function . Therefore, the function is moderately increasing for real and moderately decreasing for with the Taylor series
and, clearly, possesses the same zeros as the function .
We now give a graphical representation of the zeros of the Taylor series approximations
of the function up to in Figure 12. Apparently, the computer calculated correctly up to this high approximation that we judged only from the optical impression of the figure in comparison to figures of such kind for smaller values N. To join neighbored points of each Taylor series approximation, at least, for such high maximal N becomes unfavorable.
In Figure 12 we see the first three pairs of zeros as some accumulation points.
Figure 12. Zeros of unorthodox though entire function in its first 120 Taylor approximations with . The first pairs of zeros at , , appear as accumulation points in this scheme.
Up to now we calculated only the first 4 pairs of zeros with sufficient accuracy  but S. Skorokhodov from the “Computing Centre of the Russian Academy of Sciences” calculated much more pairs of zeros with an essentially higher accuracy as he informed me in a nice email with the calculated zeros in the appended file  (see also  ). These zeros agreed with my few in lower accuracy calculated zeros5.
In this article we illustrated the behavior of the zeros for low-order Taylor series approximations of the Xi function as equivalent to the Riemann zeta function (in view of their nontrivial zeros) and the same of the modified Bessel functions . By scaling of the variable z in the modified Bessel functions we derived by limiting transition the Gaussian function proportional to (with additional scaling of variable z). The considerations do not pretend to be a proof of the Riemann hypothesis but support the belief to its validity.
There are some bothering problems connected with the zeros of the Riemann zeta function and of functions of the kind (2.8) with symmetrical functions which are monotonically decreasing up to their vanishing on the u-axis and which are the differences between strictly monotonically and discontinuously vanishing of these functions. The application of the second mean-value theorem leads to the Bessel function of the kind with only zeros on the imaginary axis. One may imagine how the mapping onto these function by the second mean-value theorem transforms continuously functions with zeros only on the imaginary axis to the mentioned Bessel function but it is difficult to imagine how this is implemented in this Bessel function for functions of the basic integral representation (2.8) with additional zeros outside the imaginary axis but in Appendix A it is shown how this apparent contradiction can be solved by a limiting transition. It seems that the presence of zeros off the imaginary axis is restricted to step-wise discontinuous monotonically decreasing functions with periodic steps where the function to itself is step-wise discontinuous and this vanishes in the limiting transition to step lengths zero. Is the second mean-value theore applicable to include such cases? We think that it is applicable.
The time is now mature for a final solution of all aspects of the Riemann hypothesis concerning the nontrivial zeros of the Riemann zeta function.
Figure 3, Figure 8, Figure 9 and Figure A1 were made with “Mathematica 10.4” in preparation of this article. All other Figures were made years ago (about 2010) by “Mathematica 6” and were now after their translation to “Mathematica 10.4” a little improved in few details.
Appendix A: Stepwise Constant Functions and a Further Building Stone for the Proof of the Riemann Hypothesis
In this Appendix we deal with the general case of monotonically decreasing nonnegative real-valued ladder functions with equal step lengths which satisfy the condition
with small improvements in comparison to  , in particular, also in the notations.
Thus we consider now the following stepwise constant functions with equal step lengths and with monotonically decreasing discrete amplitude
where denotes the Heaviside step function. The function possesses then the form
Due to monotonicity all amplitudes in front
of are nonnegative. We calculate now these functions more in detail. Using the relation
to the Chebyshev polynomial we find from (A.3)
The stable zeros on the imaginary axis under variation of u0 are determined by but due to symmetry we discuss sometimes so as if we take into account only positive y. We now bring the factor which determines these zeros in (A.5) to the left-hand side and consider the zeros of the second factor by its separation in real and imaginary part
If we insert for y any zero on the imaginary axis determined by with then the imaginary part in (A.6) vanishes and we find
In dependence on the chosen m and the amplitudes
one may solve the equations
and all solutions provide zeros of for the Omega function in (A.2) off the imaginary axis. If we here formally insert on the
left-hand side of (A.8) we see using that
this left-hand side becomes positive with no solution for x at all. If we make the steps smaller then the first zero on the imaginary axis grows
according to . For the limiting transition and it is
not important that we choose the amplitude differences in the middle between
two steps as and this could be, e.g. also
In Figure A1 we illustrate this by a stepwise function of the kind (A.2) which approximates the Omega function (2.14) to the Riemann zeta function. If we now make the length of the steps smaller then the first imaginary zero
of the ladder approximation becomes larger according to and
goes for to infinity and since any zero off the imaginary axis must possess an imaginary value which is the same as from zeros of the imaginary axis  (see also (6.12)) such values also go to infinity in the complex plane. Therefore, in the limiting case when and the ladder curve approaches the function exactly their zeros off the imaginary also have to go to infinity and they vanish from the pictures of finite-valued zeros in
Figure A1. Omega function to Riemann Xi function (2.14) approximated by stepwise constant functions. The steps of equal length are here chosen and 0.025 and the first zero of the approximation curve of on the imaginary axis are here at and . With step length the imaginary parts of possible zeros off the imaginary axis for the ladder approximations of to the functions go to infinity when they approaches to the exact function .
the complex plane. This concerns many functions among them also the Xi function to the Riemann hypothesis.
For the modified Bessel function where this is known from other more direct proofs, e.g.   this is only an affirmation. This also means that for any function with an Omega function which is different from the described stepwise constant ladder functions with equal step length (A.2) by the same limiting procedure follows that all zeros off the imaginary axis y are absent in the complex z-plane, for example:
For this function it is directly seen that possesses zeros only on the imaginary axis. Only for the monotonically decreasing ladder functions (A.2) as mentioned the described limiting transition does not lead to other functions than (A.2) and only these functions may possess zeros outside the imaginary axis which y values in addition have to agree with one of the zeros on the imaginary axis.
We could not find a reason why the Second mean-value theorem should not be applicable to the problem of zeros in the Riemann hypothesis and for other described functions.
1Riemann denotes the complex variable by that is for the Xi function and for the zeta function.
2We made the calculations two times with “Mathematica 3” (up to 2M = 80 and 18 digits) and with “Mathematica 6” (up to 2M = 100 and 20 digits) with a time difference of some years. We did not get full agreement mainly in the last 5 digits in the higher coefficients and came already near to the limit of capabilities of our PC.
3From about on the values given in the third and fourth column on the right-hand side did not become fully stable in our calculations in dependence on the number of sum terms taken into account in (2.14) and the chosen upper limit of integration that is rather due to the limits of our computer capabilities. The genuine value for the third root is near to 25.0109 instead of the stabilized 25.0101 seen in the table. The fourth root is at 30.4249. Therefore, the table reflects one such calculation. Onto the following graphical illustrations of zeros this does not have a visible influence. From the series (4.3) one may guess how difficult such calculations are even for a computer.
4The full text in the Email from 06.01.2017 with the subject line “Riemannsche Vermutung” was the following: “The result is wrong. Counterexample:
for , but small enough.
5It seems to be possible that with other now derived representations we can calculate more pairs of zeros than given here.
 Neuberger, J.W., Feiler, C., Maier, H. and Schleich, W.P. (2014) Newton Flow of the Riemann Zeta Function: Separatrices Control the Appearance of Zeros. New Journal of Physics, 16, Article ID: 103023.
 Wünsche, A. (2016) Approach to a Proof of the Riemann Hypothesis by the Second Mean-Value Theorem of Calculus, General Mathematics. Advances in Pure Mathematics, 6, 972-1021.
 Wünsche, A. (2017) Correction and Supplement to Approach for a Proof of Riemann Hypothesis by Second Mean-Value Theorem. Advances in Pure Mathematics, 7, 263-276.