read HIMURIM:

N:=800: p:=3/5:
f:=2*p-1:
h:=1/80:


print(`Exit capital: N=`,N):
print(`Prob. of winning a round:`, p):

for x from N/2 by 50 to N-100 do
 print(`If the inital capital is`, x, `then the Breiman profile is`):
 print(evalf(BreimanProfile(p,N,x/N,f,h))):
od:

print(`This took`, time(), `seconds.`):

quit: