Home | About | Journals | Submit | Contact Us | Français |

**|**Springer Open Choice**|**PMC5539268

Formats

Article sections

- Abstract
- Introduction
- The case ≠ 1: existence of the solution
- The case = 1
- The case 0 < < 1: existence and uniqueness of the solution
- References

Authors

Related links

Journal of Inequalities and Applications

J Inequal Appl. 2017; 2017(1): 178.

Published online 2017 August 1. doi: 10.1186/s13660-017-1454-4

PMCID: PMC5539268

Anca Ciurte, Email: or.julctu.sc@etruiC.acnA.

Received 2017 March 29; Accepted 2017 July 17.

Copyright © The Author(s) 2017

We are concerned with the positive solutions of an algebraic system depending on a parameter *α* > 0 and arising in economics. For *α* > 1 we prove that the system has at least a solution. For 0 < *α* < 1 we give three proofs of the existence and a proof of the uniqueness of the solution. Brouwer’s theorem and inequalities involving convex functions are essential tools in our proofs.

Algebraic systems with positive solutions appear in a large variety of applications. Nonlinear systems of the form *F*(*x*) = *A**x* + *p*, $x\in {\mathbb{R}}_{+}^{n}$, where *A* is a positive matrix and *p* a non-negative vector, are investigated in [1]. Generalizing some results from [2], the existence of a positive solution is proved with Brouwer’s theorem, and the uniqueness is a consequence of suitable inequalities. Algorithms for finding the solution are also presented in [1]. Several problems can be converted into systems of this form (see also [2]): second order Dirichlet problems, Dirichlet problems for partial difference equations, third and fourth order difference equations, three point boundary value problems, steady states of complex dynamical networks, etc.

Extending the results of [1], systems of the more general form ${\gamma}_{i}({x}_{i})={\sum}_{j=1}^{n}{g}_{ij}({x}_{j})$, 1 ≤ *i* ≤ *n*, are studied in [3], together with a supplementary list of applications. The existence of a positive solution is proved by using a monotone iterative method; the proof of the uniqueness is based on an extension of the method used in [1].

Several classes of other systems and several methods to investigate the existence/uniqueness of their positive solutions are described, together with applications, in [3–13] and the references therein.

In this paper we consider an algebraic system which appears in some problems from economics, for example in establishing uniqueness of equilibrium of some models of trade with increasing returns: see [14].

Our proofs use Brouwer’s fixed point theorem and properties of minimum points of convex functions. Several inequalities, in particular inequalities related to convexity, are instrumental in these proofs.

The system addressed in this paper is described in what follows. Let *N* ≥ 2 be an integer. Given the real numbers *α* > 0, *a*_{ni} > 0, *b*_{n} > 0, *i*, *n* ∈ {1, …, *N*}, consider the system of equations

$$\sum _{n=1}^{N}{b}_{n}{a}_{ni}{t}_{i}^{\alpha -1}{\left(\sum _{j=1}^{N}{a}_{nj}{t}_{j}^{\alpha}\right)}^{-1}=1,\phantom{\rule{1em}{0ex}}i=1,\dots ,N.$$

1

We are interested in solutions *t* = (*t*_{1}, …, *t*_{N}) ∈ *D* where *D*: = {*t* ∈ ℝ^{N}|*t*_{1} > 0, …, *t*_{N} > 0}.

The existence and uniqueness of such a solution depend essentially on *α*. We shall see that for *α* ≠ 1 there exists at least one solution; moreover, for 0 < *α* < 1 the solution is unique. For *α* = 1 there exist simple examples of systems (1) with no solution in *D*, or with exactly one solution, or with several solutions in *D*.

*If*
*α* > 0, *α* ≠ 1, *the system* (1) *has at least one solution in*
*D*.

For *j* ∈ {2, …, *N*}, let *t*_{j} = *s*_{j}*t*_{1}. Then (1) is equivalent to the system

