On the Motion of Substance in a Channel of a Network: Extended Model and New Classes of Probability Distributions

Abstract

We discuss the motion of substance in a channel containing nodes of a network. Each node of the channel can exchange substance with: (i) neighboring nodes of the channel, (ii) network nodes which do not belong to the channel, and (iii) environment of the network. The new point in this study is that we assume possibility for exchange of substance among flows of substance between nodes of the channel and: (i) nodes that belong to the network but do not belong to the channel and (ii) environment of the network. This leads to an extension of the model of motion of substance and the extended model contains previous models as particular cases. We use a discrete-time model of motion of substance and consider a stationary regime of motion of substance in a channel containing a finite number of nodes. As results of the study, we obtain a class of probability distributions connected to the amount of substance in nodes of the channel. We prove that the obtained class of distributions contains all truncated discrete probability distributions of discrete random variable $\omega$ which can take values $0,1,\cdots ,N$. Theory for the case of a channel containing infinite number of nodes is presented in Appendix A. The continuous version of the discussed discrete probability distributions is described in Appendix B. The discussed extended model and obtained results can be used for the study of phenomena that can be modeled by flows in networks: motion of resources, traffic flows, motion of migrants, etc.

1. Introduction

Networks of different degrees of complexity arise often in research on complex systems. Flows of substances occur frequently in these systems and because of this, network flow models are used, for example, to study flows in computer networks [1], flows in electrical and communication networks [2], to detect network structure [3], to study flows in financial networks [4], flows connected to transportation problems [5,6,7], etc. [8,9,10,11,12]. Below, we discuss flow of substance in a channel of a network containing a chain of network’s nodes which are connected by edges. The kind of substance can be different (e.g., some resource which flows through the channel) and we consider a discrete-time model of motion of the substance through the studied channel. There are exchanges between nodes of the channel and (i) other nodes of the channel; (ii) network nodes which are not part of the channel; and (iii) environment of the network. The new point in this study is that in addition to the above exchanges we assume also possibility for exchange of substance among the flows between the nodes of the channel and (i) the network; (ii) the environment of the network. Thus, the discussed model extends models studied up to now, for example, in [13,14], and has different possible applications such as: (i) to model flow of a substance through a channel and use of part of this substance in some industrial process happening in the nodes of the channel or (ii) to model human migration flows. The last application is important because of the frequent use of probability and deterministic models of human migration [15,16,17,18,19] and as the migration flows are important for taking decisions about economic development of regions of a country [20,21,22,23,24] for analysis of migration networks [13,25,26,27,28,29,30,31], and ideological struggles [32,33]; for study of waves and probability distributions in population systems [34,35,36,37], etc. Migration-like models are also used in other research areas [38,39]. The model described below can be connected to an appropriate urn model. This is a useful connection as urn models are greatly applied in the research of various problems, e.g., genetics, clinical trials, biology, social dynamics, military theory, etc. [40,41,42,43].

The text below is organized as follows. In Section 2, we formulate a general model of motion of a substance in a channel of a network containing a finite number of nodes. In Section 3, we describe the particular case of the general model which will be discussed in the article. in Section 4 we obtain a class of probability distributions connected to the stationary motion of substance through the channel and show that this class of distributions contains all possible truncated discrete distributions of a random variable which can take values $0,\cdots ,N$. As a result of this, the obtained class of distributions contains, in particular, cases of the classes of (long-tail) distributions obtained in our previous research. Short discussion of the obtained results is presented in Section 5. Theory for the case of the channel containing infinitely many nodes is presented in Appendix A. The class of probability distributions obtained there contains, in particular, cases of the classes of distributions of Katz, Ord, Kemp, etc. In Appendix B, we obtain a probability distribution for the case of a continuous random variable taking values between 0 and 1. This distribution is a continuous analogue of the discrete distributions discussed in main text and in Appendix A.

2. Mathematical Formulation of the Model

We consider a channel in a network consisting of nodes connected by edges—Figure 1. We assume that the channel consists of a chain of $N+1$ nodes (labeled from 0 to N) connected by corresponding edges (ways for motion of the substance). Each edge connects two nodes and each node is connected to two edges except for 0-th node and N-th node which are connected to one edge. Some kind of substance moves through the nodes and edges of the channel.

This motion is accompanied by processes of the exchange of substance between the channel, the rest of the network, and the environment of the network. The possible processes of exchange of substance for the i-th node of the channel are shown in Figure 2. One can image that an urn is placed in any of the nodes of the channel and in any of other nodes of the network and exchange of substance is among these urns. The i-th node can exchange substance with the $(i-1)$-th and with $(i+1)$-th nodes. The i-th node can also exchange substance with network nodes outside the channel and with the environment of the network. Let us denote as: (i) “leakage”—the process of motion of substance from a node of the channel to a node of the network or from a node of the channel to environment of the network; (ii) “inflow”—the process of motion of substance from a node of the network that does not belong to the channel or from the environment of the network to a node of the channel. The additional exchanges of substance connected to the flows between the i-th node and other nodes are shown in Figure 2 by arrows having dashed or dot-dashed lines. These exchanges are not accounted for by the previous models.

We consider a discrete-time and assume that intervals between moments of time have equal length. At each time interval, the substance in a node of the channel can participate in the following processes:

(1)

substance remains in the same node;

(2)

substance moves to previous or to the next node (e.g., substance may move from node m to node $m+1$ or from node m to node $m-1$);

(3)

substance leaks from the node m. Leaked substance does not belong anymore to the channel. Such substance may spread through the network or through the environment of the network. Thus, we have two kinds of leakage: (i) leakage from node m of the channel to nodes of the network and (ii) leakage from node m of the channel to the environment of the network;

(4)

the substance “inflows” to node m. Two kinds of “inflow” are possible: (i) inflow from nodes of the network to the node m of the channel and (ii) inflow from environment of the network to the node m of the channel.

In addition, the substance from any of the flows between nodes of the channel can participate in the following processes:

(5)

substance can move from the flow to a node of the network which does not belong to the channel (arrows with dashed line pointing from the flow in Figure 2);

(6)

substance can move from a node of the network which does not belong to the channel to the flow between two nodes of the channel (arrows with dashed line pointing to the flow in Figure 2);

(7)

substance can move from the flow to the environment of the network (arrows with dot-dashed line pointing from the flow in Figure 2);

(8)

substance can move from the environment of the network to the flow between two nodes of the channel ( arrows with dot-dashed line pointing to the flow in Figure 2).

