matematikk.net

Kode: Velg alt

import numpy as np
import matplotlib.pyplot as plt

def Trekk(N):
    b = N/2
    r = N/2
    
    blåTrukket = 0
    rødTrukket = 0

    for i in range(15):
        trekk = np.random.randint(1,N+1)
        if trekk <= b:
            blåTrukket +=1
            b -= 1
            N -= 1
        else:
            rødTrukket += 1
            N -=1
            r -=1

    return blåTrukket, rødTrukket

def simmulering(antallSim,antallBaller):
    frekvens  =  0
    for i in range(antallSim):
        b, r = Trekk(antallBaller)
        if b == 6 and r == 9:
            frekvens += 1
    sannsynlighet = frekvens/antallSim
    return sannsynlighet

def simmulering2(antallSim, slutt):
    x = []
    y = []
    for i in range(18, slutt, 2):
        x.append(i)
        y.append(simmulering(antallSim, i))
        print(i) #fremgangsvariabel for en lang simmulering
    return x, y

x, y = simmulering2(50000,50)
plt.plot(x,y,'*')
plt.show()

Kode: Velg alt

from math import factorial as fact
def bin(a , b):
  return (fact(a)//fact(b)//fact(a - b))
for i in range( -5 , 6 ):
 p = bin(15 + i,6)*bin(15 + i , 9)/bin(30 + 2*i , 15)
 print(15 + i, p)
 

Kode: Velg alt

0 0.13544891640866874
11 0.1489938080495356
12 0.15547179970386324
13 0.15880333826894602
14 0.16056781980526766
15 0.1614906233673669
16 0.16193714121999556
17 0.16210443165514019
18 0.16210443165514019
19 0.16200332702001222
20 0.16184148553447775