Ask The Wizard #15
From my section on the house edge we find the standard deviation in blackjack to be 1.17. You won't understand this if you haven't studied statistics, but the probability of being at a loss in your example will be the Z statistic of 45000*0.005/(450001/2*1.17) =~ 0.91.
Any basic statistics book should have a standard normal table which will give the Z statistic of 0.8186. So the probability of being ahead in your example is about 18%.
Let ai denote the probability if the player has $i he will reach $1,200 before losing everything. Let p the probability of winning any given bet = 49.3212%.
a0 = 0
a1 = p*a2
a2 = p*a3 + (1-p)*a1
a3 = p*a4 + (1-p)*a2
.
a1197 = p*a1198 + (1-p)*a1196
a1198 = p*a1199 + (1-p)*a1197
a1199 = p*a1200 + (1-p)*a1198
a1200 = 1
Divide the left side into two parts:
p*a1 + (1-p)*a1 = p*a2
p*a2 + (1-p)*a2 = p*a3 + (1-p)*a1
p*a3 + (1-p)*a3 = p*a4 + (1-p)*a2
.
.
.
p*a1197 + (1-p)*a1197 = p*a1198 + (1-p)*a1196
p*a1198 + (1-p)*a1198 = p*a1199 + (1-p)*a1197
p*a1199 + (1-p)*a1199 = p*a1200 + (1-p)*a1198
Rearange with (1-p) terms on the left side and p terms on the right:
(1-p)*(a1) = p*(a2 - a1)
(1-p)*(a2 - a1) = p*(a3 - a2)
(1-p)*(a3 - a2) = p*(a4 - a3)
.
.
.
(1-p)*(a1197 - a1196) = p*(a1198 - a1197)
(1-p)*(a1198 - a1197) = p*(a1199 - a1198)
Next multiply both sides by 1/p:
(1-p)/p*(a1) = (a2 - a1)
(1-p)/p*(a2 - a1) = (a3 - a2)
(1-p)/p*(a3 - a2) = (a4 - a3)
.
.
.
(1-p)/p*(a1197 - a1196) = (a1198 - a1197)
(1-p)/p*(a1198 - a1197) = (a1199 - a1198)
Next telescope sums:
(a2 - a1) = (1-p)/p*(a1)
(a3 - a2) = ((1-p)/p)2*(a1)
(a4 - a3) = ((1-p)/p)3*(a1)
.
.
.
(a1199 - a1198) = ((1-p)/p)1198*(a1)
(a1200 - a1199) = ((1-p)/p)1199*(a1)
Next add the above equations:
(a1200 - a1) = a1 * (((1-p)/p) + ((1-p)/p)2 + ((1-p)/p)3 + ... + ((1-p)/p)1199)
1 = a1 * (1 + ((1-p)/p) + ((1-p)/p)2 + ((1-p)/p)3 + ... + ((1-p)/p)1199)
a1 = 1 / (1 + ((1-p)/p) + ((1-p)/p)2 + ((1-p)/p)3 + ... + ((1-p)/p)1199)
a1 = ((1-p)/p - 1) / (((1-p)/p)1200 - 1)
Now that we know a1 we can find a1000:
(a2 - a1) = (1-p)/p*(a1)
(a3 - a2) = ((1-p)/p)2*(a1)
(a4 - a3) = ((1-p)/p)3*(a1)
.
.
.
(a999 - a18) = ((1-p)/p)9998*(a1)
(a1000 - a19) = ((1-p)/p)9999*(a1)
Add the above equations together:
(a1000 - a1) = a1 * (((1-p)/p) + ((1-p)/p)2 + ((1-p)/p)3 + ... + ((1-p)/p)999)
a1000 = a1 * (((1-p)/p)1000 - 1)) / ((1-p)/p - 1))
a1000 = [ ((1-p)/p - 1) / (((1-p)/p)1200 - 1) ] * [ (((1-p)/p)1000 - 1) / ((1-p)/p - 1) ]
a1000 = (((1-p)/p)1000 - 1) / (((1-p)/p)1200 - 1) =~ 0.004378132.
Given enough time, the odds are likely to catch up to the player in any game of luck and the bankroll will keep going down gradually. However, if you were to bet larger amounts your odds would be much better. The following are the odds of winning 20% before losing 100% at various units of bet size.
$5: 0.336507
$10: 0.564184
$25: 0.731927
$50: 0.785049
$100: 0.809914
For more on the math of this kind of problem, please see my MathProblems.info site, problem 116.