Theoretical Economics Letters
Vol.08 No.05(2018), Article ID:83788,12 pages
10.4236/tel.2018.85073

Environmental Policy for Non-Point Source Pollutions in a Bertrand Duopoly

Akio Matsumoto1, Keiko Nakayama2, Ferenc Szidarovszky3

1Department of Economics, Senior Researcher, International Center for further Development of Dynamic Economic Research, Chuo University, Tokyo, Japan

2Department of Economics, Chukyo University, Nagoya, Japan

3Department of Applied Mathematics, University of Pécs, Pécs, Hungary

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: February 15, 2018; Accepted: April 14, 2018; Published: April 17, 2018

ABSTRACT

This study investigates the effectiveness of ambient charges under non-point source (NPS) pollutions in an imperfect competition framework. To this end, following Ganguli and Raju [1] , it constructs a one-stage game and a two-stage game in which Bertrand duopolistic firms choose their best prices and abatement technology, respectively. It is demonstrated in both games that an increase in the ambient charge can lead to a decrease in pollution. This finding indicates that the ambient charge can efficiently control pollution in a Bertrand duopoly.

Keywords:

Non-Point Source Pollution, Abatement Technology, Ambient Charge, Bertrand Competition, Two-Stage Game

1. Introduction

It is now well-known that most of today’s pollutions such as air pollution, river pollution and lake pollution have multiple sources although any pollution originates from a single point source. It is, however, not known how to control these NPS (non-point source) pollutions. Segarson [2] suggests an ambient charge control: the government adopts an environmental policy to establish a cut-off level of the whole pollution and make the following rule, regardless of each firm’s specific emission level, if the actual level of the total level exceeds the cut-off level, then all firms levy the same penalty while if the actual level falls short of the cut-off level, the all firms are awarded the same subsidy. The main purpose of this study is to demonstrate that the ambient charges can control the total amount of NPS pollutions under ála Bertrand imperfect competition. Based on the analysis of Raju and Ganguli [3] in a Cournot duopoly setting, Matsumoto et al. [4] show a “good-natured” effect of ambient charges in an N-firm Cournot framework that an increase in the ambient charge leads to decreases of the total level of industry pollutions. On the other hand, Ganguli and Raju [1] consider the same subject in a Bertrand duopolistic market and numerically exhibit a “perverse” effect that an increase of ambient charge may lead to an increase in pollution in two distinct settings, one in the one-stage game and the other in the two-stage game. Ishikawa et al. [5] consider an ambient charge effect in an N-firm Bertrand framework. This study shows the good-natured effect of ambient charges in Bertrand duopoly, analytically in one-stage game and numerically in two-stage game.

The rest of the paper is organized as follows. In section 2, the optimal price strategies of Bertrand duopolistic firms are derived in one-stage game in which all actions take place simultaneously. In Section 3, the optimal choices of abatement technology at the first stage and prices at the second stage are considered in two-stage game in which the actions take place in sequence. Concluding remarks and further extension of this study are given in Section 4.

2. One-Stage Game

In this section we consider the effect of the ambient charge in one stage game in which the regulator has announced the ambient charge and a cut-off ambient standard while two firms have fixed their pollution abatement technologies. Under this circumstance the firms choose their optimal prices to maximize their profits. Each firm produces a differentiated product. Market demand function for firm i for i , j = 1 , 2 and i j is

q i = a p i + γ p j (1)

where q i denotes good i produced by firm i , p i is the price of q i , p j is the price for the good j and g is a parameter with 0 < γ < 1 measuring the substitutability between two goods.1 We exclude two extreme cases, one with γ = 1 where the two goods are homogenous and the other with γ = 0 where they are independent. The total amount of pollution E generated by the two firms is given by

E = ϕ i q i + ϕ j q j (2)

where ϕ i and ϕ j represent pollution abatement technologies of firms i and j. ϕ i = 1 means the worst technology with 100% pollution while ϕ i = 0 means the best technology with no (0%) pollution. Accordingly, it is assumed that 0 ϕ i 1 .

The profit function of firm i is

π i = p i q i c q i t ( ϕ i q i + ϕ j q j E ¯ ) (3)

