# Equiprobability, Entropy, Gamma Distributions and Other Geometrical Questions in Multi-Agent Systems

^{1}

^{2}

^{3}

^{*}

Next Article in Journal

Next Article in Special Issue

Next Article in Special Issue

Previous Article in Journal

Previous Article in Special Issue

Previous Article in Special Issue

Comment published on 22 December 2009,
see
*Entropy* **2009**, *11*(4), 1121-1122.

Department of Computer Science, Universidad de Zaragoza, Pedro Cerbuna 12, Zaragoza 50009, Spain

Department of Physics, Universidad de Extremadura, Avda. de Elvas s/n, Badajoz 06071, Spain

BIFI, Universidad de Zaragoza, Corona de Aragón 42, Zaragoza 50009, Spain

Author to whom correspondence should be addressed.

Received: 3 November 2009 / Accepted: 30 November 2009 / Published: 2 December 2009

(This article belongs to the Special Issue Information and Entropy)

A set of many identical interacting agents obeying a global additive constraint is considered. Under the hypothesis of equiprobability in the high-dimensional volume delimited in phase space by the constraint, the statistical behavior of a generic agent over the ensemble is worked out. The asymptotic distribution of that statistical behavior is derived from geometrical arguments. This distribution is related with the Gamma distributions found in several multi-agent economy models. The parallelism with all these systems is established. Also, as a collateral result, a formula for the volume of high-dimensional symmetrical bodies is proposed.

Nowadays different approaches to study economic systems can be taken within a multi-agent framework. Econophysics applies the techniques and tools of statistical physics to the understanding of this type of systems [1]. From this point of view, an economic system is regarded as a set of many agents exchanging money or trading with products of the most diverse origins. This exchange, which can be done in a deterministic or in a random way, finally generates an asymptotic wealth distribution. One of the goals of econophysics is to propose different models for the agents’ interactions and, by making the respective simulation experiments, to try to reproduce in a statistical way the real asymptotic wealth distribution in society [2].

In this paper, a general multi-agent system evolving under an additive constraint is considered. The geometrical properties of this system in phase space are exploited in order to obtain its statistical behavior. A striking coincidence is observed. Some dynamical mechanisms that have been proposed in the literature to model the interaction among agents in economic systems [3,4,5,6,7,8,9] provoke the same statistical results in the asymptotic wealth distribution than those derived from the geometrical properties in our multi-agent system under the assumption of equiprobability. This fact bring us to speculate on the possible close relationship that there exists between the local interactions among the agents in the former economic systems and the global geometrical conformation in phase space of our general system. Hence, with an adequate change of coordinates, we suggest that those economic systems asymptotically evolve in an equiprobable way over the volume of accessible states in the transformed phase space.

We start in Section 2 by recalling the derivation done for some particular geometries for the cases in which the constraint has a linear or quadratic dependence on the variables defining the agents. Then, in Section 3, the statistical behavior for a more general constraint is obtained. In Section 4, we speculate on the possible relationship with some economic systems [3,4,5,6,7,8,9] in which the Gamma distributions are also obtained. A formula for the volume of high-dimensional symmetrical bodies is proposed in Section 5. The last Section 6 contains our conclusions.

Let us assume N agents interacting in an open economy, each one with coordinate ${x}_{i}$, $i=1,\dots ,N$, with ${x}_{i}\ge 0$ representing the wealth or money of the agent i, and a total available amount of money E. The additive linear constraint reads:
The result here explained for an open economy is also obtained when constraint (1) is an equality, i.e., when the economy is closed [10]. Under random evolution rules for the exchanging of money among agents [3], let us suppose that this system evolves in the interior of the N-dimensional pyramid given by Equation (1). We can suppose that the state or the bank system of western societies plays in this model the role of a heat reservoir that supplies money instead of energy. The formula for the volume ${V}_{N}(E)$ of an equilateral N-dimensional pyramid formed by $N+1$ vertices linked by N perpendicular sides of length E is
If each point on the N-dimensional pyramid is equiprobable, then the probability $f({x}_{i})d{x}_{i}$ of finding the agent i with money ${x}_{i}$, with normalization condition ${\int}_{0}^{E}f({x}_{i})d{x}_{i}=1$, is proportional to the volume formed by all the points into the $(N-1)$-dimensional pyramid having the ith-coordinate equal to ${x}_{i}$. Then $f({x}_{i})$ verifies
If we call ϵ the mean wealth per agent, $E=N\u03f5$, then in the limit of large N ($N\to \infty $), we have
where the index i has been removed because the distribution is the same for each agent, and thus the wealth distribution can be obtained by averaging over all the agents,
This Boltzmann-Gibbs distribution has been found to fit the real distribution of incomes in western societies [3].

