## Mathematical Applications to Reliability and Maintenance Problems in Engineering Systems

View this Special IssueResearch Article | Open Access

# Reliability Analysis of the Proportional Mean Residual Life Order

**Academic Editor:**Shaomin Wu

#### Abstract

The concept of mean residual life plays an important role in reliability and life testing. In this paper, we introduce and study a new stochastic order called proportional mean residual life order. Several characterizations and preservation properties of the new order under some reliability operations are discussed. As a consequence, a new class of life distributions is introduced on the basis of the anti-star-shaped property of the mean residual life function. We study some reliability properties and some characterizations of this class and provide some examples of interest in reliability.

#### 1. Introduction

Stochastic orders have shown that they are very useful in applied probability, statistics, reliability, operation research, economics, and other related fields. Various types of stochastic orders and associated properties have been developed rapidly over the years. Let be a nonnegative random variable which denotes the lifetime of a system with distribution function , survival function , and density function . The conditional random variable , , is known as the residual life of the system after given that it has already survived up to . The mean residual life (MRL) function of is the expectation of , which is given by

The MRL function is an important characteristic in various fields such as reliability engineering, survival analysis, and actuarial studies. It has been extensively studied in the literature especially for binary systems, that is, when there are only two possible states for the system as either working or failed. Another useful reliability measure is the hazard rate (HR) function of which is given by

The HR function is particularly useful in determining the appropriate failure distributions utilizing qualitative information about the mechanism of failure and for describing the way in which the chance of experiencing the event changes with time. In replacement and repair strategies, although the shape of the HR function plays an important role, the MRL function is found to be more relevant than the HR function because the former summarizes the entire residual life function whereas the latter involves only the risk of instantaneous failure at some time . For an exhaustive monograph on the MRL and HR functions and their reliability analysis, we refer the readers to Ramos-Romero and Sordo-Díaz [1], Belzunce et al. [2], and Lai and Xie [3]. Based on the MRL function, a well-known MRL order has been introduced and studied in the literature. Gupta and Kirmani [4] and Alzaid [5] were among the first who proposed the MRL order. Over the years, many authors have investigated reliability properties and applications of the MRL order in reliability and survival analysis (cf. Shaked and Shanthikumar [6] and Müller and Stoyan [7]). On the other hand, the proportional stochastic orders are considered in the literature to generalize some existing notions of stochastic comparisons of random variables. Proportional stochastic orders as extended versions of the existing common stochastic orders in the literature were studied by some researchers such as Ramos-Romero and Sordo-Díaz [1] and Belzunce et al. [2]. Recently, Nanda et al. [8] gave an effective review of the different partial ordering results related to the MRL order and studied some reliability models in terms of the MRL function.

The purpose of this paper is to propose a new stochastic order called proportional mean residual life (PMRL) order which extends the MRL order to a more general setting. Some implications, characterization properties, and preservation results under weighted distributions of this new order including its relationships with other well-known orders are derived. In addition, two characterizations of this order based on residual life at random time and the excess lifetime in renewal processes are obtained. As a consequence, a new class of lifetime distributions, namely, anti-star-shaped mean residual life (ASMRL) class of life distribution, which is closely related to the concept of the PMRL order, is introduced and studied. A number of useful implications, characterizations, and examples for this class of life distributions are discussed along with some reliability applications. The paper is organized as follows. The precise definitions of some stochastic orders as well as some classes of life distributions which will be used in the sequel are given in Section 2. In that section, the PMRL order is introduced and studied. Several characterizations and preservation properties of this new order under some reliability operations are discussed. In addition, to illustrate the concepts, some applications in the context of reliability theory are included. In Section 3, the ASMRL class of life distributions is introduced and studied. Finally, in Section 4, we give a brief conclusion and some remarks of the current research and its future.

Throughout this paper, the term increasing is used instead of monotone nondecreasing and the term decreasing is used instead of monotone nonincreasing. Let us consider two random variables and having distribution functions and , respectively, and denote by and their respective survival (density) functions. We also assume that all random variables under consideration are absolutely continuous and have 0 as the common left endpoint of their supports, and all expectations are implicitly assumed to be finite whenever they appear. In addition, we use the notations , denotes the equality in distribution, and is the weighted version of according to the weight .

