n = 4
m = 5
c = [-2 -1 0 1 2]
P = transpose([
  .1 0.4 0 0 .5;
  .3 0 0 0 .7;
  .3 0 0 .4 .3;
  0 0 .75 .25 0
])
beta = 0.1
L = 0