$${x}_{1}+{x}_{2}+\cdots +{x}_{N-1}+{x}_{N}\le E$$

$${V}_{N}(E)=\frac{{E}^{N}}{N!}$$

$$f({x}_{i})=\frac{{V}_{N-1}(E-{x}_{i})}{{V}_{N}(E)}$$

$$\underset{N\gg 1}{lim}\frac{{V}_{N-1}(E-x)}{{V}_{N}(E)}=\frac{1}{\u03f5}\phantom{\rule{0.166667em}{0ex}}{e}^{-x/\u03f5}$$

$$f(x)dx={\u03f5}^{-1}\phantom{\rule{0.166667em}{0ex}}{e}^{-x/\u03f5}dx$$

Now let us suppose a one-dimensional ideal gas of N non-identical classical particles with masses ${m}_{i}$, with $i=1,\dots ,N$, and total maximum energy E. If particle i has a momentum ${m}_{i}{v}_{i}$, we define a kinetic energy:
where ${p}_{i}$ is the square root of the kinetic energy. Then the quadratic constraint reads:
The distribution for ${p}_{i}$ here derived for an open system is also obtained when constraint (7) is an equality, i.e., when the energy is fixed [11]. When a finite number of particles are present in the system the asymptotic distribution maximizes the Tsallis entropy [12]. In the case with constraint (7) the system has accessible states with different energy, which is supposed to be supplied by a heat reservoir. These states are all those enclosed into the volume of the N-sphere given by Equation (7), with radius ${E}^{1/2}$. The formula for the volume ${V}_{N}(R)$ of an N-sphere of radius R is
where $\Gamma (\xb7)$ is the Gamma function. If we suppose that each point into the N-sphere is equiprobable, then the probability $f({p}_{i})d{p}_{i}$ of finding the particle i, with coordinate ${p}_{i}$ (energy ${p}_{i}^{2}$) and normalization condition ${\int}_{-R}^{R}f({p}_{i})d{p}_{i}=1$, is proportional to the volume formed by all the points into the $(N-1)$-sphere having the ith-coordinate equal to ${p}_{i}$. Then $f({p}_{i})$ verifies

$$K\equiv {p}_{i}^{2}\equiv \frac{1}{2}{m}_{i}{v}_{i}^{2}$$

$${p}_{1}^{2}+{p}_{2}^{2}+\cdots +{p}_{N-1}^{2}+{p}_{N}^{2}\le E$$

$${V}_{N}(R)=\frac{{\pi}^{\frac{N}{2}}}{\Gamma (\frac{N}{2}+1)}{R}^{N}$$

$$f({p}_{i})=\frac{{V}_{N-1}({(E-{p}_{i}^{2})}^{1/2})}{{V}_{N}({E}^{1/2})}$$

If we call ϵ the mean energy per particle, $E=N\u03f5$, then in the limit of large N ($N\to \infty $), we have
where the index i has been removed because the distribution is the same for each particle. Thus the asymptotic distribution
can be obtained by averaging over all the particles. If the change of variables $p=\sqrt{\frac{m}{2}}\phantom{\rule{0.166667em}{0ex}}v$ is performed, with v the generic velocity of a particle, then the Maxwellian distribution is just derived from geometrical arguments.

$$\underset{N\gg 1}{lim}\frac{{V}_{N-1}({(E-{p}^{2})}^{1/2})}{{V}_{N}({E}^{1/2})}=\sqrt{\frac{1}{2\pi}}\phantom{\rule{0.166667em}{0ex}}{\u03f5}^{-1/2}\phantom{\rule{0.166667em}{0ex}}{e}^{-{p}^{2}/2\u03f5}$$