where E ¯ > 0 denotes ambient standard specified by the regulator, c q i is the production cost where c > 0 is the common marginal cost of production and t is an ambient charge or tax with 0 t 1 . According to the spirit of the ambient charge, although the two firms’ contributions to pollutions might be different, each firm will pay the identical fine t ( E E ¯ ) if E E ¯ > 0 and receive the identical subsidy t ( E E ¯ ) if E E ¯ < 0 . Substituting (1) into (3) and differentiating the resultant profit function with respect to p i give the first order condition for an interior solution maximizing profit of firm i ,

d π i d p i = a p i + γ p j + ( p i c ) ( 1 ) t [ ϕ i ( 1 ) + ϕ j γ ] = 0

or

2 p i γ p j = a + c + t ( ϕ i γ ϕ j ) . (4)

Maximizing π j with respect to p j presents a similar first-oder condition for firm j. Hence solving the following simultaneous system, which is obtained from first order conditions for firms i and j,

( 2 γ γ 2 ) ( p i p j ) = ( a + c + t ( ϕ i γ ϕ j ) a + c + t ( ϕ j γ ϕ i ) )

yields the Bertrand equilibrium prices,

( p i p j ) = 1 4 γ 2 ( 2 γ γ 2 ) ( a + c + t ( ϕ i γ ϕ j ) a + c + t ( ϕ j γ ϕ i ) )

that are, after arranging the terms,

p i B ( γ , t ) = 1 4 γ 2 { ( 2 + γ ) ( a + c ) + t [ ( 2 γ 2 ) ϕ i γ ϕ j ] } , p j B ( γ , t ) = 1 4 γ 2 { ( 2 + γ ) ( a + c ) + t [ ( 2 γ 2 ) ϕ j γ ϕ i ] } . (5)

Concerning the positivity of the Bertrand price, we have the following results.

Theorem 1 If ϕ i ϕ j or if ϕ i < ϕ j and a + c 1 / 3 , then p i B > 0 .

Proof. If ϕ i ϕ j holds, then

( 2 γ 2 ) ϕ i γ ϕ j ( 2 γ 2 γ ) ϕ i = ( 1 γ ) ( 2 + γ ) ϕ i > 0.

Then the first equation of (5) implies p i B > 0 . Now suppose that ϕ i < ϕ j . If the right hand side of the first equation of (5) is equal to zero, then solving it for ϕ j gives the form of

ϕ j = 2 γ 2 γ ϕ i + ( 2 + γ ) ( a + c ) γ t .

Assumption t 1 implies

( 2 + γ ) ( a + c ) γ t ( 2 + γ ) ( a + c ) γ .

Since γ / ( 2 + γ ) < 1 / 3 for γ < 1 , the term on the right hand side is greater than unity if a + c 1 / 3 under which, for all ϕ j 1 ,

ϕ j < 2 γ 2 γ ϕ i + ( 2 + γ ) ( a + c ) γ t .

The inequality implies that p i B > 0 . Q.E.D.

Differentiating the Bertrand price of firm i with respect to t reveals that the sign of the derivative is the same as the sign of the terms in the square brackets in (5),

s i g n [ d p i B d t ] = s i g n [ ( 2 γ 2 ) ϕ i γ ϕ j ] .

Hence the effect caused by a change in the ambient charge on the Bertrand prices can summarized as follows.

Theorem 2 A firm with a larger or equal abatement technology positively responds to a change in the ambient charge whereas the response of a firm with a smaller abatement technology is ambiguous,

If ϕ i > ϕ j , then d p i B d t > 0 and d p j B d t 0, If ϕ i = ϕ j , then d p i B d t > 0 and d p j B d t > 0, If ϕ i < ϕ j , then d p i B d t 0 and d p j B d t > 0.

Proof. If ϕ i ϕ j , then for firm i , the bracketed terms in (5) are

( 2 γ 2 ) ϕ i γ ϕ j > ( 2 + γ ) ( 1 γ ) ϕ i > 0 implying that d p i B d t > 0

and for firm j, from the bracketed terms in the second equation of (5) to be equal to zero, we can define the ratio of the abatement technologies

( ϕ j ϕ i ) = γ 2 γ 2 .

This ratio is less than unity implying

d p j B d t 0 according to ϕ j ϕ i ( ϕ j ϕ i ) when ϕ j ϕ i 1.

The same procedure can be applied for the case of ϕ i < ϕ j . Q.E.D

