Set product could be produced as:: dp = list((xi,yi) for xi in x for yi in y) dpasarray = np.array(dp) but unless you really need the array, you can use a generator. Alan Isaac