kode kode kode

[konduktans]

#Exercise 8.8
#Emil Johansen
from numpy import random
import sys

#Function
def game(N,s,n,q): #N=nmbr games, s=max sum, n=nmbr dice, q=game cost
wins = 0
losses = 0
sum = 0

for o in range(N):
	for i in random.randint(1,7,n):
		sum += i
	if (sum < s):
		wins += 1
	else:
		losses += 1
	sum = 0

   #Probabilities and stuff
if wins == 0:
    print 'zero wins!'
    sys.exit()
else:
    p = wins/float(N)     #Probability of winnning
    reward = q/float(p) #game cost / probability
    mag = wins*reward - N*q #money after game = income - outcome

  	win = p
lose = 1-p
exp_pay = -q*(1-p) + p*(reward-q)
return wins, losses, p, reward, mag

#Variables
N = 10000
q = 1  #units of money per game
s = 9  #set the max sum
n = 4  #number of dice

#Init and print
#probability = game(N,s,n,q)

print game(N,s,n,q)