$$f(p)=\sqrt{\frac{1}{2\pi}}\phantom{\rule{0.166667em}{0ex}}{\u03f5}^{-1/2}\phantom{\rule{0.166667em}{0ex}}{e}^{-{p}^{2}/2\u03f5}$$

In this section, we address the same problem above presented but in a general way. Let b be a positive real constant (cases $b=1,2$ have been indicated in the former section). If we have a set of positive variables $({x}_{1},{x}_{2},\dots ,{x}_{N})$ verifying the constraint
with an adequate mechanism assuring the equiprobability of all the possible states $({x}_{1},{x}_{2},\dots ,{x}_{N})$ into the volume given by expression (12), will we have for the generic variable x the distribution
when we average over the ensemble in the limit $N,E\to \infty $, with $E=N\u03f5$, and constant ϵ?. Now it is shown that the answer is affirmative. Similarly, we claim that if the weak inequality (12) is transformed in equality the result will be the same, as it has been proved for the cases $b=1,2$ in References [10,11].

$${x}_{1}^{b}+{x}_{2}^{b}+\cdots +{x}_{N-1}^{b}+{x}_{N}^{b}\le E$$

$$f(x)dx\sim {\u03f5}^{-1/b}\phantom{\rule{0.166667em}{0ex}}{e}^{-{x}^{b}/b\u03f5}dx$$

From the cases $b=1,2$, (see Equations (3) and (9)), we can extrapolate the general formula that will give us the statistical behavior $f(x)$ of the generic variable x, when the system runs equiprobably into the volume defined by a constraint of type (12). The probability $f(x)dx$ of finding an agent with generic coordinate x is proportional to the volume ${V}_{N-1}({(E-{x}^{b})}^{1/b})$ formed by all the points into the $(N-1)$-dimensional symmetrical body limited by the constraint $(E-{x}^{b})$. Thus, the N-dimensional volume can be written as
Taking into account the normalization condition ${\int}_{0}^{{E}^{1/b}}f(x)dx=1$, the expression for $f(x)$ is obtained:

$${V}_{N}({E}^{1/b})={\int}_{0}^{{E}^{1/b}}\phantom{\rule{0.166667em}{0ex}}{V}_{N-1}({(E-{x}^{b})}^{1/b})\phantom{\rule{0.166667em}{0ex}}dx$$

$$f(x)=\frac{{V}_{N-1}({(E-{x}^{b})}^{1/b})}{{V}_{N}({E}^{1/b})}$$

The N-dimensional volume, ${V}_{N}(b,\rho )$, of a b-symmetrical body with side of length ρ is proportional to the term ${\rho}^{N}$ and to a coefficient ${g}_{b}(N)$ that depends on N:
The parameter b indicates the original Equation (12) that defines the boundaries of the volume ${V}_{N}(b,\rho )$. Thus, for instance, from Equation (2), we have ${g}_{b=1}(N)=1/N!$.

$${V}_{N}(b,\rho )={g}_{b}(N)\phantom{\rule{0.166667em}{0ex}}{\rho}^{N}$$

Coming back to Equation (15), we can manipulate ${V}_{N}({(E-{x}^{b})}^{1/b})$ to obtain (the index b is omitted in the formula of ${V}_{N}$):
If we suppose $E=N\u03f5$, then ϵ represents the mean value of ${x}^{b}$ in the collectivity, that is, $\u03f5=<{x}^{b}>$. If N tends toward infinity, it results:
Thus,
Substituting this last expression in formula (15), the exact form for $f(x)$ is found in the thermodynamic limit ($N,E\to \infty $):
with ${c}_{b}$ given by
Hence, the conjecture (13) is proved.

$${V}_{N}({(E-{x}^{b})}^{1/b})={g}_{b}(N)\phantom{\rule{0.166667em}{0ex}}{\left[{(E-{x}^{b})}^{1/b}\right]}^{N}={g}_{b}(N)\phantom{\rule{0.166667em}{0ex}}{E}^{\frac{N}{b}}\phantom{\rule{4pt}{0ex}}{\left(1-\frac{{x}^{b}}{E}\right)}^{\frac{N}{b}}$$

$$\underset{N\gg 1}{lim}{\left(1-\frac{{x}^{b}}{E}\right)}^{\frac{N}{b}}\phantom{\rule{0.166667em}{0ex}}=\phantom{\rule{0.166667em}{0ex}}{e}^{-{x}^{b}/b\u03f5}$$

