np.fromfunction(lambda x, y: (x+1)*(y+1) , (9, 9), dtype=np.int)