|
|
Line 1: |
Line 1: |
| {{quote box
| | Oscar is what my wife loves to contact me and I completely dig that title. North Dakota is her beginning place but she will have to move 1 working day or another. He utilized to be unemployed but now he is a meter reader. What I adore performing is performing ceramics but I haven't produced a dime with it.<br><br>Look at my blog :: home std test kit [[http://www.hotporn123.com/blog/150819 click through the following website]] |
| |align=right
| |
| |width=33%
| |
| |quote= [[Émile Borel|Borel]], then an unknown young man, discovered that his summation method gave the 'right' answer for many classical divergent series. He decided to make a pilgrimage to Stockholm to see [[Gösta Mittag-Leffler|Mittag-Leffler]], who was the recognized lord of complex analysis. Mittag-Leffler listened politely to what Borel had to say and then, placing his hand upon the complete works by [[Karl Weierstrass|Weierstrass]], his teacher, he said in Latin, 'The Master forbids it'.
| |
| |source=[[Mark Kac]], quoted by {{harvtxt|Reed|Simon|1978|p=38}}
| |
| }}
| |
| | |
| In mathematics, '''Borel summation''' is a [[summation method]] for [[divergent series]], introduced by {{harvs|txt|first=Émile|last=Borel|year=1899|authorlink=Émile Borel}}. It is particularly useful for summing divergent asymptotic series, and in some sense gives the best possible sum for such series. There are several variations of this method that are also called Borel summation, and a generalization of it called [[Mittag-Leffler summation]].
| |
| | |
| ==Definition==
| |
| | |
| There are (at least) three slightly different methods called Borel summation. They differ in which series they can sum, but are consistent, meaning that if two of the methods sum the same series they give the same answer.
| |
| | |
| Throughout let ''A(z)'' denote a formal power series
| |
| | |
| :<math>A(z) = \sum_{k = 0}^\infty a_kz^k</math>,
| |
| | |
| and define the Borel transform of ''A'' to be its equivalent exponential series
| |
| | |
| :<math>\mathcal{B}A(t) \equiv \sum_{k=0}^\infty \frac{a_k}{k!}t^k.</math>
| |
| | |
| ===Borel summation===
| |
| Suppose that the Borel transform converges to an [[analytic function]] of ''t'' near 0 that can be [[analytic continuation|analytically continued]] along the positive real axis to a function growing sufficiently slowly that the following integral is well defined (as an improper integral), the '''Borel sum''' of ''A'' is given by
| |
| | |
| :<math>\int_0^\infty e^{-t} \mathcal{B}A(tz) \, dt. </math>
| |
| | |
| If the integral converges at ''z'' ∈ '''C''' to some ''a(z)'', we say that the Borel sum of ''A'' converges at ''z'', and write <math> {\textstyle \sum} a_kz^k = a(z) \,(\boldsymbol B) </math>. (Some authors use the term "Borel sum" to mean in addition that the Borel transform converges for all positive ''t''.)
| |
| | |
| ===Weak Borel summation===
| |
| Let ''A<sub>n</sub>(z)'' denote the partial sum
| |
| | |
| :<math>A_n(z) = \sum_{k=0}^n a_k z^k.</math>
| |
| | |
| A weaker form of Borel's summation method defines the Borel sum of ''A'' to be
| |
| | |
| :<math> \lim_{t\rightarrow\infty} e^{-t}\sum_{n=0}^\infty \frac{t^n}{n!}A_n(z). </math>
| |
| | |
| If this converges at ''z'' ∈ '''C''' to some ''a(z)'', we say that the weak Borel sum of ''A'' converges at ''z'', and write <math> {\textstyle \sum} a_kz^k = a(z) \, (\boldsymbol{wB}) </math>.
| |
| | |
| ==Basic properties==
| |
| | |
| ===Regularity===
| |
| The methods ('''B''') and ('''wB''') are both [[Divergent_series#Properties_of_summation_methods|regular]] summation methods, meaning that whenever ''A(z)'' converges (in the standard sense), then the Borel sum and weak Borel sum also converge, and do so to the same value. i.e.
| |
| | |
| :<math> \sum_{k=0}^\infty a_k z^k = A(z) < \infty \quad \Rightarrow \quad {\textstyle \sum} a_kz^k = A(z) \,\, (\boldsymbol{B},\,\boldsymbol{wB}). </math>
| |
| | |
| Regularity of ('''B''') is easily seen by a change in order of integration: if ''A(z)'' is convergent at ''z'',then
| |
| | |
| :<math> A(z) = \sum_{k=0}^\infty a_k z^k = \sum_{k=0}^\infty a^k \left( \int_{0}^\infty e^{-t}t^k dt \right) \frac{z^k}{k!} = \int_{0}^\infty e^{-t} \sum_{k=0}^\infty a_k \frac{(tz)^k}{k!}dt, </math>
| |
| | |
| where the rightmost expression is exactly the Borel sum at ''z''.
| |
| | |
| Regularity of ('''B''') and ('''wB''') imply that these methods provide analytic extensions to ''A(z)''.
| |
| | |
| ===Equivalence of Borel and weak Borel summation===
| |
| Any series ''A(z)'' that is weak Borel summable at ''z'' ∈ ''C'' is also Borel summable at ''z''. However, one can construct [[#An_example_in_which_equivalence_fails|examples]] of series which are divergent under weak Borel summation, but which are Borel summable. The following theorem characterises the equivalence of the two methods.
| |
| | |
| :'''Theorem''' ({{harv|Hardy|1992|loc=8.5}}).
| |
| :Let ''A(z)'' be a formal power series, and fix ''z'' ∈ ''C'', then:
| |
| :# If <math> {\textstyle \sum} a_kz^k = a(z) \, (\boldsymbol{wB}) </math>, then <math> {\textstyle \sum}a_kz^k = a(z) \, (\boldsymbol{B})</math>.
| |
| :# If <math> {\textstyle \sum} a_kz^k = a(z) \, (\boldsymbol{B}) </math>, and <math> \lim_{t \rightarrow \infty} e^{-t}\mathcal B A(zt) = 0, </math> then <math> {\textstyle \sum} a_kz^k = a(z) \, (\boldsymbol{wB}) </math>.
| |
| | |
| ===Relationship to other summation methods===
| |
| * ('''B''') is the special case of [[Mittag-Leffler summation]] with α = 1.
| |
| * ('''wB''') can be seen as the limiting case of generalized [[Euler summation|Euler summation method]] ('''E''',q) in the sense that as ''q'' → ∞ the domain of convergence of the ('''E''',q) method converges up to the domain of convergence for ('''B''').<ref name="Hardy1992">Hardy, G. H. (1992). ''Divergent Series''. AMS Chelsea, Rhode Island.</ref>
| |
| | |
| ==Uniqueness theorems==
| |
| | |
| There are always many different functions with any given asymptotic expansion. However there is sometimes a best possible function, in the sense that the errors in the finite dimensional epproximations are as small as possible in some region. Watson's theorem and Carleman's theorem show that Borel summation produces such a best possible sum of the series.
| |
| | |
| ===Watson's theorem===
| |
| | |
| Watson's theorem gives conditions for a function to be the Borel sum of its asymptotic series. Suppose that ''f'' is a function satisfying the following conditions:
| |
| *''f'' is holomorphic in some region |''z''|<''R'', |arg(''z'')| < π/2 + ε for some positive ''R'' and ε.
| |
| *In this region ''f'' has an asymptotic series ''a''<sub>0</sub> + ''a''<sub>1</sub>''z''+... with the property that the error
| |
| :<math>|f(z)-a_0 -a_1z -\cdots -a_{n-1}z^{n-1}|</math>
| |
| is bounded by | |
| :<math>C^{n+1}n!z^n</math>
| |
| for all ''z'' in the region (for some positive constant ''C'').
| |
| | |
| Then Watson's theorem says that in this region ''f'' is given by the Borel sum of its asymptotic series. More precisely, the series for the Borel transform converges in a neighborhood of the origin, and can be analytically continued to the positive real axis, and the integral defining the Borel sum converges to ''f''(''z'') for ''z'' in the region above.
| |
| | |
| Slightly more generally, ''f'' is given by the Borel sum of its asymptotic series if the ''n''! in the error estimate above is replaced by ''kn''! provided the condition |arg(''z'')| < π/2 + ε is replaced by |arg(''z'')| < ''k''π/2 + ε. This is in some sense best possible, as there are counterexamples if the number ''k''π/2 is replaced by any smaller number.
| |
| | |
| ===Carleman's theorem===
| |
| | |
| Carleman's theorem shows that a function is uniquely determined by an asymptotic series in a sector provided the errors in the finite order approximations do not grow too fast. More precisely it states that if ''f'' is analytic in the interior of the sector |''z''|<''C'', Re(''z'')>0 and |''f''(''z'')|<|''b''<sub>''n''</sub>''z''|<sup>''n''</sup> in this region for all ''n'', then ''f'' is zero provided that the series 1/''b''<sub>0</sub> + 1/''b''<sub>1</sub> + ... diverges.
| |
| | |
| Carleman's theorem gives a summation method for any asymptotic series whose terms do not grow too fast, as the sum can be defined to be the unique function with this asymptotic series in a suitable sector if it exists. Borel summation is slightly weaker than special case of this when ''b''<sub>''n''</sub> =''c''/''n'' for some constant ''c''. More generally one can define summation methods slightly stronger than Borel's by taking the numbers ''b''<sub>''n''</sub> to be slightly larger, for example ''b''<sub>''n''</sub> =''c''/''n''log ''n'' or ''b''<sub>''n''</sub> =''c''/''n''log ''n'' log log ''n''. In practice this generalization is of little use, as there are almost no natural examples of series summable by this method that cannot also be summed by Borel's method.
| |
| | |
| ===Example===
| |
| | |
| The function ''f''(''z'')=exp(–1/''z'') has the asymptotic series 0 + 0''z''+ ... with an error bound of the form above in the region |arg(''z'')|<θ for any θ<π/2, but is not given by the Borel sum of its asymptotic series. This shows that the number π/2 in Watsons theorem cannot be replaced by any smaller number (unless the bound on the error is made smaller).
| |
| | |
| ==Examples==
| |
| | |
| ===The geometric series===
| |
| Consider the [[geometric series]]
| |
| | |
| :<math>A(z) = \sum_{k = 0}^\infty z^k,</math>
| |
| | |
| which converges (in the standard sense) to 1/(1 − ''z'') for |''z''| < 1. The Borel transform is
| |
| | |
| :<math>\mathcal{B}A(t) \equiv \sum_{k=0}^\infty \frac{1}{k!}t^k = e^t,</math>
| |
| | |
| from which we obtain the Borel sum
| |
| | |
| :<math>\int_0^\infty e^{-t}\mathcal{B}A(tz) \, dt = \int_0^\infty e^{-t} e^{tz} \, dt =\frac{1}{1-z}</math>
| |
| | |
| which converges in the larger region Re(''z'') < 1, giving an [[analytic continuation]] of the original series.
| |
| | |
| Considering instead the weak Borel transform, the partial sums are given by ''A<sub>N</sub>(z) = (1-z<sup>N+1</sup>)/(1-z)'', and so the weak Borel sum is
| |
| | |
| :<math> \lim_{t \rightarrow \infty}e^{-t} \sum_{n=0}^\infty \frac{1 -z^{n+1}}{1-z} \frac{t^n}{n!} = \lim_{t \rightarrow \infty} \frac{e^{-t}}{1-z} \big( e^t - z e^{tz} \big) = \frac{1}{1-z}, </math>
| |
| | |
| where, again, convergence is on Re(''z'') < 1. Alternatively this can be seen by appealing to part 2 of the equivalence theorem, since for Re(''z'') < 1
| |
| | |
| :<math> \lim_{t \rightarrow \infty} e^{-t} (\mathcal{B} A)(zt) = e^{t(z-1)} = 0. </math>
| |
| | |
| ===An alternating factorial series===
| |
| Consider the series
| |
| | |
| :<math>A(z) = \sum_{k = 0}^\infty k!\left(-1 \cdot z\right)^k,</math>
| |
| | |
| then ''A(z)'' does not converge for any nonzero ''z'' ∈ '''C'''. The Borel transform is
| |
| | |
| :<math>\mathcal{B}A(t) \equiv \sum_{k=0}^\infty \left(-1 \cdot t\right)^k = \frac{1}{1+t} </math>
| |
| | |
| for |''t''| < 1, which can be analytically continued to all ''t'' ≥ 0. So the Borel sum is
| |
| | |
| :<math>\int_0^\infty e^{-t}\mathcal{B}A(tz) \, dt = \int_0^\infty \frac{e^{-t}} {1+tz} \, dt = \frac 1 z \cdot e^\frac 1 z \cdot \Gamma\left(0,\frac 1 z \right)</math>
| |
| | |
| (where Γ is the [[incomplete Gamma function]]).
| |
| | |
| This integral converges for all ''z'' ≥ 0, so the original divergent series is Borel summable for all such ''z''. This function has an [[asymptotic expansion]] as ''z'' tends to 0 that is given by the original divergent series. This is a typical example of the fact that Borel summation will sometimes "correctly" sum divergent asymptotic expansions.
| |
| | |
| Again, since
| |
| :<math> \lim_{t \rightarrow \infty}e^{-t} (\mathcal B A)(zt) = \lim_{t \rightarrow \infty} \frac{e^{-t}}{1 + zt} = 0, </math>
| |
| | |
| for all ''z'', the equivalence theorem ensures that weak Borel summation has the same domain of convergence, ''z'' ≥ 0.
| |
| | |
| ===An example in which equivalence fails===
| |
| | |
| The following example extends on that given in {{harv|Hardy|1992|loc=8.5}}. Consider
| |
| | |
| :<math>A(z) = \sum_{k = 0}^\infty \left(
| |
| \sum_{l=0}^\infty \frac{(-1)^l(2l + 2)^k}{(2l+1)!} \right)
| |
| z^k. </math>
| |
| | |
| After changing the order of summation, the Borel transform is given by
| |
| | |
| :<math>
| |
| \begin{align}
| |
| \mathcal B A(t)&= \sum_{l = 0}^\infty
| |
| \left( \sum_{k=0}^\infty \frac{\big((2l+2) t\big)^k}{k!} \right) \frac{(-1)^l}{(2l+1)!} \\
| |
| &= \sum_{l=0}^\infty e^{(2l+2)t}\frac{(-1)^l}{(2l+1)!} \\
| |
| &= e^t \sum_{l=0}^\infty \big(e^t\big)^{2l+1} \frac{(-1)^l}{(2l+1)!} \\
| |
| & = e^t \sin\left( e^t \right).
| |
| \end{align}
| |
| </math>
| |
| | |
| At ''z'' = 2 the Borel sum is given by
| |
| | |
| :<math> \int_{0}^\infty e^t \sin(e^{2t})dt = \int_{1}^\infty \sin(u^2)du = \frac{\sqrt{\pi}}{8} - S(1) < \infty,
| |
| </math>
| |
| | |
| where ''S(x)'' is the [[Fresnel integral]]. Via the [[#Convergence_Properties|convergence theorem]] along chords, the Borel integral converges for all ''z'' ≤ 2 (clearly the integral diverges for ''z'' > 2).
| |
| | |
| For the weak Borel sum we note that
| |
| | |
| :<math> \lim_{t \rightarrow \infty} e^{(z-1)t}\sin \left( e^{zt} \right) = 0 </math>
| |
| | |
| holds only for ''z'' < 1, and so the weak Borel sum converges on this smaller domain.
| |
| | |
| ==Existence results and the domain of convergence==
| |
| | |
| ===Summability on chords===
| |
| If a formal series ''A(z)'' is Borel summable at ''z<sub>0</sub>'' ∈ '''C''', then it is also Borel summable at all points on the chord ''Oz<sub>0</sub>'' connecting ''z<sub>0</sub>'' to the origin. Moreover, there exists a function ''a(z)'' analytic throughout the disk with radius ''Oz<sub>0</sub>'' such that
| |
| | |
| : <math> {\textstyle \sum} a_kz^k = a(z) \, (\boldsymbol B), </math>
| |
| | |
| for all ''z'' = θ''z''<sub>0</sub>, θ ∈ [0,1].
| |
| | |
| An immediate consequence is that the domain of convergence of the Borel sum is a [[star domain]] in '''C'''. More can be said about the domain of convergence of the Borel sum, than that it is a star domain, which is referred to as the Borel polygon, and is determined by the singularities of the series ''A(z)''.
| |
| | |
| ===The Borel polygon===
| |
| Suppose that ''A(z)'' has strictly positive radius of convergence, so that it is analytic in a non-trivial region containing the origin, and let ''S<sub>A</sub>'' denote the set of singularities of ''A''. This means that ''P'' ∈ ''S<sub>A</sub>'' if and only if ''A'' can be continued analytically along the open chord from ''0'' to ''P'', but not to ''P'' itself. For ''P'' ∈ ''S<sub>A</sub>'', let ''L<sub>P</sub>'' denote the line passing through ''P'' which is perpendicular to the chord ''OP''. Define the sets
| |
| | |
| :<math> \Pi_P = \{z \in \mathbb{C} \, \colon \, Oz \cap L_P = \varnothing \}, </math>
| |
| | |
| the set of points which lie on the same side of ''L<sub>P</sub>'' as the origin. The Borel polygon of ''A'' is the set
| |
| | |
| :<math> \Pi_A = \text{cl}\Big( \bigcap_{P \in S_A} \Pi_P \Big). </math>
| |
| | |
| An alternative definition was used by Borel and Phragmén {{harv|Sansone|Gerretsen|1960|loc=8.3}}. Let <math> S \subset \mathbb{C} </math> denote the largest star domain on which there is an analytic extension of ''A'', then <math> \Pi_A </math> is the largest subset of <math> S </math> such that for all <math> P \in \Pi_A </math> the interior of the circle with diameter ''OP'' is contained in <math>S </math>. Referring to the set <math> \Pi_A </math> as a polygon is somewhat of a misnomer, since the set need not be polygonal at all; if, however, A(z) has only finitely many singularities then <math> \Pi_A </math> will in fact be a polygon.
| |
| | |
| The following theorem, due to Borel and [[Lars Edvard Phragmén|Phragmén]] provides convergence criteria for Borel summation.
| |
| | |
| :'''Theorem''' {{harv|Hardy|1992|loc=8.8}}. | |
| :The series ''A(z)'' is ('''B''') summable at all <math>z \in \text{int}(\Pi_A)</math>, and is ('''B''') divergent at all <math>z \in \mathbb{C}\backslash \Pi_A</math>. | |
| | |
| Note that ('''B''') summability for <math> z \in \partial \Pi_A</math> depends on the nature of the point.
| |
| | |
| ====Example 1====
| |
| Let ω<sub>i</sub> ∈ '''C''' denote the ''m''-th roots of unity, ''i'' =1,...''m'', and consider
| |
| | |
| : <math>\begin{align}
| |
| A(z) & = \sum_{k=0}^\infty (\omega_1^k + \ldots + \omega_m^k)z^k \\
| |
| & = \sum_{i=1}^m \frac{1}{1-\omega_iz},
| |
| \end{align} </math>
| |
| | |
| which converges on ''B''(0,1) ⊂ '''C'''. Seen as a function on '''C''', ''A(z)'' has singularities at ''S<sub>A</sub>'' = {ω<sub>i</sub> : ''i'' = 1,...''m''}, and consequently the Borel polygon <math> \Pi_A</math> is given by the regular [[Regular polygon|m-gon]] centred at the origin, and such that 1 ∈ '''C''' is a midpoint of an edge.
| |
| | |
| ====Example 2====
| |
| The formal series
| |
| | |
| :<math> A(z) = \sum_{k=0}^\infty z^{2^k}, </math>
| |
| | |
| converges for all <math> |z| < 1 </math> (for instance, by the [[Direct comparison test|comparison test]] with the geometric series). It can however be shown<ref>http://mathworld.wolfram.com/NaturalBoundary.html</ref> that ''A'' does not converge for all points ''z'' ∈ '''C''' with ''z''<sup>2<sup>n</sup></sup> = 1. Since the set of such ''z'' is dense in the unit circle, there can be no analytic extension of ''A'' outside of ''B(0,1)''. Subsequently the largest star domain to which ''A'' can be analytically extended is ''S'' = ''B(0,1)'' from which (via the second definition) one obtains <math> \Pi_A = B(0,1) </math>. In particular we see that the Borel polygon is in fact not polygonal.
| |
| | |
| ===A Tauberian Theorem===
| |
| A [[Abelian_and_tauberian_theorems#Tauberian_theorems|Tauberian theorem]] provides conditions under which convergence of one summation method, implies convergence under another method. The principal Tauberian theorem<ref name="Hardy1992" /> for Borel summation provides conditions under which the weak Borel method implies convergence of the series.
| |
| | |
| :'''Theorem''' {{harv|Hardy|1992|loc=9.13}}. :If ''A'' is ('''wB''') summable at ''z<sub>0</sub>'' ∈ '''C''', <math>{\textstyle \sum}a_kz_0^k = a(z_0) \, (\boldsymbol{wB}) </math>, and
| |
| | |
| :: <math> a_kz_0^k = O \left( k^{-{\textstyle \frac{1}{2}}} \right), \qquad \forall k \geq 0, </math>
| |
| | |
| :then <math> \sum_{k=0}^\infty a_kz_0^k = a(z_0) </math>, and the series converges for all ''|z|<|z<sub>0</sub>|''.
| |
| | |
| ==Applications==
| |
| Borel summation finds application in [[perturbation theory (quantum mechanics)|perturbation expansions]] in quantum field theory. In particular in 2-dimensional Euclidean field theory the Schwinger functions can often be recovered from their perturbation series using Borel summation {{harv|Glimm|Jaffe|1987|p=461}}. Some of the singularities of the Borel transform are related to [[instanton]]s and [[renormalon]]s in quantum field theory {{harv|Weinberg|2005|loc=20.7}}.
| |
| | |
| ==Generalizations==
| |
| | |
| Borel summation requires that the coefficients do not grow too fast: more precisely, ''a''<sub>''n''</sub> has to be bounded by ''n''!''C''<sup>''n''+1</sup> for some ''C''. There is a variation of Borel summation that replaces factorials ''n''! with ''(kn)''! for some positive integer ''k'', which allows the summation of some series with ''a''<sub>''n''</sub> bounded by ''(kn)''!''C''<sup>''n''+1</sup> for some ''C''. This generalization is similar to ordinary Borel summation.
| |
| | |
| | |
| ==See also==
| |
| * [[Abel summation]]
| |
| * [[Abel's theorem]]
| |
| * [[Abel–Plana formula]]
| |
| * [[Euler summation]]
| |
| * [[Cesàro summation]]
| |
| * [[Lambert summation]]
| |
| * [[Nachbin resummation]]
| |
| * [[Abelian and tauberian theorems]]
| |
| * [[Van Wijngaarden transformation]]
| |
| | |
| ==Notes==
| |
| <references/>
| |
| | |
| ==References==
| |
| *{{citation |last=Borel |first=E. |title=Mémoire sur les séries divergentes |journal=Ann. Sci. École Norm. Sup. (3) |volume=16 |year=1899 |pages=9–131 |url=http://www.numdam.org/item?id=ASENS_1899_3_16__9_0}}
| |
| *{{Citation |last1=Glimm |first1=James |last2=Jaffe |first2=Arthur |title=Quantum physics |publisher=[[Springer-Verlag]] |location=Berlin, New York |edition=2nd |isbn=978-0-387-96476-8 |mr=887102 |year=1987}}
| |
| *{{Citation |last=Hardy |first=Godfrey Harold |authorlink=G. H. Hardy |title=Divergent Series |origyear=1949 | url=http://books.google.com/books?isbn=0821826492 |publisher=Chelsea |location=New York |isbn=978-0-8218-2649-2 | mr=0030620 |year=1992}}
| |
| *{{Citation |last1=Reed |first1=Michael |last2=Simon |first2=Barry |authorlink2=Barry Simon |title=Methods of modern mathematical physics. IV. Analysis of operators |publisher=Academic Press [Harcourt Brace Jovanovich Publishers] |location=New York |isbn=978-0-12-585004-9 |mr=0493421 |year=1978}}
| |
| *{{Citation |last1=Sansone |first1=Giovanni |last2=Gerretsen |first2=Johan |title=Lectures on the theory of functions of a complex variable. I. Holomorphic functions |publisher=P. Noordhoff, Groningen |mr=0113988 |year=1960}}
| |
| *{{Citation |last=Weinberg | first=Steven |authorlink=Steven Weinberg |title=The quantum theory of fields. |volume=II |publisher=[[Cambridge University Press]] |isbn=978-0-521-55002-4 |mr=2148467 |year=2005}}
| |
| *{{eom |id=B/b017170 |first=A. A. |last=Zakharov |title=Borel summation method}}
| |
| | |
| [[Category:Mathematical series]]
| |
| [[Category:Summability methods]]
| |
| [[Category:Quantum chromodynamics]]
| |