Let us formalize mathematically the above considerations. We consider discrete-time $t_k$, $k=0,1,2,\cdots$ and denote by $x_i\left(t_k\right)$ the amount of substance in the i-th node of the channel at the beginning of time interval $[t_k,t_k+\mathsf{\Delta }t]$. For processes happening in this time interval in the n-th node of the channel, we use the following notations:

1.

$i_n^e\left(t_k\right)$ and $o_n^e\left(t_k\right)$ are amounts of inflow and outflow of substance from the environment to the n-th node of the channel (upper index e means that the quantities are for the environment);

2.

$o_n^c\left(t_k\right)$ is the amount of outflow of substance from the n-th node of the channel to the $(n+1)$-th node of the channel (upper index c means that the quantity is for the channel);

3.

$i_n^c\left(t_k\right)$ is the amount of inflow of substance from the $(n+1)$ node of the channel to the n-th node of the channel;

4.

$o_n^n\left(t_k\right)$ and $i_n^n\left(t_k\right)$ are the amounts of outflow and inflow of substance between the n-th node of the channel and the network (upper index n means that quantities are for the network).

In addition, there are the following exchanges for the flows between nodes of the channel, nodes of the network, and environment of the network. We note that in general we have two flows of substance between the nodes of the channel: (a) flow directed from i-th node of the channel to $i+1$-th node of the channel, and (b) flow directed from $i+1$-th node of the channel to i-th node of the channel. We denote the corresponding flows as follows:

5.

$f_n^e\left(t_k\right)$: amount of the outflow of substance to the environment of network from the flow between the n-th and $n+1$-th nodes of the channel;

6.

$f_n^n\left(t_k\right)$: amount of the outflow of substance to the network from the flow between the n-th and $n+1$-th nodes of the channel;

7.

$g_n^e\left(t_k\right)$: amount of the inflow of substance from the environment of the network to the flow between the n-th and $n+1$-th nodes of the channel;

8.

$g_n^n\left(t_k\right)$: amount of the inflow of substance from the network to the flow between the n-th and $n+1$-th nodes of the channel;

9.

$h_n^e\left(t_k\right)$: amount of the outflow of substance to the environment of network from the flow between the $n+1$-th and n-th nodes of the channel;

10.

$h_n^n\left(t_k\right)$: amount of the outflow of substance to the network from the flow between the $n+1$-th and n-th node of the channel;

11.

$p_n^e\left(t_k\right)$: amount of the inflow of substance from the environment of the network to the flow between the $n+1$-th and n-th nodes of the channel;

12.

$p_n^n\left(t_k\right)$: amount of the inflow of substance from the network to the flow between the $n+1$-th and n-th nodes of the channel.

We are interested in the amounts of substance available in the nodes of the channel. For the 0-th node there are: exchange of substance with the environment (inflow and outflow); exchange of substance with next node of channel (inflow and outflow); and exchange of substance with the network (inflow and outflow). In addition to the inflow from the neighbor node of the channel, there are exchanges of substance (inflow and outflow) with the network and with the environment of the network. Similarly, 4 exchanges are available also for the outflow from the 0-th node of the channel to the next node of the channel. From all of these exchanges, what are significant for the 0-th node are the inflows and outflows from the network and from the environment of the network towards the inflow of substance in the 0-th node as these exchanges of substance contribute to the change of substance in the 0-th node of the channel. Thus the change of amount of substance in 0-th node of channel is described by the relationship

$$\begin{array}{ccc}\hfill x_0\left(t_{k+1}\right)& =& x_0\left(t_k\right)+i_0^e\left(t_k\right)-o_0^e\left(t_k\right)-o_0^c\left(t_k\right)+i_0^c\left(t_k\right)-o_0^n\left(t_k\right)+i_0^n\left(t_k\right)-h_0^e\left(t_k\right)-\hfill \\ & & h_0^n\left(t_k\right)+p_0^e\left(t_k\right)+p_0^n\left(t_k\right).\hfill \end{array}$$

(1)

For the nodes of the channel numbered by $i=1,\cdots ,N-1$, there is an exchange of substance with the environment, exchange of substance with the network, and exchange of substance with the $(i-1)$-th and $(i+1)$-th node of the channel. Thus, the change of amount of substance in the i-th node of channel is described by relationship

$$\begin{array}{ccc}\hfill x_i\left(t_{k+1}\right)& =& x_i\left(t_k\right)+i_i^e\left(t_k\right)-o_i^e\left(t_k\right)+o_{i-1}^c\left(t_k\right)-i_{i-1}^c\left(t_k\right)-o_i^c\left(t_k\right)+i_i^c\left(t_k\right)-o_i^n\left(t_k\right)+i_i^n\left(t_k\right)-\hfill \\ & & h_i^e\left(t_k\right)-h_i^n\left(t_k\right)+p_i^e\left(t_k\right)+p_i^n\left(t_k\right)-f_{i-1}^e\left(t_k\right)-f_{i-1}^n\left(t_k\right)+g_{i-1}^e\left(t_k\right)+g_{i-1}^n\left(t_k\right),\hfill \\ & & i=1,\cdots ,N-1\hfill \end{array}$$

(2)

For the N-th node of the channel, there is an exchange of substance with the environment, exchange of substance with the network, and exchange of substance with the $(N-1)$-th node of the channel. Thus, the change of the amount of substance in the N-th node of channel is described by relationship

$$\begin{array}{ccc}\hfill x_N\left(t_{k+1}\right)& =& x_N\left(t_k\right)+i_N^e\left(t_k\right)-o_N^e\left(t_k\right)+o_{N-1}^c\left(t_k\right)-i_{N-1}^c\left(t_k\right)-o_N^n\left(t_k\right)+i_N^n\left(t_k\right)-f_{N-1}^e\left(t_k\right)-\hfill \\ & & f_{N-1}^n\left(t_k\right)+g_{N-1}^e\left(t_k\right)+g_{N-1}^n\left(t_k\right).\hfill \end{array}$$

(3)

Equations (1)–(3) describe the general case of motion of substance along a channel of nodes connected to a network and in direction to the environment of this network.

3. Studied Particular Case of the General Model (1)–(3)

We continue our study by consideration of a particular case for the relationships for the quantities from the system of Equations (1)–(3). In this particular case, we have linear relationships connecting the exchanges between the nodes of the channel and the amounts of substance in the nodes. The relationships are as follows.