$${V}_{N}({(E-{x}^{b})}^{1/b})={V}_{N}({E}^{1/b})\phantom{\rule{0.166667em}{0ex}}{e}^{-{x}^{b}/b\u03f5}$$

$$f(x)dx={c}_{b}\phantom{\rule{0.166667em}{0ex}}{\u03f5}^{-1/b}\phantom{\rule{0.166667em}{0ex}}{e}^{-{x}^{b}/b\u03f5}dx$$

$${c}_{b}\phantom{\rule{0.166667em}{0ex}}=\phantom{\rule{0.166667em}{0ex}}\frac{{g}_{b}(N-1)}{{g}_{b}(N)\phantom{\rule{0.166667em}{0ex}}{N}^{1/b}}$$

Doing a thermodinamical simile, we can calculate the dependence of ϵ on the temperature by differentiating the entropy with respect to the energy. The entropy can be written as $S=-kN\phantom{\rule{-0.166667em}{0ex}}{\int}_{0}^{\infty}f(x)lnf(x)\phantom{\rule{0.166667em}{0ex}}dx$, where $f(x)$ is given by Equation (20) and k is the Boltzmann constant. If we recall that $\u03f5=E/N$, we obtain
where it has been used that $\u03f5=<{x}^{b}>=\phantom{\rule{-0.166667em}{0ex}}{\int}_{0}^{\infty}{x}^{b}f(x)dx$

$$S(E)=\frac{kN}{b}ln\left(\frac{E}{N}\right)+\frac{kN}{b}(1-bln{c}_{b})$$

The calculation of the temperature T gives
Thus $\u03f5=kT/b$, a result that recovers the theorem of equipartition of energy for the quadratic case $b=2$. The distribution for all b is finally obtained:

$${T}^{-1}={\left(\frac{\partial S}{\partial E}\right)}_{N}=\frac{kN}{bE}=\frac{k}{b\u03f5}$$

$$f(x)dx={c}_{b}{\left(\frac{b}{kT}\right)}^{1/b}\phantom{\rule{0.166667em}{0ex}}{e}^{-{x}^{b}/kT}dx$$

If we perform the change of variables $y={\u03f5}^{-1/b}x$ in the normalization condition of $f(x)$, $\phantom{\rule{-0.166667em}{0ex}}{\int}_{0}^{\infty}f(x)dx=1$, where $f(x)$ is expressed in (20), we find that
If we introduce the new variable $z={y}^{b}/b$, the distribution $f(x)$ as function of z reads:
Let us observe that the Gamma function appears in the normalization condition,
This implies that
By using Mathematica the positive constant ${c}_{b}$ is plotted versus b in Figure 1. We see that ${lim}_{b\to 0}{c}_{b}=\infty $, and that ${lim}_{b\to \infty}{c}_{b}=1$. The minimum of ${c}_{b}$ is reached for $b=3.1605$, taking the value ${c}_{b}=0.7762$. Still further, we can calculate from Equation (28) the asymptotic dependence of ${c}_{b}$ on b:
where γ is the Euler constant, $\gamma =0.5772$. The asymptotic function (29) is obtained after substituting in (28) the value of $\Gamma (1/b)$ by $(1/b-1)!$, and performing the Stirling approximation on this last expression, knowing that $1/b\to \infty $. The function (30) is found after looking for the first Taylor expansion terms of the Gamma function around the origin $x=0$. They can be derived from the Euler’s reflection formula, $\Gamma (x)\Gamma (1-x)=\pi /sin(\pi x)$. We obtain $\Gamma (x\to 0)={x}^{-1}+{\Gamma}^{\prime}(1)+\cdots $. From here, recalling that ${\Gamma}^{\prime}(1)=-\gamma $, we get $\Gamma (1/b)=b-\gamma +\cdots $, when $b\to \infty $. Although this last term of the Taylor expansion, $-\gamma $, is negligible we maintain it in expression (30). The only minimum of ${c}_{b}$ is reached for the solution $b=3.1605$ of the equation $\psi (1/b)+logb+b-1=0$, where $\psi (\xb7)$ is the digamma function (see Figure 1).

