consList = []
for y in range(4):
    newConstraint = constraint(y)
    consList.append(NonlinearConstraint(newConstraint.constraintFunction(), 0, 1, keep_feasible = True))
result = scipy.differential_evolution(objFunction, myBounds, constraint = consList )