• Exchange between the nodes of the channel and the environment of the network

  $$\begin{array}{ccc}\hfill i_0^e\left(t_k\right)& =& {\sigma }_0\left(t_k\right)x_0\left(t_k\right);o_0^e\left(t_k\right)={\mu }_0\left(t_k\right)x_0\left(t_k\right)\hfill \\ \hfill i_i^e\left(t_k\right)& =& {\sigma }_i\left(t_k\right)x_i\left(t_k\right);o_i^e\left(t_k\right)={\mu }_i\left(t_k\right)x_i\left(t_k\right),i=1,\cdots ,N-1\hfill \\ \hfill i_N^e\left(t_k\right)& =& {\sigma }_N\left(t_k\right)x_N\left(t_k\right);o_N^e\left(t_k\right)={\mu }_N\left(t_k\right)x_N\left(t_k\right)\hfill \end{array}$$

  (4)
• Exchange between the channel and the network

  $$\begin{array}{ccc}\hfill i_0^n\left(t_k\right)& =& {ϵ}_0\left(t_k\right)x_0\left(t_k\right);o_0^n\left(t_k\right)={\gamma }_0\left(t_k\right)x_0\left(t_k\right)\hfill \\ \hfill i_i^n\left(t_k\right)& =& {ϵ}_i\left(t_k\right)x_i\left(t_k\right);o_i^n\left(t_k\right)={\gamma }_i\left(t_k\right)x_i\left(t_k\right),i=1,\cdots ,N-1\hfill \\ \hfill i_N^n\left(t_k\right)& =& {ϵ}_N\left(t_k\right)x_N\left(t_k\right);o_N^n\left(t_k\right)={\gamma }_N\left(t_k\right)x_N\left(t_k\right)\hfill \end{array}$$

  (5)
• Exchange within the channel

  $$\begin{array}{ccc}\hfill o_0^c\left(t_k\right)& =& {\phi }_0\left(t_k\right)x_0\left(t_k\right);i_0^c\left(t_k\right)={\delta }_1\left(t_k\right)x_1\left(t_k\right)\hfill \\ \hfill o_i^c\left(t_k\right)& =& {\phi }_i\left(t_k\right)x_i\left(t_k\right);i_i^c\left(t_k\right)={\delta }_{i+1}\left(t_k\right)x_{i+1}\left(t_k\right),i=1,\cdots ,N-2\hfill \\ \hfill o_{N-1}^c\left(t_k\right)& =& {\phi }_{N-1}\left(t_k\right)x_{N-1}\left(t_k\right);\hfill \\ \hfill i_{N-1}^c\left(t_k\right)& =& {\delta }_N\left(t_k\right)x_N\left(t_k\right)\hfill \end{array}$$

  (6)
• Exchanges between flows among the nodes of the channel and the network/environment of the network

  $$\begin{array}{ccc}\hfill f_n^e\left(t_k\right)& =& {\zeta }_n\left(t_k\right);f_n^n\left(t_k\right)={\theta }_n\left(t_k\right);\hfill \\ \hfill g_n^e\left(t_k\right)& =& {\kappa }_n\left(t_k\right);g_n^n\left(t_k\right)={\lambda }_n\left(t_k\right)\hfill \\ \hfill h_n^e\left(t_k\right)& =& {\nu }_n\left(t_k\right);h_n^n\left(t_k\right)={\pi }_n\left(t_k\right)\hfill \\ \hfill p_n^e\left(t_k\right)& =& {\rho }_n\left(t_k\right);p_n^n\left(t_k\right)={\tau }_n\left(t_k\right)\hfill \end{array}$$

  (7)

For this particular case, the system of Equations (1)–(3) becomes

$$\begin{array}{ccc}\hfill x_0\left(t_{k+1}\right)& =& x_0\left(t_k\right)+{\sigma }_0\left(t_k\right)x_0\left(t_k\right)-{\mu }_0\left(t_k\right)x_0\left(t_k\right)-{\phi }_0\left(t_k\right)x_0\left(t_k\right)+{\delta }_1\left(t_k\right)x_1\left(t_k\right)-\hfill \\ & & {\gamma }_0\left(t_k\right)x_0\left(t_k\right)+{ϵ}_0\left(t_k\right)x_0\left(t_k\right)-{\nu }_0\left(t_k\right)-{\pi }_0\left(t_k\right)+{\rho }_0\left(t_k\right)+{\tau }_0\left(t_k\right),\hfill \end{array}$$

(8)

$$\begin{array}{ccc}\hfill x_i\left(t_{k+1}\right)& =& x_i\left(t_k\right)+{\sigma }_i\left(t_k\right)x_i\left(t_k\right)-{\mu }_i\left(t_k\right)x_i\left(t_k\right)+{\phi }_{i-1}\left(t_k\right)x_{i-1}\left(t_k\right)-{\delta }_i\left(t_k\right)x_i\left(t_k\right)-\hfill \\ & & {\phi }_i\left(t_k\right)x_i\left(t_k\right)+{\delta }_{i+1}\left(t_k\right)x_{i+1}\left(t_k\right)-{\gamma }_i\left(t_k\right)x_i\left(t_k\right)+{ϵ}_i\left(t_k\right)x_i\left(t_k\right)-{\nu }_i\left(t_k\right)-\hfill \\ & & {\pi }_i\left(t_k\right)+{\rho }_i\left(t_k\right)+{\tau }_i\left(t_k\right)-{\zeta }_{i-1}\left(t_k\right)-{\theta }_{i-1}\left(t_k\right)+{\kappa }_{i-1}\left(t_k\right)+{\lambda }_{i-1}\left(t_k\right),\hfill \\ & & i=1,\cdots ,N-1,\hfill \end{array}$$

(9)

$$\begin{array}{ccc}\hfill x_N\left(t_{k+1}\right)& =& x_N\left(t_k\right)+{\sigma }_N\left(t_k\right)x_N\left(t_k\right)-{\mu }_N\left(t_k\right)x_i\left(t_k\right)+{\phi }_{N-1}\left(t_k\right)x_{N-1}\left(t_k\right)-{\delta }_N\left(t_k\right)x_N\left(t_k\right)-\hfill \\ & & {\gamma }_N\left(t_k\right)x_N\left(t_k\right)+{ϵ}_N\left(t_k\right)x_N\left(t_k\right)-{\zeta }_{N-1}\left(t_k\right)-{\theta }_{N-1}\left(t_k\right)+{\kappa }_{N-1}\left(t_k\right)+{\lambda }_{N-1}\left(t_k\right).\hfill \end{array}$$

Below, we consider the case of absence of an inflow from the $(i+1)$-th node to the i-th node of the channel (no flow of substance in the direction opposite to the direction from the 0-th node to the N-th node of the channel). In this case, the system of Equations (8)–(10) becomes