$${c}_{b}={\left[\phantom{\rule{-0.166667em}{0ex}}{\int}_{0}^{\infty}{e}^{-{y}^{b}/b}\phantom{\rule{0.166667em}{0ex}}dy\right]}^{-1}$$

$$f(z)dz=\frac{{c}_{b}}{{b}^{1-\frac{1}{b}}}\phantom{\rule{0.166667em}{0ex}}{z}^{\frac{1}{b}-1}\phantom{\rule{0.166667em}{0ex}}{e}^{-z}\phantom{\rule{0.166667em}{0ex}}dz$$

$${\int}_{0}^{\infty}f(z)dz=\frac{{c}_{b}}{{b}^{1-\frac{1}{b}}}\phantom{\rule{0.166667em}{0ex}}{\int}_{0}^{\infty}\phantom{\rule{0.166667em}{0ex}}{z}^{\frac{1}{b}-1}\phantom{\rule{0.166667em}{0ex}}{e}^{-z}\phantom{\rule{0.166667em}{0ex}}dz=\frac{{c}_{b}}{{b}^{1-\frac{1}{b}}}\phantom{\rule{0.166667em}{0ex}}\Gamma \left(\frac{1}{b}\right)=1$$

$${c}_{b}=\frac{{b}^{1-\frac{1}{b}}}{\Gamma \left(\frac{1}{b}\right)}$$

$$\begin{array}{ccc}\hfill \underset{b\to 0}{lim}{c}_{b}& =& \sqrt{\frac{1}{2\pi}}\phantom{\rule{0.166667em}{0ex}}\sqrt{b}\phantom{\rule{0.166667em}{0ex}}{e}^{1/b}\left(1-\frac{b}{12}+\cdots \right)\hfill \end{array}$$

$$\begin{array}{ccc}\hfill \underset{b\to \infty}{lim}{c}_{b}& =& {b}^{-1/b}\left(1+\frac{\gamma}{b}+\cdots \right)\hfill \end{array}$$

Let us now recall two interesting statistical economic models that display a statistical behavior given by distributions of the form (26), that is, the standard Gamma distributions with shape parameter $1/b$,

$$f(z)dz=\frac{1}{\Gamma (\frac{1}{b})}\phantom{\rule{0.166667em}{0ex}}{z}^{\frac{1}{b}-1}\phantom{\rule{0.166667em}{0ex}}{e}^{-z}\phantom{\rule{0.166667em}{0ex}}dz$$

$$\begin{array}{ccc}\hfill {u}_{i}^{\prime}& =& \lambda {u}_{i}+\u03f5(1-\lambda )({u}_{i}+{u}_{j})\hfill \end{array}$$

$$\begin{array}{ccc}\hfill {u}_{j}^{\prime}& =& \lambda {u}_{j}+\overline{\u03f5}(1-\lambda )({u}_{i}+{u}_{j})\hfill \end{array}$$

$$n(\lambda )=\frac{1+2\lambda}{1-\lambda}$$

$$f(\overline{z})d\overline{z}=\frac{1}{\Gamma (n)}\phantom{\rule{0.166667em}{0ex}}{\overline{z}}^{n-1}\phantom{\rule{0.166667em}{0ex}}{e}^{-\overline{z}}\phantom{\rule{0.166667em}{0ex}}d\overline{z}$$

$$\begin{array}{ccc}\hfill n& =& 1/b\hfill \end{array}$$

$$\begin{array}{ccc}\hfill \overline{z}& =& z\hfill \end{array}$$

$$b=\frac{1-\lambda}{1+2\lambda}$$

$$x={\left[\frac{\u03f5}{<u>}\phantom{\rule{0.277778em}{0ex}}u\phantom{\rule{0.277778em}{0ex}}\right]}^{1/b}$$

$$\begin{array}{ccc}\hfill {u}_{i}^{\prime}& =& {u}_{i}-\Delta u\hfill \end{array}$$

$$\begin{array}{ccc}\hfill {u}_{j}^{\prime}& =& {u}_{j}+\Delta u\hfill \end{array}$$

$$\Delta u=\eta ({x}_{i}-{x}_{j})\phantom{\rule{0.166667em}{0ex}}\u03f5\omega {x}_{i}-[1-\eta ({x}_{i}-{x}_{j})]\phantom{\rule{0.166667em}{0ex}}\u03f5\omega {x}_{j}$$