Substituting the Bertrand prices into the demand functions in (1) presents the Bertrand outputs of firm i and j,

q i B = a p i B + γ p j B , q j B = a p j B + γ p i B . (6)

To check whether q i B is positive, we subtract the second equation of (6) from the first equation to obtain

q i B q j B = ( 1 + γ ) ( p j B p i B ) (7)

where, in the same way, from (5)

p i B p j B = 1 + γ 2 + γ t ( ϕ i ϕ j ) . (8)

Concerning the positivity of the Bertrand output, we have the following:

Theorem 3 If a c + 1 , then q i B > 0 and q j B > 0 for 0 ϕ k 1 for k = 1 , 2 .

Proof. If ϕ i ϕ j holds, then (8) leads to p i B p j B with which (7) implies q i B q j B . On the other hand, the first equation of (6) with the forms of p i B and p j B given in (5) is reduced to

q i B = a 3 γ 3 + a 2 γ 2 + a 1 γ + a 0 ( 4 γ 2 ) (9)

where

a 3 = t ϕ j > 0 ,

a 2 = c < 0 ,

a 1 = ( a + c + 3 t ϕ j ) < 0 ,

a 0 = 2 ( a c t ϕ i ) 0

The direction of the last inequality is due to a c + 1 and t ϕ i 1 . Let the numerator of (9) be f ( γ ) . Then due to Descartes’ rule of sign, f ( γ ) = 0 has only one positive root, γ 0 > 0 . Substituting γ = 1 gives

f ( 1 ) = a 3 + a 2 + a 1 + a 0 = 2 ( a + t ϕ j ) ( a t ϕ i ) < 0.

The last inequality means that f ( γ ) < 0 = f ( γ 0 ) for γ < 1 , with which then (9) leads to q i B > 0 , implying that q j B > 0 as well. If ϕ j > ϕ i , then interchanging the two firms generates the same result. Q.E.D.

Now consider the effect of a change in the ambient charge on each output level.

d q i B d t = d p i B d t + γ d p j B d t = 1 4 γ 2 [ 2 ϕ i + γ ( 3 γ 2 ) ϕ j ] , d q j B d t = d p j B d t + γ d p i B d t = 1 4 γ 2 [ 2 ϕ j + γ ( 3 γ 2 ) ϕ i ] . (10)

Concerning the ambient charge effect on output, we have the following results.

Theorem 4 A firm with a larger or equal abatement technology negatively responds to a change in the ambient charge whereas the response of a firm with a smaller abatement technology is ambiguous:

If ϕ i > ϕ j , then d q i B d t < 0 and d q j B d t 0, If ϕ i = ϕ j , then d q i B d t < 0 and d q j B d t < 0, If ϕ i < ϕ j , then d q i B d t 0 and d q j B d t < 0.

Proof. If ϕ i = ϕ j , then γ ( 3 γ 2 ) < 2 for γ < 1 immediately implies the results. If ϕ i > ϕ j , then γ ( 3 γ 2 ) > 0 for γ > 0 indicates that for firm i,

2 ϕ i + γ ( 3 γ 2 ) ϕ j < [ 2 + γ ( 3 γ 2 ) ] ϕ i < 0 implying d q i B d t < 0

and for firm j, there is a threshold value of the parameter ratio such as

γ ( 3 γ 2 ) 2 = ( ϕ j ϕ i ) < 1 leading to d q j B d t 0 according to ϕ j ϕ i ( ϕ j ϕ i ) and ϕ j ϕ i < 1.

The same procedure can be applied for the case of ϕ i < ϕ j . Q.E.D.

It should be noticed first that firm i with a larger ϕ i has an inefficient technology because it generates larger emission. Theorems 2 and 4 imply that, the firm with inefficient technology exhibits natural response to the change in the ambient charge, that is, it increases price and decreases output. On the other hand the firm with efficient technology responds ambiguously. It should be noticed second that these firm-specific responses are non-observable for the regulator which can see only the total amount in the case of NPS pollution.

The total level of pollution at the Bertrand equilibrium is obtained by substituting q i B and q j B into (2)

E B = ϕ i q i B + ϕ j q j B .

Concerning the effect of a change in the ambient charge on the total pollution, we have the following result.