$$\begin{array}{ccc}\hfill x_0\left(t_{k+1}\right)& =& x_0\left(t_k\right)+{\sigma }_0\left(t_k\right)x_0\left(t_k\right)-{\mu }_0\left(t_k\right)x_0\left(t_k\right)-{\phi }_0\left(t_k\right)x_0\left(t_k\right)-{\gamma }_0\left(t_k\right)x_0\left(t_k\right)\hfill \\ & & +{ϵ}_0\left(t_k\right)x_0\left(t_k\right)\hfill \end{array}$$

(10)

$$\begin{array}{ccc}\hfill x_i\left(t_{k+1}\right)& =& x_i\left(t_k\right)+{\sigma }_i\left(t_k\right)x_i\left(t_k\right)-{\mu }_i\left(t_k\right)x_i\left(t_k\right)+{\phi }_{i-1}\left(t_k\right)x_{i-1}\left(t_k\right)-{\phi }_i\left(t_k\right)x_i\left(t_k\right)-\hfill \\ & & {\gamma }_i\left(t_k\right)x_i\left(t_k\right)+{ϵ}_i\left(t_k\right)x_i\left(t_k\right)-{\zeta }_{i-1}\left(t_k\right)-{\theta }_{i-1}\left(t_k\right)+{\kappa }_{i-1}\left(t_k\right)+{\lambda }_{i-1}\left(t_k\right)\hfill \\ & & i=1,\cdots ,N-1\hfill \end{array}$$

(11)

$$\begin{array}{ccc}\hfill x_N\left(t_{k+1}\right)& =& x_N\left(t_k\right)+{\sigma }_N\left(t_k\right)x_N\left(t_k\right)-{\mu }_N\left(t_k\right)x_i\left(t_k\right)+{\phi }_{N-1}\left(t_k\right)x_{N-1}\left(t_k\right)-{\gamma }_N\left(t_k\right)x_N\left(t_k\right)+\hfill \\ & & {ϵ}_N\left(t_k\right)x_N\left(t_k\right)-{\zeta }_{N-1}\left(t_k\right)-{\theta }_{N-1}\left(t_k\right)+{\kappa }_{N-1}\left(t_k\right)+{\lambda }_{N-1}\left(t_k\right)\hfill \end{array}$$

(12)

We shall study the stationary case of the model Equations (10)–(12) in more detail below.

4. Results

We consider the case of stationary motion of the substance through the channel of nodes. In this case, $x_i\left(t_{k+1}\right)=x_i\left(t_k\right)$, $i=0,\cdots ,N$ and the system of Equations (10)–(12) becomes

$$\begin{array}{c}\hfill \left[{\sigma }_0\left(t_k\right)-{\mu }_0\left(t_k\right)-{\phi }_0\left(t_k\right)-{\gamma }_0\left(t_k\right)+{ϵ}_0\left(t_k\right)\right]x_0\left(t_k\right)=0\end{array}$$

(13)

$$\begin{array}{c}\hfill \left[{\mu }_i\left(t_k\right)+{\phi }_i\left(t_k\right)+{\gamma }_i\left(t_k\right)-{\sigma }_i\left(t_k\right)-{ϵ}_i\left(t_k\right)\right]x_i\left(t_k\right)={\phi }_{i-1}\left(t_k\right)x_{i-1}\left(t_k\right)+(-{\zeta }_{i-1}\left(t_k\right)-{\theta }_{i-1}\left(t_k\right)+\\ \hfill {\kappa }_{i-1}\left(t_k\right)+{\lambda }_{i-1}\left(t_k\right)),i=1,\cdots ,N-1\end{array}$$

(14)

$$\begin{array}{c}\hfill \left[{\mu }_N\left(t_k\right)+{\gamma }_N\left(t_k\right)-{\sigma }_N\left(t_k\right)-{ϵ}_N\left(t_k\right)\right]x_N\left(t_k\right)={\phi }_{N-1}\left(t_k\right)x_{N-1}\left(t_k\right)+(-{\zeta }_{N-1}\left(t_k\right)-{\theta }_{N-1}\left(t_k\right)+\\ \hfill {\kappa }_{N-1}\left(t_k\right)+{\lambda }_{N-1}\left(t_k\right))\end{array}$$

(15)

