In [1]:

    

import numpy as np
from game_tools import NormalFormGame
from logitdyn import LogitDynamics


    

In [2]:

    

g_coo = NormalFormGame([[4, 0], [3, 2]])


    

In [3]:

    

ld = LogitDynamics(g_coo)


    

In [4]:

    

ld.beta = 4.0
T = 10**3
num_reps = 500
init_actions = (0, 0)
out = ld.replicate(T=T, num_reps=num_reps, init_actions=init_actions)
np.histogram2d(out[:, 0], out[:, 1], bins=np.arange(3))


    

    
Out[4]:





(array([[  12.,    0.],
        [   0.,  488.]]), array([ 0.,  1.,  2.]), array([ 0.,  1.,  2.]))

In [4]: