gistfile1.py

import numpy as np
from scipy.optimize import minimize
import pylab as plt
import itertools

num_sides1 = 6
num_sides2 = 10
mean1 = 4.5
mean2 = mean1


def f(x):
    return  - x*np.log(x)


def sum_form(p):
    s = 0
    for pi in p:
        # as we will to minimize, add minus
        s += - f(pi)
    return s