Below we discuss the model described by Equations (13)–(15) for the case when the parameters of the model are time independent (i.e., when ${\sigma }_i\left(t_k\right)=\sigma$; ${\mu }_i\left(t_k\right)={\mu }_i$; ${\gamma }_i\left(t_k\right)={\gamma }_i$; ${ϵ}_i\left(t_k\right)={ϵ}_i$; $f_0\left(t_k\right)=f_0$; $i=0,\cdots ,N$. Note that these parameters do not depend on time but they may depend on i and also on other parameters connected to the network and to the environment of the network. In this case, the system of equations becomes

$$\begin{array}{c}\hfill \left[{\sigma }_0-{\mu }_0-{\phi }_0-{\gamma }_0+{ϵ}_0\right]x_0=0\end{array}$$

(16)

$$\begin{array}{c}\hfill \left[{\mu }_i+{\phi }_i+{\gamma }_i-{\sigma }_i-{ϵ}_i\right]x_i={\phi }_{i-1}{x}_{i-1}+(-{\zeta }_{i-1}-{\theta }_{i-1}+{\kappa }_{i-1}+{\lambda }_{i-1}),i=1,\cdots ,N-1\end{array}$$

(17)

$$\begin{array}{c}\hfill \left[{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N\right]x_N={\phi }_{N-1}{x}_{N-1}+(-{\zeta }_{N-1}-{\theta }_{N-1}+{\kappa }_{N-1}+{\lambda }_{N-1})\end{array}$$

(18)

From the system of Equations (16)–(18), we obtain the following relationships (note that ${\delta }_{i,j}$ below is the Kronecker delta symbol):

$$\begin{array}{ccc}\hfill {\phi }_0& =& {\sigma }_0-{\mu }_0-{\gamma }_0+{ϵ}_0,\hfill \end{array}$$

$$\begin{array}{ccc}\hfill x_k& =& x_0\prod _{i=1}^k\frac{{\phi }_{i-1}}{{\mu }_i+{\phi }_i+{\gamma }_i-{\sigma }_i-{ϵ}_i}+\hfill \\ & & (1-{\delta }_{k1})\sum _{j=1}^{k-1}\left(\frac{{\kappa }_{j-1}+{\lambda }_{j-1}-{\zeta }_{j-1}-{\theta }_{j-1}}{{\mu }_j+{\phi }_j+{\gamma }_j-{\sigma }_j-{ϵ}_j}\prod _{l=j+1}^k\frac{{\phi }_{l-1}}{{\mu }_l+{\phi }_l+{\gamma }_l-{\sigma }_l-{ϵ}_l}\right)+\hfill \\ & & \frac{{\kappa }_{k-1}+{\lambda }_{k-1}-{\zeta }_{k-1}-{\theta }_{k-1}}{{\mu }_k+{\phi }_k+{\gamma }_k-{\sigma }_k-{ϵ}_k},k=1,\cdots ,N-1,\hfill \end{array}$$

$$\begin{array}{ccc}\hfill x_N& =& x_0\frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\prod _{i=1}^{N-1}\frac{{\phi }_{i-1}}{{\mu }_i+{\phi }_i+{\gamma }_i-{\sigma }_i-{ϵ}_i}+\hfill \\ & & \frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}(1-{\delta }_{N,2})\sum _{j=1}^{N-2}\left(\frac{{\kappa }_{j-1}+{\lambda }_{j-1}-{\zeta }_{j-1}-{\theta }_{j-1}}{{\mu }_j+{\phi }_j+{\gamma }_j-{\sigma }_j-{ϵ}_j}\prod _{l=j+1}^{N-1}\frac{{\phi }_{l-1}}{{\mu }_l+{\phi }_l+{\gamma }_l-{\sigma }_l-{ϵ}_l}\right)+\hfill \\ & & \frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\frac{{\kappa }_{N-2}+{\lambda }_{N-2}-{\zeta }_{N-2}-{\theta }_{N-2}}{{\mu }_{N-1}+{\phi }_{N-1}+{\gamma }_{N-1}-{\sigma }_{N-1}-{ϵ}_{N-1}}+\frac{{\kappa }_{N-1}+{\lambda }_{N-1}-{\zeta }_{N-1}-{\theta }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N},\hfill \end{array}$$

(19)

and

$$\begin{array}{ccc}\hfill x& =& \sum _{k=0}^N{x}_k=x_0\left[1+\sum _{k=1}^{N-1}\prod _{i=1}^k\frac{{\phi }_{i-1}}{{\mu }_i+{\phi }_i+{\gamma }_i-{\sigma }_i-{ϵ}_i}+\frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\prod _{i=1}^{N-1}\frac{{\phi }_{i-1}}{{\mu }_i+{\phi }_i+{\gamma }_i-{\sigma }_i-{ϵ}_i}\right]+\hfill \\ & & \sum _{k=1}^{N-1}[(1-{\delta }_{k,1})\sum _{j=1}^{k-1}\left(\frac{{\kappa }_{j-1}+{\lambda }_{j-1}-{\zeta }_{j-1}-{\theta }_{j-1}}{{\mu }_j+{\phi }_j+{\gamma }_j-{\sigma }_j-{ϵ}_j}\prod _{l=j+1}^k\frac{{\phi }_{l-1}}{{\mu }_l+{\phi }_l+{\gamma }_l-{\sigma }_l-{ϵ}_l}\right)+\hfill \\ & & \frac{{\kappa }_{k-1}+{\lambda }_{k-1}-{\zeta }_{k-1}-{\theta }_{k-1}}{{\mu }_k+{\phi }_k+{\gamma }_k-{\sigma }_k-{ϵ}_k}]+\frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}(1-{\delta }_{N,2})\sum _{j=1}^{N-2}(\frac{{\kappa }_{j-1}+{\lambda }_{j-1}-{\zeta }_{j-1}-{\theta }_{j-1}}{{\mu }_j+{\phi }_j+{\gamma }_j-{\sigma }_j-{ϵ}_j}\times \hfill \\ & & \prod _{l=j+1}^{N-1}\frac{{\phi }_{l-1}}{{\mu }_l+{\phi }_l+{\gamma }_l-{\sigma }_l-{ϵ}_l})+\frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\frac{{\kappa }_{N-2}+{\lambda }_{N-2}-{\zeta }_{N-2}-{\theta }_{N-2}}{{\mu }_{N-1}+{\phi }_{N-1}+{\gamma }_{N-1}-{\sigma }_{N-1}-{ϵ}_{N-1}}+\hfill \\ & & \frac{{\kappa }_{N-1}+{\lambda }_{N-1}-{\zeta }_{N-1}-{\theta }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}.\hfill \end{array}$$

(20)

Equations (19) and (20) lead to a class of probability distributions as follows. We have $x_i$ and x and we can consider probability distribution $y_i=x_i/x$ connected to amount of substance in nodes of the channel. $y_i$ can be considered as probability values of distribution of a discrete random variable $\omega$: $y_i=p(\omega =i)$, $i=0,\cdots ,N$. For this distribution, we obtain

$$\begin{array}{ccc}\hfill y_0& =& x_0/\{x_0\left[1+\sum _{k=1}^{N-1}\prod _{i=1}^k\frac{{\phi }_{i-1}}{{\mu }_i+{\phi }_i+{\gamma }_i-{\sigma }_i-{ϵ}_i}+\frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\prod _{i=1}^{N-1}\frac{{\phi }_{i-1}}{{\mu }_i+{\phi }_i+{\gamma }_i-{\sigma }_i-{ϵ}_i}\right]+\hfill \\ & & \sum _{k=1}^{N-1}[(1-{\delta }_{k,1})\sum _{j=1}^{k-1}\left(\frac{{\kappa }_{j-1}+{\lambda }_{j-1}-{\zeta }_{j-1}-{\theta }_{j-1}}{{\mu }_j+{\phi }_j+{\gamma }_j-{\sigma }_j-{ϵ}_j}\prod _{l=j+1}^k\frac{{\phi }_{l-1}}{{\mu }_l+{\phi }_l+{\gamma }_l-{\sigma }_l-{ϵ}_l}\right)+\hfill \\ & & \frac{{\kappa }_{k-1}+{\lambda }_{k-1}-{\zeta }_{k-1}-{\theta }_{k-1}}{{\mu }_k+{\phi }_k+{\gamma }_k-{\sigma }_k-{ϵ}_k}]+\frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}(1-{\delta }_{N,2})\sum _{j=1}^{N-2}(\frac{{\kappa }_{j-1}+{\lambda }_{j-1}-{\zeta }_{j-1}-{\theta }_{j-1}}{{\mu }_j+{\phi }_j+{\gamma }_j-{\sigma }_j-{ϵ}_j}\times \hfill \\ & & \prod _{l=j+1}^{N-1}\frac{{\phi }_{l-1}}{{\mu }_l+{\phi }_l+{\gamma }_l-{\sigma }_l-{ϵ}_l})+\frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\frac{{\kappa }_{N-2}+{\lambda }_{N-2}-{\zeta }_{N-2}-{\theta }_{N-2}}{{\mu }_{N-1}+{\phi }_{N-1}+{\gamma }_{N-1}-{\sigma }_{N-1}-{ϵ}_{N-1}}+\hfill \\ & & \frac{{\kappa }_{N-1}+{\lambda }_{N-1}-{\zeta }_{N-1}-{\theta }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\},\hfill \end{array}$$

