you could just do this
import numpy as np a = np.linspace(-2.5, 2.5, 6, endpoint=True) c = np.array([2, 1]) result = (a[:, np.newaxis] * c) print(result)