Theorem 5 An increase in the ambient charge decreases the total level of pollution,

d E B d t < 0.

Proof. Differentiating E B with respect t gives

d E B d t = ϕ i d q i B d t + ϕ j d q j B d t = 2 4 γ 2 { ϕ i 2 + ϕ j 2 γ ( 3 γ 2 ) ϕ i ϕ j } .

Notice that γ ( 3 γ 2 ) 2 and the equality holds when γ = 1 . Hence

ϕ i 2 + ϕ j 2 γ ( 3 γ 2 ) ϕ i ϕ j > ϕ i 2 + ϕ j 2 2 ϕ i ϕ j = ( ϕ i ϕ j ) 2 0.

Therefore we arrive at the result where the strict inequality is due to the assumption, γ < 1 . Q.E.D.

Although Theorem 4 implies a possibility of the perverse effect on emission of the individual firm with the efficient abetment technology, Theorem 5 implies that the total effect is always negative, implying that the negative effect of the inefficient firm dominates the positive effect of the efficient firm.

3. Two-Stage Game

In this section the firms and the regulator take actions in two stages. At the first stage, each firm determines its optimal abatement technology whereas the regulator announces the ambient charge and the cut-off level of total pollution. Then at the second stage, the firms choose their prices to maximize their profits, given the ambient charge, the cut-of level and their abatement technologies.

The decision-making at the second stage have been already considered in the one-stage game. Given the Bertrand prices and outputs in (5) and (6), we consider the actions of choosing abatement technology at the first stage. The Bertrand profit function of firm i under Bertrand prices and Bertrand output is defined as

π i B = p i B q i B [ c q i B + ( 1 ϕ i ) 2 ] t ( ϕ i q i B + ϕ j q j B E ¯ ) (11)

for i , j = 1 , 2 and i j . Notice that there is a small difference between the definitions of (3) and (11). There is a term ( 1 ϕ i ) 2 in (11) and no such a term in (3). This term reflects the cost associated with selecting the abatement technology.2 At the second stage, the abatement technology has been already selected somehow and thus it does not affect the determination of the optimal price whether this cost is included or not. However this cost function is effective for choosing technology at the first stage where the firm i determines its abatement technology, ϕ i , so as to maximize its Bertrand profit. Differentiating (11) with respect to ϕ i gives the first-order condition for the optimal level of ϕ i ,

π i B ϕ i = π i B p i p i B ϕ i + π i B p j p j B ϕ i + π i B ϕ i | p i B , p j B ; given = 0

where

π i B p i = 0 bythefirstorderconditionatthesecondstage, π i B p j = γ p i B c γ t ( γ ϕ i ϕ j ) , π i B ϕ i | p i B , p j B ; given = 2 ( 1 ϕ i ) t ( a p i B + γ p j B ) , p j B ϕ i = γ t 4 γ 2 .

Arranging the terms in π i B / ϕ i = 0 with taking account of the forms of the Bertrand prices in (5) yield the modified FOC at the first stage,

2 [ ( 2 t ) 2 ( 4 γ 2 ) 2 ] ϕ i ( γ 4 9 γ 2 + 16 ) γ t 2 ϕ j = 2 ( 4 γ 2 ) 2 + t d (12)

where

d = 4 ( 2 + γ ) [ a ( 1 γ ) c ]

In the same way, arranging the terms in π j B / ϕ j = 0 yields the modified FOC for firm j

( γ 4 9 γ 2 +16 )γ t 2 ϕ i +2[ ( 2t ) 2 ( 4 γ 2 ) 2 ] ϕ j =2 ( 4 γ 2 ) 2 +td. (13)

Solving (12) and (13) simultaneously for ϕ i and ϕ j presents the optimal abatement technology for both firms,

ϕ i = ϕ j = ϕ ( γ , t ) = 2 ( 4 γ 2 ) 2 + t d 2 [ ( 2 t ) 2 ( 4 γ 2 ) 2 ] ( γ 4 9 γ 2 + 16 ) γ t 2 . (14)

Although the form of the solution seems to be highly complicated, the following result is obtained.

Proposition 1 If 2 / 3 a > c , then the optimal abatement technology ϕ ( γ , t ) is positive for 0 < γ < 1 and 0 t 1 .

