Let x,y, and z denote the relative frequencies of individuals playing rock, paper, and scissors, respectively. Then x + y + z = 1 or z = 1 − x − y. By eliminating z in this fashion, one can capture the dynamics of the three strategies by studying x and y alone: x'=x*(f_x-p)+u*(-2*x+y+z), y'=y*(f_y-p)+u*(-2*y+x+z), f_x and f_y denote the expected fitness of individuals playing rock and paper, respectively, and p = x*f_x + y*f_y + z*f_z is the average fitness in the whole population.
Then the question is how can I plot the phase diagram like this:
I'm sorry for missing some key informations: the payoffs of the three strategies are:
f_z=-(1.0+e)*x+y，
f_x=1-x-(e+2)*y,
f_y=(e+2)*x+(e+1)*(y-1).
p = x*f_x + y*f_y + z*f_z
is the average fitness in the whole population.
According to the restriction: x+y+z=1.0, the ODEs of the system could be simplified as:
x'=x*(f_x-p)+u*(-3*x+1.0)
y'=y*(f_y-p)+u*(-3*y+1.0)
And the initial conditions are
x_0=1/3. y_0=1/3, z_0=1/3.
The values of the parameters are：
(1）u=0.4 and e=2,
(2) u=0.05 and e=5
respectively.