(21)

$$\begin{array}{ccc}\hfill y_k& =& \{x_0\prod _{i=1}^k\frac{{\phi }_{i-1}}{{\mu }_i+{\phi }_i+{\gamma }_i-{\sigma }_i-{ϵ}_i}+(1-{\delta }_{k1})\sum _{j=1}^{k-1}\left(\frac{{\kappa }_{j-1}+{\lambda }_{j-1}-{\zeta }_{j-1}-{\theta }_{j-1}}{{\mu }_j+{\phi }_j+{\gamma }_j-{\sigma }_j-{ϵ}_j}\prod _{l=j+1}^k\frac{{\phi }_{l-1}}{{\mu }_l+{\phi }_l+{\gamma }_l-{\sigma }_l-{ϵ}_l}\right)+\hfill \\ & & \frac{{\kappa }_{k-1}+{\lambda }_{k-1}-{\zeta }_{k-1}-{\theta }_{k-1}}{{\mu }_k+{\phi }_k+{\gamma }_k-{\sigma }_k-{ϵ}_k}\}/\{x_0[1+\sum _{k=1}^{N-1}\prod _{i=1}^k\frac{{\phi }_{i-1}}{{\mu }_i+{\phi }_i+{\gamma }_i-{\sigma }_i-{ϵ}_i}+\hfill \\ & & \frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\prod _{i=1}^{N-1}\frac{{\phi }_{i-1}}{{\mu }_i+{\phi }_i+{\gamma }_i-{\sigma }_i-{ϵ}_i}]+\hfill \\ & & \sum _{k=1}^{N-1}\left[(1-{\delta }_{k,1})\sum _{j=1}^{k-1}\left(\frac{{\kappa }_{j-1}+{\lambda }_{j-1}-{\zeta }_{j-1}-{\theta }_{j-1}}{{\mu }_j+{\phi }_j+{\gamma }_j-{\sigma }_j-{ϵ}_j}\prod _{l=j+1}^k\frac{{\phi }_{l-1}}{{\mu }_l+{\phi }_l+{\gamma }_l-{\sigma }_l-{ϵ}_l}\right)+\frac{{\kappa }_{k-1}+{\lambda }_{k-1}-{\zeta }_{k-1}-{\theta }_{k-1}}{{\mu }_k+{\phi }_k+{\gamma }_k-{\sigma }_k-{ϵ}_k}\right]+\hfill \\ & & \frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}(1-{\delta }_{N,2})\sum _{j=1}^{N-2}\left(\frac{{\kappa }_{j-1}+{\lambda }_{j-1}-{\zeta }_{j-1}-{\theta }_{j-1}}{{\mu }_j+{\phi }_j+{\gamma }_j-{\sigma }_j-{ϵ}_j}\prod _{l=j+1}^{N-1}\frac{{\phi }_{l-1}}{{\mu }_l+{\phi }_l+{\gamma }_l-{\sigma }_l-{ϵ}_l}\right)+\hfill \\ & & \frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\frac{{\kappa }_{N-2}+{\lambda }_{N-2}-{\zeta }_{N-2}-{\theta }_{N-2}}{{\mu }_{N-1}+{\phi }_{N-1}+{\gamma }_{N-1}-{\sigma }_{N-1}-{ϵ}_{N-1}}+\frac{{\kappa }_{N-1}+{\lambda }_{N-1}-{\zeta }_{N-1}-{\theta }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\},\hfill \\ & & k=1,\cdots ,N-1,\hfill \end{array}$$

(22)

$$\begin{array}{ccc}\hfill y_N& =& \{x_0\frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\prod _{i=1}^{N-1}\frac{{\phi }_{i-1}}{{\mu }_i+{\phi }_i+{\gamma }_i-{\sigma }_i-{ϵ}_i}+\hfill \\ & & \frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}(1-{\delta }_{N,2})\sum _{j=1}^{N-2}\left(\frac{{\kappa }_{j-1}+{\lambda }_{j-1}-{\zeta }_{j-1}-{\theta }_{j-1}}{{\mu }_j+{\phi }_j+{\gamma }_j-{\sigma }_j-{ϵ}_j}\times \prod _{l=j+1}^{N-1}\frac{{\phi }_{l-1}}{{\mu }_l+{\phi }_l+{\gamma }_l-{\sigma }_l-{ϵ}_l}\right)+\hfill \\ & & \frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\frac{{\kappa }_{N-2}+{\lambda }_{N-2}-{\zeta }_{N-2}-{\theta }_{N-2}}{{\mu }_{N-1}+{\phi }_{N-1}+{\gamma }_{N-1}-{\sigma }_{N-1}-{ϵ}_{N-1}}+\frac{{\kappa }_{N-1}+{\lambda }_{N-1}-{\zeta }_{N-1}-{\theta }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\}/\hfill \\ & & \{x_0\left[1+\sum _{k=1}^{N-1}\prod _{i=1}^k\frac{{\phi }_{i-1}}{{\mu }_i+{\phi }_i+{\gamma }_i-{\sigma }_i-{ϵ}_i}+\frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\prod _{i=1}^{N-1}\frac{{\phi }_{i-1}}{{\mu }_i+{\phi }_i+{\gamma }_i-{\sigma }_i-{ϵ}_i}\right]+\hfill \\ & & \sum _{k=1}^{N-1}[(1-{\delta }_{k,1})\sum _{j=1}^{k-1}\left(\frac{{\kappa }_{j-1}+{\lambda }_{j-1}-{\zeta }_{j-1}-{\theta }_{j-1}}{{\mu }_j+{\phi }_j+{\gamma }_j-{\sigma }_j-{ϵ}_j}\prod _{l=j+1}^k\frac{{\phi }_{l-1}}{{\mu }_l+{\phi }_l+{\gamma }_l-{\sigma }_l-{ϵ}_l}\right)+\hfill \\ & & \frac{{\kappa }_{k-1}+{\lambda }_{k-1}-{\zeta }_{k-1}-{\theta }_{k-1}}{{\mu }_k+{\phi }_k+{\gamma }_k-{\sigma }_k-{ϵ}_k}]+\frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}(1-{\delta }_{N,2})\sum _{j=1}^{N-2}(\frac{{\kappa }_{j-1}+{\lambda }_{j-1}-{\zeta }_{j-1}-{\theta }_{j-1}}{{\mu }_j+{\phi }_j+{\gamma }_j-{\sigma }_j-{ϵ}_j}\times \hfill \\ & & \prod _{l=j+1}^{N-1}\frac{{\phi }_{l-1}}{{\mu }_l+{\phi }_l+{\gamma }_l-{\sigma }_l-{ϵ}_l})+\frac{{\phi }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\frac{{\kappa }_{N-2}+{\lambda }_{N-2}-{\zeta }_{N-2}-{\theta }_{N-2}}{{\mu }_{N-1}+{\phi }_{N-1}+{\gamma }_{N-1}-{\sigma }_{N-1}-{ϵ}_{N-1}}+\hfill \\ & & \frac{{\kappa }_{N-1}+{\lambda }_{N-1}-{\zeta }_{N-1}-{\theta }_{N-1}}{{\mu }_N+{\gamma }_N-{\sigma }_N-{ϵ}_N}\}.\hfill \end{array}$$