$$\{\begin{array}{ll}{\sum}_{n=1}^{N}{b}_{n}{a}_{n1}{({a}_{n1}+{\sum}_{j=2}^{N}{a}_{nj}{s}_{j}^{\alpha})}^{-1}={t}_{1},& \\ {\sum}_{n=1}^{N}{b}_{n}{a}_{ni}{s}_{i}^{\alpha -1}{({a}_{n1}+{\sum}_{j=2}^{N}{a}_{nj}{s}_{j}^{\alpha})}^{-1}={t}_{1},& i=2,\dots ,N,\end{array}$$

2

with respect to the unknowns *t*_{1} > 0, *s*_{2} > 0, …, *s*_{N} > 0.

From (2) we infer that *s*_{2}, …, *s*_{N} satisfy

$$\begin{array}{rl}& {s}_{i}^{\alpha -1}\sum _{n=1}^{N}{b}_{n}{a}_{ni}{({a}_{n1}+\sum _{j=2}^{N}{a}_{nj}{s}_{j}^{\alpha})}^{-1}\\ & \phantom{\rule{1em}{0ex}}=\sum _{n=1}^{N}{b}_{n}{a}_{n1}{({a}_{n1}+\sum _{j=2}^{N}{a}_{nj}{s}_{j}^{\alpha})}^{-1},\phantom{\rule{1em}{0ex}}i=2,\dots ,N.\end{array}$$

3

Let *s*: = (*s*_{2}, …, *s*_{N}) and, for *i* = 2, …, *N*,

$$\begin{array}{rl}{F}_{i}(s):=& {\left(\sum _{n=1}^{N}{b}_{n}{a}_{ni}{({a}_{n1}+\sum _{j=2}^{N}{a}_{nj}{s}_{j}^{\alpha})}^{-1}\frac{{a}_{n1}}{{a}_{ni}}\right)}^{1/(\alpha -1)}\\ & \times {\left(\sum _{n=1}^{N}{b}_{n}{a}_{ni}{({a}_{n1}+\sum _{j=2}^{N}{a}_{nj}{s}_{j}^{\alpha})}^{-1}\right)}^{-1/(\alpha -1)}.\end{array}$$

Let ${U}_{i}:=\{\frac{{a}_{n1}}{{a}_{ni}}|n=1,\dots ,N\}$, *m*_{i}: = min *U*_{i}, *M*_{i} = max *U*_{i}, *i* = 2, …, *N*.

Then ${F}_{i}^{\alpha -1}(s)$ is a weighted mean of the numbers from *U*_{i}, and so

$${m}_{i}\le {F}_{i}^{\alpha -1}(s)\le {M}_{i},\phantom{\rule{1em}{0ex}}i=2,\dots ,N,$$

5

for all *s* = (*s*_{2}, …, *s*_{N}) such that *s*_{2} > 0, …, *s*_{N} > 0.

If *α* > 1, let $V:={\prod}_{i=2}^{N}[{m}_{i}^{1/(\alpha -1)},{M}_{i}^{1/(\alpha -1)}]$; if 0 < *α* < 1, let $V:={\prod}_{i=2}^{N}[{M}_{i}^{1/(\alpha -1)},{m}_{i}^{1/(\alpha -1)}]$.

In any case, due to (5) we can consider the continuous function *F*:*V* → *V*. Since *V* is compact and convex, Brouwer’s theorem guarantees the existence of a solution *s*^{0} ∈ *V* to (4). Then ${s}_{2}^{0},\dots ,{s}_{N}^{0}$ will satisfy (3). Now let

$${t}_{1}^{0}:=\sum _{n=1}^{N}{b}_{n}{a}_{n1}{({a}_{n1}+\sum _{j=2}^{N}{a}_{nj}{\left({s}_{j}^{0}\right)}^{\alpha})}^{-1}.$$

Then $({t}_{1}^{0},{s}_{2}^{0},\dots ,{s}_{N}^{0})$ is a solution to (2), and consequently ${t}^{0}:=({t}_{1}^{0},{s}_{2}^{0}{t}_{1}^{0},\dots ,{s}_{N}^{0}{t}_{1}^{0})\in D$ is a solution to (1).□

In the next example we present a system (1) with three solutions in *D*.

Let *N* = 2, *α* = 2, *b*_{1} = *b*_{2} = 1, *a*_{11} = *a*_{22} = 6, *a*_{12} = *a*_{21} = 1. Then (1) becomes

$$\{\begin{array}{l}\frac{6{t}_{1}}{6{t}_{1}^{2}+{t}_{2}^{2}}+\frac{{t}_{1}}{{t}_{1}^{2}+6{t}_{2}^{2}}=1,\\ \frac{{t}_{2}}{6{t}_{1}^{2}+{t}_{2}^{2}}+\frac{6{t}_{2}}{{t}_{1}^{2}+6{t}_{2}^{2}}=1,\end{array}$$

6

and has the solutions (1, 1), $(\frac{8}{7},\frac{6}{7}),(\frac{6}{7},\frac{8}{7})\in D$.

In the next three examples we provide systems (1) with (i) one solution in *D*, (ii) no solution in *D*, (iii) infinitely many solutions in *D*. In each example we take *N* = 2, *α* = 1.

Let *a*_{11} = *a*_{21} = ½, *a*_{12} = ⅓, *a*_{22} = ⅔, *b*_{1} = 2, *b*_{2} = 3. The corresponding system (1) has a unique solution in *D*, namely (2, 3).

Let again *a*_{11} = *a*_{21} = ½, *a*_{12} = ⅓, *a*_{22} = ⅔, but now *b*_{1} = 2, *b*_{2} = 1. Then the system has no solution in *D*.

Let *a*_{11} = *a*_{21} = *a*_{12} = *a*_{22} = *b*_{1} = *b*_{2} = 1. Then (*t*, 2 − *t*) is a solution of (1) for all 0 < *t* < 2.

We begin with another proof of the existence of solution in *D* to the system (1) if 0 < *α* < 1.

Fix an *ϵ* > 0 sufficiently small, such that *N**ϵ* ≤ *b*_{1} + ⋯ + *b*_{N} and

$$\frac{{\u03f5}^{\alpha -1}}{{({b}_{1}+\cdots +{b}_{N}-(N-1)\u03f5)}^{\alpha}}\ge max\left\{{\left(\sum _{n=1}^{N}{b}_{n}{a}_{ni}{\left(\sum _{j=1}^{N}{a}_{nj}\right)}^{-1}\right)}^{-1}\right|i=1,\dots ,N\}.$$

Let $W:=\{t\in {\mathbb{R}}^{N}\phantom{\rule{0.25em}{0ex}}|\phantom{\rule{0.25em}{0ex}}{\sum}_{j=1}^{N}{t}_{j}={\sum}_{j=1}^{N}{b}_{j},\phantom{\rule{0.25em}{0ex}}{t}_{j}\ge \u03f5,\phantom{\rule{0.25em}{0ex}}j=1,\dots ,N\}$. For *t* ∈ *W* and *i* ∈ {1, …, *N*} let ${H}_{i}(t):={\sum}_{n=1}^{N}{b}_{n}{a}_{ni}{t}_{i}^{\alpha}{({\sum}_{j=1}^{N}{a}_{nj}{t}_{j}^{\alpha})}^{-1}$.

It is not difficult to verify that *H*_{i}(*t*) ∈ *W*, so that we can consider the continuous function *H* = (*H*_{1}, …, *H*_{N}):*W* → *W*. Since *W* is compact and convex, from Brouwer’s theorem we deduce the existence of a *t* ∈ *W* such that *H*(*t*) = *t*. Obviously *t* ∈ *D* and *t* is a solution to (1).

The next theorem offers a third proof of the existence and also a proof of the uniqueness of the solution to (1) in the case 0 < *α* < 1.

*If*
0 < *α* < 1, *the system* (1) *has a unique solution in*
*D*.

Let *K*: = {*t* ∈ ℝ^{N}|*t* ≠ 0, *t*_{1} ≥ 0, …, *t*_{N} ≥ 0} and *f*:*K* → ℝ, $f(t):=\alpha {\sum}_{n=1}^{N}{t}_{n}-{\sum}_{n=1}^{N}{b}_{n}log{\sum}_{j=1}^{N}{a}_{nj}{t}_{j}^{\alpha}$.

The proof of the theorem is divided into the following steps.

*(1)*
*f*
*is strictly convex.*

Indeed, let *s*, *t* ∈ *K*, *s* ≠ *t*. Since the function *u* → *u*^{α} (*u* ≥ 0) is strictly concave, we have

$${\left(\frac{{s}_{j}+{t}_{j}}{2}\right)}^{\alpha}\ge \frac{{s}_{j}^{\alpha}+{t}_{j}^{\alpha}}{2},\phantom{\rule{1em}{0ex}}j=1,\dots ,N,$$

with at least one strict inequality.

Thus

$$\sum _{j=1}^{N}{a}_{nj}{\left(\frac{{s}_{j}+{t}_{j}}{2}\right)}^{\alpha}>\frac{1}{2}(\sum _{j=1}^{N}{a}_{nj}{s}_{j}^{\alpha}+\sum _{j=1}^{N}{a}_{nj}{t}_{j}^{\alpha}),\phantom{\rule{1em}{0ex}}n=1,\dots ,N.$$

The function log is strictly increasing and strictly concave, so that

$$\begin{array}{rl}\sum _{n=1}^{N}{b}_{n}log\sum _{j=1}^{N}{a}_{nj}{\left(\frac{{s}_{j}+{t}_{j}}{2}\right)}^{\alpha}& >\sum _{n=1}^{N}{b}_{n}log\frac{1}{2}(\sum _{j=1}^{N}{a}_{nj}{s}_{j}^{\alpha}+\sum _{j=1}^{N}{a}_{nj}{t}_{j}^{\alpha})\\ & >\frac{1}{2}(\sum _{n=1}^{N}{b}_{n}log\sum _{j=1}^{N}{a}_{nj}{s}_{j}^{\alpha}+\sum _{n=1}^{N}{b}_{n}log\sum _{j=1}^{N}{a}_{nj}{t}_{j}^{\alpha}).\end{array}$$

It follows immediately that

$$f\left(\frac{s+t}{2}\right)<\frac{1}{2}(f(s)+f(t)),$$

which means that *f* is strictly convex.

*(2)*
*f*:*K* → ℝ
*has a global minimum point*
*x, and*
*x* ∈ *D*
*.*

Let $\parallel t\parallel ={({t}_{1}^{2}+\cdots +{t}_{N}^{2})}^{1/2}$. Obviously

$$\underset{\parallel t\parallel \to 0}{lim}f(t)=+\mathrm{\infty}.$$

7

Writing

$$f(t)=\sum _{j=1}^{N}(\alpha {t}_{j}-\left(\sum _{n=1}^{N}{b}_{n}{a}_{nj}\right){t}_{j}^{\alpha})+\sum _{n=1}^{N}{b}_{n}(\sum _{j=1}^{N}{a}_{nj}{t}_{j}^{\alpha}-log\sum _{j=1}^{N}{a}_{nj}{t}_{j}^{\alpha}),$$

we see that

$$\underset{\parallel t\parallel \to +\mathrm{\infty}}{lim}f(t)=+\mathrm{\infty}.$$

8

For 0 < *r* < *R*, let *K*_{r,R}: = {*t* ∈ *K* | *r* ≤ ∥ *t* ∥ ≤ *R*}.

Let *s* = (1, 1, …, 1) ∈ *K*. Due to (7) and (8), there exist *r* and *R* such that, for all *t* ∈ *K*∖*K*_{r,R},

9

Since *K*_{r,R} is compact, the continuous function *f* restricted to *K*_{r,R} has a global minimum point in *K*_{r,R}; let us denote it by *x* ∈ *K*_{r,R}. Due to (9), *s* ∈ *K*_{r,R}, and so *f*(*x*) ≤ *f*(*s*) < *f*(*t*) for all *t* ∈ *K*∖*K*_{r,R}. It follows that *x* is a global minimum point of *f*:*K* → ℝ.

Suppose that *x* ∉ *D*. Then, for example, *x* = (0, *x*_{2}, …, *x*_{N}) and *x* ≠ 0. Consider the function *g*(*u*): = *f*(*u*, *x*_{2}, …, *x*_{N}), *u* ≥ 0. We get *g*^{′}(0) = −∞, and so *x* cannot be a minimum point of *f*. This concludes the proof of (2).

*(3) The critical points*
*t* ∈ *D*
*of*
*f*
*are the solutions in*
*D*
*of the system* (1).

Indeed, for *t* ∈ *D*, $\frac{\partial f}{\partial {t}_{i}}(t)=0$, *i* = 1, …, *N*, is equivalent to (1).

*(4) The system* (1) *has a unique solution in*
*D, namely the global minimum point*
*x.*

As a minimum point of *f* in the open set *D*, *x* is a critical point, hence a solution of (1). Let *y* ∈ *D* be another solution, *y* ≠ *x*. Then *y* is a critical point of *f*, and *f* is convex; it follows that *y* is a global minimum point of *f* on *D*. (See, *e.g.*, [15], p. 14, Theorem 1.17.) Thus *f*(*y*) = *f*(*x*); since *f* is *strictly* convex, we have

$$f\left(\frac{x+y}{2}\right)<\frac{f(x)+f(y)}{2}=f(x).$$

This contradicts the global minimality property of *x* and the proof of Theorem 2 is finished.□

The authors are grateful to Prof. Andrés Rodriguez-Clareb from University of California - Berkeley for suggesting the problem discussed in this article.

**Competing interests**

The authors declare that they have no competing interests.

**Authors’ contributions**

All authors contributed equally to this work. The authors read and approved the final manuscript.

**Publisher’s Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Anca Ciurte, Email: or.julctu.sc@etruiC.acnA.

Sergiu Nedevschi, Email: or.julctu.sc@ihcsvedeN.uigreS.

Ioan Rasa, Email: or.julctu.htam@asaR.naoI.

1. Ciurte A, Nedevschi S, Rasa I. Systems of nonlinear algebraic equations with unique solution. Numer. Algorithms. 2015;68:367–376. doi: 10.1007/s11075-014-9849-5. [Cross Ref]

2. Zhang G, Feng W. On the number of positive solutions of a nonlinear algebraic system. Linear Algebra Appl. 2007;422:404–421. doi: 10.1016/j.laa.2006.10.026. [Cross Ref]

3. Győri I, Hartung F, Mohamady NA. Existence and uniqueness of positive solutions of a system of nonlinear algebraic equations. Period. Math. Hung. 2016

4. Ciurte A, Nedevschi S, Rasa I. An algorithm for solving some nonlinear systems with applications to extremum problems. Taiwan. J. Math. 2012;16:1137–1150. doi: 10.11650/twjm/1500406683. [Cross Ref]

5. Bonanno G, Candito P, D’Agui G. Positive solutions for a nonlinear parameter-depending algebraic system. Electron. J. Differ. Equ. 2015;2015 doi: 10.1186/s13662-014-0347-9. [Cross Ref]

6. Farina L, Rinaldi S. Positive Linear Systems: Theory and Application. New York: Wiley-lnterscience; 2000.

7. Kaykobad M. Positive solutions of positive linear systems. Linear Algebra Appl. 1985;64:133–140. doi: 10.1016/0024-3795(85)90271-X. [Cross Ref]

8. Nelson P. Positive solutions of positive linear equations. Proc. Am. Math. Soc. 1972;31:453–457. doi: 10.1090/S0002-9939-1972-0288605-5. [Cross Ref]

9. Yang Y, Zhang JH. Existence results for a nonlinear system with a parameter. J. Math. Anal. Appl. 2008;340:658–668. doi: 10.1016/j.jmaa.2007.06.012. [Cross Ref]

10. Yang Y, Zhang J. Existence and multiple solutions for a nonlinear system with a parameter. Nonlinear Anal. 2009;70:2542–2548. doi: 10.1016/j.na.2008.03.040. [Cross Ref]

11. Zhang G. Existence of non-zero solutions for a nonlinear system with a parameter. Nonlinear Anal. 2007;66:1410–1416. doi: 10.1016/j.na.2006.01.024. [Cross Ref]

12. Zhang G, Bai L. Existence of solutions for a nonlinear algebraic system. Discrete Dyn. Nat. Soc. 2009;2009

13. Zhang G, Cheng SS. Existence of solutions for a nonlinear system with a parameter. J. Math. Anal. Appl. 2006;314:311–319. doi: 10.1016/j.jmaa.2005.03.098. [Cross Ref]

14. Costinota A, Rodriguez-Clareb A. Handbook of International Economics. Amterdam: North-Holland; 2014. Trade theory with numbers: quantifying the consequences of globalization; pp. 197–261.

15. Horst R, Pardalos PM, Thoai NV. Introduction to Global Optimization. Dordrecht: Kluwer Academic; 1995.

Articles from Springer Open Choice are provided here courtesy of **Springer**