dx = [3, 4, 5, 9]
dy = [5, 3, 2, 1]
output_matrix = []
 
for dx_n in dx:
    dx_row = []
    for dy_n in dy:
        dx_row.append(dx_n * dy_n)
    output_matrix.append(dx_row)
 
# Output is [[15, 9, 6, 3], [20, 12, 8, 4], [25, 15, 10, 5], [45, 27, 18, 9]]