Proof. It is numerically confirmed that the denominator of (14) is negative for 0 < γ < 1 and 0 t 1 . 3 The numerator is rewritten as

2 ( 2 + γ ) { 2 [ a ( 1 γ ) c ] t f ( γ ) } (15)

where

f ( γ ) = ( 2 γ ) 2 ( 2 + γ ) and 3 < f ( γ ) < 8 for 0 < γ < 1.

Assumption a > c implies a ( 1 γ ) c > 0 . The sign of the terms in the braces of (15) is negative for t = 0 (i.e., f ( γ ) < 0 ) and is also negative for t = 1 if 2 a < 3 . As the terms satisfies the inequality

2 [ a ( 1 γ ) c ] f ( γ ) < 2 a f ( γ ) .

and the lower bound of f ( γ ) is 3, then the right hand side is negative for 0 < γ < 1 if 2 a 3 . Thus the sign of the left hand side is also negative. Hence 2 [ a ( 1 γ ) c ] t f ( γ ) < 0 for 0 t 1 . Therefore (14) with the negative denominator and the negative numerator implies ϕ ( γ , t ) > 0 . Q.E.D.

Notice that the following set is not empty,

{ ( a , c ) | 0 < a 3 / 2 and 0 < c a 1 } ,

implying that the conditions imposed on the values of a and c given in Theorem 3 and Proposition 1 are compatible. Substituting ϕ into the Bertrand prices in (5) gives the optimal Bertrand price

p i = p j = p ( γ , t ) = 1 2 γ [ ( a + c ) + ( 1 γ ) t ϕ ( t ) ] (16)

that is clearly positive for 0 < γ < 1 and 0 t 1 , which is summarized as follows.

Proposition 2 If 3 / 2 a > c , then the optimal Bertrand price p ( γ , t ) is positive for 0 < γ < 1 and 0 t 1 .

The optimal Bertrand output is obtained by substituting p ( γ , t ) into the demand function, (1),

q ( γ , t ) = a ( 1 γ ) p ( γ , t ) . (17)

Concerning the optimal output, we can have the following,

Proposition 3 q ( γ , t ) < 0 for γ = 0 and t = 1 .

Proof. Substituting γ = 0 into (14) and (16) present

p ( 0 , t ) = 1 2 [ ( a + c ) + t t ( a c ) 4 2 ( t 2 1 ) ]

that is, in turn, substituted into (17) to obtain

q ( 0 , t ) = 1 4 ( t 2 1 ) [ 3 ( a c ) t 2 + 4 t 2 ( a c ) ]

then

lim t 1 , t < 1 q ( 0 , t ) = lim t 1 , t < 1 a c + 4 4 ( t 2 1 ) < 0

where the numerator is positive and the denominator is negative. This completes the proof. Q.E.D.

4 a = 3 / 2 and c = 1 / 2 are taken for this and the following examples. The vertical and horisontal intersecpts of the black curve are

γ 0 0.187 and t 0 0.563.

Proposition 3 implies that q ( γ , t ) is inevitably negative for a small neighborhood of point ( 0,1 ) . Numerical confirmation of Proposition 3 is given in Figure 1 in which the Bertrand output is negative in the shaded region and positive otherwise.4 The shaded region should be eliminated from further considerations.

Figure 1. q* > 0 in the white region and q* < 0 in the shaded region.

Finally, the optimal level of total pollution is

E ( γ , t ) = 2 ϕ ( γ , t ) q ( γ , t ) (18)

that is apparently negative if q ( γ , t ) < 0 , that is, in the shaded region in Figure 1. Its derivative with respect to t is

d E d t = 2 q ϕ t [ t ϕ ϕ t + t q q t ] (19)

where ϕ / t and q / t can be of either sign and thus the sign of the terms in the square bracket seems ambiguous. Numerically, as as shown in Figure 2(a) and Figure 2(b), each derivative is confirmed to be positive in the corresponding shaded region in which

γ ^ 0 0.160 and t ^ 0 0.53 in Figure 2(a)

and

γ ¯ 0 0.594 and t ¯ 0 = 0.5 in Figure 2(b).

