Eu tenho uma matriz:
U = np.array([3, 5, 7, 9, 11])
Quero obter um resultado como:
result = np.array([
[ np.nan, np.nan, np.nan, np.nan, np.nan],
[U[0] - U[1], np.nan, np.nan, np.nan, np.nan],
[U[0] - U[2], U[1] - U[2], np.nan, np.nan, np.nan],
[U[0] - U[3], U[1] - U[3], U[2] - U[3], np.nan, np.nan],
[U[0] - U[4], U[1] - U[4], U[2] - U[4], U[3] - U[4], np.nan]
])
Posso usar np.tril_indices(4, k=-1)
para obter índices do triângulo inferior sem diagonal, mas o que vem depois?