Here is a syntax which works:
from scipy import optimize
from math import log
bounds = ((0, 1), (0, 1), (0, 1), (0, 1))
def some_function(x):
return log(1+x[0])/(1+x[1]) + log(1+x[2])/(1+x[3])
results = dict()
print(optimize.shgo(some_function, bounds))