The dotted upward-sloping curve in Figure 2(a) is described by the q = 0 curve. The shaded region is included in the region with q < 0 so that the derivative of the optimal technology is negative in the feasible region with q > 0 . Hence the sign of the tax-derivative of the total emission level is definitely negative in the white region in Figure 2(b) as both derivatives are negative. On the other hand, it is sensitive to the relative magnitude between the elasticities of the optimal ambient technology and the optimal output with respect to t in the shaded region in Figure 2(b).

Given the parameter values, Figure 3 illustrates the values of d E / d t against a pair of ( γ , t ) with 0 < γ < 1 and 0 t 1 . Since it is seen that any point on the 3D surface in Figure 3 is negative, it is numerically shown that d E / d t < 0 for 0 < γ < 1 and 0 t 1 . According to equation (19), the sign of d E / d t is

(a) (b)

Figure 2. Divisions of the nonnegative (γ, t) region.

Figure 3. Effective ambient charge under a = 3/2 and c = 1/2.

5We obtain this results with various values of a and c . However, we are unable to prove it analytically so this is an numerically-shown result.

ambiguous in the shaded region of Figure 2(b). However, Figure 3 implies that the elasticity of ϕ in absolute value is larger than that of q , leading to the negative sign of d E / d t . Therefore we have our main result that an increase in ambient charge always decreases the total level of optimal pollution.5 We summarize this result.

Proposition 4 It is numerically confirmed that a change in the ambient charge has the good-natured effect,

d E d t < 0 for 0 < γ < 1 and 0 t 1.

4. Concluding Remark

In this paper, we reconsider the “perverse” effect caused by a change in ambient charges shown by Ganguli and Raju [1] . To this end, following their basic framework, we first re-examine the effect in one-stage game in which the Bertrand firms determine their prices so as to maximize their profits, given the abatement technology. Our first result analytically demonstrates that an increase of ambient charges decreases the total level of NPS pollutions. We then turn attention to the effect in two-stage game in which the optimal abatement technology is selected at the first stage and the optimal prices are determined at the second stage. Our second result numerically shows the good-natured effect on the total level of pollution. With these results, we conclude that the ambient charge might be an efficient method to control NPS pollutions even in a duopoly Bertrand market.

Funding

The authors highly acknowledge the financial supports from the MEXT-Supported Program for the Strategic Research Foundation at Private Universities 2013-2017, the Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research (C), 16K03556, 18K01631 and 17H01931) and Chuo University (Joint Research Grant). The usual disclaimer applies.

Cite this paper

Matsumoto, A., Nakayama, K. and Szidarovszky, F. (2018) Environmental Policy for Non-Point Source Pollutions in a Bertrand Duopoly. Theoretical Economics Letters, 8, 1058-1069. https://doi.org/10.4236/tel.2018.85073

References

  1. 1. Ganguli, S. and Raju, S. (2012) Perverse Environmental Effects of Ambient Charges in a Bertrand Duopoly. Journal of Environmental Economics and Policy, 1, 289-296. https://doi.org/10.1080/21606544.2012.714972

  2. 2. Segerson, K. (1988) Uncertainty and Incentives for Non-Point Pollution Control. Journal of Environmental Economics and Management, 15, 87-98. https://doi.org/10.1016/0095-0696(88)90030-7

  3. 3. Raju, S. and Ganguli, S. (2013) Strategic Firm Interaction, Return to Scale, Environmental Regulation and Ambient Charges in a Cournot Duopoly. Technology and Investment, 4, 113-122. https://doi.org/10.4236/ti.2013.42014

  4. 4. Matsumoto, A., Szidarovszky, F. and Yabuta, M. (2018) Environmental Effects of Ambient Charge in Cournot Oligopoly. Journal of Environmental Economic Policy, 7, 41-56. https://doi.org/10.1080/21606544.2017.1347527

  5. 5. Ishikawa, T., Matsumoto, A. and Szidarovszky, F. (2017) Regulation of Non-Point Source Pollution under n-Firm Bertrand Competition. IERCU DP #284, Institute of Economic Research, Chuo University.

NOTES

1Intuitively, x i and x j are substitutes in the following sence that Δ p j > 0 implies Δ x j < 0 and Δ p j > 0 with γ > 0 implies Δ x i > 0. Hence to a change in price j , the quantity response of good i runs in an opposite direction of the response of good j .

2See Raju and Ganguli [3] for this cost.

3Numerical calculations of this and any others that follow are done with Mathematica, ver. 11.