(23)

To the best of our knowledge, the class of distributions (21)–(23) was not discussed by other authors. The corresponding class of distributions for the case $N=\infty$ is discussed in Appendix A. We note that the system of Equations (16)–(18) is connected to the system of equations

$$x_i=F_i{x}_{i-1}+{\alpha }_{i-1},i=1,\cdots ,N,$$

(24)

where $x_0$ and ${\alpha }_i$ are parameters and $F_i$ is a function of i and eventually also a function of other variables and parameters. This connection can be easily verified. We just have to set

$$\begin{array}{ccc}\hfill F_i& =& \frac{{\phi }_{i-1}}{{\mu }_i+{\phi }_i-{\gamma }_i-{\sigma }_i-{ϵ}_i};\hfill \\ \hfill {\alpha }_{i-1}& =& \frac{{\kappa }_{i-1}+{\lambda }_{i-1}+{\zeta }_{i-1}-{\theta }_{i-1}}{{\mu }_i+{\phi }_i-{\gamma }_i-{\sigma }_i-{ϵ}_i}\hfill \end{array}$$

(25)

in Equation (17) and

$$\begin{array}{ccc}\hfill F_N& =& \frac{{\phi }_{N-1}}{{\mu }_N-{\gamma }_N-{\sigma }_N-{ϵ}_N};\hfill \\ \hfill {\alpha }_{N-1}& =& \frac{{\kappa }_{N-1}+{\lambda }_{N-1}+{\zeta }_{N-1}-{\theta }_{N-1}}{{\mu }_N-{\gamma }_N-{\sigma }_N-{ϵ}_N}\hfill \end{array}$$

(26)

in Equation (18). In addition, we stress that $x_0$ is a free parameter.

Equation (24) leads to a class of probability distributions as follows. From Equation (24) we obtain

$$\begin{array}{c}\hfill x_k=x_0\prod _{i=1}^k{F}_i+(1-{\delta }_{k1})\sum _{j=1}^{k-1}\left({\alpha }_{j-1}\prod _{l=j+1}^k{F}_l\right)+{\alpha }_{k-1},k=1,\cdots ,N\end{array}$$

(27)

Then the amount of substance in the channel will be

$$\begin{array}{c}\hfill x=x_0+\sum _{k=1}^N{x}_k=x_0\left[1+\sum _{k=1}^N\prod _{i=1}^k{F}_i\right]+\sum _{k=1}^N\left[(1-{\delta }_{k1})\sum _{j=1}^{k-1}\left({\alpha }_{j-1}\prod _{l=j+1}^k{F}_l\right)+{\alpha }_{k-1}\right]\end{array}$$

(28)

We can consider probability distribution $y_i=x_i/x$ connected to the amount of substance in the nodes of the studied channel. $y_i$ can be considered as probability values of a distribution of a discrete random variable $\omega$: $y_i=p(\omega =i)$, $i=0,\cdots ,N$. For this distribution, we obtain

$$\begin{array}{ccc}\hfill y_0& =& \frac{x_0}{x_0\left[1+{\displaystyle \sum _{k=1}^N}{\displaystyle \prod _{i=1}^k}{F}_i\right]+{\displaystyle \sum _{k=1}^N}\left[(1-{\delta }_{k1}){\displaystyle \sum _{j=1}^{k-1}}\left({\alpha }_{j-1}{\displaystyle \prod _{l=j+1}^k}{F}_l\right)+{\alpha }_{k-1}\right]},\hfill \\ \hfill y_k& =& \frac{x_0{\displaystyle \prod _{i=1}^k}{F}_i+(1-{\delta }_{k1}){\displaystyle \sum _{j=1}^{k-1}}\left({\alpha }_{j-1}{\displaystyle \prod _{l=j+1}^k}{F}_l\right)+{\alpha }_{k-1}}{x_0\left[1+{\displaystyle \sum _{k=1}^N}{\displaystyle \prod _{i=1}^k}{F}_i\right]+{\displaystyle \sum _{k=1}^N}\left[(1-{\delta }_{k1}){\displaystyle \sum _{j=1}^{k-1}}\left({\alpha }_{j-1}{\displaystyle \prod _{l=j+1}^k}{F}_l\right)+{\alpha }_{k-1}\right]},k=1,\cdots ,N.\hfill \end{array}$$

(29)

As a result of the presence of functions $F_k$, the shapes of the distributions from the class (29) (and the shapes of distributions from class (21), respectively) can be quite different from one another— Figure 3 and Figure 4.

Let us consider Equation (24). We note that we consider relationships connected to probability distributions and our interest will be for the values of $F_i$ and $x_i$, which are nonnegative real numbers. We shall prove the following lemma.

Lemma 1.

Both quantities $F_i$ and ${\alpha }_i$ form Equation (24) (given by Equations (25) and (26)) can have arbitrary values at the same time.

Proof. 

Let us show first that $F_i$ can have an arbitrary (appropriate positive) real value. Indeed, in Equations (25) and (26), we have many parameters and we can choose one of them in such a way that $F_i$ has an appropriate arbitrary value. We obtain from Equation (25).

$${\gamma }_i={\mu }_i+{\phi }_i-{\sigma }_i-{ϵ}_i-\frac{{\phi }_{i-1}}{F_i}$$

(30)

and we can choose ${\gamma }_i$ in such a way that $F_i$ will be an arbitrary nonzero real number (we have just to put the desired value of $F_i$ in Equation (30) and we obtain the corresponding value of ${\gamma }_i$ for any combination of the other parameters and for ${\phi }_{i-1}\ne 0$). There is no problem to obtain also $F_i=0$. For this, we just have to set ${\phi }_{i-1}=0$ in Equation (25). Thus, $F_i$, $i=1,\cdots ,N-1$ can have arbitrary value. The same conclusion holds for $F_N$ and it can be obtained on the basis of Equation (26).