$$n(\omega )=\frac{3-2\omega}{2\omega}$$

$$\begin{array}{ccc}\hfill n& =& 1/b\hfill \end{array}$$

$$\begin{array}{ccc}\hfill \overline{z}& =& z\hfill \end{array}$$

$$b=\frac{2\omega}{3-2\omega}$$

We shall proceed now to derive an asymptotic formula ($N\to \infty $) for the volume of the N-dimensional symmetrical body enclosed by the surface
The linear dimension ρ of this volume, i.e., the length of one of its sides verifies $\rho \sim {E}^{1/b}$. As argued in Equation (16), the N-dimensional volume, ${V}_{N}(b,\rho )$, is proportional to the term ${\rho}^{N}$ and to a coefficient ${g}_{b}(N)$ that depends on N. Thus,
where the characteristic b indicates the particular boundary given by Equation (49).

$${x}_{1}^{b}+{x}_{2}^{b}+\cdots +{x}_{N-1}^{b}+{x}_{N}^{b}=E$$

$${V}_{N}(b,\rho )={g}_{b}(N)\phantom{\rule{0.166667em}{0ex}}{\rho}^{N}$$

For instance, from Equation (2), we can write in a formal way:
From Equation (8), if we take the diameter, $\rho =2R$, as the linear dimension of the N-sphere, we obtain:
These expressions (51) and (52) suggest a possible general formula for the factor ${g}_{b}(N)$, let us say
where a is a b-dependent constant to be determined. For example, $a=1$ for $b=1$ and $a=\pi /4$ for $b=2$.

$${g}_{b=1}(N)=\frac{{1}^{\frac{N}{1}}}{\Gamma (\frac{N}{1}+1)}$$

$${g}_{b=2}(N)=\frac{{\left(\frac{\pi}{4}\right)}^{\frac{N}{2}}}{\Gamma \left(\frac{N}{2}+1\right)}$$

$${g}_{b}(N)=\frac{{a}^{\frac{N}{b}}}{\Gamma \left(\frac{N}{b}+1\right)}$$

In order to find the dependence of a on the parameter b, the regime $N\to \infty $ is supposed. Applying Stirling approximation for the factorial $(\frac{N}{b})!$ in the denominator of expression (53), and inserting it in expression (21), it is straightforward to find out the relationship:
From here and formula (28), we get:
that recovers the exact results for $b=1,2$. The behavior of a is monotonous decreasing when b is varied from $b=0$, where a diverges as $a\sim 1/b+\cdots $, up to the limit $b\to \infty $, where a decays asymptotically toward the value ${a}_{\infty}={e}^{-\gamma}=0.5614$.

$${c}_{b}\phantom{\rule{0.166667em}{0ex}}=\phantom{\rule{0.166667em}{0ex}}{(ab)}^{-1/b}$$

$$a\phantom{\rule{0.166667em}{0ex}}=\phantom{\rule{0.166667em}{0ex}}{\left[\Gamma \left(\frac{1}{b}+1\right)\right]}^{b}$$

Hence, the formula for ${g}_{b}(N)$ is obtained:
It would be also possible to multiply this last expression (56) by a general polynomial $K(N)$ in the variable N, and all the derivation done from Equation (53) would continue to be correct. We omit this possibility in our calculations. For a fixed N, we have that ${g}_{b}(N)$ increases monotonously from ${g}_{b}(N)=0$, for $b=0$, up to ${g}_{b}(N)=1$, in the limit $b\to \infty $ (see Figure 2). For a fixed b, we have that ${g}_{b}(N)$ decreases monotonously from ${g}_{b}(N)=1$, for $N=1$, up to ${g}_{b}(N)=0$, in the limit $N\to \infty $ (see Figure 3).

$${g}_{b}(N)\phantom{\rule{0.166667em}{0ex}}=\phantom{\rule{0.166667em}{0ex}}\frac{\phantom{\rule{0.166667em}{0ex}}\phantom{\rule{0.166667em}{0ex}}\Gamma {\left(\frac{1}{b}+1\right)}^{N}}{\Gamma \left(\frac{N}{b}+1\right)}$$