#### 2. Proportional Mean Residual Life Order

For ease of reference, before stating our main results, let us recall some stochastic orders, classes of life distributions, and dependence concepts which will be used in the sequel.

*Definition 1. *The random variable is said to be smaller than in the (i)HR order (denoted as if
(ii) reversed hazard (RH) order (denoted as if
which denotes the reversed hazard (RH) rate order,(iii) MRL order (denoted as ) if

*Definition 2 (Lai and Xie [3]). *The nonnegative random variable is said to have a decreasing mean residual life (DMRL) whenever the MRL of is decreasing.

*Definition 3 (Lariviere and Porteus [9]). *The nonnegative random variable is said to have an increasing generalized failure rate (IGFR) whenever the generalized failure rate function of which is given by is increasing in .

Note that, in view of a result in Lariviere [10], has IGFR property if and only if , for all , or equivalently if , for any .

*Definition 4 (Karlin [11]). *A nonnegative measurable function is said to be totally positive of order 2 () in and , whenever

*Definition 5 (Shaked and Shanthikumar [6]). *A nonnegative function is said to be anti-star-shaped on a set if , for all on and for every . Equivalently, is anti-star-shaped on if is nonincreasing in .

Below, we present the definition of the proportional hazard rate (PHR) order and its related proportional aging class.

*Definition 6 (see Belzunce et al. [2]). *Let and be two nonnegative continuous random variables. It is said that(i) is smaller than in the PHR order (denoted as ) if , for all ,(ii) is increasing proportional hazard rate (IPHR), if , for all .

Consider the situation wherein denotes the risk that the direct insurer faces and the corresponding reinsurance contract. One important reinsurance agreement is quota-share treaty defined as , for . The random variable denotes the risk that an independent insurer faces. Insurers sometimes seek a quota-share treaty when they require financial support from their reinsurers, thus maintaining an adequate relation between net income and capital reserves. Motivated by this, we propose the following new stochastic order.

*Definition 7. *Let and be two nonnegative random variables. The random variable is smaller than in the PMRL order (denoted as ), if , for all .

*Remark 8. *Note that , for all , thus; with we have , that is; the PMRL order is stronger than the MRL order.

The first results of this section provide an equivalent condition for the PMRL order.

Theorem 9. *The following assertions are equivalent:*(i)*;
*(ii)*, for all , and each ;*(iii)

*is increasing in , for all .*

*Proof. *First, we prove that (i) and (ii) are equivalent. Note that the MRL of as a function of is given by , for all and for any . Now, we have if, for all , it holds that
To prove that (ii) and (iii) are equivalent, we have
It is obvious that the last term is nonnegative if and only if , for all and for any

In the context of reliability engineering and survival analysis, weighted distributions are of tremendous practical importance (cf. Jain et al. [12], Bartoszewicz and Skolimowska [13], Misra et al. [14], Izadkhah and Kayid [15], and Kayid et al. [16]). In renewal theory the residual lifetime has a limiting distribution that is a weighted distribution with the weight function equal to the reciprocal of the HR function. Some of the well-known and important distributions in statistics and applied probability may be expressed as weighted distributions such as truncated distributions, the equilibrium renewal distribution, distributions of order statistics, and distributions arisen in proportional hazards and proportional reversed hazards models. Recently, Izadkhah et al. [17] have considered the preservation property of the MRL order under weighted distributions. Here we develop a similar preservation property for the PMRL order under weighted distributions. For two weight functions and , assume that and denote the weighted versions of the random variables of and , respectively, with respective density functions where and . Let and . Then survival functions of and are, respectively, given by

First, we consider the following useful lemma which is straightforward and hence the proof is omitted.

Lemma 10. *Let be a nonnegative absolutely continuous random variable. Then, for any weight function ,
**
where is a weight function of the form .*

Theorem 11. *Let be an increasing function and let increase in , for all . Then
*

*Proof. *Let be fixed. Then, gives . We know by assumption that is increasing and the ratio
is increasing in , when . In view of Theorem 2 in Izadkhah et al. [17], we conclude that . Because of Lemma 10 and because the equality in distribution of and implies the equality in their MRL functions, it follows that . So, for all we have which means that

On the other hand, in many reliability engineering problems, it is interesting to study , the residual life of with a random age . The residual life at random time (RLRT) represents the actual working time of the standby unit if is regarded as the total random life of a warm standby unit with its age . For more details about RLRT we refer the readers to Yue and Cao [18], Li and Zuo [19], and Misra et al. [20], among others. Suppose that and are independent. Then, the survival function of , for any , is given by

Theorem 12. *Let and be two nonnegative random variables. for any which is independent of , if and only if
*

*Proof. *To prove the “if” part, let for all . It then follows that, for all and ,
By integrating both sides of (16) with respect to through the measure , we have
which is equivalent to saying that , for all ’s that are independent of . For the “only if” part, suppose that holds for any nonnegative random variable . Then , for all , follows by taking as a degenerate random variable.

Let be a sequence of mutually independent and identically distributed (i.i.d.) nonnegative random variables with common distribution function . For , denote which is the time of the th arrival and , and let Sup represent the number of arrivals during the interval . Then, is a renewal process with underlying distribution (see Ross [21]). Let be the excess lifetime at time ; that is, . In this context we denote the renewal function by which satisfies the following well-known fundamental renewal equation: According to Barlow and Proschan [22], it holds for all and that

In the literature, several results have been given to characterize the stochastic orders by the excess lifetime in a renewal process. Next, we investigate the behavior of the excess lifetime of a renewal process with respect to the PMRL order.

Theorem 13. *If , for all , then for all .*

*Proof. *First note that , for all , if and only if for any , , and
In view of the identity of (19) and the inequality in (20) we can get
Hence, it holds that, for all , and for any ,
which means for all .

#### 3. Anti-Star-Shaped Mean Residual Life Class

Statisticians and reliability analysts have shown a growing interest in modeling survival data using classifications of life distributions by means of various stochastic orders. These categories are useful for modeling situations, maintenance, inventory theory, and biometry. In this section, we propose a new class of life distributions which is related to the MRL function. We study some characterizations, preservations, and applications of this new class. Some examples of interest in the context of reliability engineering and survival analysis are also presented.

*Definition 14. *The lifetime variable is said to have an anti-star-shaped mean residual life (ASMRL), if the MRL function of is anti-star-shaped.

It is simply derived that ASMRL whenever is decreasing in . Useful description and motivation for the definition of the ASMRL class which is due to Nanda et al. [8] are the following. Consider a situation in which represents the risk that the direct insurer faces and the corresponding reinsurance contract. The ASMRL class provides that the quota-share treaty related to a risk is less than risk itself in the sense of the MRL order. In what follows, we focus on the ASMRL class as a weaker class than the DMRL class to get some basic results. First, consider the following characterization property which can be immediately obtained by Theorem 9(ii).

Theorem 15. *The lifetime random variable is ASMRL if and only if .*

Theorem 16. *The lifetime random variable is ASMRL if and only if
*

*Proof. *Denote , for . The MRL function of is then given by , for all and .In view of the fact that , for all , if and only if
By taking and the above inequality is equivalent to saying that , for all and for any . This means that is ASMRL.

*Remark 17. *The result of Theorem 16 indicates that the family of distributions , is stochastically increasing in with respect to the MRL order if and only if the distribution has an anti-star-shaped MRL function. Another conclusion of Theorem 16 is to say that if and only if , for all .

Theorem 18. *If and if either or has an anti-star-shaped MRL function, then .*

*Proof. *Let and let be ASMRL. Then, we have
Hence it holds that , for all , which means . The proof of the result when is ASMRL is similar by taking the fact that is ASMRL if and only if , for all , into account. Note also that if and only if , for any

The following counterexample shows that the MRL order does not generally imply the ASMRL order and hence the sufficient condition in Theorem 18 cannot be removed.

*Counterexample 1. *Let have MRL , for , and let have MRL , for , and , for . These MRL functions are readily shown not to be ASMRL. We can also see that , for all ; that is, . It can be easily checked now that , for , which means that .

As an obvious conclusion of Theorem 18 above and Theorem 2.9 in Nanda et al. [8], if is DMRL, then is ASMRL. The next result presents another characterization of the ASMRL class.

Theorem 19. *A lifetime random variable is ASMRL if and only if , for each random variable with , which is independent of .*

*Proof. *To prove the “if” part, note that , for each and any . Take , for each one at a time, as a degenerate random variable implying , for all , which means ASMRL. For the “only if” part, assume that has distribution function . From the assumption and the well-known Fubini theorem, for all , it follows that
That is,

The following result presents a sufficient condition for a probability distribution to be ASMRL.

Theorem 20. *Let the lifetime random variable be IGFR. Then, is ASMRL.*

*Proof. *Recall that is IGFR if and only if is increasing in . Because of the identity
we can write, for all ,
Thus, is decreasing in if and only if the ratio
Let
as a function of and of , where

Note that the ratio given in (29) is increasing in if and only if is in . From the assumption, since is increasing, then is in . Also it is easy to see that is in . By applying the general composition theorem of Karlin [11] to the equality of (30), the proof is complete.

To demonstrate the usefulness of the ASMRL class in reliability engineering problems, we consider the following examples.

*Example 21. *The Weibull distribution is one of the most widely used lifetime distributions in reliability engineering. It is a versatile distribution that can take on the characteristics of other types of distributions, based on the value of the shape parameter. Let have the Weibull distribution with survival function
The HR function is given by . Thus we have which is increasing in for all parameter values and hence according to Theorem 20 is ASMRL.

*Example 22. *The generalized Pareto distribution has been extensively used in reliability studies when robustness is required against heavier tailed or lighter tailed alternatives to an exponential distribution. Let have generalized Pareto distribution with survival function
The HR function is given by . Thus we get which is increasing in for all parameter values and so Theorem 20 concludes that is ASMRL.

*Example 23. *Let be a lifetime variable having survival function given by , where is a nonnegative random variable and is the survival function of a lifetime variable , for each . This is called scale change random effects model in Ling et al. [23]. Noting the fact that , for all , is equivalent to saying that is ASMRL, according to Theorem 3.10 of Ling et al. [23] if and is ASMRL, then .

In the context of reliability theory, shock models are of great interest. The system is assumed to have an ability to withstand a random number of these shocks, and it is commonly assumed that the number of shocks and the interarrival times of shocks are* s-*independent. Let denote the number of shocks survived by the system, and let denote the random interarrival time between the th and th shocks. Then the lifetime of the system is given by . Therefore, shock models are particular cases of random sums. In particular, if the interarrivals are assumed to be* s-*independent and exponentially distributed (with common parameter ), then the distribution function of can be written as
where for all (and ). Shock models of this kind, called Poisson shock models, have been studied extensively. For more details, we refer to Fagiuoli and Pellerey [24], Shaked and Wong [25], Belzunce et al. [26], and Kayid and Izadkhah [27].

In the following, we make conditions on the random number of shocks under which has ASMRL property. First, let us define the discrete version of the ASMRL class.

*Definition 24. *A discrete distribution is said to have discrete anti-star-shaped mean residual life (D-ASMRL) property if is nonincreasing in .

Theorem 25. *If , in (34) is D-ASMRL, then with the sf as given in (34) is ASMRL.*

*Proof. *We may note that, for all ,
Hence is ASMRL if and only if
is increasing in , or equivalently if
for and , where
By the assumption, is in , for and . It is also evident that is in , for and . The result now follows from the general composition theorem of Karlin [11].

Lemma 26. *Let be an i.i.d. sample from and let be an i.i.d. sample from . Then implies .*

*Example 27. *Reliability engineers often need to work with systems having elements connected in series. Let be i.i.d. random lifetimes such that has the ASMRL property. Then, according to Theorem 15, , for all . This means that
By appealing to Lemma 26, it follows that , for all . That is, , is ASMRL. Hence, the ASMRL property passes from the lifetime of the series system to the lifetime of its i.i.d. components.

Accelerated life models relate the lifetime distribution to the explanatory variables (stress, covariates, and regressor). This distribution can be defined by the survival, cumulative distribution, or probability density functions. Nevertheless, the sense of accelerated life models is best seen if they are formulated in terms of the hazard rate function. In the following example, we state an application of Theorem 16 in accelerated life models.

*Example 28. *Consider units (not necessarily independent) with lifetimes . Suppose that the units are working in a common operating environment, which is represented by a random vector , independent of , and has an effect on the units of the form
If has support on , then the components are working in a harsh environment, and if they have support on , then the components are working in a gentler environment (see Ma [28]). In a harsh environment let ASMRL for some . Then, Theorem 19 states that, for each with support on , we must have . Thus, by taking , we must have . Hence, by (40) it stands that . With a similar discussion, in a gentler environment if ASMRL for some , then we must have .

In the following we state the preservation property of the ASMRL class under weighted distribution. Let have density function and survival function . The following result states the preservation of the ASMRL class under weighted distributions. The proof is quite similar to that of Theorem 11 and hence omitted.

Theorem 29. *Let be an increasing function and let increase in , for all . Then is ASMRL implying that is ASMRL.*

#### 4. Conclusion

Due to economic consequences and safety issues, it is necessary for the industry to perform systematic studies using reliability concepts. There exist plenty of scenarios where a statistical comparison of reliability measures is required in both reliability engineering and biomedical fields. In this paper, we have proposed a new stochastic order based on the MRL function called proportional mean residual life (PMRL) order. The relationships of this new stochastic order with other well-known stochastic orders are discussed. It was shown that the PMRL order enjoys several reliability properties which provide several applications in reliability and survival analysis. We discussed several characterization and preservation properties of this new order under some reliability operations. To enhance the study, we proposed a new class of life distributions called anti-star-shaped mean residual life (ASMRL) class. Several reliability properties of the new class as well as a number of applications in the context of reliability and survival analysis are included. Our results provide new concepts and applications in reliability, statistics, and risk theory. Further properties and applications of the new stochastic order and the new proposed class can be considered in the future of this research. In particular, the following topics are interesting and still remain as open problems:(i)closure properties of the PMRL order and the ASMRL class under convolution and coherent structures,(ii)discrete version of the PMRL order and enhancing the obtained results related to the D-ASMRL class,(iii)testing exponentiality against the ASMRL class.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors would like to thank two reviewers for their valuable comments and suggestions, which were helpful in improving the paper. The authors would also like to extend their sincere appreciation to the Deanship of Scientific Research at King Saud University for funding this Research Group (no. RG-1435-036).

#### References

- H. M. Ramos-Romero and M. A. Sordo-Díaz, “The proportional likelihood ratio order and applications,”
*Questiio*, vol. 25, no. 2, pp. 211–223, 2001. View at: Google Scholar | MathSciNet - F. Belzunce, J. M. Ruiz, and M. C. Ruiz, “On preservation of some shifted and proportional orders by systems,”
*Statistics and Probability Letters*, vol. 60, no. 2, pp. 141–154, 2002. View at: Publisher Site | Google Scholar | MathSciNet - C. D. Lai and M. Xie,
*Stochastic Ageing and Dependence for Reliability*, Springer, New York, NY, USA, 2006. View at: MathSciNet - R. C. Gupta and S. N. U. Kirmani, “On order relations between reliability measures,”
*Communications in Statistics. Stochastic Models*, vol. 3, no. 1, pp. 149–156, 1987. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - A. A. Alzaid, “Mean residual life ordering,”
*Statistical Papers*, vol. 29, no. 1, pp. 35–43, 1988. View at: Publisher Site | Google Scholar | MathSciNet - M. Shaked and J. G. Shanthikumar,
*Stochastic Orders*, Springer Series in Statistics, Springer, New York, NY, USA, 2007. View at: Publisher Site | MathSciNet - A. Müller and D. Stoyan,
*Comparison Methods for Stochastic Models and Risks*, John Wiley, New York, NY, USA, 2002. View at: MathSciNet - A. K. Nanda, S. Bhattacharjee, and N. Balakrishnan, “Mean residual life function, associated orderings and properties,”
*IEEE Transactions on Reliability*, vol. 59, no. 1, pp. 55–65, 2010. View at: Publisher Site | Google Scholar - M. A. Lariviere and E. L. Porteus, “Selling to a news vendor: an analysis of price-only contracts,”
*Manufacturing and Service Operations Management*, vol. 3, no. 4, pp. 293–305, 2001. View at: Google Scholar - M. A. Lariviere, “A note on probability distributions with increasing generalized failure rates,”
*Operations Research*, vol. 54, no. 3, pp. 602–604, 2006. View at: Publisher Site | Google Scholar | MathSciNet - S. Karlin,
*Total Positivity*, Stanford University Press, Stanford, Calif, USA, 1968. View at: MathSciNet - K. Jain, H. Singh, and I. Bagai, “Relations for reliability measures of weighted distributions,”
*Communications in Statistics, Theory and Methods*, vol. 18, no. 12, pp. 4393–4412, 1990. View at: Publisher Site | Google Scholar | MathSciNet - J. Bartoszewicz and M. Skolimowska, “Preservation of classes of life distributions and stochastic orders under weighting,”
*Statistics & Probability Letters*, vol. 76, no. 6, pp. 587–596, 2006. View at: Publisher Site | Google Scholar | MathSciNet - N. Misra, N. Gupta, and I. D. Dhariyal, “Preservation of some aging properties and stochastic orders by weighted distributions,”
*Communications in Statistics, Theory and Methods*, vol. 37, no. 5, pp. 627–644, 2008. View at: Publisher Site | Google Scholar | MathSciNet - S. Izadkhah and M. Kayid, “Reliability analysis of the harmonic mean inactivity time order,”
*IEEE Transactions on Reliability*, vol. 62, no. 2, pp. 329–337, 2013. View at: Publisher Site | Google Scholar - M. Kayid, I. A. Ahmad, S. Izadkhah, and A. M. Abouammoh, “Further results involving the mean time to failure order, and the decreasing mean time to failure class,”
*IEEE Transactions on Reliability*, vol. 62, no. 3, pp. 670–678, 2013. View at: Publisher Site | Google Scholar - S. Izadkhah, A. H. Rezaei Roknabadi, and G. R. M. Borzadaran, “Aspects of the mean residual life order for weighted distributions,”
*Statistics*, vol. 48, no. 4, pp. 851–861, 2014. View at: Publisher Site | Google Scholar | MathSciNet - D. Yue and J. Cao, “Some results on the residual life at random time,”
*Acta Mathematicae Applicatae Sinica*, vol. 16, no. 4, pp. 435–443, 2000. View at: Publisher Site | Google Scholar | MathSciNet - X. Li and M. J. Zuo, “Stochastic comparison of residual life and inactivity time at a random time,”
*Stochastic Models*, vol. 20, no. 2, pp. 229–235, 2004. View at: Publisher Site | Google Scholar | MathSciNet - N. Misra, N. Gupta, and I. D. Dhariyal, “Stochastic properties of residual life and inactivity time at a random time,”
*Stochastic Models*, vol. 24, no. 1, pp. 89–102, 2008. View at: Publisher Site | Google Scholar | MathSciNet - S. M. Ross,
*Stochastic Processes*, John Wiley, New York, NY, USA, 1996. View at: MathSciNet - R. E. Barlow and F. Proschan,
*Statistical Theory of Reliability and Life Testing*, Silver Spring, Md, USA, 1981. - X. Ling, P. Zhao, and P. Li, “A note on the stochastic properties of a scale change random effects model,”
*Statistics & Probability Letters*, vol. 83, no. 10, pp. 2407–2414, 2013. View at: Publisher Site | Google Scholar | MathSciNet - E. Fagiuoli and F. Pellerey, “Mean residual life and increasing convex comparison of shock models,”
*Statistics & Probability Letters*, vol. 20, no. 5, pp. 337–345, 1994. View at: Publisher Site | Google Scholar | MathSciNet - M. Shaked and T. Wong, “Preservation of stochastic orderings under random mapping by point processes,”
*Probability in the Engineering and Informational Sciences*, vol. 9, no. 4, pp. 563–580, 1995. View at: Publisher Site | Google Scholar | MathSciNet - F. Belzunce, E. Ortega, and J. M. Ruiz, “The Laplace order and ordering of residual lives,”
*Statistics & Probability Letters*, vol. 42, no. 2, pp. 145–156, 1999. View at: Publisher Site | Google Scholar | MathSciNet - M. Kayid and S. Izadkhah, “Mean inactivity time function, associated orderings and classes of life distributions,”
*IEEE Transactions on Reliability*, vol. 63, no. 2, pp. 593–602, 2014. View at: Google Scholar - C. Ma, “Convex orders for linear combinations of random variables,”
*Journal of Statistical Planning and Inference*, vol. 84, no. 1-2, pp. 11–25, 2000. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet

#### Copyright

Copyright © 2014 M. Kayid et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.