Let us show that the value of ${\alpha }_{i-1}$ in Equation (25) and the value of ${\alpha }_{N-1}$ in Equation (26) can be arbitrary for any value of $F_i$. From Equation (25), we obtain

$${\theta }_{i-1}={\kappa }_{i-1}+{\lambda }_{i-1}+{\zeta }_{i-1}-{\alpha }_{i-1}({\mu }_i+{\phi }_i-{\gamma }_i-{\sigma }_i-{ϵ}_i)$$

(31)

From Equation (31), we obtain the needed value of ${\theta }_{i-1}$ for desired arbitrary value of ${\alpha }_{i-1}$ and fixed values of the other parameters. Note that we can obtain an arbitrary value of ${\alpha }_{i-1}$ by the use of ${\theta }_{i-1}$ and we can do this for an arbitrary value of parameter ${\gamma }_i$.

It follows from the above that $F_i$ and ${\alpha }_{i-1}$ can have arbitrary values at the same time. □

Lemma 1 has interesting consequences. Let us prove:

Lemma 2.

The relationship $x_i=F_i{x}_{i-1}+{\alpha }_{i-1}$ from Equation (24) connected to distribution (21)–(23) describes any sequence $\left\{X_i\right\}$ of nonnegative real numbers.

Proof. 

Let $\left\{X_i\right\}$ be an arbitrary sequence of nonnegative real numbers. We note that $x_0$ connected to distribution (21)–(23) can have arbitrary (nonnegative real) value. Let us set $x_0=X_0$. According to Lemma 1, $F_1$ and ${\alpha }_0$ can have arbitrary values at the same time. Let us consider the following equation for $X_1$: $X_1=F_1{X}_0+{\alpha }_0$. As $F_1$ and ${\alpha }_0$ can have arbitrary values, we can choose these two parameters in such a way that the equation for $X_1$ is satisfied for an arbitrary value of $X_1$ given the value of $X_0$, e.g., ${\alpha }_0=X_1-F_1{X}_0$. We can do the same for $X_2$. We can write the equation $X_2=F_2{X}_2+{\alpha }_1$ and we can choose free parameters $F_2$ and ${\alpha }_1$ in such a way that the equation is satisfied for arbitrary $X_2$ given the value of $X_1$. We can proceed in this way with respect to $X_3$, ⋯ and as we have enough free parameters, we can choose them in such a way that we can obtain any sequence of numbers $\left\{X_i\right\}$. Thus, the relationship $x_i=F_i{x}_{i-1}+{\alpha }_{i-1}$ from Equation (24) connected to distribution (21)–(23) describes any sequence $\left\{X_i\right\}$ of nonnegative real numbers. □

Lemma 2 leads us to the main result:

Theorem 1.

Any truncated discrete probability distribution of the random variable ω that can take values $0,1,\cdots ,N$ is a particular case of the distribution (21)–(23).

Proof. 

Any discrete probability distribution $\left\{Z_i\right\}$, $i=0,\cdots ,N$ of the random variable $\omega$ can be written in the from

$$Z_i=F_i{Z}_{i-1}+{\beta }_{i-1},i=1,\cdots ,N$$

(32)

where $F_i$ has a nonnegative real value and ${\beta }_{i-1}$ has a real value. If the above is not true then there exists a sequence $\left\{Z_i^{*}\right\}$, $i=0,\cdots ,N$ such that it can not be represented in the form (32). We can take $Z_0=Z_0^{*}$ and then we can construct all of the other $Z_i^{*}$ by relationship of the kind (32) fixing the parameters $F_i$ and ${\beta }_{i-1}$ (as it has been done in the Proof of Lemma 2). Thus, we can represent any truncated discrete probability distribution by a sequence of the kind (32). However, the sequence of kind (32) is a particular case of the sequence $x_i=F_i{x}_{i-1}+{\alpha }_{i-1}$ from Lemma 2 for the case when x varies between 0 and 1. Thus, any truncated discrete probability distribution can be described by a sequence of the kind $x_i=F_i{x}_{i-1}+{\alpha }_{i-1}$ from Lemma 2. The sequence $x_i=F_i{x}_{i-1}+{\alpha }_{i-1}$ from Lemma 2 describes probability distributions from the kind (21)–(23). Then the class of probability distributions described by (32) belongs to the class of probability distributions described by (21)–(23). This means that any truncated discrete probability distribution of the random variable $\omega$ is a particular case of the distribution (21)–(23). □

5. Discussion

In this text, we propose a discrete model of motion of substance in a channel of a network which accounts for: (i) exchanges of substance between the nodes of the channel and the nodes of the network, (ii) exchanges of substance between the nodes of the channel and the environment of the network, (iii) exchanges of substance between the flows in the channel and the nodes of the network, and (iv) exchanges of substance between the flows in the channel and the environment of the network. The general model is complicated but nevertheless one can obtain analytical results for the distribution of substance in the nodes of the channel for the particular case of stationary flow of substance through the channel. We obtain the analytical relationship for the general distribution of substance under these conditions of flow for the case of a channel containing a finite number of nodes (in Appendix A we present the corresponding distribution for the case of a channel containing an infinite number of nodes and in Appendix B we discuss the continuous version of the distribution). We show that the obtained distribution contains as particular cases all possible discrete distributions of a random variable which can take values $0,1,\cdots ,N$. This class of distributions contains famous named distributions and entire families of discrete distributions such as families of distributions of Katz, Ord, Kemp, etc. [44,45,46,47,48,49,50,51,52,53,54].

The new point in this study is connected with the presence of the processes (5)–(8) described in Section 2. The processes are connected to exchange of substance among flows between the nodes and the environment of the studied channel. If the values of the parameters listed in points 5–12 of Section 2 are set to 0 then the model is reduced to the model discussed in [14]. Thus, the results obtained in [14] are a particular case of the results obtained in this text. In [17], we discuss a continuous model for flow of substance in a channel of a network. The probability distributions for the case of stationary flow of substance along the channel obtained in [17] are particular cases of the corresponding probability distributions obtained in this text and this is a consequence of Theorem 2 from Appendix A.

The theory discussed above leads to interesting results and can have numerous applications. We have started our research on flows in channels of the network by a study of problems of migration [27,28,29,30,31]. The theory presented above can be applied for many other practical situations connected to the flow of substances in networks. The flows can be, for example, traffic flows and the substances can be resources, goods, or even humans. The discussed model has potential to describe various situations connected to different kinds of flows and one can obtain analytical results for cases of simple flows or can make numerical simulations for the case of more complicated flows.