The final result for the volume of an N-dimensional symmetrical body of characteristic b given by the boundary (49) reads:
with $\rho \sim {E}^{1/b}$.

$${V}_{N}(b,\rho )\phantom{\rule{0.166667em}{0ex}}=\phantom{\rule{0.166667em}{0ex}}\frac{\phantom{\rule{0.166667em}{0ex}}\phantom{\rule{0.166667em}{0ex}}\Gamma {\left(\frac{1}{b}+1\right)}^{N}}{\Gamma \left(\frac{N}{b}+1\right)}\phantom{\rule{0.166667em}{0ex}}{\rho}^{N}$$

In this work, we have considered a general multi-agent open system verifying an additive constraint. Its statistical behavior has been derived from geometrical arguments. The Maxwellian and the Boltzmann-Gibbs distributions are particular cases of this type of systems. Also, other multi-agent economy models, such as the Dragalescu and Yakovenko’s model [3], the Chakraborti and Chakrabarti’s model [4] and the modified Angle’s model [6], show similar statistical behaviors compared to our general system. This fact fosters a geometrical interpretation of all those models. This geometrical speculation allows us to suggest the equivalence with the Chakraborti and Chakrabarti’s model when the geometrical characteristic b of our model runs in the interval $(0,1)$. The equivalence with the modified Angle’s model is suggested when b varies in the interval $(0,2)$. As a particular case of both types of models, the Dragulescu and Yakovenko’s model is obtained for $b=1$.

We have not found in the literature other multi-agent models to establish an equivalence with our system in the range $b\in (2,\infty )$. This point remains an open question and a challenge that will probably trigger other works in this direction.

The authors acknowledge some financial support from Spanish DGICYT Projects FIS2006-12781-C02-01.

- Mantegna, R.; Stanley, H.E. An Introduction to Econophysics: Correlations and Complexity in Finance; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
- Yakovenko, V.M. Econophysics, Statistical Mechanics Approach to. In Encyclopedia of Complexity and System Science; Meyers, R.A., Ed.; Springer: Berlin, Germany, 2009; pp. 2800–2826. [Google Scholar]
- Dragulescu, A.; Yakovenko, V.M. Statistical mechanics of money. Eur. Phys. J.
**2000**, B17, 723–729. [Google Scholar] [CrossRef] - Chakraborti, A.; Chakrabarti, B.K. Statistical mechanics of money: How saving propensity affects its distribution. Eur. Phys. J. B
**2000**, 17, 167–170. [Google Scholar] [CrossRef] - Patriarca, M.; Chakraborti, A.; Kaski, K. Statistical model with a standard Gamma distribution. Phys. Rev. E
**2004**, 70, 016104–016105. [Google Scholar] [CrossRef] - Patriarca, M.; Heinsalu, E.; Chakraborti, A. The ABCD’s of statistical many-agent economy models. 2006; arXiv:physics/0611245. [Google Scholar]
- Angle, J. The inequality process as a wealth maximizing process. Physica A
**2006**, 367, 388–414, (references therein). [Google Scholar] [CrossRef] - Gonzalez-Estevez, J.; Cosenza, M.G.; Lopez-Ruiz, R.; Sanchez, J.R. Pareto and Boltzmann-Gibbs behaviors in a deterministic multi-agent system. Physica A
**2008**, 387, 4637–4642. [Google Scholar] [CrossRef] - Pellicer-Lostao, C.; Lopez-Ruiz, R. Economic models with chaotic money exchange. In Proceedings of the ICCS 2009, Baton Rouge, LA, USA, 2009; Part I. pp. 43–52.
- Lopez-Ruiz, R.; Sañudo, J.; Calbet, X. Geometrical derivation of the Boltzmann factor. Am. J. Phys.
**2008**, 76, 780–781. [Google Scholar] [CrossRef] - Lopez-Ruiz, R.; Calbet, X. Derivation of the Maxwellian distribution from the microcanonical ensemble. Am. J. Phys.
**2007**, 75, 752–753. [Google Scholar] [CrossRef] - Harremoës, P. Projections maximizing Tsallis entropy. AIP Conf. Proc.
**2007**, 965, 90–95. [Google Scholar]

© 2009 by the authors; licensee Molecular Diversity Preservation International, Basel, Switzerland. This article is an open-access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).