Plane Thermoelastic Waves in Infinite Half-Space Caused Decision Making: Applications in Management and Engineering Vol. 2, Issue 1, 2019, pp. 1-12. ISSN: 2560-6018 eISSN: 2620-0104 DOI: https://doi.org/10.31181/dmame1901001a * Corresponding author. E-mail addresses: krish.math23@gmail.com (K. Adhikary), jaga.math23@gmail.com (J. Roy), kar_s_k@yahoo.com (S. Kar) Newsboy problem with birandom demand Krishnendu Adhikary 1, Jagannath Roy 1* and Samarjit Kar 1 1 Department of Mathematics, National Institute of Technology Durgapur, Durgapur, India Received: 15 October 2018; Accepted: 12 December 2018; Available online: 18 December 2018. Original scientific paper Abstract: Estimation of accurate product demand in a single period inventory model (SPIM) is an essential prerequisite for successfully managing the supply chain in large and medium merchandise. Managers/ decision makers (DMs) often find it difficult to forecast the exact inventory level of a product due to complex market situations and its volatility caused by several factors like customers uncertain behavior, natural disasters, and uncertain demand information. In order to make fruitful decisions under such complicated environment, managers seek applicable models that can be implemented in profit maximization problems. Many authors studied SPIM (also known as newsboy problem) considering the demand as a normal random variable with fixed mean and variance. But for more practical situations the mean demand also varies time to time yielding two-folded randomness in demand distribution. Thus, it becomes more difficult for DMs to apprehend the actual demand having two-folded random/birandom distribution. A blend of birandom theory and newsboy model has been employed to propose birandom newsboy model (BNM) in this research to find out the optimal order quantity as well as maximize the expected profit. The practicality of the projected BNM is illustrated by a numerical example followed by a real case study of SPIM. The results will help DMs to know how much they should order in order to maximize the expected profit and avoid potential loss from excess ordering. Finally, the BNM will enhance the ability of the managers to keep parity of product demand and supply satisfying customersโ€™ needs effectively under uncertain environment. Key words: Newsboy problem, Uncertain variable, Birandom variable, Expectation. 1. Introduction The classical newsboy problem (CNP) aims at determining the optimal order quantity of products, which minimize the expected total cost and / or maximize the mailto:krish.math23@gmail.com mailto:jaga.math23@gmail.com Adhikary et al./Decis. Mak. Appl. Manag. Eng. 2 (1) (2019) 1-12 2 expected profit in a SPIM. Thus, the characteristics of a SPIM become more complicated due to complex market situations and its volatility caused by several factors like customers uncertain behavior. In response to handle such problems, the CNP is extended in many directions and several researchers have proposed modified versions of CNP. Some of them made proper and thoughtful use of probability theory to solve newsboy problems where product demands follow Poisson and Normal distribution (Hadley and Whitin, 1963), Weibull distribution (Tadikamalla, 1978), Erlang distribution (Mahoney and Sivazlian, 1980), compound Poisson distribution (Dominey and Hill, 2004). Gallego and Moon (1993, 1994) analyzed distribution-free newsboy problems where DMs have no idea about the demand distribution. The only information they have are mean and variance of demand. Their newsboy models can be used as strategic tools in deciding the stock of products that have a limited selling period. This paragraph is dedicated to articulate the recent developments of SPIM and find the literature gap. Agarwal and Seshadri (2000) worked in CNP where they assumed the demand distribution as a function of selling price and the objective of the risk- averse retailers. To maximize DMs expected utility they presented two models for comparing the risk-neutral retailers (who charge a higher price for less order) with a risk-averse retailer (who charge a low price). The distribution-free newsboy problem under the worst-case and best-case scenario was revealed by Kamburowski (2014). Further, Kamburowski (2015) studied a newsboy problem where the distribution of the random variable is only known when to be non-skewed with given support, mean and variance. For the distribution-free newsboy problem, Gler (2014) extended the model developed by Lee and Hsu (2011). Here, the authors showed the expected profit increases with a proper advertisement policy while an unorganized advertising policy can have its backfire effect or make a very small improvement of the optimal profit value. Ding (2013) proposed a chance constraint multi-product newsboy problem with uncertain demand and uncertain storage capacity. Abdel-Aal et al. (2017) studied a multi-product newsboy problem assuming the service level as a constraint to offer the DMs to select the market to serve. Watt and Vzquez (2017) considered newsboy problem under two new assumptions. First, they assumed that the wholesaler is an expert who sets the wholesale price optimally and a newsboy can return the unsold item with some salvage value. In the second one, the salvage value acts as a standard insurance demand. Sun and Guo (2017) built a newsboy model with fuzzy random demand based on fuzzy random expected value model. Vipin and Amit (2017) proposed a loss aversion SPIM under alternative option and proved the rationality of the decision maker to predict the order quantity by imposing loss aversion in the newsboy model with the change of selling price and purchase cost factors. Additionally, they showed the models based on utility functions perform better in forecasting the rational behavior due to loss aversion. Natarajan et al. (2017) allowed asymmetry and ambiguity in newsvendor models. The effects of decision makers emergency order in SPIM are discussed and analyzed by Pando et al. (2013) and Zhang et al. (2017). Zhang et al. (2017) compared two ways to treat the excess demand and came up with the better one. In the aforementioned works, the product demand is assumed to be either normally distributed with ๐‘(๐œ‡, ๐œŽ 2) or exponentially distributed with constant mean (๐œ†) or somewhere distribution free. But the DMs face difficulties to forecast the exact demands of products in many practical problems. The demand distribution changes dynamically from time to time, which yields randomness in the mean demand. For example, the demand (๐ท) is normally distributed with ๐ท ~ ๐‘(๐œ‡ ; 400) where ๐œ‡ ~ ๐‘ˆ(3000; 4000) or ๐ท ~ ๐‘(๐œ‡ ; 400) with ๐œ‡ ~ ๐‘’๐‘ฅ๐‘(0.0003). However, in newsboy Newsboy problem with birandom demand 3 problems, it is largely appreciated to consider demand variable having the standard normal distribution. From the probabilistic viewpoint and the above arguments, it would be more realistic to assume the values of ๐œ‡ and ๐œŽ 2 should also be treated as random variables. For such cases, it is more convincing and practical to represent the product demand as birandom variable effectively captures two-folded randomness. Traditional probabilistic approaches cannot handle such complicated real world problems. In response, Peng and Liu (2007) developed a birandom theory to tackle such problems. Zhaojun and Liu (2013) showed the common formula on birandom variable and further, several researchers used this theory to solve inventory problems in the birandom environment. In the recent years, birandom theory is widely accepted as the mathematical language of uncertainty. Some more notable extensions and applications of newsboy problem can be found in the recent literature (Abdel-Malek and Otegbeye 2013; Chen and Ho 2013). From its inception until today, birandom theory has progressed and been applied to different areas. Xu and Zhou (2009) introduced a class of multiple objective decision making problems using birandom variables and by transforming the birandom uncertain problem into its crisp equivalent form through expected value operator and used it in the flow shop scheduling problem. A Portfolio selection problem is analysed by Yan (2009) assumed the security returns as birandom variables. Xu and Ding (2011) developed the general chance constrained multi objective linear programming model with birandom parameters for solving a vendor selection problem. They presented a crisp equivalent model for a special case and gave a traditional method to solve the crisp model. Wang et al. (2012) established a class of job search problem with birandom variables, where the job searcher examined job offers from a finite set jobs having equivalent probability. A multi-mode resource constrained project scheduling problem (Zhang and Xu, 2013) of drilling grounding construction projects considering the uncertain parameters as birandom variables. In the earlier year, Xu et al. (2012) used birandom theory to develop the nonlinear multi objective bi-level models for finding the minimum cost in a network flow problem dealing a large scale construction project. Tavana et al. (2013) measured the efficiencies of decision making units after developing a data envelopment analysis (DEA) model with birandom input and output data. Nevertheless, many more real life applications can be found in the following literature: a multi-objective birandom inventory problem (Tao and Xu, 2013), a birandom multi-objective scheduling problem (Xu et al., 2013) in ship transportation, optimal portfolio selection with birandom returns(Cao and Shan, 2013), a modified genetic algorithm (Maity et al., 2015), chance-constrained programming model for municipal waste management with birandom variables (Zhou et al., 2015), and the CCUS (carbon capture, utilization, and storage) management system in birandom environment (Wang et al., 2017). To the best of our knowledge, no researcher has investigated the SPIM with birandom demand till date. With these considerations, we discuss an inventory model with single period considering the demand for the birandom variable. We solve a real problem using the presented model in searching the optimum order quantity for maximum profit by using the expected value model. The principal aim of this paper is to deliver basic knowledge and suggest precise results of a complex inventory practical problem for the management. The remaining part of our paper is presented in the following way. Section 2 presents some basic knowledge of birandom variable and related theorems. Section 3 introduces the birandom simulation for finding the expected value of the birandom http://www.dmame.org/index.php/dmame/article/view/24 Adhikary et al./Decis. Mak. Appl. Manag. Eng. 2 (1) (2019) 1-12 4 variable. In section 4, we provide a mathematical model for newsboy problem with birandom demand. A numerical example is discussed in section 5. To validate the applicability of the proposed model we discuss a real case study in section 6. Finally, the conclusion and future research directions are presented in section 7. 2. Preliminaries In this section we discuss the basic notations of birandom variables. 2.1. Birandom Variable Roughly speaking a birandom variable is a random variable of a random variable, i.e., a function defined from a probability space to a collection of random variables is said to be a birandom variable. The formal definition of birandom variable and related theorems are defined in the following way. Definition 1 (Peng and Liu, 2007). A birandom variable ๐œ‰ is a mapping from a probability space (๐›บ, ๐ด, ๐‘ƒ๐‘Ÿ) to a collection ๐‘† of a random variable such that for any Borel subset ๐ต of the real line โ„œ the induced function ๐‘ƒ๐‘Ÿ {๐œ‰(๐œ”) โˆˆ ๐ต} is a measurable function with respect to ๐œ”. For each given Borel subset ๐ต of the real line โ„œ, the function ๐‘ƒ๐‘Ÿ {๐œ‰(๐œ”) โˆˆ ๐ต} is a random variable defined on the probability space (๐œ”, ๐ด, ๐‘ƒ๐‘Ÿ). Lemma 1 (Peng and Liu 2007). Let an n dimensional birandom vector ๐œ‰ = (๐œ‰1, ๐œ‰2, โ€ฆ , ๐œ‰๐‘› ) and ๐‘“: โ„œ ๐‘› โ†’ โ„œ be a measurable function. Then ๐‘“(๐œ‰) is a birandom variable. Let two probability spaces (ฮฉ1, ๐ด1, Pr1) and (ฮฉ2, ๐ด2, Pr2), ๐œ‰1 and ๐œ‰2 be two birandom variables respectively taken from that probability spaces. Then ๐œ‰ = ๐œ‰1 + ๐œ‰2 is a birandom variable on (ฮฉ1 ร— ฮฉ2, ๐ด1 ร— ๐ด2, Pr1 ร— Pr2) defined by ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ1 2 1 1 2 2 1 2 1 2, , , ฮฉ ฮฉ๏ธ ๏ท ๏ท ๏ธ ๏ท ๏ธ ๏ท ๏ท ๏ท๏€ฝ ๏€ซ ๏€ข ๏ƒŽ ๏‚ด Widely, for the n-tuple operation on birandom variables defined as follows. Let a Borel measurable function defined as ๐‘“: โ„œ๐‘› โ†’ โ„œ and ๐œ‰๐‘– be birandom variable defined on (ฮฉ๐‘– , ๐ด๐‘– , Pri), ๐‘– = 1,2, โ€ฆ , ๐‘› respectively. Then ๐œ‰ = ๐‘“(๐œ‰1, ๐œ‰2, โ€ฆ , ๐œ‰๐‘› ) is birandom variable on (ฮฉ1 ร— ฮฉ2 ร— โ€ฆ ร— ฮฉ๐‘› , ๐ด1 ร— ๐ด2 ร— โ€ฆ ร— ๐ด๐‘›, Pr1 ร— Pr2 ร— โ€ฆ ร— Prn), defined by ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ๏€จ ๏€ฉ ๏€จ ๏€ฉ 1 2 1 1 2 2 1 2 1 2 , , , , , ., , , ., ฮฉ ฮฉ ฮฉ n n n n n f๏ธ ๏ท ๏ท ๏ท ๏ธ ๏ท ๏ธ ๏ท ๏ธ ๏ท ๏ท ๏ท ๏ท ๏‚ผ ๏€ฝ ๏‚ผ ๏€ข ๏‚ผ ๏ƒŽ ๏‚ด ๏‚ด๏‚ผ๏‚ด 2.2. Expected Value of birandom variables We can transform the complex uncertain problems into their equivalent crisp models, which will be easier to solve. Generally, expected value operator is applied to transform the birandom problem into its deterministic value model for calculating the objective functional value. First, we present the definition of the expected value operator of a birandom variable and then the expected value model of SPIM. The effective tool of an uncertain variable is expectation, which is applied in a different field of applications. Therefore, the idea of expected value of birandom variable is useful. The expected value operator of birandom variable is defined as follows. Newsboy problem with birandom demand 5 Definition 2 (Peng and Liu, 2007). Let ๐œ‰ be a birandom variable defined on the probability (ฮฉ, ๐ด, ๐‘ƒ๐‘Ÿ). Then the expected value of birandom variable ๐œ‰ is defined as ๏› ๏ ๏€จ ๏€ฉ๏ป ๏ฝ ๏€จ ๏€ฉ๏ป ๏ฝ 0 0 Pr ฮฉ| Pr ฮฉ|E E t dt E t dt๏ธ ๏ท ๏ธ ๏ท ๏ท ๏ธ ๏ท ๏‚ฅ ๏€ญ๏‚ฅ ๏€ฝ ๏ƒŽ ๏ƒฉ ๏ƒน ๏‚ณ ๏€ญ ๏ƒŽ ๏ƒฉ ๏ƒน ๏‚ฃ๏ƒซ ๏ƒป ๏ƒซ ๏ƒป๏ƒฒ ๏ƒฒ (1) Provided that at least one of the above two integrals is finite. Lemma 2 (Peng and Liu, 2007). Let ๐œ‰ be a birandom variable defined on the probability (๐œ”, ๐ด, ๐‘ƒ๐‘Ÿ). If the expected value ๐ธ[๐œ‰(๐œ”)] of the random variable ๐œ‰(๐œ”) is finite for each ๐œ”, then ๐ธ[๐œ‰(๐œ”)] is a random variable on (๐œ”, ๐ด, ๐‘ƒ๐‘Ÿ). Lemma 3 (Peng and Liu, 2007). Let us consider two birandom variable ๐œ‰ and ๐œ‚ with finite expected value, then for any two real numbers a and b, we have ๏› ๏ ๏› ๏ ๏› ๏E a b aE bE๏ธ ๏จ ๏ธ ๏จ๏€ซ ๏€ฝ ๏€ซ We know that a function from a probability space (ฮฉ, ๐ด, ๐‘ƒ๐‘Ÿ) to a collection of random variables is called a birandom variable, from the definition of birandom variable. Birandom variables are two type. They are either discrete birandom variable or continuous birandom variable. Expectation theory of birandom variables will be discussed in the current subsection. Definition 3 (Xu et al., 2009). For a birandom variable ๐œ‰, in the probability space (๐›บ, ๐ด, ๐‘ƒ๐‘Ÿ), if ๐œ‰(๐œ”) is a random variable with a continuous distribution function when ๐œ” โˆˆ ๐›บ and its expected value ๐ธ[๐œ‰(๐œ”)] is a birandom variable. Then we call ๐œ‰ continuous birandom variable. Definition 4 (Xu et al., 2009). Suppose ๐œ‰ is a birandom variable, then ๐œ‰(๐œ”) is a random variable. If ๐‘“(๐‘ฅ, ๐œ‰) be the density function of ๐œ‰(๐œ”), and ๏€จ ๏€ฉ ๏€จ ๏€ฉ ฮฉx E xf x dx๏ธ ๏ท ๏ƒŽ ๏ƒฉ ๏ƒน ๏€ฝ๏ƒซ ๏ƒป ๏ƒฒ (2) Then the density function of birandom variable ๐œ‰ is ๐‘“(๐‘ฅ, ๐œ‰). Definition 5 (Xu et al., 2009). For the continuous birandom variable ๐œ‰, if its density function is ๐‘“(๐‘ฅ), we can define the expected value of ๐œ‰ as follows ๏› ๏ ๏€จ ๏€ฉ ๏€จ ๏€ฉ 0 0 ฮฉ ฮฉx x E Pr xf x r dr Pr xf x r dr๏ธ ๏‚ฅ ๏ƒŽ ๏€ญ๏‚ฅ ๏ƒŽ ๏ƒฌ ๏ƒผ ๏ƒฌ ๏ƒผ๏ƒฏ ๏ƒฏ ๏ƒฏ ๏ƒฏ ๏€ฝ ๏‚ณ ๏€ญ ๏‚ฃ๏ƒญ ๏ƒฝ ๏ƒญ ๏ƒฝ ๏ƒฏ ๏ƒฏ ๏ƒฏ ๏ƒฏ๏ƒฎ ๏ƒพ ๏ƒฎ ๏ƒพ ๏ƒฒ ๏ƒฒ ๏ƒฒ ๏ƒฒ (3) Definition 6 (Xu et al., 2009). If the density function of a birandom variable ๐œ‰ is ๐‘“(๐‘ฅ, ๐œ‰) and ๐‘”(๐‘ฅ) is a continuous function. Then expectation for the birandom variable ๐‘”(๐œ‰) is defined as ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ 0 0 ฮฉ ฮฉx x E g Pr g x f x r dr Pr g x f x r dr๏ธ ๏‚ฅ ๏ƒŽ ๏€ญ๏‚ฅ ๏ƒŽ ๏ƒฌ ๏ƒผ ๏ƒฌ ๏ƒผ๏ƒฏ ๏ƒฏ ๏ƒฏ ๏ƒฏ ๏ƒฉ ๏ƒน ๏€ฝ ๏‚ณ ๏€ญ ๏‚ฃ๏ƒญ ๏ƒฝ ๏ƒญ ๏ƒฝ๏ƒซ ๏ƒป ๏ƒฏ ๏ƒฏ ๏ƒฏ ๏ƒฏ๏ƒฎ ๏ƒพ ๏ƒฎ ๏ƒพ ๏ƒฒ ๏ƒฒ ๏ƒฒ ๏ƒฒ (4) Theorem 1 (Xu et al., 2009). Let ๐‘“(๐‘ฅ, ๐œ‰) is the density function of the birandom variable ๐œ‰. Then the expected value of ๐œ‰ exists if only if the expected value of random variable ๐œ‰(๐œ”) exists. Theorem 2 (Xu et al., 2009). The expectation of a birandom variable. ๐œ‰~๐‘(๐œ‡, ๐œŽ 2), where, ๐œ‡~U(a, b) is ๐‘Ž+๐‘ 2 . http://www.dmame.org/index.php/dmame/article/view/24 Adhikary et al./Decis. Mak. Appl. Manag. Eng. 2 (1) (2019) 1-12 6 Theorem 3. The expectation of a birandom variable ๐œ‰~๐‘(๐œ‡, ๐œŽ 2) where ๐œ‡~exp (๐œ†) is 1 ๐œ† . Proof: By definition (2), we know ๏› ๏ ๏€จ ๏€ฉ๏ป ๏ฝ ๏€จ ๏€ฉ๏ป ๏ฝ 0 0 Pr | Pr | E E t dt E t dt๏ธ ๏ท ๏ธ ๏ท ๏ท ๏ธ ๏ท ๏‚ฅ ๏€ญ๏‚ฅ ๏€ฝ ๏ƒŽ๏— ๏ƒฉ ๏ƒน ๏‚ณ ๏€ญ ๏ƒŽ๏— ๏ƒฉ ๏ƒน ๏‚ฃ๏ƒซ ๏ƒป ๏ƒซ ๏ƒป๏ƒฒ ๏ƒฒ Since ๐œ‡~ exp(๐œ†), and obviously ๐ธ[๐œ‰(๐œ”)] = ๐œ‡, by definition 4 and 5, the above function can be transformed as follows, ๏› ๏ ๏ป ๏ฝ ๏ป ๏ฝ 0 0 Pr PrE t dt t dt๏ธ ๏ญ ๏ญ ๏‚ฅ ๏€ญ๏‚ฅ ๏€ฝ ๏‚ณ ๏€ญ ๏‚ฃ๏ƒฒ ๏ƒฒ Since, ๐œ‡~ exp(๐œ†), and we know that the density function and the distribution function of exponential distribution are as follows, ๏€จ ๏€ฉ ๏› ๏€ฉ, 0, xf x e x and๏ฌ๏ฌ ๏€ญ๏€ฝ ๏ƒŽ ๏‚ฅ ๏€จ ๏€ฉ ๏› ๏€ฉ1 , 0, .xF x e x๏ฌ๏€ญ๏€ฝ ๏€ญ ๏ƒŽ ๏‚ฅ According to the definition of the distribution function we can obtain the following two functions from the distribution function Pr(๐œ‡ โ‰ค ๐‘ฅ) = 1 โˆ’ ๐‘’ โˆ’๐œ†๐‘ฅ , 0 โ‰ค ๐‘ฅ < โˆž ๐‘Ž๐‘›๐‘‘ Pr(๐œ‡ โ‰ฅ ๐‘ฅ) = ๐‘’ โˆ’๐œ†๐‘ฅ , 0 โ‰ค ๐‘ฅ < โˆž Obviously, โˆซ Pr{๐œ‡ โ‰ค ๐‘ก} ๐‘‘๐‘ก = 0 0 โˆ’โˆž Therefore ๐ธ[๐œ‰] = โˆซ Pr{๐œ‡ โ‰ฅ ๐‘ก} ๐‘‘๐‘ก = โˆซ ๐‘’ โˆ’๐œ†๐‘ก ๐‘‘๐‘ก = [ ๐‘’ โˆ’๐œ†๐‘ก ๐œ† ] 0 โˆž = 1 ๐œ† โˆž 0 โˆž 0 . However, it is very hard to accomplish the mathematical expression of expected value for all types of birandom variables. But using birandom simulation, we could calculate the expected value of birandom variables, with the help of Strong Number Law. 3. Birandom Simulation Let (๐›บ, ๐ด, ๐‘ƒ๐‘Ÿ), be a probability space and a ๐‘“: โ„œ๐‘› โ†’ โ„œ be a measurable function. Consider that ๐œ‰ is an ๐‘› โˆ’ dimension birandom vector on the given probability space. Now we have to find the expectation ๐ธ[๐‘“(๐œ‰)] of birandom variable. Using stochastic simulation, we can find the expected value for every ฯ‰ โˆˆ โ„ฆ. Here, we have used an algorithm for birandom simulation to find the expectation of ๐ธ[๐‘“(๐œ‰(๐œ”))], which is defined as follows. Algorithm: Step 1. Start Step 2. Set l=0 and N= Number of iteration Step 3. Sample ๐œ” from ๐›บ according to the probability measure ๐‘ƒ๐‘Ÿ Newsboy problem with birandom demand 7 Step 4. ๐ธ[๐‘“(๐œ‰(๐œ”))] is find by the stochastic simulation. Step 5. Then ๐‘™ โ† ๐‘™ + ๐ธ[๐‘“(๐œ‰(๐œ”))] Step 6. Repeat the steps from second to fifth steps ๐‘ ttimes. Step 7. ๐ธ[๐‘“(๐œ‰(๐œ”))] = ๐‘™ ๐‘› . Step 8. Stop 4. Mathematical Formulation We are assuming a single period inventory problem with single product. Here all the costs (buying cost and selling cost) are deterministic. Salvage value is taken which is also deterministic. But the demand is birandom variable. The mathematical notation of a birandom newsboy problem is defined as follows: ๐œ‰ ฬ… : The demand of market, a birandom variable ๐‘ฅ : The quantity which to be order, a decision variable ๐‘ = ๐‘(1 + ๐‘š) : Selling price per unit ๐‘  = ๐‘(1 โˆ’ ๐‘‘) : Salvage value per unit ๐‘ : Purchasing cost per unit ๐‘” (๐‘ฅ, ๐œ‰ ฬ…) : The profit for the order quantity ๐‘ฅ and demand ๐œ‰ ฬ… ๐œ‡ : Expected value of birandom demand ๐œ‰ ฬ… ๐œŽ 2 : Variance of the birandom demand ๐œ‰ฬ… ๐‘š : Mark-up, i.e., return per dollar on unit sold ๐‘‘ : Discount rate, i.e., loss per dollar on unit unsold ๐‘ฅ + = max{๐‘ฅ, 0} : The positive part of x Then the profit can be expressed as ๏€จ ๏€ฉ ๏ป ๏ฝ ๏€จ ๏€ฉ, min ,g x p x s x cx๏ธ ๏ธ ๏ธ ๏€ซ ๏ƒฉ ๏ƒน๏€ฝ ๏€ซ ๏€ญ ๏€ญ ๏ƒช ๏ƒบ๏ƒซ ๏ƒป (5) Now min (๐‘ฅ, ๐œ‰ ฬ… ) = ๐œ‰ฬ… โˆ’ (๐œ‰ฬ… โˆ’ x) + Where ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉx x x๏ธ ๏ธ ๏ธ ๏€ซ ๏€ซ ๏€ญ ๏€ฝ ๏€ญ ๏€ซ ๏€ญ (6) ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ๏€จ ๏€ฉ, g x p s x cx p s x๏ธ ๏ธ ๏ธ ๏ธ ๏€ซ ๏€ฝ ๏€ซ ๏€ญ ๏€ญ ๏€ญ ๏€ญ ๏€ญ (7) Since the demand for the product is birandom, the profit function ๐‘” (๐‘ฅ, ๐œ‰ฬ… ) is also consists of birandom variable. Hence, the expectation criteria is used for handling the birandom variable. Therefore, to find the optimal quantity, the decision maker will maximize the total expected value. We can write the expected profit as ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉฮ  x p s s c x p s E x๏ญ ๏ธ ๏€ซ ๏€ฝ ๏€ญ ๏€ซ ๏€ญ ๏€ญ ๏€ญ ๏€ญ (8) http://www.dmame.org/index.php/dmame/article/view/24 Adhikary et al./Decis. Mak. Appl. Manag. Eng. 2 (1) (2019) 1-12 8 Or using the definition of m and d, as ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉฮ  [x c m d xd m d E x๏ญ ๏ธ ๏€ซ ๏€ฝ ๏€ซ ๏€ญ ๏€ญ ๏€ซ ๏€ญ (9) The information of ๐œ‰ ฬ…is known. To maximize the profit function, we need the following lemma. Lemma 4. For given a birandom variable ๐œ‰,ฬ… we have the following inequality, ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ 1 22 2 2 x x E x ๏ณ ๏ญ ๏ญ ๏ธ ๏€ซ ๏ƒฉ ๏ƒน๏€ซ ๏€ญ ๏€ญ ๏€ญ ๏ƒซ ๏ƒป ๏€ญ ๏‚ฃ (10) Proof: Notice that (๐œ‰ฬ… โˆ’ ๐‘ฅ) + = |๏ฟฝฬ…ฬƒ๏ฟฝโˆ’๐‘ฅ|+(๏ฟฝฬ…ฬƒ๏ฟฝโˆ’๐‘ฅ) 2 The result follows by taking expectations and by using the Cauchy-Schwarz inequality ๐ธ [๐œ‰ฬ… โˆ’ ๐‘ฅ] โ‰ค [๐ธ |๐œ‰ฬ… โˆ’ ๐‘ฅ| 2 ] 1 2 = [๐œŽ 2 + (๐‘ฅ โˆ’ ๐œ‡)2] 1 2 By using the lemma (4) the equation (9) will be rewritten as ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ ๏€จ ๏€ฉ 1 22 2 ฮ  2 x x x c m d xd m d ๏ณ ๏ญ ๏ญ ๏ญ ๏ƒฉ ๏ƒน ๏ƒฉ ๏ƒน๏€ซ ๏€ญ ๏€ญ ๏€ญ๏ƒช ๏ƒบ ๏ƒซ ๏ƒป ๏€ฝ ๏€ซ ๏€ญ ๏€ญ ๏€ซ๏ƒช ๏ƒบ ๏ƒช ๏ƒบ ๏ƒช ๏ƒบ๏ƒซ ๏ƒป (11) It is easy to validate that equation no (11) is strictly convex in ๐‘ฅ. Upon setting the derivative to zero and solving for ๐‘ฅ we obtain the ordering rule 1 1 2 2 * 2 m d x d m ๏ณ ๏ญ ๏ƒฆ ๏ƒถ ๏ƒฉ ๏ƒน ๏ƒฉ ๏ƒน๏ƒง ๏ƒท๏€ฝ ๏€ซ ๏€ญ๏ƒช ๏ƒบ ๏ƒช ๏ƒบ๏ƒง ๏ƒท๏ƒซ ๏ƒป ๏ƒซ ๏ƒป ๏ƒจ ๏ƒธ (12) 5. Numerical Example Assume that the unit purchase price of a perishable product is ๐‘ = $40, the unit selling price is ๐‘ = $60, and there is no salvage value (๐‘  = 0). Thus, ๐‘š = ๐‘ ๐‘ โˆ’ 1 = 60 40 โˆ’ 1 = 1 2 . Discount rate ๐‘‘ = 1 โˆ’ ๐‘  ๐‘ = 1. Further assume that the product demand is a birandom variable with normal distribution ๐‘(ยต1, 400), and ยต1 โˆผ ๐‘’๐‘ฅ๐‘(0.0003). From theorem (3), we have ๐œ‰ โˆผ ๐‘(ยต1, 400), where, ๐œ† = 0.0003. Hence, by theorem (3) we can say that the mean (ยต) of the birandom variable (๐œ‰) is 1 ๐œ† = 1 0.0003 = 3333.33, and ๐œŽ 2 = 400. Now it remains to calculate the optimal order quantity and expected profit. For this purpose, we apply equation (12) and obtained the optimal order quantity, ๐‘ฅ โˆ— = 3326. Hence, the expected profit is, ๐›ฑ = $66100. Newsboy problem with birandom demand 9 6. A Case Study To endorse the model developed in this study, we sent our projected framework to five leading fish merchants in West Bengal, India. They sell only the freshest and best quality fishes, and maintains quality control at every stage of packaging and delivery in many different parts or areas of West Bengal. Among them two firms positively responded to explore this research proposal and we conducted necessary preliminary tasks on these companies. We selected a reputed fish merchant (โ€œLakshmi fish enterpriseโ€, name changed), situated in the โ€southernโ€ West Bengal, which has several operational units nationwide. Our objective is to incorporate the perceptions of all participants (customer/retailer/company mangers) in the fish industry and to achieve its comprehensive outcomes since this research is purely grounded on birandom product demand information obtained from experts in the business. In this paper, we consider the perspectives of a wholesale fish merchant who buys fishes from the company, having a large market share in Kolkata zone. In West Bengal, fish merchants generally sell a special fish in monsoon. The name of this fish is Hilsa. The business of this fish is a good example of SPIM. The business of this fish totally depends on its demand and supply in the monsoon season. Merchants have to decide how many fishes should be purchased from his or her supplier depending on the customerโ€™s demand. Buying more amount of Hilsa may not bring him more profit. Rather it can cause him a great loss since it cannot be preserved for long periods and the expired fish has no market value. If he/she buys too few amount of Hilsa he/she will lose the opportunity of making a higher profit. Thus, the actual inventory level cannot be determined precisely in such complex situation. It may be assumed for simplicity that the fish demand follows normal distribution. But under such circumferences, the manager looks after of some previous data, and finds the mean demand of โ€Hilsaโ€ is also a random variable. This leads us to consider the โ€Hilsaโ€ demand as a birandom variable. Each such fish sells for $60 and costs for the shop owner $40. Investigating the previous year data, the decision maker decides the demand follows the two types two- folded random variable (birandom). Scenario 1: Normal distribution ๐‘~(๐œ‡1, 400), with ๐œ‡1~๐‘ˆ(3000, 4000). Scenario 2: Normal distribution ๐‘~(๐œ‡1, 400), with ๐œ‡1~๐‘(3500, 500). Therefore, according to our proposed model we have ๐‘ = 40, ๐‘ = 60, and there is no salvage value i.e., ๐‘  = 0. In scenario 1, by theorem (2) the mean of the birandom variable is ยต = 3500, and ๐œŽ 2 = 400. Therefore the optimal order quantity ๐‘ฅ โˆ— = 3492. Expected profit ๐›ฑ = $69434. And in scenario 2, using birandom simulation we get the mean of the birandom simulation ยต = 3508, and ๐œŽ 2 = 400. Therefore, the optimal order quantity ๐‘ฅ = 3500, expected profit = $69594. Finally, we shared the outcomes of this research work with the managers of our case enterprise, they are satisfied with the outcomes and willing to accept this result for their monsoon business of โ€Hilsaโ€. 7. Conclusion In this paper we have proposed a newsboy problem where the demand is considered as birandom variable. The market volatility and uncertainty in customersโ€™ behavior make the demand of the product a birandom variable. We use the expected value model for handling this birandom variable and to convert the BNM into its equivalent deterministic model. We discuss a case study of fish merchant to validate the usefulness and applicability of the proposed model. In this proposed model http://www.dmame.org/index.php/dmame/article/view/24 Adhikary et al./Decis. Mak. Appl. Manag. Eng. 2 (1) (2019) 1-12 10 demand is the only parameter considered as birandom variable. However, in reality cost may be assume as a birandom variable. This is the one aspects to be considered for further investigation in the current model. Also, for future research work, one can develop chance constraint technique to convert the birandom model into deterministic one. References Hadley, G. W. (1963). Analysis of inventory systems, Prentice-Hall, New Jersey, USA. Tadikamalla, P. R. (1978). Applications of the Weibull distribution in inventory control. Journal of the Operational Research Society, 29(1), 77-83. Mahoney, J. F., & Sivazlian, B. D. (1980). Probability of Shortage for Erlang Distributed Demand in a (๐œŽ, ๐‘†) Inventory Problem. SIAM Journal on Applied Mathematics, 38(1), 156-162. Gallego, G. (1992). A minmax distribution free procedure for the (๐‘„, ๐‘…) inventory model. Operations Research Letters, 11(1), 55-60. Gallego, G. & Moon, I. (1993). The distribution free newsboy problem: review and extensions. Journal of Operational Research Society, 44(8), 825-834. Abdel-Aal, M. A., Syed, M. N., & Selim, S. Z. (2017). Multi-product selective newsvendor problem with service level constraints and market selection flexibility. International Journal of Production Research, 55(1), 96-117. Watt, R., & Vzquez, F. J. (2017). An analysis of insurance demand in the newsboy problem. European Journal of Operational Research, 259(3), 1064-1072. Sun, Y., & Guo, Z. (2017). Optimal Inventory Policy for Newsboy Problem Based on Fuzzy Random Demand. In Proceedings of the Fourth International Forum on Decision Sciences. Uncertainty and Operations Research. Springer, Singapore. Vipin, B., & Amit, R. K. (2017). Loss aversion and rationality in the newsvendor problem under recourse option. European Journal of Operational Research, 261(2), 563-571. Natarajan, K., Sim, M., & Uichanco, J. (2017). Asymmetry and ambiguity in newsvendor models. Management Science, 64(7), 2973-3468. Pando, V., San-Jos, L. A., Garca-Laguna, J., & Sicilia, J. (2013). A newsboy problem with an emergency order under a general backorder rate function. Omega, 41(6), 1020- 1028. Zhang, W., Li, J., Zhang, R., & Chen, Y. (2017). Impact of Emergency Order in Price- Dependent Newsvendor Problems. Asia-Pacific Journal of Operational Research, 34(02), 1750001. Kamburowski, J. (2015). On the distribution-free newsboy problem with some non- skewed demands. Operations Research Letters, 43(2), 165-171. Kamburowski, J. (2014). The distribution-free newsboy problem under the worst-case and best-case scenarios. European Journal of Operational Research, 237(1), 106-112. Newsboy problem with birandom demand 11 Gler, M. G. (2014). A note on: The effect of optimal advertising on the distribution-free newsboy problem. International Journal of Production Economics, 148, 90-92. Lee, C. M., & Hsu, S. L. (2011). The effect of advertising on the distribution-free newsboy problem. International Journal of Production Economics, 129(1), 217-224. Ding, S. (2013). Uncertain multi-product newsboy problem with chance constraint. Applied mathematics and computation, 223(15), 139-146. Abdel-Malek, L. L., & Otegbeye, M. (2013). Separable programming/duality approach to solving the multi-product Newsboy/Gardener Problem with linear constraints. Applied Mathematical Modelling, 37(6), 4497-4508. Chen, S. P., & Ho, Y. H. (2013). Optimal inventory policy for the fuzzy newsboy problem with quantity discounts. Information Sciences, 228, 75-89. Peng, J., & Liu, B. (2007). Birandom variables and birandom programming. Computers & Industrial Engineering, 53(3), 433-453. Liu, Z. (2013). The Common Formula on Birandom Variable Variance. International Conference on Information Technology and Applications (ITA), IEEE, 486-488. Xu, J., & Zhou, X. (2009). A class of multi-objective expected value decision-making model with birandom coefficients and its application to flow shop scheduling problem. Information Sciences, 179(17), 2997-3017. Yan, L. (2009). One type of optimal portfolio selection in birandom environments. Modern Applied Science, 3(6), 126. Xu, J., & Ding, C. (2011). A class of chance constrained multiobjective linear programming with birandom coefficients and its application to vendors selection. International Journal of Production Economics, 131(2), 709-720. Wang, G., Tang, W., & Zhao, R. (2012). A birandom job search problem with risk tolerance. Journal of Applied Mathematics, Article ID 161573, 1-12. Zhang, Z., & Xu, J. (2013). A multi-mode resource-constrained project scheduling model with birandom coefficients for drilling grouting construction project. International Journal of Civil Engineering, 11(1), 1-13. Xu, J., Tu, Y., & Zeng, Z. (2012). A nonlinear multiobjective bilevel model for minimum cost network flow problem in a large-scale construction project. Mathematical Problems in Engineering. Article ID 463976, 1-40. Tavana, M., Shiraz, R. K., & Hatami-Marbini, A. (2014). A new chance-constrained DEA model with birandom input and output data. Journal of the operational research society, 65(12), 1824-1839. Tao, Z., & Xu, J. (2013). An equilibrium chance-constrained multiobjective programming model with birandom parameters and its application to inventory problem. Journal of Applied Mathematics. Article ID 382620, 1-12. Xu, J., Zhang, Z., & Mookerjee, V. S. (2013). Applying birandom MODM model to navigation coordinated scheduling: a case study of Three Gorges Project. Transport, 28(2), 140-157. http://www.dmame.org/index.php/dmame/article/view/24 Adhikary et al./Decis. Mak. Appl. Manag. Eng. 2 (1) (2019) 1-12 12 Cao, G., & Shan, D. (2013). The effect of exit strategy on optimal portfolio selection with birandom returns. Journal of Applied Mathematics. Article ID 236579, 1-6. Maity, S., Roy, A., & Maiti, M. (2015). A modified genetic algorithm for solving uncertain constrained solid travelling salesman problems. Computers & Industrial Engineering, 83, 273-296. Zhou, M., Lu, S., Tan, S., Yan, D., Ou, G., Liu, D.,& Zhu, X. (2017). A stochastic equilibrium chance-constrained programming model for municipal solid waste management of the City of Dalian, China. Quality & Quantity, 51(1), 199-218 Wang, Y., Fan, J., Xu, Y., Sun, W., & Chen, D. (2017). Designing management strategies for carbon dioxide storage and utilization under uncertainty using inexact modelling. Engineering Optimization, 49(6), 1024-1039. Agrawal, V., & Seshadri, S. (2000). Impact of uncertainty and risk aversion on price and order quantity in the newsvendor problem. Manufacturing & Service Operations Management, 2(4), 410-423. ยฉ 2018 by the authors. Submitted for possible open access publication under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).