Newsvendor Model Chapter 11
1
Learning Goals
Determine the optimal level of product availability – Demand forecasting – Profit maximization
Service measures such as a fill rate
Learning Goals
Determine the optimal level of product availability – Demand forecasting – Profit maximization
Service measures such as a fill rate
Motivation
Determining optimal levels (purchase orders) – Single order (purchase) in a season – Short lifecycle items
1 month: Printed Calendars, Rediform
6 months: Seasonal Camera, Panasonic
18 months, Cell phone, Nokia
Motivating Newspaper Article for toy manufacturer Mattel
Mattel [who introduced Barbie in 1959 and run a stock out for several years then on] was hurt last year by inventory cutbacks at Toys “R” Us, and officials are also eager to avoid a repeat of the 1998 199 8 Thanksgiving weekend. Mattel had expected to ship a lot of merchandise after the weekend, but retailers, wary of excess inventory, stopped ordering from Mattel. That led the company to report a $500 million sales shortfall in the last weeks of the year y ear ... For the crucial holiday selling season this year, Mattel said it will require retailers to place their full orders before Thanksgiving. And, for the first time, the company will no longer take reorders in December, December, Ms. Barad said. This will enable Mattel Mattel to tailor production more closely to demand and avoid building inventory for orders that - Wall Street Journal, Journal, Feb. 18, 18, 1999 don't come.
For tax (in accounting), accounting), option pricing (in finance) and revenue revenue management management
O’Neill’s Hammer 3/2 wetsuit
Hammer 3/2 timeline and economics Generate forecast of demand and submit an order to TEC, supplier
Spring selling season
Nov Dec Jan
Feb
Mar Apr
Receive order from TEC at the end of the month
Economics:
May
Jun
Jul
Aug
Leftover units are discounted
Each suit sells for p = $180 TEC charges c = $110/suit Discounted suits sell for v = $90
The “too much/too little problem”: – Order too much and inventory is left over at the end of the season – Order too little and sales are lost.
Newsvendor model implementation steps 1.
Gather economic inputs: a) selling price, b) production/procurement cost, c) salvage value of inventory
2.
Generate a demand model to represent demand a)
Use empirical demand distribution
b)
Choose a standard distribution function: the normal distribution and the Poisson distribution – for discrete items
3.
4.
Choose an aim: a)
maximize the objective of expected profit
b)
satisfy a fill rate constraint.
Choose a quantity to order.
The Newsvendor Model: Develop a Forecast
7
Historical forecast performance at O’Neill 7000 6000 .
5000 d n
4000 a m e d l
3000 a ut c A
2000 1000 0 0
1000
2000
3000
4000
5000
6000
7000
Forecast
Forecasts and actual demand for surf wet-suits from the previous season
How do we know the actual when the actual demand > forecast demand Are the number of stockout units (= unmet demand=demand-stock) observable, i.e., known to the store manager?
Yes, if the store manager issues rain checks to customers. No, if the stockout demand disappears silently.
– A vicious cycle Underestimate the demand Stock less than necessary. Stocking less than the demand Stockouts and lower sales. Lower sales Underestimate the demand.
– Demand filtering: Demand known exactly only when below the stock. – Shall we order more than optimal to learn about demand?
Yes and no, if some customers complain about a stockout; see next page.
Observing a Portion of Unmet Demand
Unmet demand are reported by partners (sales associates)
Reported lost sales are based on customer complaints Not everybody complains of a stock out, Not every sales associate records complaints, Not every complaint is reported, Only a portion of complaints are observed by IM
??
Empirical distribution of forecast accuracy Order by A/F ratio Product description Forecast JR ZEN FL 3/2 90 EPIC 5/3 W/HD 120 JR ZEN 3/2 140 WMS ZEN-ZIP 4/3 170 HEATWAVE 3/2 170 JR EPIC 3/2 180 WMS ZEN 3/2 180 ZEN-ZIP 5/4/3 W/HOOD 270 WMS EPIC 5/3 W/HD 320 EVO 3/2 380 JR EPIC 4/3 380 WMS EPIC 2MM FULL 390 HEATWAVE 4/3 430 ZEN 4/3 430 EVO 4/3 440 ZEN FL 3/2 450 HEAT 4/3 460 ZEN-ZIP 2MM FULL 470 HEAT 3/2 500 WMS EPIC 3/2 610 WMS ELITE 3/2 650 ZEN-ZIP 3/2 660 ZEN 2MM S/S FULL 680 EPIC 2MM S/S FULL 740 EPIC 4/3 1020 WMS EPIC 4/3 1060 JR HAMMER 3/2 1220 HAMMER 3/2 1300 HAMMER S/S FULL 1490 EPIC 3/2 2190 ZEN 3/2 3190 ZEN-ZIP 4/3 3810 WMS HAMMER 3/2 FULL 6490 * Error = Forecast - Actual demand
Actual demand 140 83 143 163 212 175 195 317 369 587 571 311 274 239 623 365 450 116 635 830 364 788 453 607 732 1552 721 1696 1832 3504 1195 3289 3673
Error* A/F Ratio** -50 1.56 37 0.69 -3 1.02 7 0.96 -42 1.25 5 0.97 -15 1.08 -47 1.17 -49 1.15 -207 1.54 -191 1.50 79 0.80 156 0.64 191 0.56 -183 1.42 85 0.81 10 0.98 354 0.25 -135 1.27 -220 1.36 286 0.56 -128 1.19 227 0.67 133 0.82 288 0.72 -492 1.46 499 0.59 -396 1.30 -342 1.23 -1314 1.60 1995 0.37 521 0.86 2817 0.57
33 products, so increment probability by 3%. 100% 90% 80% 70% t i l i b a b o r P
60% 50% 40% 30% 20% 10% 0% 0.00
0.25
0.50
0.75
1.00
1.25
1.50
A/F ratio Empirical distribution function for the historical A/F ratios.
1.75
Normal distribution tutorial
All normal distributions are specified by 2 parameters, mean = m and st_dev = s.
Each normal distribution is related to the standard normal that has mean = 0 and st_dev = 1.
For example: – Let Q be the order quantity, and (m , s ) the parameters of the normal demand forecast. – Prob{demand is Q or lower} = Prob{the outcome of a standard normal is z or lower}, where
z
Q m
or
Q m z s
s
– The above are two ways to write the same equation, the first allows you to calculate z from Q and the second lets you calculate Q from z . – Look up Prob{the outcome of a standard normal is z or lower} in the
Standard Normal Distribution Function Table.
Using historical A/F ratios to choose a Normal distribution for the demand forecast
Start with an initial forecast generated from hunches, guesses, etc. – O’Neill’s initial forecast for the Hammer 3/2 = 3200 units.
Evaluate the A/F ratios of the historical data: A/F ratio
Actual demand
Forecast
Set the mean of the normal distribution to Expected actual demand Expected A/F ratio Forecast
Set the standard deviation of the normal distribution to Standard deviation of actual demand Standard deviation of A/F ratios Forecast
O’Neill’s Hammer 3/2 normal distribution forecast Product description JR ZEN FL 3/2 EPIC 5/3 W/HD JR ZEN 3/2 WMS ZEN-ZIP 4/3 …
Forecast Actual demand 90 140 120 83 140 143 170 156 …
ZEN 3/2 ZEN-ZIP 4/3 WMS HAMMER 3/2 FULL Average Standard deviation
…
3190 3810 6490
1195 3289 3673
Error -50 37 -3 14 …
1995 521 2817
A/F Ratio 1.5556 0.6917 1.0214 0.9176 …
0.3746 0.8633 0.5659 0.9975 0.3690
Expected actual demand 0.9975 3200 3192 Standard deviation of actual demand
0.369
3200 1181
Choose a normal distribution with mean 3192 and st_dev 1181 to represent demand for the Hammer 3/2 during the Spring season.
Why not a mean of 3200?
Fitting Demand Distributions: Empirical vs normal demand distribution 1.00 0.90 0.80 0.70 0.60 y ti l i
0.50 b a b
0.40 o r P
0.30 0.20 0.10 0.00 0
1000
2000
3000
4000
5000
6000
Quantity
Empirical distribution function (diamonds) and normal distribution function with
An Example of Empirical Demand: Demand for Candy (in the Office Candy Jar)
An OPRE 6302 instructor believes that passing out candies (candies, chocolate, cookies) in a late evening class builds morale and spirit. This belief is shared by office workers as well. For example, secretaries keep office candy jars, which are irresistible: “… 4-week study involved the chocolate candy consumption of 40 adult secretaries.
The study utilized a 2x2 within-subject design where candy proximity was crossed with visibility. Proximity was manipulated by placing the chocolates on the desk of the participant or 2 m from the desk. Visibility was manipulated by placing the chocolates in covered bowls that were either clear or opaque. Chocolates were replenished each evening. “ People ate an average of 2.2 more candies each day when they were visible, and 1.8 candies more when they were proximately placed on their desk vs 2 m away.” They ate 3.1 candies/day when candies were in an opaque container.
Candy demand is fueled by the proximity and visibility. What fuels the candy demand in the OPRE 6302 class?
What undercuts the demand?
Hint: The aforementioned study is titled “The office candy dish: proximity's influence on estimated and actual consumption” and published in International Journal of Obesity (2006) 30: 871 – 875.
The Newsvendor Model: The order quantity that maximizes expected profit
17
“Too much” and “too little” costs
C o = overage cost – The cost of ordering one more unit than what you would have ordered had you known demand. – In other words, suppose you had left over inventory (i.e., you over ordered). C o is the increase in profit you would have enjoyed had you ordered one fewer unit. – For the Hammer C o = Cost – Salvage value = c – v = 110 – 90 = 20
C u = underage cost – The cost of ordering one fewer unit than what you would have ordered had you known demand. – In other words, suppose you had lost sales (i.e., you under ordered). C u is the increase in profit you would have enjoyed had you ordered one more unit. – For the Hammer C u = Price – Cost = p – c = 180 – 110 = 70
Balancing the risk and benefit of ordering a unit
Ordering one more unit increases the chance of overage – Probability of overage F (Q) =Prob{ Demand ≤ Q) – Expected loss on the Qth unit = C o x F (Q) = “Marginal cost of overstocking”
The benefit of ordering one more unit is the reduction in the chance of underage – Probability of underage 1- F (Q) – Expected benefit on the Qth unit = C u x (1- F (Q)) = “Marginal benefit of understocking” 80
.
70
s s o l r o n i a g d e t c e p x E
60
Expected marginal benefit of an extra unit in reducing understocking
As more units are ordered, the expected marginal benefit from ordering 1 more unit decreases while the expected marginal cost of ordering 1 more unit increases.
50 40 Expected marginal overstocking cost of an extra unit
30 20 10
0 0
800
1600
2400
3200
4000
4800
5600
6400
Expected profit maximizing order quantity
To minimize the expected total cost of underage and overage, order Q units so that the expected marginal cost with the Qth unit equals the expected marginal benefit with the Qth unit: C o F (Q) C u 1 F Q
Rearrange terms in the above equation F (Q)
C u C o
C u
The ratio C u / (C o + C u) is called the critical ratio.
Hence, to minimize the expected total cost of underage and overage, choose Q such that we do not have lost sales (i.e., demand is Q or lower) with a probability that equals to the critical ratio
Expected cost minimizing order quantity with the empirical distribution function
Inputs: Empirical distribution function table; p = 180; c = 110; v = 90; C u = 180-110 = 70; C o = 110-90 =20
Evaluate the critical ratio:
C u C o
C u
70
20 70
0.7778
Look up 0.7778 in the empirical distribution function graph
Or, look up 0.7778 among the ratios: – If the critical ratio falls between two values in the table, choose the one that leads to the greater order quantity Product description …
HEATWAVE 3/2 HEAT 3/2 HAMMER 3/2 …
Forecast Actual demand A/F Ratio Rank …
170 500 1300
…
– Convert A/F ratio into the order quantity
…
212 635 1696
…
…
1.25 1.27 1.30
…
…
24 25 26 …
Percentile
…
72.7% 75.8% 78.8%
…
Expected cost minimizing order quantity with the normal distribution
Inputs: p = 180; c = 110; v = 90; C u = 180-110 = 70; C o = 110-90 =20; critical ratio = 0.7778; mean = m = 3192; standard deviation = s = 1181
Look up critical ratio in the Standard Normal Distribution Function Table: z
0
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
0.09
0.5 0.6915 0.6950 0.6985 0.7019 0.7054 0.7088
0.7123
0.7157 0.7190 0.7224
0.6 0.7257 0.7291 0.7324 0.7357 0.7389 0.7422
0.7454
0.7486 0.7517 0.7549
0.7 0.7580 0.7611 0.7642 0.7673 0.7704 0.7734
0.7764
0.7794 0.7823 0.7852
0.8 0.7881 0.7910 0.7939 0.7967 0.7995 0.8023
0.8051
0.8078 0.8106 0.8133
0.9 0.8159 0.8186 0.8212 0.8238 0.8264 0.8289
0.8315
0.8340 0.8365 0.8389
– If the critical ratio falls between two values in the table, choose the greater z -statistic – Choose z = 0.77
Convert the z-statistic into an order quantity: Q m z s
3192 0.77 1181 4101
Or, Q = norminv(0.778,3192,1181) =3192+1181norminv(0.778,0,1) = 4096
Another Example: Apparel Industry How many L.L. Bean Parkas to order? Demand data / distribution Demand D i
Probability of demand being this size
Cumulative Probability of demand being this size or less, F (.)
Cost/Profit data Probability of demand greater than this size, 1-F (.)
Cost per parka = c = $45 Sale price per parka = p = $100 Discount price per parka = $50
4
.01
.01
.99
Holding and transportation cost = $10
5
.02
.03
.97
Salvage value per parka = v = 50-10=$40
6
.04
.07
.93
7
.08
.15
.85
8
.09
.24
.76
Profit from selling parka = p-c = 100-45 = $55
9
.11
.35
.65
Cost of overstocking = c-v = 45-40 = $5
10
.16
.51
.49
11
.20
.71
.29
12
.11
.82
.18
13
.10
.92
.08
14
.04
.96
.04
15
.02
.98
.02
16
.01
.99
.01
Cost of understocking=$55
17
.01
1.00
.00
Cost of overstocking=$5
Expected demand is 1,026 parkas, order 1026 parkas regardless of costs?
Had the costs and demand been symmetric, we would have ordered the average demand.
Costs are almost always antisymmetric.
Optimal Order Q* p = sale price; v = outlet or salvage price; c = purchase price F (Q )=CSL
= In-stock probability = Cycle Service Level = Probability that demand will be at or below reorder point
Raising the order size if the order size is already optimal
Expected Marginal Benefit = =P(Demand is above stock)*(Profit from sales)=(1-CSL)(p - c) Expected Marginal Cost = =P(Demand is below stock)*(Loss from discounting)=CSL(c - v)
Define Co= c-v=overstocking cost; Cu=p-c=understocking cost (1-CSL)Cu = CSL Co CSL= Cu / (Cu + Co)
CSL
*
C u
*
F (Q ) P(Demand Q )
C
C
55 55 5
0.917
Optimal Order Quantity 1.2 0.917
1 0.8 Cumulative Probability
0.6 0.4 0.2 0 4 5 6 7 8 9 10 11 12 13 14 15 16 87 Optimal Order Quantity = 13(‘00)
Marginal Profits at L.L. Bean Approximate additional (marginal) expected profit from ordering 1(‘00) extra parkas if 10(’00) are already ordered =(55.P(D>1000) - 5.P(D≤1000)) 100=(55.(0.49) - 5.(0.51)) 100 =2440 Approximate additional (marginal) expected profit from ordering 1(‘00) extra parkas if 11(’00) are already ordered =(55.P(D>1100) - 5.P(D≤1100)) 100=(55.(0.29) - 5.(0.71)) 100 =1240 Additional Expected Expected Expected Marginal 100s Marginal Benefit Marginal Cost Contribution 1011 2695-255 = 2440 5500.49 = 2695 500 .51 = 255 1112
5500.29 = 1595 500 .71 = 355
1595-355 = 1240
13
5500 .18 = 990 500 .82 = 410
990-410 = 580
1314
5500.08 = 440 500 .92 = 460
440-460 = -20
1415
5500.04 = 220 500 .96 = 480
220-480 = -260
1516
5500.02 = 110 500 .98 = 490
110-490 = -380
1617
5500.01 = 55 500.99 = 495
55-495 = -440
12
Revisit Newsvendor Problem with Calculus
Total cost by ordering Q units: C(Q) = overstocking cost C (Q) C o dC (Q) dQ
Q
0
+
understocking cost
(Q x) f ( x)dx C u
Q
( x Q) f ( x)dx
C o F (Q) C u (1 F (Q)) F (Q)(C o C u ) C u
Marginal cost of raising Q* - Marginal cost of decreasing Q* = 0
*
*
F (Q ) P ( D Q )
C u C o C u
0
Safety Stock Inventory held in addition to the expected demand is called the safety stock The expected demand is 1026 parkas but we order 1300 parkas. So the safety stock is 1300-1026=274 parkas.
The Newsvendor Model: Performance measures
29
Newsvendor model performance measures
For any order quantity we would like to evaluate the following performance measures: – Expected lost sales » The average number of demand units that exceed the order quantity
– Expected sales » The average number of units sold.
– Expected left over inventory » The average number of inventory units that exceed the demand
– Expected profit – Fill rate » The fraction of demand that is satisfied immediately from the stock (no backorder)
– In-stock probability » Probability all demand is satisfied
– Stockout probability » Probability that some demand is lost (unmet)
Expected (lost sales=shortage)
ESC is the expected shortage in a season (cycle)
ESC is not a percentage, it is the number of units ,
also see next page
Demand - Q if Demand Q Shortage 0 if Demand Q ESC E(max{Dema nd in a season Q,0})
ESC
(x - Q)f(x)dx , x Q
f is the probabilit y density of demand.
Inventory and Demand during a season Leftover inventory
0
Q
Leftover Inventory
0 Season
Demand During a Season
Upside down
Q
Leftover Inventory=Q-D
0
D, Demand During A Season
Inventory and Demand during a season Shortage
0
Q
Upside down
0
Q
Shortage =D-Q
Shortage Season Demand
0
n o s a e S a g n i r u D d n a m e D : D
Expected shortage during a season Expected shortage E (max{ D Q,0})
( D Q) P ( D d ) d Q 1
Ex:
d 1 9 with prob p1 1/4 Q 10, D d 2 10 with prob p2 2/4, Expected Shortage? d 11 with prob p 1/4 3 3 Expected shortage
3
11
max{0, (d Q)} p (d 10)} P ( D d ) i
i
i 1
d 11
1
2
1
1
4
4
4
4
max{0, (9 - 10)} max{0, (10 - 10)} max{0, (11 - 10)}
Expected shortage during a season Expected shortage E (max{ D Q,0})
( D Q) f ( D)dD
where f is pdf of Demand.
D Q
Ex:
Q 10, D Uniform (6,12), Expected Shortage?
12
Expected shortage
1
1 D 2
D 12
( D 10) dD 10 D 6 6 2 D 10 D 10
2/6
1 12 2
1 102 10(12) 10(10) 6 2 6 2
Expected lost sales of Hammer 3/2s with Q = 3500 Normal demand with mean 3192, standard deviation=1181 – Step 1: normalize the order quantity to find its z -statistic. z
Q
3500
m
3192
1181
s
0.26
– Step 2: Look up in the Standard Normal Loss Function Table the expected lost sales for a standard normal distribution with that z -statistic: L(0.26)=0.2824 see textbook Appendix B Loss Function Table. » or, in Excel L(z)=normdist(z,0,1,0)-z*(1-normdist(z,0,1,1)) see textbook Appendix D.
– Step 3: Evaluate lost sales for the actual normal distribution:
Expected lost sales
s
L( z ) 1181 0.2824 334
Keep 334 units in mind, we shall repeatedly use it
The Newsvendor Model: Cycle service level and fill rate
37
Type I service measure: Instock probability = CSL Cycle service level Instock probability: percentage of seasons without a stock out For example consider 10 seasons : Instock Probabilit y
11 0 111 0 1 0 1
10 Write 0 if a season has stockout,1 otherwise
Instock Probabilit y 0.7 Instock Probabilit y 0.7 Probabilit y that a single season has sufficient inventory [Sufficien t inventory] [Demand during a season Q]
InstockPro bability P(Demand Q)
Instock Probability with Normal Demand N(μ,σ) denotes a normal demand with mean μ and standard deviation σ P N ( m , s ) Q
Normdist(Q, m , s ,1) P N ( m , s ) m Q m Taking out the mean N ( m , s ) m Q m P Dividing by theStDev s s Q m P N (0,1) Obtaining standard normal distribution s Q m Normdist ,0,1,1 s
Example: Finding Instock probability for given Q μ
= 2,500; = 500;
Q =
3,000;
Instock probability if demand is Normal?
Instock probability = Normdist((3,000-2,500)/500,0,1,1)
Example: Finding Q for given Instock probability μ
= 2,500/week; = 500;
To achieve Instock Probability=0.95, what should Q be?
Q = Norminv(0.95, 2500, 500)
Type II Service measure Fill rate Recall: Expected sales = m - Expected lost sales = 3192 – 334 = 2858
Expected fill rate
Expected sales
Expected demand 1
m
Expected lost sales
2858
m
Expected sales
3192
89.6%
Is this fill rate too low? Well, lost sales of 334 is with Q=3500, which is less than optimal.
Service measures of performance 100% 90% 80% Expected fill rate
70% 60% 50%
In-stock probability CSL
40% 30% 20% 10% 0% 0
1000
2000
3000
4000
Order quantity
5000
6000
7000
Service measures: CSL and fill rate are different inventory
CSL is 0%, fill rate is almost 100%
0
time inventory
0
CSL is 0%, fill rate is almost 0% time
The Newsvendor Model: Measures that follow from lost sales
45
Measures that follow from expected lost sales
Demand=Sales+Lost Sales D=min{D,Q}+max{D-Q,0} or min{D,Q}=D- max{D-Q,0} Expected sales = m - Expected lost sales = 3192 – 334 = 2858
Inventory=Sales+Leftover Inventory Q=min{D,Q}+max{Q-D,0} or max{Q-D,0}=Q-min{D,Q} Expected Leftover Inventory = Q - Expected Sales = 3500 – 2858 = 642
Measures that follow from expected lost sales Economics: • • •
Each suit sells for p = $180 TEC charges c = $110/suit Discounted suits sell for v = $90
Expected total underage and overage cost with (Q=3500) =70*334 + 20*642 Expected profit Price-Cost Expected sales
Cost-Salvage value Expected left over inventory $70 2858 $20 642 $187, 221 What is the relevant objective? Minimize the cost or maximize the profit? Hint: What is profit + cost? It is 70*(3192=334+2858)=C u*μ, which is a constant.
Profit or [Underage+Overage] Cost; Does it matter? (p : price; v : salvage value; c : cost) per unit. D : demand; Q : order quantity.
(p - c)D - (c - v)(Q - D) if [D Q] Overage Profit(D,Q) if [D Q] Underage (p - c)Q (c - v)(Q - D) if [D Q] Overage Cost(D,Q) (p c)(D Q) if [D Q] Underage (p - c)D if D Q Profit(D,Q) Cost(D,Q) (p c)D (p - c)D if D Q E[Profit(D, Q)] E[Cost(D,Q)] (p - c)E(Demand ) Constant in Q M ax E[Profit(D, Q)] and M in E[Cost(D,Q)] are equivalent ; Q
Q
Because they yield the same optimal order quantity.
Computing the Expected Profit with Normal Demands
Expected Profit
Profit(D, Q) f(D) dD
Suppose that the demand is Normal with mean μ and standard deviation σ Expected Profit (p - v) μ normdist(Q, μ, σ,1) - (p - v) σ normdist(Q, μ, σ,0) - (c - v) Q normdist(Q, μ, σ,1) (p c) Q (1 normdist(Q, μ, σ,1)) Example: Follett Higher Education Group (FHEG) won the contract to operate the UTD
bookstore. On average, the bookstore buys textbooks at $100, sells them at $150 and unsold books are salvaged at $50. Suppose that the annual demand for textbooks has mean 8000 and standard deviation 2000. What is the annual expected profit of FHEG from ordering 10000 books? What happens to the profit when standard deviation drops to 20 and order drops to 8000? Expected Profit is $331,706 with order of 10,000 and standard deviation of 2000: = (150-50) *8000*normdist(10000,8000,2000,1)-(150-50)*2000*normdist(10000,8000,2000,0) -(100-50)*10000*normdist(10000,8000,2000,1)+(150-100)*10000*(1-normdist(10000,8000,2000,1))
Expected Profit is $399,960 with order of 8000 and standard deviation of 20: = (150-50)*8000*normdist(8000,8000,20,1)-(150-50)*20* normdist(8000,8000,20,0) -(100-50)*8000*normdist(8000,8000,20,1)+(150-100)*8000*(1-normdist(8000,8000,20,1))