### STAKE & NUMBER_OF_BETS ###
stake_1 = 0.1
stake_2 = 0.2
stake_3 = 0.3
stake_4 = 0.4
stake_5 = 0.5

### Please uncomment wanted amount of bets:
number_of_bets = 100
#number_of_bets = 250
#number_of_bets = 500
#number_of_bets = 1000
#number_of_bets = 2500


### Please uncomment wanted value for interval on x axis:
inter = 17 #for number_of_bets == 100
#inter = 23; #for number_of_bets = 250
#inter = 34; #for number_of_bets == 500
#inter = 50; #for number_of_bets == 1000
#inter = 100; #for number_of_bets == 2500


base = 100 # any value
p = 0.6 # probability == 60%
probability = dbinom(0:number_of_bets,number_of_bets, p)
s = (0:number_of_bets)

### Please feel free to change stake percent
final_bankroll_1 = base*(1+stake_1)^s*(1-stake_1)^(number_of_bets-s) #calculation for final_bankroll_1 for stake == 10%
final_bankroll_2 = base*(1+stake_2)^s*(1-stake_2)^(number_of_bets-s) #calculation for final_bankroll_2 for stake == 20%
final_bankroll_3 = base*(1+stake_3)^s*(1-stake_3)^(number_of_bets-s) #calculation for final_bankroll_3 for stake == 30%
final_bankroll_4 = base*(1+stake_4)^s*(1-stake_4)^(number_of_bets-s) #calculation for final_bankroll_4 for stake == 40%
final_bankroll_5 = base*(1+stake_5)^s*(1-stake_5)^(number_of_bets-s) #calculation for final_bankroll_5 for stake == 50%

### Please uncomment if you want to see data.frame of final_bankroll's and probability
#print (data.frame(final_bankroll_1, probability)) #printing the aktual values for final_bankroll_1
#print (data.frame(final_bankroll_2, probability)) #printing the aktual values for final_bankroll_2
#print (data.frame(final_bankroll_3, probability)) #printing the aktual values for final_bankroll_3
#print (data.frame(final_bankroll_4, probability)) #printing the aktual values for final_bankroll_4
#print (data.frame(final_bankroll_5, probability)) #printing the aktual values for final_bankroll_5


### visualization of both bankrolls using log() function
plot(log(final_bankroll_1), probability, type="h",  xlim = c(0, inter))
par(new = TRUE)
#plot(log(final_bankroll_2), probability, type="h", col = "red", xlim = c(0, inter))
#plot(log(final_bankroll_3), probability, type="h", col = "red", xlim = c(0, inter))
#plot(log(final_bankroll_4), probability, type="h", col = "red", xlim = c(0, inter))
plot(log(final_bankroll_5), probability, type="h", col = "red", xlim = c(0, inter))